Dipolar Bosons Research Articles

Erratum: Ground states of two-dimensional tilted dipolar bosons with density-induced hopping [Phys. Rev. A 103 , 043333 (2021)

Received 2 December 2021DOI:https://doi.org/10.1103/PhysRevA.105.019901©2022 American Physical SocietyPhysics Subject Headings (PhySH)Research AreasBose gasesCold gases in optical latticesDipolar gasesAtomic, Molecular & Optical

Ground states of two-dimensional tilted dipolar bosons with density-induced hopping

Motivated by recent experiments with ultracold magnetic atoms trapped in optical lattices where the orientation of atomic dipoles can be fully controlled by external fields, we study, by means of quantum Monte Carlo, the ground-state properties of dipolar bosons trapped in a two-dimensional lattice with density-induced hopping and where the dipoles are tilted along the $xz$ plane. We present ground-state phase diagrams of the above system at different tilt angles. We find that as the dipolar interaction increases, the superfluid phase at half-filling factor is destroyed in favor of either a checkerboard or stripe solid phase for tilt angle $\ensuremath{\theta}\ensuremath{\lesssim}{30}^{\ensuremath{\circ}}$ or $\ensuremath{\theta}\ensuremath{\gtrsim}{30}^{\ensuremath{\circ}}$, respectively. More interesting physics happens at tilt angles $\ensuremath{\theta}\ensuremath{\gtrsim}{58}^{\ensuremath{\circ}}$, where we find that as the dipolar interaction strength increases, solid phases first appear at filling factor lower than 0.5. Moreover, unlike what is observed at lower tilt angles, we find that at half filling, a stripe supersolid intervenes between the superfluid and stripe solid phase.

Finite-temperature instabilities of a two-dimensional dipolar Bose gas at arbitrary tilt angle

Advances in creating stable dipolar Bose systems, and ingenious box traps have generated tremendous interest. Theory study of dipolar bosons at finite temperature (T) has been limited. Motivated by these, we study 2D dipolar bosons at arbitrary tilt angle, $\theta$, using finite-T random phase approximation. We show that a comprehensive understanding of phases and instabilities at non-zero T can be obtained on concurrently considering dipole strength, density, temperature and $\theta$. We find the system to be in a homogeneous non-condensed phase that undergoes a collapse transition at large $\theta$, and a finite momentum instability, signaling a striped phase, at large dipolar strength; there are important differences with the T=0 case. At T = 0, BEC appears at critical dipolar strength, and at critical density. Our predictions for polar molecule system, $^{41}K^{87}Rb$, and $^{166}Er$ may provide tests of our results. Our approach may apply broadly to systems with long-range, anisotropic interactions.

Polarization angle dependence of the breathing mode in confined one-dimensional dipolar bosons

Probing the radial collective oscillation of a trapped quantum system is an accurate experimental tool to investigate interactions and dimensionality effects. We consider a fully polarized quasi-one dimensional dipolar quantum gas of bosonic dysprosium atoms in a parabolic trap at zero temperature. We model the dipolar gas with an effective quasi-one dimensional Hamiltonian in the single-mode approximation, and derive the equation of state using a variational approximation based on the Lieb-Liniger gas Bethe Ansatz wavefunction or perturbation theory. We calculate the breathing mode frequencies while varying polarization angles by a sum-rule approach, and find them in good agreement with recent experimental findings.

Staggered superfluid phases of dipolar bosons in two-dimensional square lattices

We study the quantum ground state of ultracold bosons in a two-dimensional square lattice. The bosons interact via the repulsive dipolar interactions and s-wave scattering. The dynamics is described by the extended Bose-Hubbard model including correlated hopping due to the dipolar interactions, the coefficients are found from the second quantized Hamiltonian using the Wannier expansion with realistic parameters. We determine the phase diagram using the Gutzwiller ansatz in the regime where the coefficients of the correlated hopping terms are negative and can interfere with the tunneling due to single-particle effects. We show that this interference gives rise to staggered superfluid and supersolid phases at vanishing kinetic energy, while we identify parameter regions at finite kinetic energy where the phases are incompressible. We compare the results with the phase diagram obtained with the cluster Gutzwiller approach and with the results found in one dimension using DMRG.

Comment on “Berezinskii-Kosterlitz-Thouless transition in two-dimensional dipolar stripes”

In a recent article [R. Bombin, F. Mazzanti and J. Boronat, Phys. Rev. A100, 063614 (2019)], it is contended that a two-dimensional system of dipolar bosons, with dipole moments aligned at particular angles with respect to the direction perpendicular to the plane of motion, featuring a "striped" crystalline ground state, in turn undergoes a Berezinskii-Kosterlitz-Thouless superfluid transition at low temperature, making it a two-dimensional supersolid. We show here that the results provided therein, obtained by means of Quantum Monte Carlo simulations, do not actually support such a conclusion. Rather, they are consistent with that expounded in our work [J. Low Temp. Phys. 196, 413 (2019)], namely that the striped ground state is insulating (i.e., non-superfluid in the conventional sense), essentially behaving like a system of quasi-one-dimensional, parallel independent chains. We attribute the incorrectness of the conclusion reached by Bombin et al. to the very small sizes of their simulated system, which do not allow for a reliable extrapolation to the thermodynamic limit.

Segregated quantum phases of dipolar bosonic mixtures in two-dimensional optical lattices

We identify the quantum phases in a binary mixture of dipolar bosons in two-dimensional optical lattices. Our study is motivated by the recent experimental realization of binary dipolar condensate mixtures of Er-Dy [Phys. Rev. Lett. 121, 213601 (2018)]. We model the system by using the extended two-species Bose-Hubbard model and calculate the ground-state phase diagrams by using mean-field theory. For selected cases we also obtain analytical phase boundaries by using the site-decoupled mean-field theory. For comparison we also examine the phase diagram of two-species Bose-Hubbard model. Our results show that the quantum phases with the long-range intraspecies interaction phase separate with no phase ordering. The introduction of the long-range interspecies interaction modifies the quantum phases of the system. It leads to the emergence of phase-separated quantum phases with phase ordering. The transition from the phase-separated quantum phases without phase ordering to phase ordered ones breaks the inversion symmetry.

Detecting One-Dimensional Dipolar Bosonic Crystal Orders via Full Distribution Functions.

We explore the ground states of a few dipolar bosons in optical lattices with incommensurate filling. The competition of kinetic, potential, and interaction energies leads to the emergence of a variety of crystal state orders with characteristic one- and two-body densities. We probe the transitions between these orders and construct the emergent state diagram as a function of the dipolar interaction strength and the lattice depth. We show that the crystal state orders can be observed using the full distribution functions of the particle number extracted from simulated single-shot images.

Breakdown of the mean field for dark solitons of dipolar bosons in a one-dimensional harmonic trap

We directly compare the mean-field and the many-body approach in a one-dimensional Bose system in a harmonic trap. Both contact and dipolar interactions are considered. We propose a multi-atom version of the phase imprinting method to generate dark solitons in the system. We begin with a general analysis of system dynamics and observe the emergence of a dark soliton and a shock wave. Center of mass and soliton motion become decoupled because the shock wave oscillates with the trap frequency and soliton does not. A detailed investigation of frequencies reveals significant differences between results obtained in the mean-field and the many-body pictures.

Discrete interactions between a few interlayer excitons trapped at a MoSe2\u2013WSe2 heterointerface

Inter-layer excitons (IXs) in hetero-bilayers of transition metal dichalcogenides (TMDs) represent an exciting emergent class of long-lived dipolar composite bosons in an atomically thin, near-ideal two-dimensional (2D) system. The long-range interactions that arise from the spatial separation of electrons and holes can give rise to novel quantum, as well as classical multi-particle correlation effects. Indeed, first indications of exciton condensation have been reported recently. In order to acquire a detailed understanding of the possible many-body effects, the fundamental interactions between individual IXs have to be studied. Here, we trap a tunable number of dipolar IXs (NIX ~ 1–5) within a nanoscale confinement potential induced by placing a MoSe2–WSe2 hetero-bilayer (HBL) onto an array of SiO2 nanopillars. We control the mean occupation of the IX trap via the optical excitation level and observe discrete sharp-line emission from different configurations of interacting IXs. The intensities of these features exhibit characteristic near linear, quadratic, cubic, quartic and quintic power dependencies, which allows us to identify them as different multiparticle configurations with NIX ~ 1–5. We directly measure the hierarchy of dipolar and exchange interactions as NIX increases. The interlayer biexciton (NIX = 2) is found to be an emission doublet that is blue-shifted from the single exciton by ΔE = (8.4 ± 0.6) meV and split by 2J = (1.2 ± 0.5) meV. The blueshift is even more pronounced for triexcitons ((12.4 ± 0.4) meV), quadexcitons ((15.5 ± 0.6) meV) and quintexcitons ((18.2 ± 0.8) meV). These values are shown to be mutually consistent with numerical modelling of dipolar excitons confined to a harmonic trapping potential having a confinement lengthscale in the range ell approx 3 nm. Our results contribute to the understanding of interactions between IXs in TMD hetero-bilayers at the discrete limit of only a few excitations and represent a key step towards exploring quantum correlations between IXs in TMD hetero-bilayers.

Superfluid phases induced by dipolar interactions

We determine the quantum ground state of dipolar bosons in a quasi-one-dimensional optical lattice and interacting via $s$-wave scattering. The Hamiltonian is an extended Bose-Hubbard model which includes hopping terms due to the interactions. We identify the parameter regime for which the coefficients of the interaction-induced hopping terms become negative. For these parameters we numerically determine the phase diagram for a canonical ensemble and by means of density matrix renormalization group. We show that at sufficiently large values of the dipolar strength there is a quantum interference between the tunneling due to single-particle effects and the one due to the interactions. Because of this phenomenon, incompressible phases appear at relatively large values of the single-particle tunneling rates. This quantum interference cuts the phase diagram into two different, disconnected superfluid phases. In particular, at vanishing kinetic energy, the phase is always superfluid with a staggered superfluid order parameter. These dynamics emerge from quantum interference phenomena between quantum fluctuations and interactions and shed light into their role in determining the thermodynamic properties of quantum matter.

Few-body bound states of two-dimensional bosons

We study clusters of the type A$_N$B$_M$ with $N\leq M\leq 3$ in a two-dimensional mixture of A and B bosons, with attractive AB and equally repulsive AA and BB interactions. In order to check universal aspects of the problem, we choose two very different models: dipolar bosons in a bilayer geometry and particles interacting via separable Gaussian potentials. We find that all the considered clusters are bound and that their energies are universal functions of the scattering lengths $a_{AB}$ and $a_{AA}=a_{BB}$, for sufficiently large attraction-to-repulsion ratios $a_{AB}/a_{BB}$. When $a_{AB}/a_{BB}$ decreases below $\approx 10$, the dimer-dimer interaction changes from attractive to repulsive and the population-balanced AABB and AAABBB clusters break into AB dimers. Calculating the AAABBB hexamer energy just below this threshold, we find an effective three-dimer repulsion which may have important implications for the many-body problem, particularly for observing liquid and supersolid states of dipolar dimers in the bilayer geometry. The population-imbalanced ABB trimer, ABBB tetramer, and AABBB pentamer remain bound beyond the dimer-dimer threshold. In the dipolar model, they break up at $a_{AB}\approx 2 a_{BB}$ where the atom-dimer interaction switches to repulsion.

Relaxation of Shannon entropy for trapped interacting bosons with dipolar interactions

We study the quantum many-body dynamics and entropy production triggered by an interaction quench of few dipolar bosons in an external harmonic trap. We solve the time-dependent many-body Schrödinger equation by using an in-principle numerically exact many-body method called the multiconfigurational time-dependent Hartree method for bosons (MCTDHB). We study the dynamical measures with high level of accuracy. We monitor the time evolution of the occupation in the natural orbitals and normalized first- and second-order Glauber’s correlation functions. In particular, we focus on the relaxation dynamics of the Shannon entropy. Comparison with the corresponding results for contact interactions is presented. We observe significant effects coming from the presence of the non-local part of the dipolar interaction. The relaxation process is very fast for dipolar bosons with a clear signature of a truly saturated maximum entropy state. We also discuss the connection between the entropy production and the occurrence of correlations and loss of coherence in the system. We identify the long-time relaxed state as a many-body state retaining only diagonal correlations in the first-order correlation function and building up anti-bunching effect in the second-order correlation function.

Entangled states of dipolar bosons generated in a triple-well potential

We study the generation of entangled states using a device constructed from dipolar bosons confined to a triple-well potential. Dipolar bosons possess controllable, long-range interactions. This property permits specific choices to be made for the coupling parameters, such that the system is integrable. Integrability assists in the analysis of the system via an effective Hamiltonian constructed through a conserved operator. Through computations of fidelity we establish that this approach, to study the time-evolution of the entanglement for a class of non-entangled initial states, yields accurate approximations given by analytic formulae.

Strongly Correlated Quantum Droplets in Quasi-1D Dipolar Bose Gas.

We exploit a few- to many-body approach to study strongly interacting dipolar bosons in the quasi-one-dimensional system. The dipoles attract each other while the short range interactions are repulsive. Solving numerically the multiatom Schrödinger equation, we discover that such systems can exhibit not only the well-known bright soliton solutions but also novel quantum droplets for a strongly coupled case. For larger systems, basing on microscopic properties of the found few-body solution, we propose a new equation for a density amplitude of atoms. It accounts for fermionization for strongly repelling bosons by incorporating the Lieb-Liniger energy in a local density approximation and approaches the standard Gross-Pitaevskii equation (GPE) in the weakly interacting limit. Not only does such a framework provide an alternative mechanism of the droplet stability, but it also introduces means to further analyze this previously unexplored quantum phase. In the limiting strong repulsion case, yet another simple multiatom model is proposed. We stress that the celebrated Lee-Huang-Yang term in the GPE is not applicable in this case.

Dipolar bosons in one dimension: The case of longitudinal dipole alignment

We study by quantum Monte Carlo simulations the low-temperature phase diagram of dipolar bosons confined to one dimension, with dipole moments aligned along the direction of particle motion. A hard core repulsive potential of varying range ($\sigma$) is added to the dipolar interaction, in order to ensure stability of the system against collapse. In the $\sigma\to 0$ limit the physics of the system is dominated by the potential energy and the ground state is quasi-crystalline; as $\sigma$ is increased the attractive part of the interaction weakens and the equilibrium phase evolves from quasi-crystalline to a non-superfluid liquid. At a critical value $\sigma_c$, the kinetic energy becomes dominant and the system undergoes a quantum phase transition from a self-bound liquid to a gas. In the gaseous phase with $\sigma\to\sigma_c$, at low density attractive interactions bring the system into a "weak" superfluid regime. However, gas-liquid coexistence also occurs, as a result of which the topologically protected superfluid regime is not approached.

Variational Bethe ansatz approach for dipolar one-dimensional bosons

We propose a variational approximation to the ground state energy of a one-dimensional gas of interacting bosons on the continuum based on the Bethe Ansatz ground state wavefunction of the Lieb-Liniger model. We apply our variational approximation to a gas of dipolar bosons in the single mode approximation and obtain its ground state energy per unit length. This allows for the calculation of the Tomonaga-Luttinger exponent as a function of density and the determination of the structure factor at small momenta. Moreover, in the case of attractive dipolar interaction, an instability is predicted at a critical density, which could be accessed in lanthanide atoms.

Quantum phases of canted dipolar bosons in a two-dimensional square optical lattice

We consider a minimal model to describe the quantum phases of ultracold dipolar bosons in two-dimensional (2D) square optical lattices. The model is a variation of the extended Bose-Hubbard model and apt to study the quantum phases arising from the variation in the tilt angle $\theta$ of the dipolar bosons. At low tilt angles $0^{\circ}\leqslant\theta\apprle25^{\circ}$, the ground state of the system are phases with checkerboard order, which could be either checkerboard supersolid or checkerboard density wave. For high tilt angles $55^{\circ}\apprge\theta\apprge35^{\circ}$, phases with striped order of supersolid or density wave are preferred. In the intermediate domain $25^{\circ}\apprle\theta\apprle35^{\circ}$ an emulsion or SF phase intervenes the transition between the checkerboard and striped phases. The attractive interaction dominates for $\theta\apprge55^{\circ}$, which renders the system unstable and there is a density collapse. For our studies we use Gutzwiller mean-field theory to obtain the quantum phases and the phase boundaries. In addition, we calculate the phase boundaries between an incompressible and a compressible phase of the system by considering second order perturbation analysis of the mean-field theory. The analytical results, where applicable, are in excellent agreement with the numerical results.

Patterned Supersolids in Dipolar Bose Systems

We study by means of first principle Quantum Monte Carlo simulations the ground state phase diagram of a system of dipolar bosons with aligned dipole moments, and with the inclusion of a two-body repulsive potential of varying range. The system is shown to display a supersolid phase in a relatively broad region of the phase diagram, featuring different crystalline patterns depending on the density and on the range of the repulsive part of the interaction (scattering length). The supersolid phase is sandwiched between a classical crystal of parallel filaments and a homogeneous superfluid phase. We show that a "roton" minimum appears in the elementary excitation spectrum of the superfluid as the system approaches crystallization.

Correlation Dynamics of Dipolar Bosons in 1D Triple Well Optical Lattice

The concept of spontaneous symmetry breaking and off-diagonal long-range order (ODLRO) are associated with Bose–Einstein condensation. However, as in the system of reduced dimension the effect of quantum fluctuation is dominating, the concept of ODLRO becomes more interesting, especially for the long-range interaction. In the present manuscript, we study the correlation dynamics triggered by lattice depth quench in a system of three dipolar bosons in a 1D triple-well optical lattice from the first principle using the multiconfigurational time-dependent Hartree method for bosons (MCTDHB). Our main motivation is to explore how ODLRO develops and decays with time when the system is brought out-of-equilibrium by a sudden change in the lattice depth. We compare results of dipolar bosons with contact interaction. For forward quench ( V f > V i ) , the system exhibits the collapse–revival dynamics in the time evolution of normalized first- and second-order Glauber’s correlation function, time evolution of Shannon information entropy both for the contact as well as for the dipolar interaction which is reminiscent of the one observed in Greiner’s experiment [Nature, 415 (2002)]. We define the collapse and revival time ratio as the figure of merit ( τ ) which can uniquely distinguish the timescale of dynamics for dipolar interaction from that of contact interaction. In the reverse quench process ( V i > V f ) , for dipolar interaction, the dynamics is complex and the system does not exhibit any definite time scale of evolution, whereas the system with contact interaction exhibits collapse–revival dynamics with a definite time-scale. The long-range repulsive tail in the dipolar interaction inhibits the spreading of correlation across the lattice sites.