Search NASA⌕ Search

SEARCH · Search NASA

Results for “system sizing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

Scalable Circuits for Preparing Ground States on Digital Quantum Computers: The Schwinger Model Vacuum on 100 Qubits

The vacuum of the lattice Schwinger model is prepared on up to 100 qubits of IBM’s Eagle-processor quantum computers. A new algorithm to prepare the ground state of a gapped translationally invariant system on a quantum computer is presented, which we call “scalable circuits ADAPT-VQE” (SC-ADAPT-VQE). This algorithm uses the exponential decay of correlations between distant regions of the ground state, together with ADAPT-VQE, to construct quantum circuits for state preparation that can be scaled to arbitrarily large systems. These scalable circuits can be determined with use of classical computers, avoiding the challenging task of optimizing parameterized circuits on a quantum computer. SC-ADAPT-VQE is applied to the Schwinger model, and is shown to be systematically improvable, with an accuracy that converges exponentially with circuit depth. Both the structure of the circuits and the deviations of prepared wave functions are found to become independent of the number of spatial sites, L . This allows a controlled extrapolation of the circuits, determined with use of small or modest-sized systems, to arbitrarily large L . The circuits for the Schwinger model are determined on lattices up to L = 14 (28 qubits) with the Qiskit classical simulator, and are subsequently scaled up to prepare the L = 50 (100 qubits) vacuum on IBM’s 127-superconducting-qubit quantum computers ibm_brisbane and ibm_cusco. After introduction of an improved error-mitigation technique, which we call “operator decoherence renormalization”, the chiral condensate and charge-charge correlators obtained from the quantum computers are found to be in good agreement with classical matrix product state simulations. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Acquisitive root exploration strategies help maintain higher peak sap flux rates during summer drought, but more root biomass does not

Roots are responsible for soil water uptake, yet little is known about how variation in fine-root traits relates to whole-tree water movement, particularly during periods of drought. By combining a 3-yr dataset monitoring sap flow rates with measures of fine-root biomass, length, and morphology across 10 tree species, we addressed hypotheses concerning the role of fine-root system size and morphology in determining tree responses to summer drought as well as potential changes in their relative importance under more severe droughts. Greater fine-root biomass and length did not enable trees to maintain high rates of sap flow during drought, whereas the morphological traits, specific root length and specific root area, were linked to sustained sap flow rates during drought. We found that all species, irrespective of root traits, progressively reduced their sap flow under more extreme drought conditions, although more acquisitive root morphology was still associated with smaller reductions. Our results run counter to long-standing assumptions that larger root systems are better able to access soil water and maintain photosynthetic activity during drought. Instead, we find evidence that root morphology at least partially determines the capacity for water uptake and movement as soil moisture declines.

drought↗

Role of boundary conditions in quantum computations of scattering observables

Quantum computing may offer the opportunity to simulate strongly interacting field theories, such as quantum chromodynamics, with physical time evolution. This would give access to Minkowski-signature correlators, in contrast to the Euclidean calculations routinely performed at present. However, as with present-day calculations, quantum computation strategies still require the restriction to a finite system size, including a finite, usually periodic, spatial volume. In this work, we investigate the consequences of this in the extraction of hadronic and Compton-like scattering amplitudes. Using the framework presented in Briceño et al., we estimate the volume effects for various 1+1D Minkowski-signature quantities and show that these can be a significant source of systematic uncertainty, even for volumes that are very large by the standards of present-day Euclidean calculations. We then present an improvement strategy, based in the fact that the finite volume has a reduced symmetry. This implies that kinematic points, which yield the same Lorentz invariants, may still be physically distinct in the periodic system. As we demonstrate, both numerically and analytically, averaging over such sets can significantly suppress the unwanted volume distortions and improve the extraction of the physical scattering amplitudes. As the improvement strategy is based only in kinematics, it can be applied without detailed knowledge of the system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measurements of $\pi ^\pm $, $K^\pm $, p and $\bar{p}$ spectra in $^{40}\hbox {Ar+}^{45}\hbox {Sc}$ collisions at 13A to 150A $\text{ Ge }\hspace{-1.00006pt}\text{ V }\!/\!c$

The NA61/SHINE experiment at the CERN Super Proton Synchrotron studies the onset of deconfinement in strongly interacting matter through a beam energy scan of particle production in collisions of nuclei of varied sizes. This paper presents results on inclusive double-differential spectra, transverse momentum and rapidity distributions and mean multiplicities of $\pi ^\pm $, $K^\pm $, p and $\bar{p}$ produced in $^{40}\hbox {Ar+}^{45}\hbox {Sc}$ collisions at beam momenta of 13A, 19A, 30A, 40A, 75A and 150A $\text{ Ge }\hspace{-1.00006pt}\text{ V }\!/\!c$. The analysis uses the 10% most central collisions, where the observed forward energy defines centrality. The energy dependence of the $K^\pm $/$\pi ^\pm $ ratios as well as of inverse slope parameters of the $K^\pm $ transverse mass distributions are placed in between those found in inelastic $p+p$ and central Pb + Pb collisions. The results obtained here establish a system-size dependence of hadron production properties that so far cannot be explained either within statistical or dynamical models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Luther-Emery liquid and dominant singlet superconductivity in the hole-doped Haldane spin-1 chain

We investigate the pairing tendencies in the hole-doped Haldane spin-1 chain. To allow for doping, we extend the original spin chain Hamiltonian into a fermionic model involving a two-orbital Hubbard chain at intermediate or strong repulsive interaction strengths U and for degenerate orbitals. At half filling and large U, the ferromagnetic Hund's coupling, J H , generates effective spin-1 moments, with antiferromagnetic correlations between sites. Using large-scale density matrix renormalization group calculations, we accurately study the system's behavior under light hole-doping. For U = 1.6 in units of the noninteracting bandwidth and for J H /U ≳ 0.275, we find that singlet pairing dominates the long-distance physics, establishing this system as a promising platform for repulsively mediated superconductivity. We provide concrete examples of materials that could realize the physics described here. We also provide evidence that the system approaches a Luther-Emery liquid state at large system sizes, reminiscent of the behavior of doped one-orbital two-leg ladders at weak coupling, which also have superconducting tendencies. The numerically calculated central charge approaches one in the thermodynamic limit, indicating a single gapless mode as is expected for the Luther-Emery state. Exponents characterizing the power-law decays of singlet pair-pair and charge density-density correlations are determined, and found to approximately satisfy the Luther-Emery identity.

1-dimensional systems↗

Small-world complex network generation on a digital quantum processor

Abstract Quantum cellular automata (QCA) evolve qubits in a quantum circuit depending only on the states of their neighborhoods and model how rich physical complexity can emerge from a simple set of underlying dynamical rules. The inability of classical computers to simulate large quantum systems hinders the elucidation of quantum cellular automata, but quantum computers offer an ideal simulation platform. Here, we experimentally realize QCA on a digital quantum processor, simulating a one-dimensional Goldilocks rule on chains of up to 23 superconducting qubits. We calculate calibrated and error-mitigated population dynamics and complex network measures, which indicate the formation of small-world mutual information networks. These networks decohere at fixed circuit depth independent of system size, the largest of which corresponding to 1,056 two-qubit gates. Such computations may enable the employment of QCA in applications like the simulation of strongly-correlated matter or beyond-classical computational demonstrations.

97 MATHEMATICS AND COMPUTING↗

Dynamics of disordered mechanical systems with large connectivity, free probability theory, and quasi-Hermitian random matrices

Disordered mechanical systems with high connectivity represent a limit opposite to the more familiar case of disordered crystals. Individual ions in a crystal are subjected essentially to nearest-neighbor interactions. In contrast, the systems studied in this paper have all their degrees of freedom coupled to each other. Thus, the problem of linearized small oscillations of such systems involves two full positive-definite and non-commuting matrices, as opposed to the sparse matrices associated with disordered crystals. Consequently, the familiar methods for determining the averaged vibrational spectra of disordered crystals, introduced many years ago by Dyson and Schmidt, are inapplicable for highly connected disordered systems. In this paper we apply random matrix theory (RMT) to calculate the averaged vibrational spectra of such systems, in the limit of infinitely large system size. At the heart of our analysis lies a calculation of the average spectrum of the product of two positive definite random matrices by means of free probability theory techniques. We also show that this problem is intimately related with quasi-hermitian random matrix theory (QHRMT), which means that the ‘hamiltonian’ matrix is hermitian with respect to a non-trivial metric. This extends ordinary hermitian matrices, for which the metric is simply the unit matrix. The analytical results we obtain for the spectrum agree well with our numerical results. The latter also exhibit oscillations at the high-frequency band edge, which fit well the Airy kernel pattern. We also compute inverse participation ratios of the corresponding amplitude eigenvectors and demonstrate that they are all extended, in contrast with conventional disordered crystals. Finally, we compute the thermodynamic properties of the system from its spectrum of vibrations. In addition to matrix model analysis, we also study the vibrational spectra of various multi-segmented disordered pendula, as concrete realizations of highly connected mechanical systems. A universal feature of the density of vibration modes, common to both pendula and the matrix model, is that it tends to a non-zero constant at vanishing frequency.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Relaxation mechanisms in a disordered system with Poisson-level statistics

Here we discuss the interplay between many-body localization and spin symmetry. To this end, we study the time evolution of several observables in the anisotropic t- j model. Like the Hubbard chain, the studied model contains charge and spin degrees of freedom, yet it has smaller Hilbert space and thus allows for numerical studies of larger systems. We compare the field disorder that breaks the $\mathbb{Z}_2$ spin symmetry and a potential disorder that preserves the latter symmetry. In the former case, sufficiently strong disorder leads to localization of all studied observables, at least for the studied system sizes. However, in the case of symmetry-preserving disorder, we observe that odd operators under the $\mathbb{Z}_2$ spin transformation relax towards the equilibrium value at relatively short timescales that grow only polynomially with the disorder strength. On the other hand, the dynamics of even operators and the level statistics within each symmetry sector are consistent with localization. Our results indicate that localization exists within each symmetry sector for symmetry-preserving disorder. Odd operators' apparent relaxation is due to their time evolution between distinct symmetry sectors.

36 MATERIALS SCIENCE↗

Principles and Options for Designing Battery Energy Storage Zoning Ordinances

Deployment of battery energy storage (BESS) systems, both standalone and as part of hybrid systems paired with generation, has rapidly increased in the United States in recent years as utilities and communities have deployed storage to improve electric grid reliability and act as a cost-effective alternative to larger grid infrastructure. The modular nature of BESS technologies means systems may be built near other existing land uses, creating the potential for conflicts with neighboring landowners that can be managed and mitigated through zoning and permitting requirements established by local jurisdictions. While many cities and counties have adopted ordinances specific to BESS into their local zoning codes, these ordinances vary significantly in their requirements and level of detail. Meanwhile, many other jurisdictions, including those home to proposed or existing BESS projects, lack any specific language related to BESS in their zoning codes. Local planning and zoning officials have limited capacity and may lack the familiarity with BESS technologies needed to develop ordinances or otherwise make reasonable zoning decisions that balance safety, community impacts, and other goals. The resulting uncertainty at the local zoning level has led developers to withdraw projects in some areas and has spurred moratoria or bans on energy storage projects in others. This report intends to provide practical resources for practitioners interested in reasonable and effective local regulation of battery energy storage. It does not present a model zoning ordinance, but rather is intended to complement model ordinances developed by others by providing additional context and analysis regarding the structure of energy storage zoning ordinances and the decision points for local officials. Zoning ordinances at the city, town, and county level across the U.S. were surveyed alongside two template model ordinances to identify common elements and options for regulating the zoning and siting of BESS. Common elements identified and analyzed include definitions and general requirements, including cutoffs or tiers used to apply regulations to different system sizes and the permitted zones where jurisdictions allow BESS to be sited; visual, noise, and aesthetic requirements, including property line setbacks, fencing and visual screening, noise, and lighting requirements; and safety and planning requirements, such as site plans, decommissioning plans or funds, and requirements for access by emergency services. The report also summarizes some of the more unique regulations, including those that place additional restrictions on BESS at the local level.

25 ENERGY STORAGE↗

Power Electronics Manufacturing Improvements for Heavy-Duty Fuel Cell Vehicles

The Marel Power Solutions project, funded by the U.S. Department of Energy under Award DE-SC0023801, focused on advancing manufacturing techniques for power electronics in heavy-duty fuel cell vehicles. The research aimed to enhance system efficiency, reduce costs, and support broader adoption of hydrogen fuel cell technology. Key areas of investigation included power topology, thermal modeling, system architecture, and accessibility through software tools. Key Accomplishments: 1. Power Topology: - Developed an interleaved boost converter with optimized phase count, leveraging Marel’s proprietary Power Stacks. - Achieved reduced parasitic inductance and resistance, enabling high efficiency in DC-DC converters. 2. Thermal Modeling: - Integrated innovative cooling systems into compact Silicon Carbide (SiC) modules. - Simulations demonstrated the ability to dissipate significant heat (up to 7.5 kW), ensuring device reliability under heavy loads. 3. System Architecture: - Utilized simulation tools to analyze the impact of various fuel cell and vehicle parameters on efficiency. - Highlighted the role of smaller, modular improvements, such as enhanced DC-DC converters, in achieving system-wide gains. 4. Accessibility: - Evaluated and implemented MATLAB/Simulink code generation tools for real-world hardware applications. - Demonstrated the potential for rapid prototyping of custom power systems with reduced development costs. Impact and Benefits: - Efficiency and Cost Reduction: Marel’s cooling technology enhances SiC die performance, reducing the number of dies required and overall system size. - Scalability and Flexibility: The innovations support tailored solutions for diverse applications, from mass transit to mining vehicles. - Sustainability: The research promotes the integration of electrification technologies, helping meet rising energy demands sustainably. Conclusion: The project’s outcomes advance the state of power electronics for hydrogen fuel cell vehicles, enabling more efficient, compact, and cost-effective solutions. These developments lay a foundation for future innovation, contributing to the broader adoption of clean energy technologies in transportation and other industries.

08 HYDROGEN↗

Skewness and kurtosis of mean transverse momentum fluctuations at the LHC energies

The first measurements of skewness and kurtosis of mean transverse momentum ($\langle$$p_T\rangle$) fluctuations are reported in Pb–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV, Xe–Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV and pp collisions at $\sqrt{s}$ = 5.02 TeV using the ALICE detector. The measurements are carried out as a function of system size $\langle$dN ch /dη$\rangle$$^{1/3}_{|η|<0.5}$, using charged particles with transverse momentum ($p_T$) and pseudorapidity (η), in the range 0.2<$p_T$<3.0 GeV/c and |η|<0.8, respectively. In Pb–Pb and Xe–Xe collisions, positive skewness is observed in the fluctuations of $\langle$$p_T\rangle$ for all centralities, which is significantly larger than what would be expected in the scenario of independent particle emission. This positive skewness is considered a crucial consequence of the hydrodynamic evolution of the hot and dense nuclear matter created in heavy-ion collisions. Furthermore, similar observations of positive skewness for minimum bias pp collisions are also reported here. Kurtosis of $\langle$$p_T\rangle$ fluctuations is found to be in good agreement with the kurtosis of Gaussian distribution, for most central Pb–Pb collisions. Hydrodynamic model calculations with MUSIC using Monte Carlo Glauber initial conditions are able to explain the measurements of both skewness and kurtosis qualitatively from semicentral to central collisions in Pb–Pb system. Color reconnection mechanism in PYTHIA8 model seems to play a pivotal role in capturing the qualitative behavior of the same measurements in pp collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Accurate localization of Kosterlitz-Thouless-type quantum phase transitions for one-dimensional spinless fermions

We investigate the charge-density wave (CDW) transition for one-dimensional spinless fermions at half band filling with nearest-neighbor electron transfer amplitude t and interaction V. The model is equivalent to the anisotropic XXZ Heisenberg model for which the Bethe Ansatz provides an exact solution. For V>V c =2t, the CDW order parameter and the single-particle gap are finite but exponentially small, as is characteristic for a Kosterlitz-Thouless transition. It is notoriously difficult to locate such infinite-order phase transitions in the phase diagram using approximate analytical and numerical approaches. Second-order Hartree-Fock theory is qualitatively applicable for all interaction strengths, and predicts the CDW transition to occur at V$^{(2)}_{c,2}$≈1.5t. Second-order Hartree Fock theory is almost variational because the density of quasiparticle excitations is small. We apply the density-matrix renormalization group (DMRG) for periodic boundary conditions for system sizes up to 514 sites, which permits a reliable extrapolation of all physical quantities to the thermodynamic limit, apart from the critical region. We investigate the ground-state energy, the gap, the order parameter, the momentum distribution, the quasiparticle density, and the density-density correlation function to locate V c from the DMRG data. In conclusion, tracing the breakdown of the Luttinger liquid and the peak in the quasiparticle density at the band edge permits us to reproduce V c with an accuracy of one percent.

1-dimensional spin chains↗

Classification of dynamical Lie algebras generated by spin interactions on undirected graphs

Dynamical Lie algebras (DLAs) are a versatile tool for various topics that span from the expressibility-trainability of variational quantum algorithms (VQAs), to simulation of many body Hamiltonians. Quantum gates and most of the Hamiltonians of interest consist of local interactions; therefore, the analysis of all possible DLAs generated by 1- and 2-local operators is crucial for quantum simulation and VQAs on current hardware. Previously in [R. Wiersema et al ., npj Quantum Inf. 10 , 110 (2024)], we analyzed the DLAs on linear, circular and all-to-all topologies, and obtained results about their dimensions and algebraic structure. Here, in this work, we extend our analysis into any possible hardware topology and provide a classification of all DLAs generated by Pauli strings on any undirected interaction graph. Our results indicate that the DLAs depend solely on whether the connectivity or interaction graph is bipartite or not. In addition, we find that the non-trivial polynomially scaling DLAs appear only on 1D line or circle topologies, and all other DLAs have dimensions scaling exponentially with the system size. Together with the current VQA literature, our results imply that either the majority of VQAs are non-trainable, or we are yet to understand the role of DLAs on the trainability of VQAs.

Algebraic structures↗

Preliminary Thermal-Hydraulic Analysis of the UTA–2 Subcritical Linear Accelerator–Driven System

The National Nuclear Security Administration’s mission includes establishing a reliable supply of 99 Mo without highly enriched uranium. Oak Ridge National Laboratory (ORNL) supports this objective through collaborative research and development with industrial partners. Niowave Inc., a current partner, is currently designing a subcritical linear accelerator-driven system (ADS) and preparing for the US Nuclear Regulatory Commission’s licensing process. Niowave’s system has the potential to efficiently supply medical radioisotopes. The technology includes a superconducting electron accelerator and a pile of both natural and low-enriched uranium (LEU) targets. The process fissions uranium and many valuable isotopes can then be extracted from the targets. Niowave is currently iterating through conceptual and detailed design processes for several system sizes. This report discusses UTA–2, which is at the demonstration stage. UTA–2 will validate numerical modeling results with experimental measurements before progressing to the detailed design of UTA–3, the commercial-sized ADS. The thermal-hydraulic behavior of the UTA–2 core design was numerically investigated using STAR-CCM+, as described in this report. STAR-CCM+, a state-of-the-art computational fluid dynamics (CFD) software that was commercially developed by Siemens, has an extensive user base and a set of validation studies. It is also compliant with the American Society of Mechanical Engineers’ Nuclear Quality Assurance 1 standard. Three cases are investigated in this report. Case 1 quantified the temperature field within the UTA–2 assembly and water tank using only conduction as the method of thermal energy transport. This simplified approach was overly conservative and yielded wetted cladding temperatures above the coolant saturation temperature. In this case, the maximum temperature of the wetted cladding surface of the highest power rod exceeded the saturation temperature of water by 63.8°C. Because of the overly conservative approach taken in Case 1 and its negative subcooled margin, buoyancy-driven natural circulation flow physics were implemented in Case 2. Adding coolant motion significantly distributed the thermal energy of the system through convective heat transfer. This relatively small amount of convective heat transfer significantly reduced system temperatures and increased the subcooled margin from –63.8 to 61.3°C. This margin confirmed that no boiling was expected during normal operation of UTA–2 at 230 W. Case 3 had no additional physics models but considered an overpower event during which the power of each LEU and natural uranium rod was at its respective peak values. This resulted in a study with the total assembly power equal to 176% of the nominal power of 230 W considered in Cases 1 and 2. The resulting natural circulation flows were slightly enhanced. The subcooled margin decreased slightly to 46.3°C, which is still a significant margin to local boiling of the water within the tank. This margin confirmed that no boiling was expected during an abnormal operation of UTA–2 at 406 W.

07 ISOTOPE AND RADIATION SOURCES↗

Approximate two-body generating Hamiltonian for the particle-hole Pfaffian wave function

We present two two-body Hamiltonians that approximate the exact particle-hole Pfaffian wave function with their ground states for all the system sizes where this wave function has been numerically constructed to date. The approximate wave functions have high overlap with the original and reproduce well the low-lying entanglement spectrum and structure factor. The approximate generating Hamiltonians are obtained by an optimization procedure where three to four pseudopotentials are varied in the neighbourhood of second Landau level Coulomb interaction or of a noninteracting model. They belong to a finite region in the variational space of Hamiltonians where each point approximately generates the particle-hole Pfaffian. Here we diagonalize the identified Hamiltonians for up to 20 electrons and find that for them the particle-hole Pfaffian shift appears energetically more favorable. The possibility to interpret the data in terms of composite fermions is discussed.

36 MATERIALS SCIENCE↗

A simplified, efficient approach to hybrid wind and solar plant site optimization

Abstract. Wind plant layout optimization is a difficult, complex problem with a large number of variables and many local minima. Layout optimization only becomes more difficult with the addition of solar generation. In this paper, we propose a parameterized approach to wind and solar hybrid power plant layout optimization that greatly reduces problem dimensionality while guaranteeing that the generated layouts have a desirable regular structure. Thus far, hybrid power plant optimization research has focused on system sizing. We go beyond sizing and present a practical approach to optimizing the physical layout of a wind–solar hybrid power plant. We argue that the evolution strategy class of derivative-free optimization methods is well-suited to the parameterized hybrid layout problem, and we demonstrate how hard layout constraints (e.g., placement restrictions) can be transformed into soft constraints that are amenable to optimization using evolution strategies. Next, we present experimental results on four test sites, demonstrating the viability, reliability, and effectiveness of the parameterized evolution strategy approach for generating optimized hybrid plant layouts. Completing the tool kit for parameterized layout generation, we include a brief tutorial describing how the parameterized evolutionary approach can be inspected, understood, and debugged when applied to hybrid plant layouts.

14 SOLAR ENERGY↗

Lifetime of Almost Strong Edge-Mode Operators in One-Dimensional, Interacting, Symmetry Protected Topological Phases

Almost strong edge-mode operators arising at the boundaries of certain interacting one-dimensional symmetry protected topological phases with Z 2 symmetry have infinite temperature lifetimes that are nonperturbatively long in the integrability breaking terms, making them promising as bits for quantum information processing. We extract the lifetime of these edge-mode operators for small system sizes as well as in the thermodynamic limit. For the latter, a Lanczos scheme is employed to map the operator dynamics to a one-dimensional tight-binding model of a single particle in Krylov space. We suggest this model to be that of a spatially inhomogeneous Su-Schrieffer-Heeger model with a hopping amplitude that increases away from the boundary, and a dimerization that decreases away from the boundary. We associate this dimerized or staggered structure with the existence of the almost strong mode. Thus, the short time dynamics of the almost strong mode is that of the edge mode of the Su-Schrieffer-Heeger model, while the long time dynamics involves decay due to tunneling out of that mode, followed by chaotic operator spreading. Furthermore, we show that competing scattering processes can lead to interference effects that can significantly enhance the lifetime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Stabilizing Non-Abelian Topological Order Against Heralded Noise via Local Lindbladian Dynamics

An important open question for the current generation of highly controllable quantum devices is understanding which phases can be realized as stable steady states under local quantum dynamics. In this work, we show how robust steady-state phases with both Abelian and non-Abelian mixed-state topological order can be stabilized, in two spatial dimensions, against generic “heralded” noise using active dynamics that incorporate measurement and feedback, modeled as a fully local Lindblad master equation. These topologically ordered steady states are two-way connected to pure topologically ordered ground states using local quantum channels, and preserve quantum information for a time that is exponentially large in the system size. Specifically, we present explicit constructions of families of local Lindbladians for both Abelian (ℤ 2 ) and non-Abelian (𝐷 4 ) topological order whose steady states host mixed-state topological order when the noise is below a threshold strength. As the noise strength is increased, these models exhibit first-order transitions to intermediate mixed-state phases where they encode robust classical memories, followed by (first-order) transitions to a trivial steady state at high noise rates. When the noise is imperfectly heralded, steady-state order disappears but our active dynamics significantly enhances the lifetime of the encoded logical information. To carry out the numerical simulations for the non-Abelian 𝐷 4 case, we introduce a generalized stabilizer tableau formalism that permits efficient simulation of the non-Abelian Lindbladian dynamics.

Monte Carlo methods↗