Many-body physics of low-density dipolar bosons in box potentials

TL;DR

This paper reviews the many-body physics of low-density dipolar bosons in box potentials, focusing on crystallization phenomena, interactions, and phases, with insights into correlations and superfluidity.

Contribution

It provides a comprehensive review of crystallization and phase behavior in ultracold dipolar bosons, highlighting the interplay of interactions and density.

Findings

Crystallization occurs due to dipole interactions and density effects.

Correlations and superfluidity are key to phase characterization.

Numerical methods effectively analyze many-body phases.

Abstract

Crystallization is a generic phenomenon in classical and quantum mechanics arising in a variety of physical systems. In this work we focus on a specific platform, ultracold dipolar bosons, which can be realized in experiments with dilute gases. We review the relevant ingredients leading to crystallization, namely the interplay of contact and dipole-dipole interactions and system density, as well as the numerical algorithm employed. We characterize the many-body phases investigating correlations and superfluidity.

Figures

Equations

$$V(\mathbf{r})=\frac{4\pi {\hslash }^2}{m}{a}_s\delta (\mathbf{r})+\frac{C_{dd}}{4\pi }\frac{\left(1-3{\cos }^2\theta \right)}{r^3}$$

$${\cal V}({\bf r})=\left\{\begin{array}[]{ll}\infty&\text{if}~{}r<r_{0},\\ \\ \frac{C_{dd}}{4\pi}\,\frac{1-3\cos^{2}\theta}{r^{3}}&\text{if}~{}r\geq r_{0}.\end{array}\right.$$

$$\mathcal{H}=-\frac{1}{2}\sum _{i=1}^N{\nabla }_i^2+\sum _{i<j}^N\mathcal{V}({\mathbf{r}}_i-{\mathbf{r}}_j).$$

$$n{r}_0^3=\frac{n{a}_d^3}{(a_d/r_0{)}^3}.$$

$$n_G(r)=\frac{N_p}{\pi a_{\perp }^{2}L}{e}^{-\frac{r^2}{a_{\perp }^2}}$$

$$⟨r⟩\approx \frac{1}{n_f^{1/3}}={\left(\frac{\pi a_{\perp }^{2}L}{N_p}\right)}^{1/3}=2.1r_0$$

Full text

\pubvolume

xx \issuenum1 \articlenumber5

\historyReceived: date; Accepted: date; Published: date \TitleMany-body physics of low-density dipolar bosons in box potentials \AuthorTommaso Macrì 1 and Fabio Cinti 2 \AuthorNamesFirstname Lastname, Firstname Lastname and Firstname Lastname

1 Introduction

Cluster phases are ubiquitous in diverse systems, from biological structures Butenko et al. (2009) to soft matter Likos (2001); Díaz-Méndez et al. (2017) to ultra-cold atomic and molecular condensates Cinti et al. (2014). In many respects, particle aggregates display fascinating properties, whose features can be mainly described microscopically by means of effective two-body potentials. In a simple mean field picture, a cluster crystalline phase (or, more formally, a spontaneous breaking of translational symmetry) appears when the Fourier transform (when it exists) of the particle-particle interaction exhibits a negative region with a local minima. The cluster phase has then a density modulation with wavelength corresponding to those minima. The validity of this condition has been investigated also in the quantum regime, and it holds also if one considers fermionic or bosonic quantum gases Cinti et al. (2010). In particular, having potentials displaying a soft-shoulder-like shape at short distance is a crucial ingredient for designing novel quantum phases. As an example, alkaline atoms off-resonantly excited to Rydberg states display an effective interaction with a soft-core at short distances that causes a supersolid phase Henkel et al. (2012); Macrì et al. (2014, 2013); Cinti et al. (2014). A supersolid is a phase of matter that simultaneously accommodates diagonal as well as off-diagonal long-range order, which means that particles self-assemble into a rigid, regular crystal but at the same time can they can rotate with a corresponding reduced moment of inertia, which is an indication of the presence of a finite superfluid fraction Cinti et al. (2010); Saccani et al. (2012, 2011); Boninsegni (2012); Moroni and Boninsegni (2014); Boninsegni and Prokof’ev (2012)

Recently, ground-breaking experiments with bosonic dipolar condensates (Dysprosium and Erbium) demonstrated the existence of self-bound droplets in trapped configurations as well as in free space Kadau et al. (2016); Chomaz et al. (2016). These experiments provides a useful isolated system to probe general quantum-mechanical properties related to cluster phases. A mean-field treatment mainly based on a generalized nonlocal nonlinear Schrödinger equation reveal a good agreement with the density distribution and the excitation spectra observed in the laboratory Wächtler and Santos (2016). The same methodology has been successfully implemented to Bose-Bose mixtures Cappellaro et al. (2018, 2017). In addition, quantum Monte Carlo simulations (QMC) have investigated the zero temperature phase diagram of bosons interacting via dipolar interactions in three dimensions in free space Cinti et al. (2017). Other theoretical works have shown how the physical nature of self-bound droplets smoothly evolves from classical to quantum mechanical as the range of the repulsive two-body potential increases Cinti and Boninsegni (2017); Cidrim et al. (2018).

In the present paper, by means of first-principles numerical simulations, we characterize the many-body physics of bosons interacting via anisotropic dipole-dipole potentials with different densities and interaction strengths. In particular we carry out simulation to better understand the quantum behaviour of the system that supports recent experimental findings about the stability of droplets.

The rest of the paper is organized as follows: In section 2 we introduce the microscopic potential that describes the physics of dipolar systems in three-dimension, whereas in section 3 we present the details of the many-body properties known in literature.In section 4, we illustrate our results, outlining our conclusions in section 5.

2 Two-body physics of dipole-dipole interactions

The two body-potential describing a short-range and an anisotropic long-range dipolar interaction reads

[TABLE]

where $m$ is the particle mass, while $a_{s}$ is the s-wave scattering length and $a_{d}=mC_{dd}/4\pi\hbar^{2}$ the dipolar length representing the characteristic length of the dipole interactions $C_{dd}/4\pi$ Lahaye et al. (2009). $\theta$ denotes the angle between the vector ${\bf r}$ and the polarization axis of the dipoles, the $z$-axis. Equation (1) applies both to vertically aligned dipoles interacting via magnetic or electric dipole moments Lahaye et al. (2009). The pseudo-potential $V({\bf r})$ has been already discussed by Yi and Yu Yi and You (2000, 2001) making use of the Born approximation. More realistic model potentials, that include also a Van der Walls short range potentials, were recently considered in Ołdziejewski and Jachymski (2016) for a more detailed description of Dy two-body collisions. These potentials eventually lead to the same effective potential of Yi and You (2000, 2001) , with a dipolar length depending on the short range scattering length.

In our first-principle simulations based on a QMC technique (see section 3), $a_{d}$ as well as $a_{s}$ are defined in terms of a unit length $r_{0}$, $\varepsilon_{0}=\hbar^{2}/mr^{2}_{0}$ being the corresponding energy with the effective pairwise potential ${\cal V}({\bf r})$:

[TABLE]

${\cal V}({\bf r})$ describes a short-range hard core with cut-off at $r_{0}$ and the usual anisotropic dipolar component. The authors of Bortolotti et al. (2006); Ronen et al. (2007) have computed numerically the full low-energy scattering amplitude for the same model potential in equation (2). A comparison with equation (1) found a quite good agreement.

As a useful example, we now consider the Dysprosium’s isotopes 162Dy ($a_{d}=129.2\,a_{0}$, $a_{0}$ being the Bohr radius) and 164Dy ($a_{d}=130.8\,a_{0}$), relevant to the experimental observations by Kadau et al. Kadau et al. (2016). In addition, the background scattering lengths were recently measured for both isotopes (6-8): $a_{s}=122\,a_{0}$ and $a_{s}=92a_{0}$, for 162Dy and 164Dy, respectively.

In Bortolotti et al. (2006) the authors extracted the relation between the dipolar length $a_{d}$ and the scattering length $a_{s}$ for the model potential ${\cal V}({\bf r})$ as shown in fig. 1. For $a_{d}/r_{0}\approx 2.9$ and $a_{d}/r_{0}\approx 6.5$ we observe two resonances at which the scattering length $a_{s}$ diverges. Moreover, one also extracts the value of $\epsilon_{d}=a_{d}/a_{s}$ as a function of the corresponding ratios $a_{s}/r_{0}$ and $a_{d}/r_{0}$ (see Fig. 2). From the calculation of $\epsilon_{d}$ for the two isotopes we can also compute the value of the parameter $r_{0}$. In particular, for 162Dy we have $r_{0}=177a_{0}$, while considering 164 Dy we get $r_{0}=154a_{0}$.

3 Many-body properties

The present section is devoted to discuss many-body properties of bosons interacting via a dipolar interaction. The quantum-mechanical Hamiltonian of an ensemble of $N$ interacting identical bosons reads

[TABLE]

The Hamiltonian is reported in units of $r_{0}$ and $\varepsilon_{0}$. It is important to note that the zero-temperature physics is exclusively controlled by the dimensionless interaction strength ${\cal V}_{0}=mC_{dd}/4\pi\hbar^{2}r_{0}$ and the dimensionless density $nr^{3}_{0}$. Recently the zero-temperature phase diagram of equation (3) has been throughly investigated applying QMC techniques Cinti et al. (2017). In particular, considering low ${\cal V}_{0}$ (that is $\lesssim 2.1$ in the present units), simulations agree with the mean-field phase boundary predicted by standard Bogoliubov analysis. Upon increasing the dipole-dipole interaction we can distinguish two different phases. At low densities ($nr^{3}_{0}\lesssim 5\cdot 10^{-3}$) the system displays a cluster phase with vanishing superfluidity and characterized by droplet structures with few particles. For higher densities, one observes a phase marked by elongated filaments with an anisotropic superfluid fraction; that is, whereas the superfluidity is observed along the direction of the filaments, it results greatly suppressed on the corresponding orthogonal plane, excluding then the occurrence of a global supersolid phase. In addition, finite-temperature calculations confirmed the stability of the filament phase against thermal fluctuations and provide an estimate of the superfluid fraction in the weak coupling limit in the framework of the Landau two-fluid model. It is important to stress that both cluster phase and filament phase extend from a small positive induced contact potential ($\epsilon_{d}\gtrsim 1$), corresponding to the experimentally relevant regime of Kadau et al. (2016), to the strongly coupled limit of large dipolar interactions.

Theoretical results obtained solving Hamiltonian (3) can be compared with recent experiments with 162Dy and 164Dy making use of the following condition:

[TABLE]

As an example, setting an average density $n=5\cdot 10^{20}m^{-3}$, and with $na^{3}_{d}=2\cdot 10^{-4}$, taking $a_{d}=130a_{0}$ Kadau et al. (2016); Schmitt et al. (2016), we have shown in Ref.Cinti et al. (2017) that experimental results lie in the transition region (within the error bars) from superfluid to cluster phase by tuning the scattering length.

4 Results

4.1 Methodology

The low-temperature equilibrium properties of the systems described via the Hamiltonian (3) have been investigated employing first-principle computer simulations based on the worm algorithm in the continuous-space path integral representation Boninsegni et al. (2006a), which allows one to essentially obtain the exact thermodynamics properties of Bose systems. Over the last ten years this computational technique has been successfully tested on a large variety of bosonic systems, including, for instance, 4He Boninsegni et al. (2006b), Rydberg atoms Cinti et al. (2014) as well as dipolar systems Cinti (2015); Cinti et al. (2017); Lechner et al. (2015); Díaz-Méndez et al. (2015); Jain et al. (2011). For the details of the implementation of the worm algorithm the reader may refer to Boninsegni et al. Boninsegni et al. (2006a). Due to the intrinsic anisotropy of the potential (2), we have implemented in this work the so-called primitive approximation for the imaginary time propagator, which requires a greater number of time slices with respect to more complex propagators, such as the pair-product ansatz Ceperley (1995) or the fourth-order one Jang et al. (2001). All the results reported below are extrapolated to the limit of zero temperature. In this work we study the equilibrium properties of scaled average density $nr^{3}_{0}<0.1$ and different interaction strengths ${\cal V}_{0}$. We work with $N$ atoms in a cubic box of linear dimension $L$ (fixing the density $n=N/L^{3}$) and using periodic boundary conditions. We have performed simulations with $N$ between 100 and 400; in such a limit one can exclude finite-size issues. From these simulations we obtain density profiles, radial correlation functions, and the superfluid fraction. The limit of zero temperature is obtained lowering the temperature until observables (e.g. energy, or superfluid fraction) do not change on further decreasing $T$.

4.2 Density profile of the Filament Phase

In this section we present results applying the QMC method discussed in section 4.1. We consider for now a regime of high density such as $nr_{0}^{3}=10^{-2}$. For this peculiar case our calculations obtained a stable filament phase at temperature much lower then the transition temperature from normal to superfluid liquid. Fig. 3a shows to a filament phase with $N=100$ and $P=1000$, $P$ being the number of imaginary-time slices. This snapshot displays the projection of the single-particle imaginary-time evolution onto the real space. This representation of the condensate provides information of the delocalisation of particles in the box, and therefore of their probability distribution.

The system shapes in four different well defined filaments, that we now characterize. We fit the single-filament density with a gaussian $n_{G}(r)$

[TABLE]

and extract the average number of particles within each filament. The inset in Fig. 3b shows the width of the gaussian $a_{\perp}$ in units of $r_{0}$, as a function of the particle number $N_{p}$. The estimate of the average inter-particle distance with a gaussian wave packet in the radial direction and a cylinder of length $L$ along $z$ yields

[TABLE]

for the data in Fig. 3a, where $n_{f}$ is the density of the filament. From the fit of the parameters above we get $N_{p}=25\pm 7$ and $a_{\perp}/r_{0}=1.7\pm 0.4$.

4.3 Weakly interacting regime

With the many-body Hamiltonian (3) one can also investigate the weakly interacting regime. To do that, we consider a limit where $a_{s}$ and $a_{d}$ length scales are much smaller than the inter-particle distance. This corresponds to maintain the general conditions $n|a_{s}|^{3}\ll 1$ and $na_{d}^{3}\ll 1$ simultaneously fulfilled. Quantitatively, as already discussed in Ref.Cinti et al. (2017), both inequalities hold keeping $n|a_{s}|^{3}$ and $na_{d}^{3}\leq 5\cdot 10^{-3}$.

We present QMC data at $\epsilon_{d}=0.6,\,1.2,\,1.8$. The upper part of Fig. 4 reports snapshot configurations at a rescaled low temperature equal to $T/T_{0}=0.25$, $T_{0}$ being the critical temperature of the ideal Bose gas $k_{B}T_{0}=2\pi(nr_{0}^{3})^{2/3}/\zeta(3/2)^{2/3}$; whereas the bottom panels show their corresponding radial distribution functions $g(r)$. Also in this case, simulations have been obtained with $N=100$, $P=1000$, keeping a fixed density of $nr_{0}^{3}=1\cdot 10^{-4}$.

For $\epsilon_{d}=0.6$ (Fig. 4a) the condensate appears modified by the dipolar strength only marginally. Radial distribution functions display strong correlations at short distances, presenting a typical fluid behavior at large $r$. Upon increasing the dipolar interaction, for $\epsilon_{d}>1$ (e.g. Fig. 4b-c) we observe a stabilization of clusters. As expected, their emergence exhibits a $g(r)$ with stronger correlations at short $r$ if compared to Fig. 4a. Yet, the same observable in Fig. 4b-c shows a gentle modulation at long distances, excluding long range crystalline order at larger $r$, and therefore any cluster crystal phase. Also, contrarily to the prediction of mean field theory for the simulation parameters, this phase is not globally superfluid. In order to illustrate this point, we show the computed the frequency $P(n)$ of the occurrence of permutation cycles involving $n$ particles ($1\leq n\leq N$) for the system with $\epsilon_{d}=1.2$ (see Fig. 5). The important point here is that permutations extend only along cycles that remain confined within the single cluster, supporting the interpretation that locally particles exchange and quantum droplets are locally superfluid. Finally, increasing further $\epsilon_{d}$ (Fig. 4c), we observe nothing but the stabilization of the cluster phase. This completely excludes the existence of the filament regime discussed at higher densities.

5 Conclusion

We discussed the appearance of cluster phases in ultracold dipolar bosons. We reviewed the relevant pairwise interactions among dipoles, leading to the formation of clusters and filaments, as well as the numerical algorithm employed. We characterized the many-body phases varying the strength of the potential as well as the density of the system via extensive QMC calculations.

\funding

T.M. acknowledges for support CAPES through the Project CAPES/Nuffic n. 88887.156521/2017-00 and CNPq through Bolsa de produtividade em Pesquisa n. 311079/2015-6.

\reftitle

References