arxiv: v1 [cond-mat.quant-gas] 12 Feb 2018 Dilute Bose gas in classical environment at low temperatures

Transcription

1 arxiv: v [cond-mat.quant-gas] Feb 08 Dilute Bose gas in classical environment at low temperatures Volodymyr Pastukhov Ivan Franko National University of Lviv, Department for Theoretical Physics Drahomanov St., Lviv, Ukraine Abstract The properties of a dilute Bose gas with the non-gaussian quenched disorder are analysed. Being more specific we have considered a system of bosons immersed in the classical bath consisting of the non-interacting particles with infinite mass. Making use of perturbation theory up to second order we have studied the impact of environment on the ground-state thermodynamic and superfluid characteristics of the Bose component. Introduction The properties of Bose gas with quenched disorder was studied extensively during last two decades. This rise of an interest to such a system was stimulated by the possibility to observe Bose glass state transition [] where superfluidity disappears [] even at very low temperatures. The first attempts [3] for theoretical description within the approximate second-quantization method showed that at absolute zero the superfluid density totally depletes when the disorder strength reaches some critical value while the condensate never vanishes. The further developments [4], particularly the non-perturbative extensions on case of an arbitrary two-body coupling strength [5] and strong external potential [6, 7] generally confirm these findings. Diffusive Monte Carlo simulations [8] also agree with the Bogoliubov-like result in the dilute limit, but the increase of disorder potential makes the differences more visible. The finite-temperature phase diagram of the system was clarified extensively in [9, 0, ]. A shift of the Bose-Einstein condensation transition temperature was determined in [, 3] and the critical parameters of the Bose gas in a disordered medium was calculated [4] using quantum Monte Carlo methods. No less interesting is the structure of quasiparticle spectrum and damping for a Bose system in the weak random external potential [5, 6, 7, 8]. In particular, it was shown that the presence of disorder broadens the phonon peak of the dynamic structure factor which for a case of liquid 4 He with the randomly distributed static impurities volodyapastukhov@gmail.com

2 was computed [9] by means of Path Integral Monte Carlo method and experimentally measured in Ref. [0]. Another consequence of disorder is the necessity to correct the exact universal identities which are characteristic for many-boson systems like the Hugenholtz-Pines theorem [] and Josephson s relation []. Practically in all available theoretical studies concerning Bose systems with disorder the distribution of external potential is assumed to be the Gaussian one. However, in order to obtain some reasonable results even for a simplest model of the Bose gas with randomly located infinite-mass impurities the increase of boson-impurity interaction should be accompanied by an account for the distribution non-gaussianity. Formulation of problem We consider the system of N interacting bosons immersed in the bath formed by N non-interacting classical (infinite mass) particles. This model is described by the following Euclidean action { S = dxψ (x) τ + } +µ gρ(r) ψ(x) g dx ψ(x) 4, () m where complex field ψ(x) describe bosonic degrees of freedom, µ is the chemical potential that fixes the Bose gas density. The integration over x = (τ,r) is carried out in a (3+)-domain of volume βv (β is the inverse temperature) with periodicboundaryconditions. The quantityρ(r) = j N δ(r r j) represents the density of homogeneously distributed classical particles. Both boson-boson and boson-impurity two-body interactions are assumed to be short-ranged that are characterised by the coupling constants g and g, respectively. The latter should be related to the appropriate s-wave scattering lengths in the end of calculations g = t V k ε k, g = t V k ε k, () where t = 4π a/m, t = π ã/m and ε k = k /m is the free-particle dispersion. Introducing phase-density representation for bosonic fields ψ (x) = n(x)e iϕ(x), ψ(x) = n(x)e iϕ(x) and making use of the Fourier transformation for n(x) and ϕ(x) n(x) = n+ e ikx n K, ϕ(x) = e ikx ϕ K, (3) βv βv K where K = (ω k,k) stands for the bosonic Matsubara frequency ω k and threedimensional wave-vector k, as well as for the density of classical component ρ(r) = ρ+ e ikr ρ k, (4) V K

3 whereρ = N/V isaveragedensity ofbath particlesandρ k = V j N eikrj, one rewrites action () in the following way: The first term describes Bose gas itself [3] S B = βvµn βvgn + 3! βv S = S B +S d. (5) K +nε k ϕ K ϕ K + K+Q+P =0 + βv K,Q { ω k ϕ K n K ω k ϕ K n K [ εk n +g ] n K n K } 4n (ε k +ε q +ε p )n K n Q n P while S d takes into account the presence of environment m kqϕ Kϕ Q n K Q +..., (6) S d = βv gnρ β g K δ ωk,0ρ k n K. (7) The thermodynamic relation Ω/ µ = N for the grand potential together with explicit form of Eq. (6) fix n = N/V [4, 5] to be the density of the Bose system. This observation allows to proceed in the canonical ensemble and in order to obtain physically meaningful result the averaging over the positions of classical particles should be performed for the free energy of our system. It particularly means that first we have to calculate the free energy F of Bose gas in the presence of local external potential gρ(r) and then identify the Helmholtz potential of the system Bose gas + classical bath with F = V N V dr... V dr N F. (8) In the following we will assume that bosons are weakly-coupled to each other and study the impact of the classical component on the properties of Bose gas in the zero-temperature limit. It is well-known that presence of disorder or interaction with other quantum systems depletes the superfluid density of the Bose gas even at absolute zero [6, 7, 8, 9]. Of course, this phenomena is observed in our case too, and in order to calculate the normal density of superfluid we have to assume that the Bose subsystem moves as a whole with velocity v. Simple analysis [30] shows that account of this motion only leads to the shift ω k ω k i vk of Matsubara frequency in action S B. Then the general consideration in the spirit of perturbation theory in terms of g leads to Ē v = E B +N mv +V gnρ g ρ k ρ k n K n K ωk =0 k + g 3 ρ k ρ q ρ p β nk n Q n P ωk =ω 3! q=ω p=0 +..., (9) k+q+p=0 3

4 for the ground-state energy of system, where E B is the contribution of the Bose gas in rest alone and n K n K, n K n Q n P,... denote the irreducible density correlation functions of pure bosons moving with velocity v. The structure factors of classical non-interacting particles are fully determined by their density and can be easily evaluated by using the procedure described above: ρ k ρ k = ρ, ρ k ρ q ρ p = ρδ k+q+p,0 / V, etc. Expanding r.h.s. of Eq. (9) in powers of velocity Ē v = Ēv=0 +V mv n s +o(v ), (0) (actually in powers of dimensionless parameter v/c, where c is the sound velocity) to quadratic order one obtains the density n s of superfluid component. In addition to the depletion of superfluid density the interaction with bath also decreases the number of Bose particles with zero momentum. To calculate the condensate density of Bose gas we use the following prescription: first within variational differentiation ( ) δ Ē N k =, () δε k n we determine the distribution function of particles with non-zero momentum and then obtain the condensation fraction n 0 n = N k. () N To this stage our consideration is formally exact and the problem is actually reduced to the calculation of irreducible density correlators of a pure Bose gas. But even in the Bogoliubov approximation these calculations are very cumbersome, therefore below we restrict ourselves to the case of weak interspecies interaction, where the perturbation theory in terms of g can be used. 3 Perturbation theory By treating the Bose subsystem on the basis of Bogoliubov s theory we greatly simplify the analysis below but at that time restrict it to the consideration of dilute gases. This is exactly the situation realized in experiments with cold alkali atoms. From the point of view of further calculations in the dilute limit we are free to drop the beyond Bogoliubov corrections to the various density correlation functions. Additionally, working with the same accuracy one should treat the ground-stateenergyofpure bosonse B = E LHY on level ofthe Lee-Huang-Yang [3] formula. Following the above-mentioned approximation scheme in the first order of perturbation theory we have to neglect the last term in Eq. (9) and substitute the pair density correlation function n K n K = nε k +(ω k i kv), (3) 4

5 with E k = ε k +ngε k being Bogoliubov s spectrum. The resulting formula has to be used for obtaining the particle distribution N k and condensate density n 0 n = nb 0 n V ρ t ε k 4, (4) (n B 0 /n is the Bogoliubov result for pure bosons) in the adopted approximation, and after renormalization of the coupling constant () for the explicit evaluation of the energy correction Ē () v=0 N = ρ t [ ρ t εk V ], (5) ε k and depletion of the superfluid component n s n = 4 3V ρ t ε k 4. (6) From general principles as well as from the above formulas it is clear that in the first-order approximation the results derived for our system are identical to those obtained for a dilute Bose gas with weak disorder [8]. The differences appear in the next orders of the formulated perturbation theory. The second order calculations require the account of the last term in Eq. (9). The appropriate three-point density correlator reads n K n Q n P = δ [ K+Q+P,0 βv m kq n Kϕ K n Q ϕ Q n P n P + ε k 4n n Kn K n Q n Q n P n P +perm. ], (7) in the dilute limit, where we have used shorthand notation for the phase-density correlator ϕ K n K = ω k i kv +(ω k i kv). (8) The further strategy is the same as previously used. By calculating the next correctiontotheparticledistributionn k weareinpositiontoobtainthecondensate density up to the second-order of perturbation theory. Then the application of the coupling constant renormalization procedure () yields 6 V ρ t 3 nt ε k E k,q 0 k ε q ε k+q Eq E k+q 4 V ρ t ε k 4 V q 0 t[ εq E q ε q ], (9) for the fraction of non-condensed particles. In the same fashion we obtain the energy correction { [ ] Ē () v=0 N = V ρ t 3 ε k ε q ε [ k+q εk Eq + ][ εq ε k Eq ] }, (0) ε q k,q E k+q 5

6 and the second-order normal density fraction ( 3V ρ t 3 ) kq ε k+q m E k,q k E qe k+q { [ 8 3V ρ t 3 ε k ε q k,q E 4 k E q associated with the presence of disorder. 4 Results 4 3V ] ε k+q E k+q k,q + ε k E 4 k ρ t 3 ε3 k+q ε q E k+q 4 Eq [ εq E q ε k ε q ] }, () Interesting feature of the model with a delta-like potential is that all integrals appearing during the calculation of energy corrections, condensate and superfluid densities depletion can be evaluated analytically to the very end. Particularly for the fraction of condensed particles in the presence of classical bath we obtained the following expansion: n 0 n = nb 0 π n ρã 6πρã 3. () na After tedious integration an expression for a superfluid density with a similar structure was also derived n s n = π ρã { 4π 3 } 6 ln ρã 3. (3) 3 na 3 3 The situation with energy is more complicated. There is no problem in calculation of the first-order correction, but the second one is logarithmically divergent. Let us remind that the same type of problem originally occurs [3] during the computation of the beyond Lee-Huang-Yang ground-state energy of pure Bose gas and is totally connected with the point-like approximation of the two-body potential. Applying similar regularization procedure, i.e., cutting off the upper integration limit on a scale of order /ã we finally have (with logarithmic accuracy) Ē v=0 N = E { LHY N +ρ t +4 π naã 8πnaã ln } naã. (4) Of course, this result can be justified within more sophisticated consideration which particularly intends the explicit momentum dependence of the t-matrix (). Indeed, in the limit of k the leading asymptote is the following t /(kã) that provides the convergence of integrals and correctness of the above cut-off procedure. 6

7 5 Conclusions In summary, by means of hydrodynamic approach we have studied properties of a dilute Bose gas with the non-gaussian quenched disorder. The realization of external random potential is performed by inserting into the system a macroscopic number of non-interacting classical particles with infinite mass. Assuming that the two-body potential describing interaction between bath particles and bosons is short-ranged and weak we have perturbatively analyzed the thermodynamic and superfluid characteristics of the system. Particularly we obtained, in addition to the well-known Bogoliubov-like result, the second-order beyondmean-field corrections to the energy, condensate fraction and superfluid density of Bose gas. It is instructive to note that the presence of environment generally depletes the superfluid and condensate densities, furthermore, the second-order terms of these observables do not depend on number of bosons and are totally determined by the interaction with impurities. Acknowledgements We thank Prof. A. Rovenchak for invaluable comments. This work was partly supported by Project FF-30F (No. 06U00539) from the Ministry of Education and Science of Ukraine. References [] M. P. A. Fisher, P. B. Weichman, G. Grinstein, and D. S. Fisher, Phys. Rev. B 40, 546 (989). [] M. P. A. Fisher and D. S. Fisher, Phys. Rev. Lett. 6, 847 (988). [3] K. Huang and H.-F. Meng, Phys. Rev. Lett. 69, 644 (99). [4] C. A. Müller and C. Gaul, New J. Phys. 4, (0). [5] V. I. Yukalov and R. Graham, Phys. Rev. A 75, 0369 (007). [6] P. Navez, A. Pelster, and R. Graham, Appl. Phys. B 86, 395 (007). [7] V. I. Yukalov, E. P. Yukalova, K. V. Krutitsky, and R. Graham, Phys. Rev. A 76, (007). [8] G. E. Astrakharchik, J. Boronat, J. Casulleras, and S. Giorgini, Phys. Rev. A 66, (00). [9] M. Kobayashi and M. Tsubota, Phys. Rev. B 66, 7456 (00). [0] G. M. Falco, A. Pelster, and R. Graham, Phys. Rev. A 75, (007). [] T. Khellil, and A. Pelster, J. Stat. Mech. 06, (06); J. Stat. Mech. 07, (07). 7

8 [] O. Zobay, Phys. Rev. A 73, 0366 (006). [3] M. Timmer, A. Pelster and R. Graham, EPL (Europhys. Lett.) 76, 760 (006). [4] S. Pilati, S. Giorgini, and N. Prokof ev, Phys. Rev. Lett. 0, 5040 (009). [5] S. Giorgini, L. Pitaevskii, and S. Stringari, Phys. Rev. B 49, 938 (994). [6] C. Gaul and C. A.Müller, Phys. Rev. A 83, (0). [7] P. Lugan and L. Sanchez-Palencia, Phys. Rev. A 84, 036 (0). [8] S. Lellouch, L.-K. Lim, and L. Sanchez-Palencia, Phys. Rev. A 9, 0436 (05). [9] M. Boninsegni and H. R. Glyde, J. Low Temp. Phys., 5 (998). [0] H. R. Glyde, B. Fåk, and O. Plantevin, J. Low Temp. Phys. 3, 537 (998). [] A. V. Lopatin and V. M. Vinokur, Phys. Rev. Lett. 88, (00). [] C. A. Müller, Phys. Rev. A 9, 0360 (05). [3] V. Pastukhov, Ann. Phys., 37, 49 (06). [4] V. Pastukhov, J. Phys. A: Math. Theor. 48, (05). [5] V. Pastukhov, J. Low Temp. Phys. 86, 48 (07). [6] A. F. Andreev, E. P. Bashkin, Zh. sp. Teor. Fiz. 69, 39 (975) [Sov. Phys. JETP 4, 64 (975)]. [7] L. Viverit and S. Giorgini, Phys. Rev. A 66, (00). [8] V. Pastukhov, Phys. Rev. A 95, 0364 (07). [9] O. I. Utesov, M. I. Baglay, S. V. Andreev, arxiv: [30] P. Konietin, V. Pastukhov, J. Low Temp. Phys. 90, 56 (08). [3] T. D. Lee, K. Huang, and C. N. Yang, Phys. Rev. 06, 35 (957). [3] N. M. Hugenholtz and D. Pines, Phys. Rev. 6, 489 (959). 8