Next Article in Journal
Quantum Dynamics of Hydrogen-like Ions in a Spatially Nonuniform Magnetic Field: A Possible Application to Fusion Plasma
Next Article in Special Issue
Ultra-Dilute Gas of Polarons in a Bose–Einstein Condensate
Previous Article in Journal
Structure Calculations in Nd III and U III Relevant for Kilonovae Modelling
Previous Article in Special Issue
Pattern Formation in One-Dimensional Polaron Systems and Temporal Orthogonality Catastrophe
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Static Impurities in a Weakly Interacting Bose Gas

by
Galyna Panochko
1 and
Volodymyr Pastukhov
2,*
1
Department of Optoelectronics and Information Technologies, Ivan Franko National University of Lviv, 107 Tarnavskyj Str., 79000 Lviv, Ukraine
2
Department for Theoretical Physics, Ivan Franko National University of Lviv, 12 Drahomanov Str., 79000 Lviv, Ukraine
*
Author to whom correspondence should be addressed.
Atoms 2022, 10(1), 19; https://doi.org/10.3390/atoms10010019
Submission received: 30 September 2021 / Revised: 27 January 2022 / Accepted: 28 January 2022 / Published: 8 February 2022
(This article belongs to the Special Issue Physics of Impurities in Quantum Gases)

Abstract

:
We present a comprehensive discussion of the ground-state properties of dilute D-dimensional Bose gas interacting with a few static impurities. Assuming the short-ranged character of the boson-impurity interaction, we calculated the energy of three- and two-dimensional Bose systems with one and two impurities immersed.

1. Introduction

The problem of impurities in mediums formed by bosons is comprehensively studied in condensed matter physics. Even properties of a single atom immersed in the weakly interacting Bose gas change drastically [1,2,3,4,5]. Depending on the strength of the boson–impurity interaction, a number of physically distinct impurity phases can be realized, namely, the Bose-polaronic state [6,7,8,9] in various spatial dimensions, which is very similar to the free-particle one but with the modified, due to the presence of bath, kinematic characteristics; the molecular state [10,11], when the impurity captures one boson with the formation of a dimer; a set of the Efimov states [12,13,14,15] with the universal scaling behavior of energy levels, and higher-order conglomerates, [16,17,18,19,20,21] which involve a larger number of host atoms. Remarkably, some of these phases can be observed in experiments [22,23]. The experimental progress in the field of ultra-cold atomic gases has recently lead to the observation [24] of Bose polarons at finite temperatures. This experiment confirmed previous theoretical predictions [25,26,27,28,29,30] about the breakdown of the quasi-particle picture description of Bose polarons in a close vicinity of the Bose–Einstein condensation (BEC) point.
Recently, the problem of two impurities immersed in the dilute one and three-dimensional Bose gases has become a subject of extensive examination. Physically, this problem is substantially distinguishable from the single Bose polaron one due to the emergence of the induced effective interaction [31,32,33,34] between impurity particles. In 1D, the character of this interaction crucially depends on a sign of the boson–impurity coupling constant [35]; the effective attraction is found for positive couplings, while the induced repulsive potential is inherent for the negative interactions. While it increases, the induced attractive interaction between impurities leads to the formation of bipolarons [36] in the continuum and on the lattice [37] and even to the emergence of the two-polaron bound states [38]. In one-dimensional geometries with harmonic trapping, the induced interaction causes the clustering [39] of two initially non-interacting atoms and modifies their quench dynamics [40]. By switching the boson–impurity interaction in 3D dilute BEC with two impurities, the transition from weakly interacting through the Yukawa potential bipolarons to the Efimov trimer state was predicted in Ref. [41]. Recently, properties of a single polaron in 2D BEC have been discussed both analytically [42,43] and numerically [44,45]. The arbitrary D one-polaron case was considered in Ref. [46]. As far as we know, the problem of two Bose polarons in 2D Bose gas has never been discussed; therefore, the objective of this study was to make the first step toward the revealing of peculiarities of the bipolaron physics and the boson-induced effective interaction between impurities by considering the static limit. The absence of the impurity dynamics in this limit allows to find the exact solution of the problem in the dilute 1D Bose mediums both in one- [47] and two-particle [48,49] cases. In 3D, only a case of the ideal Bose gas [50,51] is the exactly tractable one, while the presence of a weak boson–boson interaction requires [52] a substantial numerical efforts.

2. Formulation

2.1. Model

The discussed model consists of the D-dimensional (here we focus on D = 2 , 3 cases) Bose gas loaded in volume L D (with the periodic boundary conditions imposed) with the weak interparticle interaction and microscopic number N of heavy (infinite-mass) impurities immersed in it. Heavy particles are supposed to be randomly placed in positions { r j } . In the following, we adopt the imaginary-time path-integral approach with Euclidean action
S = d x ψ * ( x ) τ ε + μ Φ ( r ) ψ ( x ) g B , Λ 2 d x | ψ ( x ) | 4 ,
where x = ( τ , r ) denotes the “position” in D + 1 -dimensional space (and consequently d x = 0 β d τ L D d r ), and the complex field ψ ( x ) is periodic in τ with period β (which is the inverse temperature of the system). We also use the shorthand notations for bosonic dispersion ε = 2 2 2 m and the chemical potential μ that fixes average density n of Bose gas and for the term
Φ ( r ) = 1 j N g I , Λ δ Λ ( r r j ) ,
that describes the interaction between Bose particles and impurities. The δ -like two-body potential is ill-defined in the higher ( D 2 ) dimensions, and therefore, in order to obtain any reasonable results one should adopt some renormalization scheme. The latter is typically realized by the implication of the ultraviolet cutoff Λ in all momentum summations and in the simultaneous rewriting of bare couplings g B , Λ and g I , Λ via the two-body vacuum binding energies ϵ B and ϵ I
g B , Λ 1 = g B 1 1 L D k 1 2 ε k ,
g I , Λ 1 = g I 1 1 L D k 1 ε k ,
respectively, (from now on, we assume that all summations over the wave-vector k are restricted from the above | k | < Λ ). Such a “regularization” is already used in the definition of the point-like boson–impurity interaction potential, δ Λ ( r ) = 1 L D | k | < Λ e i k r , in Equation (2). The ”observable” couplings g B and g I are specified as follows
g B 1 = Γ ( 2 D 2 ) ( 4 π ) D 2 m 2 D 2 | ϵ B | D 2 1 ,
g I 1 = Γ ( 2 D 2 ) ( 2 π ) D 2 m 2 D 2 | ϵ I | D 2 1 ,
where Γ ( z ) stands for the gamma function. Note that the bound states are only possible for positive g B s and g I s, but it is convenient to parameterize negative couplings by the binding energies. By careful inspection of the D 2 limit, one can conclude that Equations (3)–(6) provide a correct description of zero-range potentials even in the two-dimensional case. Moreover, the D = 2 pseudo-potential always provides the existence of one bound state.
The alternative method (see, for instance [53]) to deal with point-like interactions is to initially start from some “physical” (Gaussian, for instance) potentials and then relate the appropriate coupling constant to the s-wave scattering lengths a B and a I in the limit where the effective ranges are the smallest parameters with the dimension of length in the system. In the following, no restrictions are set on a magnitude of the boson–impurity interaction, while the Bose gas itself is expected to be extremely dilute.

2.2. Effective Field Theory Approach

The further analysis will be performed in a spirit of the effective field theory approach (see, for a review [54]), which is known to be extremely convenient for the many-boson systems. Particularly, this formulation automatically guarantees the implementation of the Hugengoltz–Pines theorem (which is a concrete manifestation of the Goldstone theorem) in every order of a loop expansion. Moreover, the effective field theory approach provides a non-perturbative predictions for the Bose gas thermodynamics. In the limit of weak boson–boson coupling, the loop expansion is identical to the perturbation theory in terms of characteristic small parameter a B D n . The main idea of the method relies on the separation of “classical” dynamics during the computations of the partition function by means of the path integral
ψ ( x ) = ψ 0 ( r ) + ψ ˜ ( x ) , ψ * ( x ) = ψ 0 * ( r ) + ψ ˜ * ( x ) ,
where the introduced classical fields are determined by the minimization of the action (1): δ S 0 = δ S [ ψ 0 * , ψ 0 ] = 0 . Note that in general | ψ 0 ( r ) | 2 should not be confused with the Bose condensate density. In the absence of impurities, Φ ( r ) = 0 , the solution ψ 0 ( r ) is real and uniform. Putting a microscopic amount of heavy particles in the Bose condensate, we cannot principally change the character of this solution provided that ψ 0 ( r ) becomes only slightly non-uniform, i.e., L D d r | ψ 0 ( r ) | 2 L D . Of course, one may argue that the localized solutions ψ 0 ( r ) decrease the total energy by N | ϵ I | , but any non-zero repulsion between bosons immediately increases the energy of the system by N 2 g B / a I D . Therefore, the collapsed BEC state [50] is not energetically preferable in the thermodynamic limit, where both the number of the repulsively interacting bosons N and the volume of the box L D infinitely increase.
Performing the shift (7), we end up with the following effective action
S eff = S 0 1 2 d x ψ ˜ * ( x ) , ψ ˜ ( x ) K ^ ψ ˜ ( x ) ψ ˜ * ( x ) ,
where only the Gaussian in the fluctuation fields part is explicitly written down. Here the 2 × 2 matrix operator K ^ with elements
K ^ 11 = ε μ + Φ ( r ) + 2 g B , Λ | ψ 0 ( r ) | 2 τ , K ^ 12 = K ^ 21 * = g B , Λ ψ 0 2 ( r ) , K ^ 22 = ε μ + Φ ( r ) + 2 g B , Λ | ψ 0 ( r ) | 2 + τ .
is introduced. Taking into account the equation for ψ 0 ( r )
ε μ + Φ ( r ) + g B , Λ | ψ 0 ( r ) | 2 ψ 0 ( r ) = 0 ,
and performing the Gaussian integration in (8), we finally obtain the grand potential of the Bose system with the impurities immersed
Ω = g B , Λ 2 L D d r | ψ 0 ( r ) | 4 + 1 2 β Sp ln K ^ const ,
where Sp denotes the trace in the D + 1 space. A constant term (counterterm) in (11) is most straightforwardly represented in the plane-wave basis const = 1 2 k k | ε μ + 2 g B , Λ | ψ 0 ( r ) | 2 + Φ ( r ) | k but cannot be obtained by the functional integration and has to be written by hand [55] in order to resolve a standard normal-ordering routine. Consequently, the calculation of thermodynamics for “Bose gas + static impurities” reduces to finding a solution of Equation (10) and then, with ψ 0 ( r ) in hand, to the evaluation of the functional determinant. Note that by taking into account S 0 only, one reproduces the mean-field [56,57,58,59,60,61] description of the system generalized to N impurities in the static limit. In comparison to other Bose-polaron theories, the adopted approach takes into account the correct short-distance boson-impurity physics from the very beginning of the celebrated Bogoliubov treatment of dilute Bose condensates.

2.3. Limit of Dilute Bose Gas

In the general case, the above program, which can be realized to the very end in 1D [48] even at finite impurity masses [56,58,62], requires considerable numerical efforts in the higher dimensions, but the limit of weak inter-boson interaction can be handled more or less easily. Indeed, the intrinsic, for the dilute Bose gas, length-scale is represented by the coherence length ξ = m c (with c = n g B / m being the sound velocity), which is large in comparison to the average distance between particles and to the boson–boson s-wave scattering length a B . The magnitude of the boson–impurity interaction, in turn, is dictated by the boson–impurity s-wave scattering length a I . So, if we additionally assume that a I ξ , the solution of Equation (10) can be immediately found ψ 0 ( r ) = μ / g B , Λ n . In all other cases, we can apply the successive expansion in terms of the ψ 0 -field “non-uniformity”
ψ 0 ( r ) = μ / g B , Λ 1 ψ ¯ 0 ( 1 ) ( r ) ψ ¯ 0 ( 2 ) ( r ) ,
where after the substitution in Equation (10) the dimensionless functions ψ ¯ 0 ( 1 ) ( r ) , ψ ¯ 0 ( 2 ) ( r ) satisfy the following equations:
ε + 2 μ + Φ ( r ) ψ ¯ 0 ( 1 ) ( r ) = Φ ( r ) ,
ε + 2 μ + Φ ( r ) ψ ¯ 0 ( 2 ) ( r ) = 3 μ ψ ¯ 0 ( 1 ) ( r ) 2 .
Note that the above approximate procedure does not require the boson–impurity interaction to be weak. Furthermore, by a naive dimensional analysis, it is easy to argue that both at the weak and strong couplings g I , the contribution of the second-order correction ψ ¯ 0 ( 2 ) ( r ) in the thermodynamics of the system is much smaller than the one originating from ψ ¯ 0 ( 1 ) ( r ) . Therefore, in our consideration below we fully focus on the first-order correction. However, even this simple approximation effectively sums up some infinite set of terms of the standard pertubation theory for a model with the uniform condensate [47]. Equation (13) with Φ ( r ) given by (2) can be solved for arbitrary N by means of the Fourier transformation
ψ ¯ 0 ( 1 ) ( r ) = 1 j N A j 1 L D k e i k ( r r j ) ε k + 2 μ ,
with ε k = 2 k 2 2 m and coefficients A j = 1 i N T j i ( 2 μ ) , where matrix T j i ( 2 μ ) is introduced in Appendix A.
We can now proceed with the calculations of the functional determinant in (11). Taking into account the extreme diluteness of the Bose subsystem, it is enough to expand Sp ln K ^ Sp ln K ^ ( 0 ) + Sp [ K ^ ( 0 ) ] 1 Δ K ^ , where K ^ ( 0 ) is given by (9) but with ψ 0 ( r ) μ / g B , Λ and Δ K ^ = K ^ K ^ ( 0 ) . Following our previous discussion, we ignore in Δ K ^ all higher-order corrections except ψ ¯ 0 ( 1 ) ( r ) . After this, the calculations are relatively simple, and at absolute zero we obtain the Ω -potential in the adopted approximation
Ω L D μ 2 2 g B , Λ + μ g B , Λ 1 j N A j + 1 2 k k | E ε μ Φ ( r ) | k + 1 L D k 1 ε k + μ / 2 E k 1 j N A j ,
where E = ( ε + Φ ( r ) ) 2 + 2 μ ( ε + Φ ( r ) ) and E k = ε k 2 + 2 μ ε k stands for the Bogoliubov spectrum of the “pure” Bose system. It should be noted that for dilute Bose systems the impact of the quantum fluctuations (terms with the summations over the wave-vector) to Ω is much smaller than the first two terms (the mean-field contributions). The last step to be performed in these calculations is to replace the bare couplings g B , Λ and g I , Λ via (3) and (4), respectively. This procedure provides the convergence of sums over the wave-vector in the last two terms of (16). Then, the trace in the third term of (16) can be computed (see Appendix A for details). With the well-defined grand potential, we can relate, by using the thermodynamic identity n = μ Ω L D , the chemical potential of the Bose system to its equilibrium density n. Performing these calculations, one must keep in mind that the presence of a microscopic number of impurities cannot principally change the properties of the system. So, if we denote (and appropriate grand potential Ω B ) the chemical potential of Bose gas without exterior particles by μ B , the difference Δ μ = μ μ B N / L D should be small. Using this fact and n = μ Ω B L D μ Δ Ω L D , we can identify a small correction Δ μ = Δ Ω μ B / 2 Ω B μ B 2 . The latter formula allows to determine the energy that the Bose system gains when N impurities are immersed
Δ E N = Ω Ω B | μ μ B ,
which is an explicit manifestation of the well-known theorem about small corrections to the thermodynamic potentials.

3. Results

Before we proceed to describing our main results, it is necessary to analyze the case of “pure” bosons. Setting Φ ( r ) = 0 in (16) and calculating integrals, we obtain for density
n = μ B g B 1 Γ ( D ) D 2 Γ 2 ( D 2 ) μ B | ϵ B | D 2 1 ,
which allows to obtain the expression for μ B iteratively. For the weakly non-ideal three-dimensional bosons we find the well-known formula ( | ϵ B | = 2 m a B 2 in 3D)
μ B = 4 π 2 a B n m 1 + 32 3 π n a B 3 + .
Similarly, in the two-dimensional case we have the transcendental equation [63]
n = m μ B 4 π 2 ln | ϵ B | μ B 1 .
Being convinced that the limit of Bose gas without impurities is correctly reproduced by the adopted approach, we are ready to present our main results concerning the binding energy of one- and two-impurity atoms in the dilute three- and two-dimensional Bose gases.

3.1. 3D Case

In the 3D case, the general structure of the two-impurity binding energy in the dilute Bose gas ( n ξ 3 1 ) can be represented as
Δ E 2 = Δ E 2 ( 0 ) ε 1 a I ξ ; R ξ + 1 n ξ 3 ε 2 a I ξ ; R ξ + ,
where Δ E 2 ( 0 ) = 2 g I n is the contribution of the ideal Bose gas; a I is the s-wave scattering length that parameterizes the (renormalized) two-body coupling g I = 2 π 2 a I m ; and R is the distance between two static particles. The first term in (21) has a simple analytic form
ε 1 a I ξ ; R ξ = ξ / a I ξ / a I 2 + e 2 R / ξ / ( R / ξ ) ,
and originates purely from the mean-field correction to the grand potential (the second term in (16)), while ε 2 a I ξ ; R ξ contains both the mean-field and purely quantum corrections. Note that in the formula for Ω only the one-loop corrections were taken into account, and a consistent consideration of the next to a leading order terms in series expansion over the small parameter 1 / ( n ξ 3 ) necessarily requires the calculation of the two-loop diagrams to the grand potential. By setting the distance between heavy particles R to infinity, one obtains from (21) the one-impurity limit. A typical behavior of functions ε 1 , 2 a I ξ ; is presented in Figure 1.
Let us recall that the problem considered here is the exactly solvable one, when the bosons are non-interacting. Therefore, it should be clearly understood that the presented results are accurate if the coherence length ξ is the largest parameter with the dimension of the length in the system. In order to reveal the interplay between regimes of very dilute a I / ξ 0 Bose gas and intermediate boson–impurity interaction, we plotted in Figure 2 the binding energy of two heavy particles for the positive and negative s-wave scattering lengths a I .
Comparing these findings to the ideal Bose gas results [50], we can observe similar patterns in the behavior of the systems at weak coupling: at positive a I , the binding energy is the monotonic function of R, while at the negative boson–impurity scattering lengths both, ε 1 , 2 0.01 ; R ξ have a simple-pole singularity. When the interaction increases (see right panel in Figure 2) the mean-field and quantum corrections to the ground-state energy of 3D Bose gas possess an infinite discontinuities independently of a sign of a I .

3.2. 2D Case

In general, the low-dimensional dilute Bose systems with static impurities are very peculiar. When the interaction between bosons is switched off, these systems are insensible to the boson–impurity interaction in their un-collapsed ground state, and therefore, the binding energy of the heavy particles requires a finite compressibility of the host system to be non-zero. This is a general result for the low-dimensional (1D and 2D) ideal Bose gases with impurities that is independent of the approximations made. Introducing the two-body s-wave scattering length a I through the boson–impurity vacuum bound state energy | ϵ I | = 2 e 2 γ 2 / ( m a I 2 ) , we can write down the energy that the 2D Bose gas gains when two heavy particles are immersed in it
Δ E 2 = 2 2 π 2 n m ε 1 a I ξ ; R ξ + 1 n ξ 2 ε 2 a I ξ ; R ξ + .
As in the 3D case, the mean-field correction casts into a simple analytic form
ε 1 a I ξ ; R ξ = 1 ln e 2 γ ξ 2 a I 2 2 K 0 2 R ξ ,
K 0 ( z ) is the modified Bessel function of the second kind [64]. Note that in contrast to a 3D case, ε 2 a I ξ ; R ξ tends to zero (at least logarithmically) in the limit of ideal Bose gas ( ξ ). At large distances R, Equation (23) gives the double binding energy of a single impurity, which is presented in Figure 3.
Particularly, these calculations clearly demonstrate the weakening of the role of quantum fluctuations in the formation of polarons in two-dimensional Bose systems. Actually, this observation [62] seems to be intrinsic for the low-dimensional systems in general.
The numerical computations of the two-impurity energies (see Figure 4).
The 2D Bose gas demonstrates qualitative similarity between the two- and three-dimensional cases. At weak boson–impurity interactions a I / ξ 1 , where our effective field-theoretical formulation is supposed to make a quantitative predictions, the mean-field term ε 1 a I ξ ; R ξ as well as the one that includes the quantum corrections ε 2 a I ξ ; R ξ behave as monotonic functions of R. The interaction-induced effective two-body potential between static particles at large a I / ξ always contains a singularity.

4. Conclusions

In summary, by means of the effective field theory formulation, we calculated the impurity-induced shifts to the ground-state energies of the two- and three-dimensional dilute Bose gases. Particularly, by taking into account the extreme diluteness of the host bosons, we proposed the approximate procedure that allows to calculate the properties of an arbitrary (microscopic) number of static impurities in terms of a characteristic small parameter 1 / ( n ξ D ) (where n and ξ are the density and the coherence length of bosons, respectively). The numerical calculations of the binding energies of two static impurities in dilute 2D and 3D Bose gases that were performed for a wide range of the boson–impurity interactions and distances between impurities has revealed the peculiarities of the medium-induced (Casimir) forces: (i) the two-body effective potential always demonstrates singular behavior at the distances between impurities comparable to the boson–impurity s-wave scattering lengths a I ; (ii) an impact of purely quantum corrections decreases with the lowering of a spatial dimensionality. Similar singularities are also intrinsic for the binding energy of a single impurity at a I ξ , which may signal [65] about the inapplicability of the adopted approximate treatment for calculations of the “classical” solution ψ 0 ( r ) in that region, where the full numerical solution to Equation (10) is required.

Author Contributions

Conceptualization, V.P.; methodology, G.P. and V.P.; software, G.P.; validation, G.P. and V.P.; formal analysis, G.P. and V.P.; investigation, G.P. and V.P.; resources, G.P. and V.P.; data curation, G.P. and V.P.; writing—original draft preparation, G.P. and V.P.; visualization, G.P.; supervision, V.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this study.

Acknowledgments

We are indebted to Iryna Pastukhova for comments on the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

For completeness, in this section we give some details of the calculations not presented in the main text. Let us first start from the equation that determines the classical field ψ 0 ( r ) . Explicitly writing down Equation (13), after the implementation of ansatz (15)
1 j N A j δ Λ ( r r j ) + 1 j N g I , Λ δ Λ ( r r j ) 1 i N A i 1 L D k e i k ( r r i ) ε k + 2 μ = 1 j N g I , Λ δ Λ ( r r j ) ,
and combining j = i terms in double sum with the first term of equation, we obtain
A j 1 g I , Λ + 1 L D k 1 ε k + 2 μ + 1 i j N 1 L D k e i k ( r j r i ) ε k + 2 μ A i = 1 .
The divergent sum in the square brackets is now regularized by the renormalization of a coupling constant (4), so the final result contains only observable g I . One can easily recognize the square brackets as the boson–impurity two-body T-matrix
t I 1 ( ω ) = g I , Λ 1 1 L D k 1 ω ε k ,
and introducing auxiliary notations
Δ i j ( ω ) = 1 L D k e i k ( r i r j ) ω ε k ,
we find the result for coefficients A j announced in the main text
A i = 1 j N T i j ( 2 μ ) , T i j 1 ( 2 μ ) = δ i j t I 1 ( 2 μ ) Δ i j ( 2 μ ) ( 1 δ i j ) .
For the calculation of the trace in the second term of (16), we used formal identity
k k | E ε μ Φ ( r ) | k = d ω D ( ω ) ω 2 + 2 μ ω ω μ , D ( ω ) = k k | δ ( ω ε Φ ( r ) | k .
The density of states D ( ω ) is easily calculated within the Green’s function method [50]
D ( ω ) = k δ ( ω ε k ) 1 π Im k | T ( ω + i 0 ) | k ( ω + i 0 ε k ) 2 ,
where the T-matrix T ( ω ) characterizes the scattering of a single boson on N impurities
q | T ( ω ) | k = 1 i , j N e i q r i T i j ( ω ) e i k r j .
The calculations of k | T ( ω + i 0 ) | k in the density of states requires the knowledge of an explicit analytic formulas for the boson–impurity two-body T-matrix
t I 1 ( ω ) = Γ ( 2 D 2 ) ( 2 π ) D 2 m 2 D 2 ( ω ) D 2 1 | ϵ I | D 2 1 ,
and a function Δ i j ( ω ) = Δ R ( ω ) of distance R = | r i r j | between two impurities in arbitrary D
Δ R ( ω ) = 1 ( 2 π ) D 2 2 m k ω D 2 2 K D 2 1 ( R k ω ) R k ω D 2 1 ,
where k ω = 2 m ( ω ) / , and K ν ( z ) is the modified Bessel function of the second kind [64].

References

  1. Tempere, J.; Casteels, W.; Oberthaler, M.K.; Knoop, S.; Timmermans, E.; Devreese, J.T. Feynman path-integral treatment of the BEC-impurity polaron. Phys. Rev. B 2009, 80, 184504. [Google Scholar] [CrossRef] [Green Version]
  2. Vlietinck, J.; Casteels, W.; Houcke, K.V.; Tempere, J.; Ryckebusch, J.; Devreese, J.T. Diagrammatic Monte Carlo study of the acoustic and the Bose–Einstein condensate polaron. New J. Phys. 2015, 17, 033023. [Google Scholar] [CrossRef] [Green Version]
  3. Ardila, L.A.P.; Giorgini, S. Impurity in a Bose-Einstein condensate: Study of the attractive and repulsive branch using quantum Monte Carlo methods. Phys. Rev. A 2015, 92, 033612. [Google Scholar] [CrossRef] [Green Version]
  4. Grusdt, F.; Shchadilova, Y.E.; Rubtsov, A.N.; Demler, E. Renormalization group approach to the Frohlich polaron model: Application to impurity-BEC problem. Sci. Rep. 2015, 5, 12124. [Google Scholar] [CrossRef] [PubMed]
  5. Panochko, G.; Pastukhov, V.; Vakarchuk, I. Behavior of the impurity atom in a weakly-interacting Bose gas. Cond. Matt. Phys. 2017, 20, 13604. [Google Scholar] [CrossRef]
  6. Astrakharchik, G.E.; Pitaevskii, L.P. Motion of a heavy impurity through a Bose-Einstein condensate. Phys. Rev. A 2004, 70, 013608. [Google Scholar] [CrossRef] [Green Version]
  7. Novikov, A.; Ovchinnikov, M. A diagrammatic calculation of the energy spectrum of quantum impurity in degenerate Bose–Einstein condensate. J. Phys. A Math. Theor. 2009, 42, 135301. [Google Scholar] [CrossRef]
  8. Christensen, R.S.; Levinsen, J.; Bruun, G.M. Quasiparticle Properties of a Mobile Impurity in a Bose-Einstein Condensate. Phys. Rev. Lett. 2015, 115, 160401. [Google Scholar] [CrossRef] [Green Version]
  9. Panochko, G.; Pastukhov, V.; Vakarchuk, I. Impurity self-energy in the strongly-correlated Bose systems. Int. J. Mod. Phys. B 2018, 32, 1850053. [Google Scholar] [CrossRef]
  10. Rath, S.P.; Schmidt, R. Field-theoretical study of the Bose polaron. Phys. Rev. A 2013, 88, 053632. [Google Scholar] [CrossRef] [Green Version]
  11. Li, W.; Sarma, S.D. Variational study of polarons in Bose-Einstein condensate. Phys. Rev. A 2014, 90, 013618. [Google Scholar] [CrossRef] [Green Version]
  12. Levinsen, J.; Massignan, P.; Parish, M.M. Efimov Trimers under Strong Confinement. Phys. Rev. X 2014, 4, 031020. [Google Scholar] [CrossRef] [Green Version]
  13. Naidon, P.; Endo, S. Efimov physics: A review. Rep. Prog. Phys. 2017, 80, 056001. [Google Scholar] [CrossRef] [Green Version]
  14. Levinsen, J.; Parish, M.M.; Bruun, G.M. Impurity in a Bose-Einstein Condensate and the Efimov Effect. Phys. Rev. Lett. 2015, 115, 125302. [Google Scholar] [CrossRef] [PubMed]
  15. Sun, M.; Zhai, H.; Cui, X. Visualizing the Efimov Correlation in Bose Polarons. Phys. Rev. Lett. 2017, 119, 013401. [Google Scholar] [CrossRef] [PubMed]
  16. Wang, Y.; Laing, W.B.; Stecher, J.; Esry, B.D. Efimov Physics in Heteronuclear Four-Body Systems. Phys. Rev. Lett. 2012, 108, 073201. [Google Scholar] [CrossRef] [PubMed]
  17. Casteels, W.; Tempere, J.; Devreese, J.T. Bipolarons and multipolarons consisting of impurity atoms in a Bose-Einstein condensate. Phys. Rev. A 2013, 88, 013613. [Google Scholar] [CrossRef] [Green Version]
  18. Blume, D.; Yan, Y. Generalized Efimov Scenario for Heavy-Light Mixtures. Phys. Rev. Lett. 2014, 113, 213201. [Google Scholar] [CrossRef] [Green Version]
  19. Shi, Z.Y.; Yoshida, S.M.; Parish, M.M.; Levinsen, J. Impurity-Induced Multibody Resonances in a Bose Gas. Phys. Rev. Lett. 2018, 121, 243401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Yoshida, S.M.; Shi, Z.Y.; Levinsen, J.; Parish, M.M. Few-body states of bosons interacting with a heavy quantum impurity. Phys. Rev. A 2018, 98, 062705. [Google Scholar] [CrossRef] [Green Version]
  21. Blume, D. Few-boson system with a single impurity: Universal bound states tied to Efimov trimers. Phys. Rev. A 2019, 99, 013613. [Google Scholar] [CrossRef] [Green Version]
  22. Jorgensen, N.B.; Wacker, L.; Skalmstang, K.T.; Parish, M.M.; Levinsen, J.; Christensen, R.S.; Bruun, G.M.; Arlt, J.J. Observation of Attractive and Repulsive Polarons in a Bose-Einstein Condensate. Phys. Rev. Lett. 2016, 117, 055302. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Hu, M.-G.; de Graaff, M.J.V.; Kedar, D.; Corson, J.P.; Cornell, E.A.; Jin, D.S. Bose Polarons in the Strongly Interacting Regime. Phys. Rev. Lett. 2016, 117, 055301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Yan, Z.Z.; Ni, Y.; Robens, C.; Zwielein, M.W. Bose polarons near quantum criticality. Science 2020, 368, 190. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Levinsen, J.; Parish, M.M.; Christensen, R.S.; Arlt, J.J.; Bruun, G.M. Finite-temperature behavior of the Bose polaron. Phys. Rev. A 2017, 96, 063622. [Google Scholar] [CrossRef] [Green Version]
  26. Guenther, N.-E.; Massignan, P.; Lewenstein, M.; Bruun, G.M. Bose Polarons at Finite Temperature and Strong Coupling. Phys. Rev. Lett. 2018, 120, 050405. [Google Scholar] [CrossRef] [Green Version]
  27. Pastukhov, V. Polaron in the dilute critical Bose condensate. J. Phys. A Math. Theor. 2018, 51, 195003. [Google Scholar] [CrossRef] [Green Version]
  28. Liu, W.E.; Levinsen, J.; Paris, M.M. Variational Approach for Impurity Dynamics at Finite Temperature. Phys. Rev. Lett. 2019, 122, 205301. [Google Scholar] [CrossRef] [Green Version]
  29. Field, B.; Levinsen, J.; Parish, M.M. Fate of the Bose polaron at finite temperature. Phys. Rev. A 2020, 101, 013623. [Google Scholar] [CrossRef] [Green Version]
  30. Pascual, G.; Boronat, J. Quasiparticle Nature of the Bose Polaron at Finite Temperature. Phys. Rev. Lett. 2021, 127, 205301. [Google Scholar] [CrossRef]
  31. Zinner, N.T. Efimov states of heavy impurities in a Bose-Einstein condensate. Europhys. Phys. Lett. 2013, 101, 60009. [Google Scholar] [CrossRef]
  32. Zinner, N.T. Spectral flow of trimer states of two heavy impurities and one light condensed boson. Europhys. Phys. J. D 2014, 68, 216. [Google Scholar] [CrossRef] [Green Version]
  33. Camacho-Guardian, A.; Ardila, L.A.P.; Pohl, T.; Bruun, G.M. Bipolarons in a Bose-Einstein Condensate. Phys. Rev. Lett. 2018, 121, 013401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Camacho-Guardian, A.; Bruun, G.M. Landau Effective Interaction between Quasiparticles in a Bose-Einstein Condensate. Phys. Rev. X 2018, 8, 031042. [Google Scholar] [CrossRef] [Green Version]
  35. Brauneis, F.; Hammer, H.-W.; Lemeshko, M.; Volosniev, A.G. Impurities in a one-dimensional Bose gas: The flow equation approach. SciPost Phys. 2021, 11, 008. [Google Scholar] [CrossRef]
  36. Petcovich, A.; Ristivojevic, Z. Mediated interaction between polarons in a one-dimensional Bose gas. arXiv 2021, arXiv:2103.08772. [Google Scholar]
  37. Pasek, M.; Orso, G. Induced pairing of fermionic impurities in a one-dimensional strongly correlated Bose gas. Phys. Rev. B 2019, 100, 245419. [Google Scholar] [CrossRef] [Green Version]
  38. Will, M.; Astrakharchik, G.E.; Fleischhauer, M. Polaron Interactions and Bipolarons in One-Dimensional Bose Gases in the Strong Coupling Regime. Phys. Rev. Lett. 2021, 127, 103401. [Google Scholar] [CrossRef]
  39. Dehkharghani, A.S.; Volosniev, A.G.; Zinner, N.T. Coalescence of Two Impurities in a Trapped One-dimensional Bose Gas. Phys. Rev. Lett. 2018, 121, 080405. [Google Scholar] [CrossRef] [Green Version]
  40. Mistakidis, S.I.; Volosniev, A.G.; Schmelcher, P. Induced correlations between impurities in a one-dimensional quenched Bose gas. Phys. Rev. Res. 2020, 2, 023154. [Google Scholar] [CrossRef]
  41. Naidon, P. Two Impurities in a Bose–Einstein Condensate: From Yukawa to Efimov Attracted Polarons. J. Phys. Soc. Jpn. 2018, 87, 043002. [Google Scholar] [CrossRef]
  42. Pastukhov, V. Polaron in dilute 2D Bose gas at low temperatures. J. Phys. B At. Mol. Opt. Phys. 2018, 51, 155203. [Google Scholar] [CrossRef] [Green Version]
  43. Isaule, F.; Morera, I.; Massignan, P.; Juliá-Díaz, B. Renormalization-group study of Bose polarons. Phys. Rev. A 2021, 104, 023317. [Google Scholar] [CrossRef]
  44. Akaturk, E.; Tanatary, B. Two-dimensional Bose polaron using diffusion Monte Carlo method. Int. J. Mod. Phys. B 2019, 33, 1950238. [Google Scholar] [CrossRef] [Green Version]
  45. Ardila, L.A.P.; Astrakharchik, G.E.; Giorgini, S. Strong coupling Bose polarons in a two-dimensional gas. Phys. Rev. Res. 2020, 2, 023405. [Google Scholar] [CrossRef]
  46. Khan, M.M.; Tercas, H.; Mendonca, J.T.; Wehr, J.; Charalambous, C.; Lewenstein, M.; Garcia-March, M.A. Quantum dynamics of a Bose polaron in a d-dimensional Bose-Einstein condensate. Phys. Rev. A 2021, 103, 023303. [Google Scholar] [CrossRef]
  47. Kain, B.; Ling, H.Y. Analytical study of static beyond-Frohlich Bose polarons in one dimension. Phys. Rev. A 2018, 98, 033610. [Google Scholar] [CrossRef] [Green Version]
  48. Reichert, B.; Ristivojevic, Z.; Petkovic, A. The Casimir-like effect in a one-dimensional Bose gas. New J. Phys. 2019, 21, 053024. [Google Scholar] [CrossRef]
  49. Reichert, B.; Ristivojevic, Z.; Petkovic, A. Field-theoretical approach to the Casimir-like interaction in a one-dimensional Bose gas. Phys. Rev. B 2019, 99, 205414. [Google Scholar] [CrossRef] [Green Version]
  50. Panochko, G.; Pastukhov, V. Two- and three-body effective potentials between impurities in ideal BEC. J. Phys A Math. Theor. 2021, 54, 085001. [Google Scholar] [CrossRef]
  51. Drescher, M.; Salmhofer, M.; Enss, T. Quench Dynamics of the Ideal Bose Polaron at Zero and Nonzero Temperatures. Phys. Rev. A 2021, 103, 033317. [Google Scholar] [CrossRef]
  52. Levinsen, J.; Ardila, L.A.P.; Yoshida, S.M.; Parish, M.M. Quantum Behavior of a Heavy Impurity Strongly Coupled to a Bose Gas. Phys. Rev. Lett. 2021, 127, 033401. [Google Scholar] [CrossRef]
  53. Volosniev, A.G.; Hammer, H.-W.; Zinner, N.T. Real-time dynamics of an impurity in an ideal Bose gas in a trap. Phys. Rev. A 2015, 92, 023623. [Google Scholar] [CrossRef] [Green Version]
  54. Andersen, J.O. Theory of the weakly interacting Bose gas. Rev. Mod. Phys. 2004, 76, 599. [Google Scholar] [CrossRef] [Green Version]
  55. Salasnich, L.; Toigo, F. Zero-point energy of ultracold atoms. Phys. Rep. 2016, 640, 1. [Google Scholar] [CrossRef] [Green Version]
  56. Volosniev, A.G.; Hammer, H.-W. Analytical approach to the Bose-polaron problem in one dimension. Phys. Rev. A 2017, 96, 031601(R). [Google Scholar] [CrossRef] [Green Version]
  57. Pastukhov, V. Mean-field properties of impurity in Bose gas with three-body forces. Phys. Lett. A 2019, 383, 2610. [Google Scholar] [CrossRef] [Green Version]
  58. Panochko, G.; Pastukhov, V. Mean-field construction for spectrum of one-dimensional Bose polaron. Ann. Phys. 2019, 409, 167933. [Google Scholar] [CrossRef] [Green Version]
  59. Hryhorchak, O.; Panochko, G.; Pastukhov, V. Mean-field study of repulsive 2D and 3D Bose polarons. J. Phys. B At. Mol. Opt. Phys. 2020, 53, 205302. [Google Scholar] [CrossRef]
  60. Hryhorchak, O.; Panochko, G.; Pastukhov, V. Impurity in a three-dimensional unitary Bose gas. Phys. Lett. A 2020, 384, 126934. [Google Scholar] [CrossRef]
  61. Massignan, P.; Yegovtsev, N.; Gurarie, V. Universal Aspects of a Strongly Interacting Impurity in a Dilute Bose Condensate. Phys. Rev. Lett. 2021, 126, 123403. [Google Scholar] [CrossRef] [PubMed]
  62. Jager, J.; Barnett, R.; Will, M.; Fleischhauer, M. Strong-coupling Bose polarons in one dimension: Condensate deformation and modified Bogoliubov phonons. Phys. Rev. Res. 2020, 2, 033142. [Google Scholar] [CrossRef]
  63. Mora, C.; Castin, Y. Ground State Energy of the Two-Dimensional Weakly Interacting Bose Gas: First Correction Beyond Bogoliubov Theory. Phys. Rev. Lett. 2009, 102, 180404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Abramowitz, M.; Stegun, I. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables; US Government Printing Office: Washington, DC, USA, 1964.
  65. Schmidt, R.; Enss, T. Self-stabilized Bose polarons. arXiv 2021, arXiv:2102.13616. [Google Scholar]
Figure 1. Dimensionless functions ε 1 , 2 a I ξ ; determining the one-impurity energy in 3D dilute Bose gas.
Figure 1. Dimensionless functions ε 1 , 2 a I ξ ; determining the one-impurity energy in 3D dilute Bose gas.
Atoms 10 00019 g001
Figure 2. Mean-field and the first-order quantum corrections ε 1 , 2 a I ξ ; R ξ to the energy of 3D dilute Bose gas generated by two impurities for a I ξ = ± 0.01 and a I ξ = ± 1 .
Figure 2. Mean-field and the first-order quantum corrections ε 1 , 2 a I ξ ; R ξ to the energy of 3D dilute Bose gas generated by two impurities for a I ξ = ± 0.01 and a I ξ = ± 1 .
Atoms 10 00019 g002
Figure 3. Dimensioless one-impurity binding energy terms ε 1 , 2 a I ξ ; (see Equation (23)) in 2D case.
Figure 3. Dimensioless one-impurity binding energy terms ε 1 , 2 a I ξ ; (see Equation (23)) in 2D case.
Atoms 10 00019 g003
Figure 4. The two-impurity dimensionless binding energy corrections ε 1 , 2 a I ξ ; R ξ in 2D dilute Bose gas.
Figure 4. The two-impurity dimensionless binding energy corrections ε 1 , 2 a I ξ ; R ξ in 2D dilute Bose gas.
Atoms 10 00019 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Panochko, G.; Pastukhov, V. Static Impurities in a Weakly Interacting Bose Gas. Atoms 2022, 10, 19. https://doi.org/10.3390/atoms10010019

AMA Style

Panochko G, Pastukhov V. Static Impurities in a Weakly Interacting Bose Gas. Atoms. 2022; 10(1):19. https://doi.org/10.3390/atoms10010019

Chicago/Turabian Style

Panochko, Galyna, and Volodymyr Pastukhov. 2022. "Static Impurities in a Weakly Interacting Bose Gas" Atoms 10, no. 1: 19. https://doi.org/10.3390/atoms10010019

APA Style

Panochko, G., & Pastukhov, V. (2022). Static Impurities in a Weakly Interacting Bose Gas. Atoms, 10(1), 19. https://doi.org/10.3390/atoms10010019

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop