Frustration, bond disorder, and long-range interactions can produce novel, highly nontrivial phases of matter. This thesis investigates a bond-disordered, dipolar-interacting XY spin model, studied through both numerical simulations and experiments on a Rydberg quantum simulator. Three major results are achieved: (i) We analyze energetic-magnetic hysteresis in two disorder configurations, and find pronounced hysteresis in the strongly disordered configuration for low energies. This is a first experimental evidence of energetic-magnetic hysteresis in putative isolated glasses. (ii) We develop an extension of the generalized Kibble-Zurek mechanism for reverse quenches, which allows to characterize putative spin glass quantum phase transitions by global magnetization measurements. We validate this extension numerically exactly for one-dimensional systems, and subsequently experimentally extract a critical exponent consistent with spin glass critical exponents measured in other systems. This is a first tentative experimental evidence of a spin glass phase in an isolated dipolar interacting spin system. (iii) We studied time reversal on a Rydberg quantum simulator and showed numerically that, in bond-disordered power-law interacting models, time-reversal–based protocols reveal a localization mechanism distinct from conventional many-body localization at finite sizes.
Frustration, bond disorder, and long-range interactions can produce novel, highly nontrivial phases of matter. This thesis investigates a bond-disordered, dipolar-interacting XY spin model, studied through both numerical simulations and experiments on a Rydberg quantum simulator. Three major results are achieved: (i) We analyze energetic-magnetic hysteresis in two disorder configurations, and find pronounced hysteresis in the strongly disordered configuration for low energies. This is a first experimental evidence of energetic-magnetic hysteresis in putative isolated glasses. (ii) We develop an extension of the generalized Kibble-Zurek mechanism for reverse quenches, which allows to characterize putative spin glass quantum phase transitions by global magnetization measurements. We validate this extension numerically exactly for one-dimensional systems, and subsequently experimentally extract a critical exponent consistent with spin glass critical exponents measured in other systems. This is a first tentative experimental evidence of a spin glass phase in an isolated dipolar interacting spin system. (iii) We studied time reversal on a Rydberg quantum simulator and showed numerically that, in bond-disordered power-law interacting models, time-reversal–based protocols reveal a localization mechanism distinct from conventional many-body localization at finite sizes.
Bosonic excitations within long-range ordered, but strongly inhomogeneous phases have been studied in quite some detail. My thesis focuses instead on understanding the insulating, localized phase of disordered bosonic systems. In particular I study localization properties of strongly interacting bosons and spin systems in a disorder potential at zero temperature. I focus on simple, prototypical spin models (Ising model and XY model) in random fields on a Cayley tree with large connectivity. Regarding the nature of the quantum phase transition in strong disorder I find the following results: i) With a uniformly distributed disorder non-extensive excitations in the disordered phase are all localized. ii)Moreover, I find that the order arises due to a collective condensation, which is qualitatively distinct from a Bose Einstein condensation of single particle excitations into a delocalized state. In particular, in non-frustrated Bose glasses, I do not find evidence for a boson mobility edge in the Bose glass. These results are qualitatively different from claims in the recent literatures . Considering that (many body) localization of bosons is a kind of quantum glass transition, it is an interesting question to ask what phenomena occur, if the ingredients for more conventional (classical) glassy physics are added to a disordered bosons system, namely: random, frustrated interactions between the bosons. One can still think about such a system as bosons in a disordered potential, where the disordered potential is, at least partly, self-generated by random frustrated interactions between the bosons. This question takes us to the study of another type of disordered systems: glassy systems. Those are typically characterized by low temperature phases with an inhomogeneous density or magnetization pattern, which is extremely long-lived due to the occurrence of non-trivial ergodicity breaking. I study a solvable model of hard core bosons (pseudospins) subject to disorder and frustrating interactions. This solvable model provides insight into the possibility of coexistence of super uidity and glassy density order, as well as into the nature of the coexistence phase (the superglass). In particular, for the considered mean field model I prove the existence of a superglass phase. This complements the numerical evidence for such phases provided by quantum Monte Carlo investigations in finite dimensions and on random graphs. Those were, however, limited to finite temperature, and could thus not fully elucidate the structure of the phases at T = 0. In contrast, my analytical approach allows one to understand the quantum phase transition between glassy superfluid and insulator, and the non-trivial role played by glassy correlations. When the frustrated interactions are strong enough, the superfluid order may be destroyed. As I will show in a mean field model, this happens within the glass phase of the system, where a disorder induced superfluid-insulator phase transition takes place to give way to a frustrated Bose glass. The glassy background on top of which this happens leads to many interesting phenomena which seem not to have been noticed before. To understand the nature of the glassy superfluid-insulator quantum phase transition at zero temperature and the transport properties on the insulating, Bose glass side of the transition is the goal of the third part of my thesis. To address the above questions, I studied an exactly solvable model of a glassy superfluid-insulator quantum phase transition on a Bethe lattice geometry with high connectivity. My main results can be summarized as follows: i) I found that the superfluid-insulator transition is shifted to stronger hopping. This is a result of the pseudo gap in the density of states of the glass state, which tends to strongly disfavor the onset of superfluidity. ii) In the glassy insulator, the discrete local energy levels become broadened due to the quantum fluctuations.The level-broadening process appears as a phase transition which has strong similarities with an Anderson localization transition, and has implications on many body localization. By using the locator expansion for bosons I found that, the glassy insulator has a finite mobility edge for the bosonic excitations, which, however, does not close upon approaching the SI quantum phase transition point. This finding helps to understand the nature of the superfluid-to-frustrated Bose glass transition: the superfluid emerges as a collective phase ordering phenomenon at zero temperature, and not as a condensation in to a single particle delocalized state, in contrast to opposite predictions in the recent literatures. The existence of a mobility edge in the insulator suggests the possibility of phononless, activated transport in the bosonic insulator, which might be a candidate explanation for the experimentally seen activated transport, which has remained a mystery for a long time.
Understanding the effects of disorder in quantum many-body systems is one of the central challenges of contemporary statistical physics. Its interplay with interactions and quantum fluctuations profoundly modifies the nature of quantum states and their dynamics, giving rise to phenomena such as localization, ergodicity breaking, and the emergence of complex collective behavior. This thesis presents new results on two classes of disordered quantum systems. The first part is devoted to the study of the many-body localization transition, a paradigmatic example of disorder-induced ergodicity breaking. The critical properties of the transition are investigated using a numerical renormalization-group approach, establishing a connection between its scaling behavior and that of the Anderson localization transition on tree-like graphs. The second part focuses on the zero-temperature quantum spin-glass phase of the two-dimensional bond-disordered Heisenberg model. Combining large-scale numerical simulations with a semiclassical expansion, evidence is provided for the existence of the spin-glass phase. The effective model governing the semiclassical fluctuations is further analyzed to characterize the low-energy excitation spectrum and its localization properties. Finally, the weak-disorder regime is investigated and its implications for the ferromagnetic phase are discussed. Overall, the results presented in this thesis contribute to the understanding of the effects of disorder in interacting quantum many-body systems, encompassing both localization phenomena and the collective behavior of disordered quantum magnets.
We theoretically investigate the many-body localization phase transition in a one-dimensional Ising spin chain with random long-range spin-spin interactions, ${V}_{ij}\ensuremath{\propto}{\left|i\ensuremath{-}j\right|}^{\ensuremath{-}\ensuremath{\alpha}}$, where the exponent of the interaction range $\ensuremath{\alpha}$ can be tuned from zero to infinitely large. By using exact diagonalization, we calculate the half-chain entanglement entropy and the energy spectral statistics and use them to characterize the phase transition towards the many-body localization phase at infinite temperature and at sufficiently large disorder strength. We perform finite-size scaling to extract the critical disorder strength and the critical exponent of the divergent localization length. With increasing $\ensuremath{\alpha}$, the critical exponent experiences a sharp increase at about ${\ensuremath{\alpha}}_{c}\ensuremath{\simeq}1.2$ and then gradually decreases to a value found earlier in a disordered short-ranged interacting spin chain. For $\ensuremath{\alpha}<{\ensuremath{\alpha}}_{c}$, we find that the system is mostly localized and the increase in the disorder strength may drive a transition between two many-body localized phases. In contrast, for $\ensuremath{\alpha}>{\ensuremath{\alpha}}_{c}$, the transition is from a thermalized phase to the many-body localization phase. Our predictions could be experimentally tested with an ion-trap quantum emulator with programmable random long-range interactions, or with randomly distributed Rydberg atoms or polar molecules in lattices.
We study quantum phase transitions of Mott insulators between spin-glass and featureless paramagnetic phases. Specifically, we consider a triangular-lattice bilayer Heisenberg model with bond disorder. The clean system has two phases, a dimer quantum paramagnet and a non-collinear antiferromagnet, separated by a quantum critical point. Bond disorder destroys the antiferromagnetic phase via the interference of dipolar textures and results in spin-glass order, such that the system features a quantum phase transition between spin-glass and dimer phases. We study the vicinity of this transition using a variant of bond-operator theory and calculate thermodynamic observables as well as excitation spectra. Bond disorder leads to strong inhomogeneities near the transition which suppresses non-collinearities and leads to anomalously weak glassiness in the near-critical quantum spin glass. We characterize the low-energy excitations which have a strong tendency towards spatial localization, and we track the behavior of the amplitude (i.e. Higgs) mode across the glass phase whose spectral weight we find to be strongly suppressed due to interference effects.
Spin glass behavior was first seen in metallic alloys with magnetic impurities dispersed randomly in the main non-magnetic component. Due to the simultaneous presence of quenched disorder and competing ferromagnetic and antiferromagnetic interactions, the nonequilibrium behaviors of spin glasses are intricate and difficult to understand. Experimental research on spin glass materials indicate extremely slow dynamics, suggesting that the direct numerical simulation using local spin flips will take impractical time scales to replicate the nonequilibrium behaviors. This thesis introduces efficient algorithms and heuristics for the numerical simulation of simple models for spin glasses and discusses significant simulation results. The results provide insights into understanding the nonequilibrium behaviors, especially aging and memory, through studying microscopic constituent configurations. The methodology and conclusion presented in this thesis can be potentially extended to address general questions about disordered materials at nonequilibrium. Chapter 1 provides a thorough overview of spin glasses. A detailed introduction on complex nonequilibrium behaviors, including slow dynamics and memory effects, seen in typical experiments of spin glass materials, is presented. A general discussion about the characterization and categorization of memory effects found in many out-of-equilibrium systems is highlighted. Two prevalent theoretical models on spin glasses, the droplet and replica symmetry breaking theory, are described. Lastly, algorithm development and simulation techniques used to study spin glasses are covered. In Chapter 2, a patchwork heuristic is applied to accelerate the relaxation of spin glass models in dimensions d=1 and d=2 with the goal of investigating the dynamics of loss and recovery of spin configurations. A well-designed simulation protocol, consists of a aging and recovery stage, is used to mimic the temperature cycling experiment, where the temperature variations are replaced by the variations of the spin coupling realization. Eventual recovery of memory upon coarsening is seen in two-dimensional Ising spin glasses and one-dimensional clock models, while one-dimensional Ising spin systems do not exhibit configuration memory. A central result is that the coarsening recovery length scale sc grows faster than the aging length scale ℓ as a power law relation sc∼ℓ2.58. In Chapter 3, the microscopic dependence of spin glass equilibration on sample history, i.e., the initial configuration and the specific noise history of stochastic dynamics, is examined. The patchwork dynamics is used to simulate the glassy dynamics over a wide range of length scales. The dependence of the nearest-neighbor spin correlations at long time on initial spin configurations and on noise histories are compared. Most of the local spin correlations are independent of the details of evolution, due to the formation of rigid domains. The correlations on a fractal set of domain walls are found to have distinct dependence on the noise history and the initial state. The initial ``nature'' of the state is retained in the domain walls under coarsening more than the ``nurture'' of the noise, while the location of the domain walls at each scale is determined by the frozen disorder (``environment''). We also provide evidence that coarsening with local dynamics at finite temperature gives persistence of overlaps roughly consistent with coarsening at zero temperature, at least over short distances, supporting the use of zero-temperature heuristics to address finite-temperature dynamics and memory.
Understanding how closed quantum systems dynamically approach thermal equilibrium presents a major unresolved problem in statistical physics. Generically, non-integrable quantum systems are expected to thermalize as they comply with the Eigenstate Thermalization Hypothesis. However, in the presence of strong disorder, the dynamics can possibly slow down to a degree that systems fail to thermalize on experimentally accessible timescales, as in spin glasses or many-body localized systems. In general, particularly in long-range interacting quantum systems, the specific nature of the disorder necessary for the emergence of a prethermal, metastable state--distinctly separating the timescales of initial relaxation and subsequent slow thermalization--remains an open question. We study an ensemble of Heisenberg spins with a tunable distribution of random coupling strengths realized by a Rydberg quantum simulator. We observe a drastic change in the late-time magnetization when increasing disorder strength. The data is well described by models based on pairs of strongly interacting spins, which are treated as thermal for weak disorder and isolated for strong disorder. Our results indicate a crossover into a pair-localized prethermal regime in a closed quantum system of thousands of spins in the critical case where the exponent of the power law interaction matches the spatial dimension.