Dissipation-Induced Order: The S=1/2 Quantum Spin Chain Coupled to an Ohmic Bath.

We consider an S=1/2 antiferromagnetic quantum Heisenberg chain where each site is coupled to an independent bosonic bath with ohmic dissipation. The coupling to the bath preserves the global SO(3) spin symmetry. Using large-scale, approximation-free quantum Monte Carlo simulations, we show that any finite coupling to the bath suffices to stabilize long-range antiferromagnetic order. This is in stark contrast to the isolated Heisenberg chain where spontaneous breaking of the SO(3) symmetry is forbidden by the Mermin-Wagner theorem. A linear spin-wave theory analysis confirms that the memory of the bath and the concomitant retarded interaction stabilize the order. For the Heisenberg chain, the ohmic bath is a marginal perturbation so that exponentially large system sizes are required to observe long-range order at small couplings. Below this length scale, our numerics is dominated by a crossover regime where spin correlations show different power-law behaviors in space and time. We discuss the experimental relevance of this crossover phenomena.

Possible Inversion Symmetry Breaking in the S=1/2 Pyrochlore Heisenberg Magnet.

We address the ground-state properties of the long-standing and much-studied three-dimensional quantum spin liquid candidate, the S=1/2 pyrochlore Heisenberg antiferromagnet. By using SU(2) density-matrix renormalization group (DMRG), we are able to access cluster sizes of up to 128 spins. Our most striking finding is a robust spontaneous inversion symmetry breaking, reflected in an energy density difference between the two sublattices of tetrahedra, familiar as a starting point of earlier perturbative treatments. We also determine the ground-state energy, E_{0}/N_{sites}=-0.490(6)J, by combining extrapolations of DMRG with those of a numerical linked cluster expansion. These findings suggest a scenario in which a finite-temperature spin liquid regime gives way to a symmetry-broken state at low temperatures.

Hierarchy of Relaxation Timescales in Local Random Liouvillians.

To characterize the generic behavior of open quantum systems, we consider random, purely dissipative Liouvillians with a notion of locality. We find that the positivity of the map implies a sharp separation of the relaxation timescales according to the locality of observables. Specifically, we analyze a spin-1/2 system of size ℓ with up to n-body Lindblad operators, which are n local in the complexity-theory sense. Without locality (n=ℓ), the complex Liouvillian spectrum densely covers a "lemon"-shaped support, in agreement with recent findings [S. Denisov et al., Phys. Rev. Lett. 123, 140403 (2019)PRLTAO0031-900710.1103/PhysRevLett.123.140403]. However, for local Liouvillians (n<ℓ), we find that the spectrum is composed of several dense clusters with random matrix spacing statistics, each featuring a lemon-shaped support wherein all eigenvectors correspond to n-body decay modes. This implies a hierarchy of relaxation timescales of n-body observables, which we verify to be robust in the thermodynamic limit. Our findings for n locality generalize immediately to the case of spatial locality, introducing further splitting of timescales due to the additional structure.

How a Small Quantum Bath Can Thermalize Long Localized Chains.

We investigate the stability of the many-body localized phase for a system in contact with a single ergodic grain modeling a Griffiths region with low disorder. Our numerical analysis provides evidence that even a small ergodic grain consisting of only three qubits can delocalize a localized chain as soon as the localization length exceeds a critical value separating localized and extended regimes of the whole system. We present a simple theory, consistent with De Roeck and Huveneers's arguments in [Phys. Rev. B 95, 155129 (2017)PRBMDO2469-995010.1103/PhysRevB.95.155129] that assumes a system to be locally ergodic unless the local relaxation time determined by Fermi's golden rule is larger than the inverse level spacing. This theory predicts a critical value for the localization length that is perfectly consistent with our numerical calculations. We analyze in detail the behavior of local operators inside and outside the ergodic grain and find excellent agreement of numerics and theory.

Anomalous Thermalization in Ergodic Systems.

It is commonly believed that quantum isolated systems satisfying the eigenstate thermalization hypothesis (ETH) are diffusive. We show that this assumption is too restrictive since there are systems that are asymptotically in a thermal state yet exhibit anomalous, subdiffusive thermalization. We show that such systems satisfy a modified version of the ETH ansatz and derive a general connection between the scaling of the variance of the off-diagonal matrix elements of local operators, written in the eigenbasis of the Hamiltonian, and the dynamical exponent. We find that for subdiffusively thermalizing systems the variance scales more slowly with system size than expected for diffusive systems. We corroborate our findings by numerically studying the distribution of the coefficients of the eigenfunctions and the off-diagonal matrix elements of local operators of the random field Heisenberg chain, which has anomalous transport in its thermal phase. Surprisingly, this system also has non-Gaussian distributions of the eigenfunctions, thus, directly violating Berry's conjecture.

Critical properties of the superfluid-bose-glass transition in two dimensions.

We investigate the superfluid (SF) to Bose-glass (BG) quantum phase transition using extensive quantum Monte Carlo simulations of two-dimensional hard-core bosons in a random box potential. T=0 critical properties are studied by thorough finite-size scaling of condensate and SF densities, both vanishing at the same critical disorder Wc=4.80(5). Our results give the following estimates for the critical exponents: z=1.85(15), ν=1.20(12), η=-0.40(15). Furthermore, the probability distribution of the SF response P(lnρSF) displays striking differences across the transition: while it narrows with increasing system sizes L in the SF phase, it broadens in the BG regime, indicating an absence of self-averaging, and at the critical point P(lnρSF+zlnL) is scale invariant. Finally, high-precision measurements of the local density rule out a percolation picture for the SF-BG transition.

Universal behavior beyond multifractality in quantum many-body systems.

How many states of a configuration space contribute to a wave function? Attempts to answer this ubiquitous question have a long history in physics and are keys to understanding, e.g., localization phenomena. Beyond single-particle physics, a quantitative study of the ground state complexity for interacting many-body quantum systems is notoriously difficult, mainly due to the exponential growth of the configuration (Hilbert) space with the number of particles. Here we develop quantum Monte Carlo schemes to overcome this issue, focusing on Shannon-Rényi entropies of ground states of large quantum many-body systems. Our simulations reveal a generic multifractal behavior while the very nature of quantum phases of matter and associated transitions is captured by universal subleading terms in these entropies.