We investigate the possibility of engineering dissipatively long-range order that is robust against heating in strongly interacting fermionic systems, relevant for atoms in cavity QED. It was previously shown [Tindall ] that it is possible to stabilize long-range order in a Hubbard model by exploiting a dissipative mechanism in the Lindblad limit, this latter being valid for spectrally unstructured baths. Here, we first show that this mechanism still holds when including additional spin-exchange interactions in the model, that is, for the tUJ model. Moreover, by means of a Redfield approach that goes beyond the Lindblad case, we show how the stability of the engineered state depends crucially on properties of the bath spectral density and discuss the feasibility of those properties in an experiment.
We investigate the possibility of engineering dissipatively long-range order that is robust against heating in strongly interacting fermionic systems, relevant for atoms in cavity QED. It was previously shown [Tindall ] that it is possible to stabilize long-range order in a Hubbard model by exploiting a dissipative mechanism in the Lindblad limit, this latter being valid for spectrally unstructured baths. Here, we first show that this mechanism still holds when including additional spin-exchange interactions in the model, that is, for the tUJ model. Moreover, by means of a Redfield approach that goes beyond the Lindblad case, we show how the stability of the engineered state depends crucially on properties of the bath spectral density and discuss the feasibility of those properties in an experiment.
Can superconducting order be engineered as a robust attractor of open-system dynamics in strongly correlated systems? We demonstrate this possibility by proposing a minimal dissipation-engineering protocol for the particle-hole symmetric Hubbard model. By applying a rotated quantum jump operator, specifically a locally transformed η-pair lowering operator, on a single lattice site only, we show that the Lindblad evolution autonomously pumps the system from the vacuum into a nonequilibrium steady state (NESS) exhibiting macroscopic η-pair off-diagonal long-range order (ODLRO). Crucially, this local-to-global synchronization stands in stark contrast to schemes reliant on spatially extensive reservoirs: here, a single local dissipative seed suffices to establish long-range coherence across the entire interacting lattice system. We elucidate the underlying mechanism via three core features: local dark-state selection, the controlled elimination of off-manifold excursions induced by hopping, and a Liouvillian invariant-subspace structure that yields an attractive fixed point with a finite dissipative gap. Furthermore, we systematically classify the stability of this NESS with respect to static disorder, and identify a wide regime in which the superconducting attractor remains robust against Hamiltonian perturbations that preserve the effective subspace structure. We also pinpoint specific perturbations that directly cause dephasing of the -pseudospin coherence and suppress ODLRO. These findings open up a disorder-tolerant pathway for stabilizing superconducting order as a non-thermal attractor through minimal local quantum-jump control.
Understanding the stability of strongly correlated phases of matter when coupled to environmental degrees of freedom is crucial for identifying the conditions under which these states may be observed. Here, we focus on the paradigmatic one-dimensional Bose-Hubbard model, and study the stability of the Luttinger-liquid and Mott insulating phases in the presence of local particle exchange with site-independent baths of noninteracting bosons. We perform a numerically exact analysis of this model by adapting the recently developed wormhole quantum Monte Carlo method for retarded interactions to a continuous-time formulation with worm updates; we show how the wormhole updates can be easily implemented in this scheme. For an Ohmic bath, our numerical findings confirm the scaling prediction that the Luttinger-liquid phase becomes unstable at infinitesimal bath coupling. We show that the ensuing phase is a long-range-ordered superfluid with spontaneously broken U(1) symmetry. While the Mott insulator remains a distinct phase for small bath coupling, it exhibits diverging compressibility and noninteger local boson occupation in the thermodynamic limit. Upon increasing the bath coupling, this phase undergoes a transition to a long-range-ordered superfluid. Finally, we discuss the effects of super-Ohmic dissipation on the Luttinger-liquid phase. Our results are compatible with a stable dissipative Luttinger-liquid phase that transitions to a long-range-ordered superfluid at a finite system-bath coupling.
Long-range interactions can protect dissipative order against fluctuations, but determining when they do so in two dimensions remains out of reach for existing many-body methods. We develop a continuous SU(2)-covariant truncated-Wigner framework that maps the problem onto stochastic dynamics with one Bloch sphere per site. Derived systematically from the exact Stratonovich-Weyl evolution of the Lindblad equation, it generates both drift and reservoir noise for arbitrary spin S. Benchmarks against exact and tensor-network calculations validate its description of non-equilibrium dynamics and steady states. Applied to 400-site two-dimensional dissipative anisotropic XY lattices, our method reveals physics missed by mean-field and perturbative approaches: vortex-antivortex proliferation destroys the predicted staggered-XY phase for nearest-neighbour coupling, whereas algebraically decaying interactions stabilize an extended Z2-ordered phase. We obtain the first beyond-mean-field steady-state phase diagram of the long-range interacting model, revealing reentrant loss of order at strong coupling and closure of the ordered region beyond α≈3.5. Our new method provides a scalable route to strongly correlated driven-dissipative systems with long-range interactions, higher-dimensional geometries, spin-boson couplings, and spins beyond one half, relevant to modern experiments with magnetic atoms and polar molecules in optical lattices and tweezer arrays.
We study the spectral and transport properties of a Mott insulator driven by a static electric field into a nonequilibrium steady state. For the dissipation, we consider two mechanisms: wide band fermion reservoirs and phonons included within the Migdal approximation. The electron correlations are treated via nonequilibrium dynamical mean field theory with an impurity solver suitable for strong correlations. We find that dissipation via phonons is limited to restricted ranges of field values around Wannier-Stark resonances. To cover the full range of field strengths, we allow for a small coupling with fermionic baths, which stabilizes the solution. When considering both dissipation mechanisms, we find that phonons enhance the current for field strengths close to half of the gap while lowering it at the gap resonance as compared to the purely electronic dissipation used by Murakami and Werner [Murakami and Werner, Phys. Rev. B 98, 075102 (2018)]. Once phonons are the only dissipation mechanism, the current in the metallic phase is almost one order of magnitude smaller than the typical values obtained by coupling to a fermionic bath. In this case, the transport regime is characterized by an accumulation of charge in the upper Hubbard band described by two effective chemical potentials.
Strongly interacting fermions underpin some of the most challenging problems in condensed matter physics, such as high-temperature superconductivity. The low-energy states of these systems encode their essential microscopic properties, yet remain largely inaccessible to classical methods. Quantum simulation offers a promising path forward, and among state-preparation strategies, engineered dissipation has emerged as a particularly compelling approach. Existing cooling protocols, however, typically rely on knowledge of the quasiparticle spectrum or mappings to free-fermion limits. In this letter, we introduce a randomized, symmetry-preserving cooling algorithm that requires no spectral information, using only local coupling operators to ancilla degrees of freedom with randomly sampled energy splittings to drive generic fermionic systems toward their low-energy manifold. We benchmark the protocol on canonical correlated fermionic models relevant to high-temperature superconductors, spanning metallic, density-wave, paired, superconducting, and phase-separated phases. Across all models, we observe universal cooling behavior: monotonic energy relaxation, concentration of spectral weight at low energies, and stabilization of correlated ground-state order. Our results establish randomized dissipative cooling as a general strategy for preparing strongly correlated fermionic states on programmable quantum devices.
Floquet engineering, the control of a quantum system by means of time-periodic driving, allows to modify the properties of a system so that it is described by an approximate effective time-independent Hamiltonian. However, in the presence of interactions the stabilization of interesting many-body ground states of such effective Hamiltonians is possible only on a certain time scale, beyond which Floquet heating sets in, as it results from unwanted driving induced resonant excitation. Moreover, the preparation of the ground state of such Floquet engineered effective Hamiltonians usually has to be accomplished adiabatically, which is another source of heating, especially when a phase transition has to be passed in a large system. Here, we propose a general dissipative strategy for the preparation and stabilization of effective ground states that are protected by an energy gap, such as topologically ordered states. It is scalable and relies on coupling the driven system to a thermal bath, the properties of which are chosen so that it both suppresses Floquet heating and guides the system into a non-equilibrium steady state with a large occupation of the effective ground-state. The small residual occupations of effective excited states are non-thermal, reflecting the non-equilibrium nature of the scheme. We use the Floquet-Born-Markov master equation to verify the proposed strategy. In particular, we discuss the preparation of both a one-dimensional Mott insulator and a two-dimensional fractional Chern insulator in strongly driven Bose-Hubbard systems.