The question of how the low-energy properties of disordered quantum systems may be connected to exotic localization phenomena at high energy is a key open question in the context of quantum glasses and many-body localization. In the preceding Letter [L. Villa, S. J. Thomson, and L. Sanchez-Palencia, preceding Letter, Phys. Rev. A 104, L021301 (2021)] we have shown that key features of the excitation spectrum of a disordered system can be efficiently probed from out-of-equilibrium dynamics following a quantum quench, providing distinctive signatures of the various phases. Here we extend this work by providing a more-in-depth study of the behavior of the quench spectral functions associated with different observables and investigating an extended parameter regime. We provide a detailed introduction to quench spectroscopy for disordered systems and show how spectral properties can be probed using both local operators and two-point correlation functions. We benchmark the technique using the one-dimensional Bose-Hubbard model in the presence of a random external potential, focusing on the low-lying excitations, and demonstrate that quench spectroscopy can distinguish the Mott insulator, superfluid, and Bose glass phases. We then explicitly reconstruct the zero-temperature phase diagram of the disordered Bose-Hubbard at fixed filling using two independent methods, experimentally accessible via both time-of-flight imaging and quantum gas microscopy, respectively, and demonstrate that quench spectroscopy can give valuable insights into the distribution of rare regions within disordered systems.
The question of how the low-energy properties of disordered quantum systems may be connected to exotic localization phenomena at high energy is a key open question in the context of quantum glasses and many-body localization. In the preceding Letter [L. Villa, S. J. Thomson, and L. Sanchez-Palencia, preceding Letter, Phys. Rev. A 104, L021301 (2021)] we have shown that key features of the excitation spectrum of a disordered system can be efficiently probed from out-of-equilibrium dynamics following a quantum quench, providing distinctive signatures of the various phases. Here we extend this work by providing a more-in-depth study of the behavior of the quench spectral functions associated with different observables and investigating an extended parameter regime. We provide a detailed introduction to quench spectroscopy for disordered systems and show how spectral properties can be probed using both local operators and two-point correlation functions. We benchmark the technique using the one-dimensional Bose-Hubbard model in the presence of a random external potential, focusing on the low-lying excitations, and demonstrate that quench spectroscopy can distinguish the Mott insulator, superfluid, and Bose glass phases. We then explicitly reconstruct the zero-temperature phase diagram of the disordered Bose-Hubbard at fixed filling using two independent methods, experimentally accessible via both time-of-flight imaging and quantum gas microscopy, respectively, and demonstrate that quench spectroscopy can give valuable insights into the distribution of rare regions within disordered systems.
The characterization of excitations in disordered quantum systems is a central issue in connection with glass physics and many-body localization. Here we show that quench spectroscopy of a disordered model, as realized from its out-of-equilibrium dynamics following a global quench, allows us to fully characterize the spectral properties of the disordered phases. In the Bose-Hubbard model, a clear signature of gapless excitations in momentum-resolved spectroscopy enables us to accurately locate the Mott insulator to Bose glass transition, while the presence or absence of a well-defined soundlike mode distinguishes the superfluid from the Bose glass phase. Moreover, spatially resolved spectroscopy provides local spectral properties and allows us to extract the typical spacing of gapless regions, giving a second independent way to uniquely identify all three phases. Our findings have far-ranging implications for a variety of experimental platforms and offer a powerful and versatile probe of the low-energy phases of disordered systems.
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.
We study the equilibrium and nonequilibrium properties of strongly interacting bosons on a lattice in the presence of a random bounded disorder potential. Using a Gutzwiller projected variational technique, we study the equilibrium phase diagram of the disordered Bose-Hubbard model and obtain the Mott insulator, Bose glass, and superfluid phases. We also study the nonequilibrium response of the system under a periodic temporal drive where, starting from the superfluid phase, the hopping parameter is ramped down linearly in time, and back to its initial value. We study the density of excitations created, the change in the superfluid order parameter, and the energy pumped into the system in this process as a function of the inverse ramp rate $\ensuremath{\tau}$. For the clean case the density of excitations goes to a constant, while the order parameter and energy relax as $1/\ensuremath{\tau}$ and $1/{\ensuremath{\tau}}^{2}$ respectively. With disorder, the excitation density decays exponentially with $\ensuremath{\tau}$, with the decay rate increasing with the disorder, to an asymptotic value independent of the disorder. The energy and change in order parameter also decrease as $\ensuremath{\tau}$ is increased.
We establish the phase diagram of the disordered three-dimensional Bose-Hubbard model at unity filling which has been controversial for many years. The theorem of inclusions, proven by Pollet et al. [Phys. Rev. Lett. 103, 140402 (2009)] states that the Bose-glass phase always intervenes between the Mott insulating and superfluid phases. Here, we note that assumptions on which the theorem is based exclude phase transitions between gapped (Mott insulator) and gapless phases (Bose glass). The apparent paradox is resolved through a unique mechanism: such transitions have to be of the Griffiths type when the vanishing of the gap at the critical point is due to a zero concentration of rare regions where extreme fluctuations of disorder mimic a regular gapless system. An exactly solvable random transverse field Ising model in one dimension is used to illustrate the point. A highly nontrivial overall shape of the phase diagram is revealed with the worm algorithm. The phase diagram features a long superfluid finger at strong disorder and on-site interaction. Moreover, bosonic superfluidity is extremely robust against disorder in a broad range of interaction parameters; it persists in random potentials nearly 50 (!) times larger than the particle half-bandwidth. Finally, we comment on the feasibility of obtaining this phase diagram in cold-atom experiments, which work with trapped systems at finite temperature.
We study the dynamic structure factor of a one-dimensional Bose gas confined in an optical lattice and modeled by the Bose–Hubbard Hamiltonian, using a variety of numerical and analytical approaches. The dynamic structure factor, experimentally measurable by Bragg spectroscopy, is studied in three relevant cases: in the clean regime, featuring either a superfluid or a Mott phase; and in the presence of two types of (quasi-)disordered external potentials: a quasi-periodic potential obtained from a bichromatic superlattice and a random box disorder—both featuring a Bose-glass phase. In the clean case, we show the emergence of a gapped doublon mode (corresponding to a repulsively bound state) for incommensurate filling, well separated from the low-energy acoustic mode. In the disordered case, we show that the dynamic structure factor provides direct insight into the spatial structure of the excitations, unveiling their localized nature, which represents a fundamental signature of the Bose-glass phase. Furthermore, it provides a clear fingerprint of the very nature of the localization mechanism which differs for the two kinds of disorder potentials we consider. In special cases, the dynamic structure factor may provide an estimate of the position of the localization transition from superfluid to Bose glass, in a complementary manner to the information deduced from the momentum distribution.
We extend the self-energy functional theory to the case of interacting lattice bosons in the presence of symmetry breaking and quenched disorder. The self-energy functional we derive depends only on the self-energies of the disorder-averaged propagators, allowing for the construction of general non-perturbative approximations. Using a simple single-site reference system with only three variational parameters, we are able to reproduce numerically exact quantum Monte Carlo (QMC) results on local observables of the Bose–Hubbard model with box disorder with high accuracy. At strong interactions, the phase boundaries are reproduced qualitatively but shifted with respect to the ones observed with QMC due to the extremely low condensate fraction in the superfluid phase. Deep in the strongly-disordered weakly-interacting regime, the simple reference system employed is insufficient and no stationary solutions can be found within its restricted variational subspace. By systematically analyzing thermodynamical observables and the spectral function, we find that the strongly interacting Bose glass is characterized by different regimes, depending on which local occupations are activated as a function of the disorder strength. We find that the particles delocalize into isolated superfluid lakes over a strongly localized background around maximally-occupied sites whenever these sites are particularly rare. Our results indicate that the transition from the Bose glass to the superfluid phase around unit filling at strong interactions is driven by the percolation of superfluid lakes which form around doubly occupied sites.