Superfluid-insulator transition in commensurate disordered bosonic systems: large-scale worm algorithm simulations

Collective Modes at a Disordered Quantum Phase Transition

We study the collective excitations, i.e., the Goldstone (phase) mode and the Higgs (amplitude) mode, near the superfluid-Mott glass quantum phase transition in a two-dimensional system of disordered bosons. Using Monte Carlo simulations as well as an inhomogeneous quantum mean-field theory with Gaussian fluctuations, we show that the Higgs mode is strongly localized for all energies, leading to a noncritical scalar response. In contrast, the lowest-energy Goldstone mode undergoes a striking delocalization transition as the system enters the superfluid phase. We discuss the generality of these findings and experimental consequences, and we point out potential relations to many-body localization.

Quantum Criticality at the Superconductor-Insulator Transition Probed by the Nernst Effect

The superconductor-insulator transition (SIT) is an excellent example of a quantum phase transition at zero temperature, dominated by quantum fluctuations. These are expected to be very prominent close to the quantum critical point. So far, most of the experimental studies of the SIT have concentrated on transport properties and tunneling experiments that provide indirect information on criticality close to the transition. Here we present an experiment uniquely designed to study the evolution of quantum fluctuations through the quantum critical point. We utilize the Nernst effect, which has been shown to be effective in probing superconducting fluctuation. We measure the Nernst coefficient in amorphous indium oxide films tuned through the SIT and find a large signal on both the superconducting and the insulating sides, which peaks close to the critical point. The transverse Peltier coefficient α_{xy}, which is the thermodynamic quantity extracted from these measurements, follows quantum critical scaling with critical exponents ν∼0.7 and z∼1. These exponents are consistent with a clean X-Y model in 2+1 dimensions.

Quantum Glass of Interacting Bosons with Off-Diagonal Disorder

We study disordered interacting bosons described by the Bose-Hubbard model with Gaussian-distributed random tunneling amplitudes. It is shown that the off-diagonal disorder induces a spin-glass-like ground state, characterized by randomly frozen quantum-mechanical U(1) phases of bosons. To access criticality, we employ the "n-replica trick," as in the spin-glass theory, and the Trotter-Suzuki method for decomposition of the statistical density operator, along with numerical calculations. The interplay between disorder, quantum, and thermal fluctuations leads to phase diagrams exhibiting a glassy state of bosons, which are studied as a function of model parameters. The considered system may be relevant for quantum simulators of optical-lattice bosons, where the randomness can be introduced in a controlled way. The latter is supported by a proposition of experimental realization of the system in question.

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.

Anomalous quantum glass of bosons in a random potential in two dimensions

We present a quantum Monte Carlo study of the "quantum glass" phase of the two-dimensional Bose-Hubbard model with random potentials at filling ρ=1. In the narrow region between the Mott and superfluid phases, the compressibility has the form κ∼exp(-b/T^{α})+c with α<1 and c vanishing or very small. Thus, at T=0 the system is either incompressible (a Mott glass) or nearly incompressible (a Mott-glass-like anomalous Bose glass). At stronger disorder, where a glass reappears from the superfluid, we find a conventional highly compressible Bose glass. On a path connecting these states, away from the superfluid at larger Hubbard repulsion, a change of the disorder strength by only 10% changes the low-temperature compressibility by more than 4 orders of magnitude, lending support to two types of glass states separated by a phase transition or a sharp crossover.

Microscopic molecular superfluid response: theory and simulations

Since its discovery in 1938, superfluidity has been the subject of much investigation because it provides a unique example of a macroscopic manifestation of quantum mechanics. About 60 years later, scientists successfully observed this phenomenon in the microscopic world though the spectroscopic Andronikashvili experiment in helium nano-droplets. This reduction of scale suggests that not only helium but also para-H2 (pH2) can be a candidate for superfluidity. This expectation is based on the fact that the smaller number of neighbours and surface effects of a finite-size cluster may hinder solidification and promote a liquid-like phase. The first prediction of superfluidity in pH2 clusters was reported in 1991 based on quantum Monte Carlo simulations. The possible superfluidity of pH2 was later indirectly observed in a spectroscopic Andronikashvili experiment in 2000. Since then, a growing number of studies have appeared, and theoretical simulations have been playing a special role because they help guide and interpret experiments. In this review, we go over the theoretical studies of pH2 superfluid clusters since the experiment of 2000. We provide a historical perspective and introduce the basic theoretical formalism along with key experimental advances. We then present illustrative results of the theoretical studies and comment on the possible future developments in the field. We include sufficient theoretical details such that the review can serve as a guide for newcomers to the field.

Bose-glass transition and spin-wave localization for 2D bosons in a random potential

A spin-wave approach of the zero temperature superfluid-insulator transition for two-dimensional hard-core bosons in a random potential μ=±W is developed. While at the classical level there is no intervening phase between the Bose-condensed superfluid (SF) and the gapped disordered insulator, the introduction of quantum fluctuations leads to a much richer physics. Upon increasing the disorder strength W, the Bose-condensed fraction disappears first, before the SF. Then a gapless Bose-glass phase emerges over a finite region until the insulator appears. Furthermore, in the strongly disordered SF regime, a mobility edge in the spin-wave excitation spectrum is found at a finite frequency Ω(c) decreasing with W, and presumably vanishing in the Bose-glass phase.

Universal superfluid transition and transport properties of two-dimensional dirty bosons

We study the phase diagram of two-dimensional, interacting bosons in the presence of a correlated disorder in continuous space, by using large-scale quantum Monte Carlo simulations at finite temperature. We show that the superfluid transition is strongly protected against disorder. It remains of the Berezinskii-Kosterlitz-Thouless type up to disorder strengths comparable to the chemical potential. Moreover, we study the transport properties in the strong disorder regime where a zero-temperature Bose-glass phase is expected. We show that the conductance exhibits a thermally activated behavior vanishing only at zero temperature. Our results point towards the existence of a Bose bad-metal phase as a precursor of the Bose-glass phase.

Phase diagram of the commensurate two-dimensional disordered Bose-Hubbard model

We establish the full ground state phase diagram of the disordered Bose-Hubbard model in two dimensions at a unity filling factor via quantum Monte Carlo simulations. Similarly to the three-dimensional case we observe extended superfluid regions persisting up to extremely large values of disorder and interaction strength which, however, have small superfluid fractions and thus low transition temperatures. In the vicinity of the superfluid-insulator transition of the pure system, we observe an unexpectedly weak--almost not resolvable--sensitivity of the critical interaction to the strength of (weak) disorder.

Absence of a direct superfluid to mott insulator transition in disordered bose systems

We prove the absence of a direct quantum phase transition between a superfluid and a Mott insulator in a bosonic system with generic, bounded disorder. We also prove the compressibility of the system on the superfluid-insulator critical line and in its neighborhood. These conclusions follow from a general theorem of inclusions, which states that for any transition in a disordered system, one can always find rare regions of the competing phase on either side of the transition line. Quantum Monte Carlo simulations for the disordered Bose-Hubbard model show an even stronger result, important for the nature of the Mott insulator to Bose glass phase transition: the critical disorder bound Delta(c) corresponding to the onset of disorder-induced superfluidity, satisfies the relation Delta(c)>Eg/2, with Eg/2 the half-width of the Mott gap in the pure system.

Superfluid transition in a bose gas with correlated disorder

The superfluid transition of a three-dimensional gas of hard-sphere bosons in a disordered medium is studied using quantum Monte Carlo methods. Simulations are performed in continuous space both in the canonical and in the grand-canonical ensemble. At fixed density we calculate the shift of the transition temperature as a function of the disorder strength, while at fixed temperature we determine both the critical chemical potential and the critical density separating normal and superfluid phases. In the regime of strong disorder the normal phase extends up to large values of the degeneracy parameter, and the critical chemical potential exhibits a linear dependence in the intensity of the random potential. The role of interactions and disorder correlations is also discussed.

Strongly interacting bosons in a disordered optical lattice

We experimentally probe the properties of the disordered Bose-Hubbard model using an atomic Bose-Einstein condensate trapped in a 3D disordered optical lattice. Controllable disorder is introduced using a fine-grained optical speckle field with features comparable in size to the lattice spacing along every lattice direction. A precision measurement of the disordering potential is used to compute the single-particle parameters of the system. To constrain theories of the disordered Bose Hubbard model, we have measured the change in condensate fraction as a function of disorder strength for several different ratios of tunneling to interaction energy. We observe disorder-induced, reversible suppression of condensate fraction for superfluid and coexisting superfluid-Mott-insulator phases.

Quantum glass phases in the disordered Bose-Hubbard model

The phase diagram of the Bose-Hubbard model in the presence of off-diagonal disorder is determined using quantum Monte Carlo simulations. A sequence of quantum glass phases intervene at the interface between the Mott insulating and the superfluid phases of the clean system. In addition to the standard Bose glass phase, the coexistence of gapless and gapped regions close to the Mott insulating phase leads to a novel Mott glass regime which is incompressible yet gapless. Numerical evidence for the properties of these phases is given in terms of global (compressibility, superfluid stiffness) and local (compressibility, momentum distribution) observables.

Mott glass in site-diluted S=1 antiferromagnets with single-ion anisotropy

The interplay between site dilution and quantum fluctuations in S=1 Heisenberg antiferromagnets on the square lattice is investigated using quantum Monte Carlo simulations. Quantum fluctuations are tuned by a single-ion anisotropy D. In the clean limit, a sufficiently large D>Dc=5.65(2)J forces each spin into its mS=0 state, and thus destabilizes antiferromagnetic order. In the presence of site dilution, quantum fluctuations are found to destroy Néel order before the percolation threshold of the lattice is reached, if D exceeds a critical value D*=2.3(2)J. This mechanism opens up an extended quantum-disordered Mott-glass phase on the percolated lattice, characterized by a gapless spectrum and vanishing uniform susceptibility.

Critical dynamics of the dirty boson problem: revisiting the equality z=d

It is shown that previous arguments, leading to the equality z=d for the dynamical exponent describing the Bose glass to superfluid transition in d dimensions, may break down, as apparently seen in recent simulations. The key observation is that the major contribution to the compressibility, which remains finite through the transition and was predicted to scale as kappa approximately |delta|((d-z)nu) (where delta is the deviation from criticality and nu is the correlation length exponent) comes from the analytic, not the singular part of the free energy, and is not restricted by any conventional scaling hypothesis.

Quantum phase transitions of hard-core bosons in background potentials

We study the zero temperature phase diagram of hard-core bosons in two dimensions subjected to three types of background potentials: staggered, uniform, and random. In all three cases there is a quantum phase transition from a superfluid (at small potential) to a normal phase (at large potential), but with different universality classes. As expected, the staggered case belongs to the XY universality, while the uniform potential induces a mean field transition. The disorder driven transition is clearly different from both; in particular, we find z approximately 1.4, nu approximately 1, and beta approximately 0.6.

Worm algorithm and diagrammatic Monte Carlo: a new approach to continuous-space path integral Monte Carlo simulations

A detailed description is provided of a new worm algorithm, enabling the accurate computation of thermodynamic properties of quantum many-body systems in continuous space, at finite temperature. The algorithm is formulated within the general path integral Monte Carlo (PIMC) scheme, but also allows one to perform quantum simulations in the grand canonical ensemble, as well as to compute off-diagonal imaginary-time correlation functions, such as the Matsubara Green function, simultaneously with diagonal observables. Another important innovation consists of the expansion of the attractive part of the pairwise potential energy into elementary (diagrammatic) contributions, which are then statistically sampled. This affords a complete microscopic account of the long-range part of the potential energy, while keeping the computational complexity of all updates independent of the size of the simulated system. The computational scheme allows for efficient calculations of the superfluid fraction and off-diagonal correlations in space-time, for system sizes which are orders of magnitude larger than those accessible to conventional PIMC. We present illustrative results for the superfluid transition in bulk liquid 4He in two and three dimensions, as well as the calculation of the chemical potential of hcp 4He.

Worm algorithm for continuous-space path integral monte carlo simulations

We present a new approach to path integral Monte Carlo (PIMC) simulations based on the worm algorithm, originally developed for lattice models and extended here to continuous-space many-body systems. The scheme allows for efficient computation of thermodynamic properties, including winding numbers and off-diagonal correlations, for systems of much greater size than that accessible to conventional PIMC simulations. As an illustrative application of the method, we simulate the superfluid transition of 4He in two dimensions.

Ultracold atoms in optical lattices with random on-site interactions

We consider the physics of lattice bosons affected by disordered on-site interparticle interactions. Characteristic qualitative changes in the zero-temperature phase diagram are observed when compared to the case of randomness in the chemical potential. The Mott-insulating regions shrink and eventually vanish for any finite disorder strength beyond a sufficiently large filling factor. Furthermore, at low values of the chemical potential both the superfluid and Mott insulator are stable towards formation of a Bose glass leading to a possibly nontrivial tricritical point. We discuss feasible experimental realizations of our scenario in the context of ultracold atoms on optical lattices.

Superfluid-insulator transition in a commensurate one-dimensional bosonic system with off-diagonal disorder

We study the nature of the superfluid-insulator quantum phase transition in a one-dimensional system of lattice bosons with off-diagonal disorder in the limit of a large integer filling factor. Monte Carlo simulations of two strongly disordered models show that the universality class of the transition in question is the same as that of the superfluid-Mott-insulator transition in a pure system. This result can be explained by disorder self-averaging in the superfluid phase and the applicability of the standard quantum hydrodynamic action. We also formulate the necessary conditions which should be satisfied by the stong-randomness universality class, if one exists.

Superfluid interfaces in quantum solids

One scenario for the nonclassical moment of inertia of solid 4He discovered by Kim and Chan [Nature (London) 427, 225 (2004)] is the superfluidity of microcrystallite interfaces. On the basis of the most simple model of a quantum crystal--the checkerboard lattice solid--we show that the superfluidity of interfaces between solid domains can exist in a wide range of parameters. At strong enough interparticle interaction, a superfluid interface becomes an insulator via a quantum phase transition. Under the conditions of particle-hole symmetry, the transition is of the standard U(1) universality class in 3D, while in 2D the onset of superfluidity is accompanied by the interface roughening, driven by fractionally charged topological excitations.

Exact results for one-dimensional disordered bosons with strong repulsion

We study one-dimensional disordered bosons with strong repulsive interactions. A Bose-Fermi mapping expresses this problem in terms of noninteracting Anderson-localized fermions, whereby known results for the distribution function of the local density of states, the spectral statistics, and density-density correlations can be transferred to this new domain of applicability. We show that disorder destroys bosonic quasi-long-range order by calculating the momentum distribution, and comment on the experimental observability of these predictions in ultracold atomic gases.

Weak first-order superfluid-solid quantum phase transitions

We study superfluid-solid zero-temperature transitions in two-dimensional lattice boson-spin models using worm-algorithm Monte Carlo simulations. We observe that such transitions are typically first order with the exception of special high-symmetry points which require fine-tuning in the Hamiltonian parameter space. We present evidence that the superfluid-checkerboard solid and superfluid-valence-bond solid transitions at half-integer filling factor are extremely weak first-order transitions and in small systems can be confused with continuous or high-symmetry points.

Superfluid-superfluid phase transitions in a two-component Bose-Einstein condensate

Depending on the Hamiltonian parameters, two-component bosons in an optical lattice can form at least three different superfluid phases in which both components participate in the superflow: a (strongly interacting) mixture of two miscible superfluids (2SF), a paired superfluid (PSF) vacuum, and (at a commensurate total filling factor) the super-counter-fluid (SCF) state. We study the universal properties of the 2SF-PSF and 2SF-SCF quantum phase transitions and show that (i) they can be mapped onto each other and (ii) their universality class is identical to the (d+1)-dimensional normal-superfluid transition in a single-component liquid. The finite-temperature 2SF-PSF(SCF) transitions and the topological properties of 2SF-PSF(SCF) interfaces are also discussed.