Superfluid-insulator Transition Research Articles

Pairing Superfluid–Insulator Transition Induced by Atom–Molecule Conversion in Bosonic Mixtures in Optical Lattice

Motivated by the recent experiment on bosonic mixtures of atoms and molecules, we investigate pairing superfluid–insulator (SI) transition for bosonic mixtures of atoms and molecules in a one-dimensional optical lattice, which is described by an extended Bose–Hubbard model with atom–molecule conservation (AMC). It is found that AMC can induce an extra pair–superfluid phase though the system does not demonstrate pair-hopping. In particular, the system may undergo several pairing SI or insulator–superfluid transitions as the detuning from the Feshbach resonance is varied from negative to positive, and the larger positive detuning can bifurcate the pair–superfluid phases into mixed superfluid phases consisting of single-atomic and pair-atomic superfluid. The calculation of the second-order Rényi entropy reveals that the discontinuity in its first-order derivative corresponds to the phase boundary of the pairing SI transition. This means that the residual entanglement in our mean-field treatment can be used to efficiently capture the signature of the pairing SI transition induced by AMC.

CT or P problem and symmetric gapped fermion solution

An analogous ``strong $CP$ problem'' is identified in a toy model in two-dimensional spacetime: a general $1+1\mathrm{d}$ Abelian U(1) anomaly free chiral fermion and chiral gauge theory with a generic theta instanton term $\frac{\ensuremath{\theta}}{2\ensuremath{\pi}}\ensuremath{\int}F$. The theta term alone violates the charge-conjugation-time-reversal $CT$ and the parity $P$ discrete symmetries. The analogous puzzle here is the $CT$ or $P$ problem in $1+1\mathrm{d}$: Why can the $\overline{\ensuremath{\theta}}$ angle (including the effect of $\ensuremath{\theta}$ and the complex phase of a mass matrix) be zero or small for a natural reason? We show that this $CT$ or $P$ problem can be solved by a symmetric mass generation mechanism (SMG, namely generating a mass or energy gap while preserving an anomaly free symmetry). This $1+1\mathrm{d}$ toy model mimics several features of the $3+1\mathrm{d}$ Standard Model: chiral matter content, confinement, and Anderson-Higgs mass by Yukawa-Higgs term. One solution replaces some chiral fermion's Higgs-induced mean-field mass with SMG-induced non-mean-field mass. Another solution enriches this toy model by introducing several new physics beyond the Standard Model: a parity-reflection $PR$ discrete symmetry maps between the chiral and mirror fermions as fermion doubling localized on two domain walls at high energy, and SMG dynamically generates mass to the mirror fermion while still preserving the anomaly free chiral symmetry at an intermediate energy scale, much before the Higgs mechanism generates mass to the chiral fermion at lower energy. Without loss of generality, an arguably simplest $1+1\mathrm{d}$ U(1) symmetric anomaly free chiral fermion/gauge theory (e.g., Weyl fermions with ${3}_{L}\text{\ensuremath{-}}{4}_{L}\text{\ensuremath{-}}{5}_{R}\text{\ensuremath{-}}{0}_{R}$ U(1) charges) is demonstrated in detail. As an analogy to the superfluid-insulator or order-disorder quantum phase transition, in contrast to the conventional Peccei-Quinn solution sitting in the (quasi-long-range-order)-superfluid ordered phase, our solution is in the SMG insulator disordered phase.

Quantum Critical Behavior of Entanglement in Lattice Bosons with Cavity-Mediated Long-Range Interactions.

We analyze the ground-state entanglement entropy of the extended Bose-Hubbard model with infinite-range interactions. This model describes the low-energy dynamics of ultracold bosons tightly bound to an optical lattice and dispersively coupled to a cavity mode. The competition between on-site repulsion and global cavity-induced interactions leads to a rich phase diagram, which exhibits superfluid, supersolid, and insulating (Mott and checkerboard) phases. We use a slave-boson treatment of harmonic quantum fluctuations around the mean-field solution and calculate the entanglement entropy across the phase transitions. At commensurate filling, the insulator-superfluid transition is signaled by a singularity in the area-law scaling coefficient of the entanglement entropy, which is similar to the one reported for the standard Bose-Hubbard model. Remarkably, at the continuous Z_{2} superfluid-to-supersolid transition we find a critical logarithmic term, regardless of the filling. This behavior originates from the appearance of a roton mode in the excitation and entanglement spectrum, becoming gapless at the critical point, and it is characteristic of collective models.

Density Matrix Renormalization Group for Continuous Quantum Systems.

We introduce a versatile and practical framework for applying matrix product state techniques to continuous quantum systems. We divide space into multiple segments and generate continuous basis functions for the many-body state in each segment. By combining this mapping with existing numerical density matrix renormalization group routines, we show how one can accurately obtain the ground-state wave function, spatial correlations, and spatial entanglement entropy directly in the continuum. For a prototypical mesoscopic system of strongly interacting bosons we demonstrate faster convergence than standard grid-based discretization. We illustrate the power of our approach by studying a superfluid-insulator transition in an external potential. We outline how one can directly apply or generalize this technique to a wide variety of experimentally relevant problems across condensed matter physics and quantum field theory.

Criticality of the excess energy cost due to the unit-flux-quantum external field for the (2 + 1)D superfluid–insulator transition

The two-dimensional (2D) spin-S = 1 XY model was investigated numerically as a realization of the (2 + 1)D superfluid–Mott-insulator (SF–MI) transition. The interaction parameters are extended so as to suppress corrections to finite-size scaling. Thereby, the external field of a unit flux quantum (Φ = 2π) is applied to the 2D cluster by incorporating the phase factor (ϕ ij : gauge angle between the i and j sites) into the hopping amplitudes. Taking the advantage in that the exact-diagonalization method allows us to treat such a complex-valued matrix element, we evaluated the excess energy cost ΔE(2π) due to the magnetic flux Φ = 2π explicitly in the SF (XY) phase. As a result, we found that the amplitude ratio ρ s/ΔE(2π) (ρ s: spin stiffness) makes sense in proximity to the critical point, exhibiting a notable plateau in the SF-phase side. The plateau height is estimated, and compared to the related studies.

Scaling of the disorder operator at (2+1)d U(1) quantum criticality

We study disorder operator, defined as a symmetry transformation applied to a finite region, across a continuous quantum phase transition in $(2+1)d$. We show analytically that, at a conformally invariant critical point with U(1) symmetry, the disorder operator with a small U(1) rotation angle defined on a rectangle region exhibits power-law scaling with the perimeter of the rectangle. The exponent is proportional to the current central charge of the critical theory. Such a universal scaling behavior is due to the sharp corners of the region and we further obtain a general formula for the exponent when the corner is nearly smooth. To probe the full parameter regime, we carry out systematic computation of the U(1) disorder parameter in the square lattice Bose-Hubbard model across the superfluid-insulator transition with large-scale quantum Monte Carlo simulations, and confirm the presence of the universal corner correction. The exponent of the corner term determined from numerical simulations agrees well with the analytical predictions.

Mixture of scalar bosons and two-color fermions in one dimension: Superfluid-insulator transitions

Superfluid-insulator transitions in a one-dimensional mixture of two-color fermions and scalar bosons are studied within the framework of the Bose-Fermi-Hubbard model. Zero-temperature phase diagrams are constructed for repulsive intraspecies interactions and attractive or repulsive interspecies couplings. In addition to the trivial Mott insulator phases, we report the emergence of new non-trivial insulator phases that depend on the sign of the boson-fermion interaction. These non-trivial insulator phases satisfy the conditions $\rho_B\pm\rho_F=n$ and $\rho_B\pm \tfrac{1}{2}\rho_F=n$, with the plus (minus) sign for repulsive (attractive) interactions and $n$ an integer. Far from fermionic half-filling, the boson-fermion interaction drives a gapless-gapped transition in the spin sector. Our findings could be observed experimentally in state-of-the-art cold-atom setups.

Interplay of interactions and disorder at the superfluid-insulator transition: A dirty two-dimensional quantum critical point

We study the stability of the Wilson-Fisher fixed point of the quantum $\mathrm{O}(2N)$ vector model to quenched disorder in the large-$N$ limit. While a random mass is strongly relevant at the Gaussian fixed point, its effect is screened by the strong interactions of the Wilson-Fisher fixed point. This enables a perturbative renormalization group study of the interplay of disorder and interactions about this fixed point. We show that, in contrast to the spiralling flows obtained in earlier double-$\epsilon$ expansions, the theory flows directly to a quantum critical point characterized by finite disorder and interactions. The critical exponents we obtain for this transition are in remarkable agreement with numerical studies of the superfluid-Mott glass transition. We additionally discuss the stability of this fixed point to scalar and vector potential disorder and use proposed boson-fermion dualities to make conjectures regarding the effects of weak disorder on dual Abelian Higgs and Chern-Simons-Dirac fermion theories when $N=1$.

Superfluid weight and polarization amplitude in the one-dimensional bosonic Hubbard model

We calculate the superfluid weight and the polarization amplitude for the one-dimensional bosonic Hubbard model focusing on the strong-coupling regime. Other than analytic calculations we apply two methods: variational Monte Carlo based on the Baeriswyl wave function and exact diagonalization. The former gives zero superfluid response at integer filling, while the latter gives a superfluid response at finite hopping. From the polarization amplitude we derive the variance and the associated size scaling exponent. Again, the variational study does not produce a finite superfluid weight at integer filling (size scaling exponent is near one), but the Fourier transform of the polarization amplitude behaves in a similar way to the result of exact diagonalization, with a peak at small hopping, and suddenly decreasing at the insulator-superfluid transition. On the other hand, exact diagonalization studies result in a finite spread of the total position which increases with the size of the system. In the superfluid phase the size scaling exponent is two as expected. Importantly, our work addresses the ambiguities that arise in the calculation of the superfluid weight in variational calculations, and we comment on the prediction of Anderson about the superfluid response of the model at integer filling.

Quantum hydrodynamics of vorticity

We formulate a quantum theory of vorticity (hydro)dynamics on a general two-dimensional bosonic lattice. In the classical limit of a bosonic condensate, it reduces to conserved plasma-like vortex-antivortex dynamics. The nonlocal topological character of the vorticity flows is reflected in the bulk-edge correspondence dictated by the Stokes theorem. This is exploited to establish physical boundary conditions that realize, in the coarse-grained thermodynamic limit, an effective chemical-potential bias of vorticity. A Kubo formula is derived for the vorticity conductivity|which could be measured in a suggested practical device|in terms of quantum vorticity-flux correlators of the original lattice model. As an illustrative example, we discuss the superfluidity of vorticity, exploiting the particle-vortex duality at a bosonic superfluid-insulator transition.

Superfluid-Insulator Transition unambiguously detected by entanglement in one-dimensional disordered superfluids

We use entanglement to track the superfluid-insulator transition (SIT) in disordered fermionic superfluids described by the one-dimensional Hubbard model. Entanglement is found to have remarkable signatures of the SIT driven by i) the disorder strength V, ii) the concentration of impurities C and iii) the particle density n. Our results reveal the absence of a critical potential intensity on the SIT driven by V, i.e. any small V suffices to decrease considerably the degree of entanglement: it drops ∼50% for V = −0.25t. We also find that entanglement is non-monotonic with the concentration C, approaching to zero for a certain critical value CC. This critical concentration is found to be related to a special type of localization, here named as fully-localized state, which can be also reached for a particular density nC. Our results show that the SIT driven by n or C has distinct nature whether it leads to the full localization or to the ordinary one: it is a first-order quantum phase transition only when leading to full localization. In contrast, the SIT driven by V is never a first-order quantum phase transition independently on the type of localization reached.

Impact of speckle disorder on a superfluid Fermi system

Optical lattice experiments which probe the effect of disorder on superfluidity often use a speckle pattern for generating the disorder. Such speckle disorder is spatially correlated. While fermionic superfluidity in the presence of uncorrelated disorder is well studied, the impact of correlated disorder, particularly on thermal properties of the superfluid, is poorly understood. We provide a detailed study of the impact of speckle disorder, for varying speckle size and disorder magnitude, on the ground state and thermal properties of a Fermi superfluid. We work in the strong coupling regime of BCS-BEC crossover in a two dimensional lattice. For a fixed disorder strength, an increase in speckle size leads to smoothening of the self-consistent background potential, increase in the critical disorder needed for a superfluid-insulator transition, and an increase in superfluid $T_c$ . Along with these hints at decrease in effective disorder, speckle correlations also suppress the superfluid gap and the gap formation temperature - effects normally associated with increasing disorder. We correlate these effects with the effective potential and the single particle localisation effects in the ground state

Strong randomness criticality in the scratched XY model

We study the finite-temperature superfluid transition in a modified two-dimensional (2D) XY model with power-law-distributed ``scratch''-like bond disorder. As its exponent decreases, the disorder grows stronger and the mechanism driving the superfluid transition changes from conventional vortex-pair unbinding to a strong randomness criticality (termed scratched XY criticality) characterized by a nonuniversal jump of the superfluid stiffness. The existence of the scratched XY criticality at finite temperature and its description by an asymptotically exact semi-renormalization group theory, previously developed for the superfluid-insulator transition in one-dimensional disordered quantum systems, is numerically proven by designing a model with minimal finite-size effects. Possible experimental implementations are discussed.

Superfluid-insulator transition and the BEC-BCS crossover in the Rashba moat band

We study the superconducting transition in a two-dimensional electron gas with strong Rashba spin-orbit coupling. We assume low electron density, such that only the majority spin band participates in the transition. We show that the superconducting transition follows either the Bose-Einstein condensation (BEC), or the Bardeen-Cooper-Schrieffer (BCS) scenarios, depending on the position of the chemical potential with respect to the bottom of the majority band, and the strength of the Coulomb repulsion between electrons. Hence, the BEC-BCS crossover in this system can be driven either by the change in the chemical potential, or the distance to a gate.

Kane-Fisher weak link physics in the clean scratched XY model

The nature of the superfluid-insulator transition in 1D has been much debated recently. In particular, to describe the strong disorder regime characterized by weak link proliferation, a scratched-XY model has been proposed [New J. Phys. \textbf{18}, 045018 (2016)], where the transport is dominated by a single anomalously weak link and is governed by Kane-Fisher weak link physics. In this article, we consider the simplest problem to which the scratched-XY model relates: a single weak link in an otherwise \textit{clean} system, with an intensity $J_W$ which decreases algebraically with the size of the system $J_W\sim L^{-\alpha}$. Using a renormalization group approach and a vortex energy argument, we describe the Kane-Fisher physics in this model and show that it leads to a transition from a transparent regime for $K>K_c$ to a perfect cut for $K<K_c$, with an adjustable $K_c=1/(1-\alpha)$ depending on $\alpha$. We check our theoretical predictions with Monte Carlo numerical simulations complemented by finite-size scaling. Our results clarify two important assumptions at the basis of the scratched-XY scenario, the behaviors of the crossover length scale from weak link physics to transparency and of the superfluid stiffness.

Critical properties of the quantum rotor model with phase frustration at incommensurate filling in two dimensions

The critical properties of the zero-temperature superfluid-insulator transition in the fullyfrustrated quantum rotor model at incommensurate filling in two dimensions are studied. We develop a spherical approximation scheme for the model with phase frustration f = 1/2 by choosing a primitive unit cell composed of two sites. The finite-size scaling behavior of various physical quantities in the spherical limit provides the correlation length critical exponent ν = 0.5 and the dynamical critical exponent z = 2, which are consistent with the scenario of the generic superfluidinsulator transition even in the presence of phase frustration.

Ultracold bosons with cavity-mediated long-range interactions: A local mean-field analysis of the phase diagram

Ultracold bosonic atoms in optical lattices self-organize into a variety of structural and quantum phases when placed into a single-mode cavity and pumped by a laser. Cavity optomechanical effects induce an atom density modulation at the cavity-mode wavelength that competes with the optical lattice arrangement. Simultaneously short-range interactions via particle hopping promote superfluid order such that a variety of structural and quantum coherent phases can occur. We analyze the emerging phase diagram in two dimensions by means of an extended Bose-Hubbard model using a local mean-field approach combined with a superfluid cluster analysis. For commensurate ratios of the cavity and external lattice wavelengths, the Mott insulator-superfluid transition is modified by the appearance of charge density wave and supersolid phases, at which the atomic density supports the buildup of a cavity field. For incommensurate ratios, the optomechanical forces induce the formation of Bose-glass and superglass phases, namely, nonsuperfluid and superfluid phases, respectively, displaying quasiperiodic density modulations, which in addition can exhibit structural and superfluid stripe formation. The onset of such structures is constrained by the on-site interaction and is favorable at fractional densities. Experimental observables are identified and discussed.

Cluster mean-field signature of entanglement entropy in bosonic superfluid-insulator transitions

Entanglement entropy (EE), a fundamental conception in quantum information for characterizing entanglement, has been extensively employed to explore quantum phase transitions (QPTs). Although the conventional single-site mean-field (MF) approach successfully predicts the emergence of QPTs, it fails to include any entanglement. Here, for the first time, in the framework of a cluster MF treatment, we extract the signature of EE in the bosonic superfluid-insulator transitions. We consider a trimerized Kagome lattice of interacting bosons, in which each trimer is treated as a cluster, and implement the cluster MF treatment by decoupling all inter-trimer hopping. In addition to superfluid and integer insulator phases, we find that fractional insulator phases appear when the tunneling is dominated by the intra-trimer part. To quantify the residual bipartite entanglement in a cluster, we calculate the second-order Renyi entropy, which can be experimentally measured by quantum interference of many-body twins. The second-order Renyi entropy itself is continuous everywhere, however, the continuousness of its first-order derivative breaks down at the phase boundary. This means that the bosonic superfluid-insulator transitions can still be efficiently captured by the residual entanglement in our cluster MF treatment. Besides to the bosonic superfluid-insulator transitions, our cluster MF treatment may also be used to capture the signature of EE for other QPTs in quantum superlattice models.

Transverse forces on vortices in superfluids in a periodic potential

The paper analyzes the transverse forces (the Magnus and the Lorentz forces) on vortices in superfluids put into periodic potentials at $T=0$. The case of weak potential and the tight-binding limit described by the Bose-Hubbard model were addressed. The analysis was based on the balance of true momentum and quasimomentum. A special attention was paid to the superfluid close to the superfluid-insulator transition. In this area of the phase diagram the theory predicts the particle-hole symmetry line where the Magnus force changes sign with respect to that expected from the sign of velocity circulation. Our analysis has shown that the magnitude of the Magnus force is a continuous function at crossing the particle-hole symmetry line. This challenges the theory connecting the magnitude of the Magnus force with topological Chern numbers and predicting a jump at crossing this line. Disagreement is explained by the role of intrinsic pinning and guided vortex motion ignored in the topological approach. It is one more evidence that in general topological arguments are not sufficient for derivation of equations of vortex motion.

Cold Atoms in U(3) Gauge Potentials

We explore the effects of artificial U ( 3 ) gauge potentials on ultracold atoms. We study background gauge fields with both non-constant and constant Wilson loops around plaquettes, obtaining the energy spectra in each case. The scenario of metal–insulator transition for irrational fluxes is also examined. Finally, we discuss the effect of such a gauge potential on the superfluid–insulator transition for bosonic ultracold atoms.