Superfluid Order Parameter Research Articles

Multisite mean-field theory for cold bosonic atoms in optical lattices

We present a detailed derivation of a multi-site mean-field theory (MSMFT) used to describe the Mott-insulator to superfluid transition of bosonic atoms in optical lattices. The approach is based on partitioning the lattice into small clusters which are decoupled by means of a mean field approximation. This approximation invokes local superfluid order parameters defined for each of the boundary sites of the cluster. The resulting MSMFT grand potential has a non-trivial topology as a function of the various order parameters. An understanding of this topology provides two different criteria for the determination of the Mott insulator superfluid phase boundaries. We apply this formalism to $d$-dimensional hypercubic lattices in one, two and three dimensions, and demonstrate the improvement in the estimation of the phase boundaries when MSMFT is utilized for increasingly larger clusters, with the best quantitative agreement found for $d=3$. The MSMFT is then used to examine a linear dimer chain in which the on-site energies within the dimer have an energy separation of $\Delta$. This system has a complicated phase diagram within the parameter space of the model, with many distinct Mott phases separated by superfluid regions.

Ground-state properties of the disordered spin-1 Bose-Hubbard model: A stochastic mean-field theory study

We study the ground state of the disordered Bose-Hubbard model for spin-1 particles by means of the stochastic mean-field theory. This approach enables the determination of the probability distributions of various physical quantities, such as the superfluid order parameter, the average site occupation number, the standard deviation of the occupation per site and the square of the spin operator per site. We show how a stochastic method, previously used in the study of localization, can be flexibly used to solve the relevant equations with great accuracy. We have determined the phase diagram, which exhibits three phases: the polar superfluid, the Mott insulating, and the Bose glass. A complete characterization of the physical properties of these phases has been established.

Tanʼs universal contact and collective oscillations of strongly interacting Fermi gases

We study strongly interacting two component Fermi gas near a Feshbach resonance. By using a ground state energy functional constructed based on asymptotic limits and Monte Carlo calculations, we calculate the contact, structure factor, and collective oscillation frequencies in the BCS–BEC crossover region. The calculated contact and structure factor show excellent agreements with recent experiments. We compare these results with a standard mean-field theory and find that the contact is proportional to the square of superfluid order parameter. Further, we present the chemical potential and the polytropic index in terms of homogeneous energy and the contact.

Quantum Shock Waves and Domain Walls in the Real-Time Dynamics of a Superfluid Unitary Fermi Gas

We show that in the collision of two superfluid fermionic atomic clouds one observes the formation of quantum shock waves as discontinuities in the number density and collective flow velocity. Domain walls, which are topological excitations of the superfluid order parameter, are also generated and exhibit abrupt phase changes by π and slower motion than the shock waves. The domain walls are distinct from the gray soliton train or number density ripples formed in the wake of the shock waves and observed in the collisions of superfluid bosonic atomic clouds. Domain walls with opposite phase jumps appear to collide elastically.

Mean field theory of superglasses

We study the interplay of superfluidity and glassy ordering of hard core bosons with random, frustrating interactions. This is motivated by bosonic systems such as amorphous supersolid, disordered superconductors with preformed pairs, and helium in porous media. We analyze the fully connected mean field version of this problem, which exhibits three low-temperature phases, separated by two continuous phase transitions: an insulating, glassy phase with an amorphous frozen density pattern, a nonglassy superfluid phase, and an intermediate phase, in which both types of order coexist. We elucidate the nature of the phase transitions, highlighting in particular the role of glassy correlations across the superfluid-insulator transition. The latter suppress superfluidity down to $T=0$, due to the depletion of the low-energy density of states, unlike in the standard BCS scenario. Further, we investigate the properties of the coexistence (superglass) phase. We find anticorrelations between the local order parameters and a nonmonotonous superfluid order parameter as a function of $T$. The latter arises due to the weakening of the glassy correlation gap with increasing temperature. Implications of the mean field phenomenology for finite dimensional bosonic glasses with frustrating Coulomb interactions are discussed.

Phase Separation in Bose–Bose Mixtures in an Optical Lattice

We study the ground-state properties of mixtures of strongly interacting bosonic atoms in an optical lattice. Applying a mean-field approximation to the Hubbard model for Bose-Bose mixtures, we calculate the densities and superfluid order parameters for both species. Due to the repulsive interaction between the two species, the system exhibits phase separation. First, in the extreme limit of the zero-hopping case, we derive analytical expressions for the phase boundaries. In particular, we derive the conditions for phase separation in the Mott insulator phase. We find that the conditions for the phase separation depend on the on-site interactions as well as the occupation numbers. In particular, we show that the coexisting state appears by varying the on-site inter-species interaction. We also show the phase diagram of the finite hopping case. Second, we calculate the spatial density profile of $^{87}$Rb-$^{41}$K mixtures in the combined potential of a parabolic trap and an optical lattice using the local density approximation. We fixed the number of $^{87}$Rb and varied the number of $^{41}$K, and used the parameters estimated by experiments. We show that the phase separated $^{87}$Rb-$^{41}$K mixtures distribute like in a parabolic trap case. Furthermore, we find that phase separated mixtures distribute a nesting structure.

BCS-BEC Crossover and Topological Phase Transition in 3D Spin-Orbit Coupled Degenerate Fermi Gases

We investigate the BCS-BEC crossover in three-dimensional degenerate Fermi gases in the presence of spin-orbit coupling (SOC) and Zeeman field. We show that the superfluid order parameter destroyed by a large Zeeman field can be restored by the SOC. With increasing strengths of the Zeeman field, there is a series of topological quantum phase transitions from a nontopological superfluid state with fully gapped fermionic spectrum to a topological superfluid state with four topologically protected Fermi points (i.e., nodes in the quasiparticle excitation gap) and then to a second topological superfluid state with only two Fermi points. The quasiparticle excitations near the Fermi points realize the long-sought low-temperature analog of Weyl fermions of particle physics. We show that the topological phase transitions can be probed using the experimentally realized momentum-resolved photoemission spectroscopy.

Schwinger-Keldysh approach to out-of-equilibrium dynamics of the Bose-Hubbard model with time-varying hopping

We study the real time dynamics of the Bose Hubbard model in the presence of time-dependent hopping allowing for a finite temperature initial state. We use the Schwinger-Keldysh technique to find the real-time strong coupling action for the problem at both zero and finite temperature. This action allows for the description of both the superfluid and Mott insulating phases. We use this action to obtain dynamical equations for the superfluid order parameter as hopping is tuned in real time so that the system crosses the superfluid phase boundary. We find that under a quench in the hopping, the system generically enters a metastable state in which the superfluid order parameter has an oscillatory time dependence with a finite magnitude, but disappears when averaged over a period. We relate our results to recent cold atom experiments.

Spectroscopy of the soliton lattice formation in quasi-one-dimensional fermionic superfluids with population imbalance

Motivated by recent experiments in low-dimensional trapped fermionic superfluids, we study a quasi-one-dimensional (quasi-1D) superfluid with a population imbalance between two hyperfine states using an exact mean-field solution for the order parameter. When an effective "magnetic field" exceeds a critical value, the superfluid order parameter develops spatial inhomogeneity in the form of a soliton lattice. The soliton lattice generates a band of quasiparticle states inside the energy gap, which originate from the Andreev bound states localized at the solitons. Emergence of the soliton lattice is accompanied by formation of a spin-density wave, with the majority fermions residing at the points in space where the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) order parameter vanishes. We discuss possibilities for experimental detection of the quasi-1D FFLO state using elastic and inelastic optical Bragg scattering and radio-frequency spectroscopy. We show that these measurements can provide necessary information for unambiguous identification of the spatially-inhomogeneous quasi-1D FFLO state and the soliton lattice formation.

Low-temperature properties of the infinite-dimensional attractive Hubbard model

We investigate the attractive Hubbard model in infinite spatial dimensions by combining dynamical mean-field theory with a strong-coupling continuous-time quantum Monte Carlo method. By calculating the superfluid order parameter and the density of states, we discuss the stability of the superfluid state. In the intermediate-coupling region above the critical temperature, the density of states exhibits a heavy fermion behavior with a quasiparticle peak in the dense system, while a dip structure appears in the dilute system. The formation of the superfluid gap is also addressed.

Dislocation-induced superfluidity in a model supersolid

Motivated by recent experiments on the supersolid behavior of $^4$He, we study the effect of an edge dislocation in promoting superfluidity in a Bose crystal. Using Landau theory, we couple the elastic strain field of the dislocation to the superfluid density, and use a linear analysis to show that superfluidity nucleates on the dislocation before occurring in the bulk of the solid. Moving beyond the linear analysis, we develop a systematic perturbation theory in the weakly nonlinear regime, and use this method to integrate out transverse degrees of freedom and derive a one-dimensional Landau equation for the superfluid order parameter. We then extend our analysis to a network of dislocation lines, and derive an XY model for the dislocation network by integrating over fluctuations in the order parameter. Our results show that the ordering temperature for the network has a sensitive dependence on the dislocation density, consistent with numerous experiments that find a clear connection between the sample quality and the supersolid response.

Conventional and charge-six superfluids from melting hexagonal Fulde-Ferrell-Larkin-Ovchinnikov phases in two dimensions

We consider defect mediated melting of Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) and pair density wave (PDW) phases in two dimensions. Examining mean-field ground states in which the spatial oscillations of the FFLO/PDW superfluid order parameter exhibit hexagonal lattice symmetry, we find that thermal melting leads to a variety of novel phases. We find that a spatially homogeneous charge six superfluid can arise from melting a hexagonal vortex-anitvortex lattice FFLO/PDW phase. The charge six superfluid has an order parameter corresponding to a bound state of six fermions. We further find that a hexagonal vortex-free FFLO/PDW phase can melt to yield a conventional (charge two) homogeneous superfluid. A key role is played by topological defects that combine fractional vortices of the superfluid order and fractional dislocations of the lattice order.

πjunction and spontaneous current state in a superfluid Fermi gas

We discuss an idea to realize a spontaneous current in a superfluid Fermi gas. When a polarized Fermi superfluid (${N}_{\ensuremath{\uparrow}}>{N}_{\ensuremath{\downarrow}}$, where ${N}_{\ensuremath{\sigma}}$ is the number of atoms in the hyperfine state described by pseudospin $\ensuremath{\sigma}=\ensuremath{\uparrow},\ensuremath{\downarrow}$) is loaded onto a ring-shaped trap with a weak potential barrier, some excess atoms ($\ensuremath{\Delta}N={N}_{\ensuremath{\uparrow}}\ensuremath{-}{N}_{\ensuremath{\downarrow}}$) are localized around the barrier. As shown in our previous paper [T. Kashimura, S. Tsuchiya, and Y. Ohashi, Phys. Rev. A 82, 033617 (2010)], this polarized potential barrier works as a $\ensuremath{\pi}$ junction in the sense that the superfluid order parameter changes its sign across the barrier. Because of this, the phase of the superfluid order parameter outside the junction is shown to be twisted by $\ensuremath{\pi}$ along the ring, which naturally leads to a circulating supercurrent. While the ordinary supercurrent state is obtained as a metastable state, this spontaneous current state is shown to be more stable than the case with no current. Our results indicate that localized excess atoms would be useful for the manipulation of the superfluid order parameter in cold Fermi gases.

Formation of magnetic impurities and pair-breaking effect in a superfluid Fermi gas

We theoretically investigate the introduction of magnetic impurities into a superfluid Fermi gas. In the presence of a population imbalance (${N}_{\ensuremath{\uparrow}}>{N}_{\ensuremath{\downarrow}}$, where ${N}_{\ensuremath{\sigma}}$ is the number of Fermi atoms with pseudospin $\ensuremath{\sigma}=\ensuremath{\uparrow},\ensuremath{\downarrow}$), we show that nonmagnetic potential scatterers embedded in the system are magnetized in the sense that some of the excess spin-$\ensuremath{\uparrow}$ atoms are localized around them. They destroy the superfluid order parameter around them, as in the case of the magnetic impurity effect discussed in the superconductivity literature. This pair-breaking effect naturally leads to localized excited states below the superfluid excitation gap. To confirm our idea in a simple manner, we consider an attractive Fermi-Hubbard model within the mean-field theory at $T=0$. We determine consistent superfluid properties around a nonmagnetic impurity, such as the superfluid order parameter, local population imbalance, as well as the single-particle density of states, in the presence of a population imbalance. Since competition between superconductivity and magnetism is one of the most fundamental problems in condensed-matter physics, our results would be useful for the study of this important issue in cold Fermi gases.

Application of a multisite mean-field theory to the disordered Bose-Hubbard model

We present a multi-site formulation of mean-field theory applied to the disordered Bose-Hubbard model. In this approach the lattice is partitioned into clusters, each isolated cluster being treated exactly, with inter-cluster hopping being treated approximately. The theory allows for the possibility of a different superfluid order parameter at every site in the lattice, such as what has been used in previously published site-decoupled mean-field theories, but a multi-site formulation also allows for the inclusion of spatial correlations allowing us, e.g., to calculate the correlation length (over the length scale of each cluster). We present our numerical results for a two-dimensional system. This theory is shown to produce a phase diagram in which the stability of the Mott insulator phase is larger than that predicted by site-decoupled single-site mean-field theory. Two different methods are given for the identification of the Bose glass-to-superfluid transition, one an approximation based on the behaviour of the condensate fraction, and one of which relies on obtaining the spatial variation of the order parameter correlation. The relation of our results to a recent proposal that both transitions are non self-averaging is discussed.

Nonequilibrium Dynamics of the Bose-Hubbard Model: A Projection-Operator Approach

We study the phase diagram and nonequilibrium dynamics involving ramp of the hopping amplitude J(t)=Jt/τ with ramp time τ of the Bose-Hubbard model at zero temperature using a projection-operator formalism which allows us to incorporate the effects of quantum fluctuations beyond mean-field approximations in the strong-coupling regime. Our formalism yields a phase diagram which provides a near exact match with quantum Monte Carlo results in three dimensions. We also compute the residual energy Q, the superfluid order parameter Δ(t), the equal-time order parameter correlation function C(t), and the wave function overlap F which yields the defect formation probability P during nonequilibrium dynamics of the model. We find that Q, F, and P do not exhibit the expected universal scaling. We explain this absence of universality and show that our results compare well with recent experiments.

Cold atoms in a rotating optical lattice with nearest-neighbor interactions

Extended Bose Hubbard models with nearest neighbour interaction describe minimally the effect of long range interaction on ultra cold atoms in deep optical lattices. Rotation of such optical lattices subject such neutral cold atoms to the effect of an artificial magnetic field. The modification of the phase boundaries of the density wave and Mott Insulator phases due to this rotation are shown to be related to the edge spectrum of spinorial and scalar Harper equation. Corresponding profiles of the checkerboard vortex states with sublattice modulated superfluid order parameter near density wave phase boundary are calculated.

Evidence for orbital superfluidity in the P-band of a bipartite optical square lattice

The successful emulation of the Hubbard model in optical lattices has stimulated world wide efforts to extend their scope to also capture more complex, incompletely understood scenarios of many-body physics. Unfortunately, for bosons, Feynmans fundamental "no-node" theorem under very general circumstances predicts a positive definite ground state wave function with limited relevance for many-body systems of interest. A promising way around Feynmans statement is to consider higher bands in optical lattices with more than one dimension, where the orbital degree of freedom with its intrinsic anisotropy due to multiple orbital orientations gives rise to a structural diversity, highly relevant, for example, in the area of strongly correlated electronic matter. In homogeneous two-dimensional optical lattices, lifetimes of excited bands on the order of a hundred milliseconds are possible but the tunneling dynamics appears not to support cross-dimensional coherence. Here we report the first observation of a superfluid in the $P$-band of a bipartite optical square lattice with $S$-orbits and $P$-orbits arranged in a chequerboard pattern. This permits us to establish full cross-dimensional coherence with a life-time of several ten milliseconds. Depending on a small adjustable anisotropy of the lattice, we can realize real-valued striped superfluid order parameters with different orientations $P_x \pm P_y$ or a complex-valued $P_x \pm i P_y$ order parameter, which breaks time reversal symmetry and resembles the $\pi$-flux model proposed in the context of high temperature superconductors. Our experiment opens up the realms of orbital superfluids to investigations with optical lattice models.

Superfluid density of states and pseudogap phenomenon in the BCS-BEC crossover regime of a superfluid Fermi gas

We investigate single-particle excitations and strong-coupling effects in the BCS-BEC crossover regime of a superfluid Fermi gas. Including phase and amplitude fluctuations of the superfluid order parameter within a $T$-matrix theory, we calculate the superfluid density of states (DOS), as well as single-particle spectral weight, over the entire BCS-BEC crossover region below the superfluid transition temperature $T_{\rm c}$. We clarify how the pseudogap in the normal state evolves into the superfluid gap, as one passes through $T_{\rm c}$. While the pseudogap in DOS continuously evolves into the superfluid gap in the weak-coupling BCS regime, the superfluid gap in the crossover region is shown to appear in DOS after the pseudogap disappears below $T_{\rm c}$. In the phase diagram with respect to the temperature and interaction strength, we determine the region where strong pairing fluctuations dominate over single-particle properties of the system. Our results would be useful for the study of strong-coupling phenomena in the BCS-BEC crossover regime of a superfluid Fermi gas.

Superfluid-ferromagnet-superfluid junction and theπphase in a superfluid Fermi gas

We investigate the possibility of a superfluid-ferromagnet-superfluid (SFS) junction in a superfluid Fermi gas. To examine this possibility in a simple manner, we consider an attractive Hubbard model at $T=0$ within the mean-field theory. When a potential barrier is embedded in a superfluid Fermi gas with population imbalance (${N}_{\ensuremath{\uparrow}}>{N}_{\ensuremath{\downarrow}}$, where ${N}_{\ensuremath{\sigma}}$ is the number of atoms with pseudospin $\ensuremath{\sigma}= \ensuremath{\uparrow}$, $\ensuremath{\downarrow}$), this barrier is shown to be magnetized in the sense that excess $\ensuremath{\uparrow}$-spin atoms are localized around it. The resulting superfluid Fermi gas is spatially divided into two by this ferromagnet, so that one obtains a junction similar to the superconductor-ferromagnet-superconductor junction discussed in superconductivity. Indeed, we show that the so-called $\ensuremath{\pi}$ phase, which is a typical phenomenon in the SFS junction, is realized, where the superfluid order parameter changes its sign across the junction. Our results would be useful for the study of magnetic effects on fermion superfluidity using an ultracold Fermi gas.