Phase diagrams of the Bose-Hubbard model at finite temperature

The phase transitions in the Bose-Hubbard model are investigated. A single-particle Green's function is calculated in the random phase approximation and the formalism of the Hubbard operators is used. The regions of existence of the superfluid and Mott insulator phases are established and the...

Повний опис

Збережено в:

Бібліографічні деталі
Дата:	2009
Автори:	Stasyuk, I.V., Mysakovych, T.S.
Формат:	Стаття
Мова:	English
Опубліковано:	Інститут фізики конденсованих систем НАН України 2009
Назва видання:	Condensed Matter Physics
Онлайн доступ:	http://dspace.nbuv.gov.ua/handle/123456789/120322
Теги:	Додати тег Немає тегів, Будьте першим, хто поставить тег для цього запису!
Назва журналу:	Digital Library of Periodicals of National Academy of Sciences of Ukraine
Цитувати:	Phase diagrams of the Bose-Hubbard model at finite temperature / I.V. Stasyuk, T.S. Mysakovych // Condensed Matter Physics. — 2009. — Т. 12, № 4. — С. 539-546. — Бібліогр.: 15 назв. — англ.

Репозитарії

Digital Library of Periodicals of National Academy of Sciences of Ukraine

id	irk-123456789-120322
record_format	dspace
spelling	irk-123456789-1203222017-06-12T03:04:56Z Phase diagrams of the Bose-Hubbard model at finite temperature Stasyuk, I.V. Mysakovych, T.S. The phase transitions in the Bose-Hubbard model are investigated. A single-particle Green's function is calculated in the random phase approximation and the formalism of the Hubbard operators is used. The regions of existence of the superfluid and Mott insulator phases are established and the (μ, t) (the chemical potential transfer parameter) phase diagrams are built. The effect of temperature change on this transition is analyzed and the phase diagram in the (T, μ) plane is constructed. The role of thermal activation of the ion hopping is investigated by taking into account the temperature dependence of the transfer parameter. The reconstruction of the Mott-insulator lobes due to this effect is analyzed. Дослiджено фазовi переходи в моделi Бозе-Хаббарда. Одночастинкова функцiя Грiна розрахована в наближеннi хаотичних фаз iз використанням формалiзму операторiв Хаббарда. Встановлено областi iснування надплинної фази та фази моттiвського дiелектрика i побудовано фазовi дiаграми в площинi (; t). Проаналiзовано вплив змiни температури на фазовi переходи та побудовано фазовi дiаграми в площинi (T; ). Дослiджено випадок термiчної активацiї перескоку частинок та проаналiзовано перебудову фазових дiаграм у цьому випадку. 2009 Article Phase diagrams of the Bose-Hubbard model at finite temperature / I.V. Stasyuk, T.S. Mysakovych // Condensed Matter Physics. — 2009. — Т. 12, № 4. — С. 539-546. — Бібліогр.: 15 назв. — англ. 1607-324X PACS: 03.75.Lm, 05.30.Jp DOI:10.5488/CMP.12.4.539 http://dspace.nbuv.gov.ua/handle/123456789/120322 en Condensed Matter Physics Інститут фізики конденсованих систем НАН України
institution	Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection	DSpace DC
language	English
description	The phase transitions in the Bose-Hubbard model are investigated. A single-particle Green's function is calculated in the random phase approximation and the formalism of the Hubbard operators is used. The regions of existence of the superfluid and Mott insulator phases are established and the (μ, t) (the chemical potential transfer parameter) phase diagrams are built. The effect of temperature change on this transition is analyzed and the phase diagram in the (T, μ) plane is constructed. The role of thermal activation of the ion hopping is investigated by taking into account the temperature dependence of the transfer parameter. The reconstruction of the Mott-insulator lobes due to this effect is analyzed.
format	Article
author	Stasyuk, I.V. Mysakovych, T.S.
spellingShingle	Stasyuk, I.V. Mysakovych, T.S. Phase diagrams of the Bose-Hubbard model at finite temperature Condensed Matter Physics
author_facet	Stasyuk, I.V. Mysakovych, T.S.
author_sort	Stasyuk, I.V.
title	Phase diagrams of the Bose-Hubbard model at finite temperature
title_short	Phase diagrams of the Bose-Hubbard model at finite temperature
title_full	Phase diagrams of the Bose-Hubbard model at finite temperature
title_fullStr	Phase diagrams of the Bose-Hubbard model at finite temperature
title_full_unstemmed	Phase diagrams of the Bose-Hubbard model at finite temperature
title_sort	phase diagrams of the bose-hubbard model at finite temperature
publisher	Інститут фізики конденсованих систем НАН України
publishDate	2009
url	http://dspace.nbuv.gov.ua/handle/123456789/120322
citation_txt	Phase diagrams of the Bose-Hubbard model at finite temperature / I.V. Stasyuk, T.S. Mysakovych // Condensed Matter Physics. — 2009. — Т. 12, № 4. — С. 539-546. — Бібліогр.: 15 назв. — англ.
series	Condensed Matter Physics
work_keys_str_mv	AT stasyukiv phasediagramsofthebosehubbardmodelatfinitetemperature AT mysakovychts phasediagramsofthebosehubbardmodelatfinitetemperature
first_indexed	2025-07-08T17:39:24Z
last_indexed	2025-07-08T17:39:24Z
_version_	1837101340702015488
fulltext	Condensed Matter Physics 2009, Vol. 12, No 4, pp. 539–546 Phase diagrams of the Bose-Hubbard model at finite temperature I.V. Stasyuk, T.S. Mysakovych Institute for Condensed Matter Physics of the National Academy of Sciences of Ukraine, 1 Svientsitskii Str., 79011 Lviv, Ukraine Received July 8, 2009 The phase transitions in the Bose-Hubbard model are investigated. A single-particle Green’s function is calcu- lated in the random phase approximation and the formalism of the Hubbard operators is used. The regions of existence of the superfluid and Mott insulator phases are established and the (µ, t) (the chemical potential – transfer parameter) phase diagrams are built. The effect of temperature change on this transition is analyzed and the phase diagram in the (T, µ) plane is constructed. The role of thermal activation of the ion hopping is investigated by taking into account the temperature dependence of the transfer parameter. The reconstruction of the Mott-insulator lobes due to this effect is analyzed. Key words: Bose-Hubbard model, phase transition PACS: 03.75.Lm, 05.30.Jp 1. Introduction The Bose-Hubbard model (BHM) has been intensively investigated in the last 15 years. This model is of great interest due to the experimental realization of Bose-Einstein condensation of ultra-cold atoms in optical lattices (see, for example, [1]). Various theoretical methods were used to study the model: mean-field theory [2], random phase approximation [3,4], strong coupling approach [5], quantum Monte-Carlo method [6], bosonic version of dynamical mean field theory [7]. The Hamiltonian of the model includes two terms. One is connected with the on-site energy U , which describes repulsion of bosons at a lattice site; another one describes the nearest-neighbor hopping with tunneling parameter t: H = − ∑ ij tijb + i bj + U 2 ∑ i ni(ni − 1) − µ ∑ i ni , (1) µ is the chemical potential. The competition between the kinetic and on-site repulsion terms defines the equilibrium state of the system. The existence of the superfluid and Mott-insulator phases is a characteristic feature of this model. When the kinetic energy dominates (t/U � 1) the ground state of the system is superfluid (SF). In the opposite case the ground state is a Mott insulator (MI). In the limit U → ∞ this model reduces to the hard-core boson model. Models of a lattice gas type with Pauli statistics of particles are often used for the description of ionic conductors and calculation of their conductivity starting from the works of Mahan [8] and others [9,10]. Recently the one-particle spectrum has been investigated in the one-dimensional limit [10]. As noted above, the Bose-Hubbard model can be directly applied to the description of optical lattices. Besides, the model of Bose-Hubbard type can be also useful not only for investigation of ionic conductivity in crystalline ionic conductors, but it can also be applied to description of the intercalation in crystals [11] and kinetics of ionic adsorption on the crystal surfaces [12]. In this work we consider the case of the finite value of the on-site interaction U and study phase transitions in the BHM at finite temperatures. Most of the previous investigations by the Bose- Hubbard model were restricted to the case of zero temperature. Recently only the low-temperature c© I.V. Stasyuk, T.S. Mysakovych 539 I.V. Stasyuk, T.S. Mysakovych case (T � U) has been investigated [3,13]. A special attention will be paid here to the investigation of the effect of thermal activation of the ion hopping on the shape of phase diagrams (we take into account the temperature dependence of the transfer parameter). The effect of thermal activation can be important in the case of ionic conductors. It was not studied in the framework of the BHM; such an investigation has been performed only in the hard-core boson limit [14] so far. 2. The random phase approximation We introduce the on-site basis \|n〈i (n is a number of bosons on a site i) and use the Hubbard operators Xnm i = \|n〉i〈m\|i. The Bose-operators of creation and annihilation can be expressed in terms of the X-operators: bi = ∑ n √ n + 1Xn,n+1 i , ni = ∑ n nXnn i . (2) Then the Hamiltonian (1) can be written as follows: H = H0 + H1, H0 = ∑ in λnXnn i ; λn = U 2 n(n − 1) − µn, H1 = − ∑ ij tijb + i bj . (3) To find the Green’s function 〈〈b\|b+〉〉 we employ the equation of motion technique: ω 〈〈 Xm,m+1 l \|Xr+1,r p 〉〉 = 1 2π δlpδmr 〈 Xmr l − Xr+1,m+1 l 〉 + (λm+1 − λm) 〈〈 Xm,m+1 l \|Xr+1,r l 〉〉 − ∑ ij tij √ m + 1δli 〈〈( Xmm l − Xm+1,m+1 l ) bj \|Xr+1,r p 〉〉 − ∑ ij tijδlj 〈〈 b+ i (√ m + 2Xm,m+2 l − √ mXm−1,m+1 l ) \|Xr+1,r p 〉〉 . (4) We use the random phase approximation (which is an analogy to the Hubbard-I approximation for the case of the fermionic Hubbard model) and perform the following decoupling: 〈〈 (Xm,m l − Xm+1,m+1 l )bj \|Xr+1,r p 〉〉 ≈ 〈 Xmm − Xm+1,m+1 〉 〈〈 bj \|Xr+1,r p 〉〉 , 〈〈 b+ i Xm,m+2 l \|Xr+1,r p 〉〉 ≈ 0. (5) It should be noted that here we consider the Mott-insulator phase, where 〈bi〉 = 0. If we want to investigate the superfluid phase we should take into account that in this phase 〈bi〉 6= 0 (this case is considered below). As a result, the following expression for the Green’s function is obtained: 〈〈 Xm,m+1 l \|Xr+1,r p 〉〉 = 1 2π δlpδmr Qm ω − ∆m − Qm ω − ∆m √ m + 1 ∑ j tlj 〈〈 bj \|Xr+1,r p 〉〉 , (6) ∆m = λm+1 − λm, 〈 Xmm − Xm+1,m+1 〉 = Qm . (7) We pass to k-representation and the Green’s function Gk(ω) ≡ 〈〈b\|b+〉〉kω can be written as Gk(ω) = 1 2π g0 1 + g0tk , g0(ω) = ∑ m Qm(m + 1) ω − ∆m . (8) 540 Phase diagrams of the Bose-Hubbard model at finite temperature We should note that in the independent subband approximation (this approximation is valid in the case U � t and in the regions of the integer value of µ/U) the Green’s function takes the form Gk = 1 2π ∑ m (m + 1)Qm ω − ∆m + tk(m + 1)Qm , (9) and the subband energy is equal to λm = ∆m − (m + 1)Qmtk. In order to study the phase with the Bose-Einstein condensate we can also apply the decoupling of the (5) type, which in this case is more complicated: 〈〈 (Xm,m l − Xm+1,m+1 l )bj \|Xr+1,r p 〉〉 ≈ Qm 〈〈 bj \|Xr+1,r p 〉〉 + 〈bj〉 〈〈 (Xmm l − Xm+1,m+1 l )\|Xr+1,r p 〉〉 . (10) The function 〈〈 (Xm,m l − Xm+1,m+1 l )bj \|Xr+1,r p 〉〉 satisfies the following equation: ω 〈〈 Xm,m l − Xm+1,m+1 l \|Xr+1,r p 〉〉 = 1 2π δlp( 〈 Xm,m−1 〉 δm,r+1 + 〈 Xm+2,m+1 〉 δm+1,r −2 〈 Xm+1,m 〉 δmr) − ∑ j tlj (√ m 〈〈 Xm,m−1 l bj − b+ j Xm−1,m l \|Xr+1,r p 〉〉 + √ m + 2 〈〈 Xm+2,m+1 l bj − b+ j Xm+1,m+2 l \|Xr+1,r p 〉〉 −2 √ m + 1 〈〈 Xm+1,m l bj − b+ j Xm,m+1 l \|Xr+1,r p 〉〉) . (11) Using decoupling (10) in the independent subband approximation (which is equivalent to taking into consideration the transition between \|m〉 and \|m+1〉 levels only) the equation takes a simpler form ω 〈〈 Xm,m l − Xm+1,m+1 l \|Xr+1,r p 〉〉 = − 1 π δlpδmr( 〈 Xm+1,m 〉 + 2 (m + 1) ∑ j tlj (〈 Xm,m+1 j 〉〈〈 Xm+1,m l \|Xr+1,r p 〉〉 − 〈 Xm,m+1 l 〉〈〈 Xm+1,m j \|Xr+1,r p 〉〉) − 2 (m + 1) ∑ j tlj (〈 Xm+1,m j 〉〈〈 Xm,m+1 l \|Xr+1,r p 〉〉 − 〈 Xm+1,m l 〉〈〈 Xm,m+1 j \|Xr+1,r p 〉〉) . (12) In a similar way, we can obtain equations for the Green’s functions 〈〈 Xm+1,m l \|Xr+1,r p 〉〉 and 〈〈 Xm,m+1 l \|Xr+1,r p 〉〉 : (ω + ∆m) 〈〈 Xm+1,m l \|Xr+1,r p 〉〉 = (m + 1) ∑ j tljQm 〈〈 Xm+1,m j \|Xr+1,r p 〉〉 + (m + 1) ∑ j tlj 〈 Xm+1,m j 〉〈〈 Xm,m l − Xm+1,m+1 l \|Xr+1,r p 〉〉 , (13) (ω − ∆m) 〈〈 Xm,m+1 l \|Xr+1,r p 〉〉 = 1 2π δlpδmrQm − (m + 1) ∑ j tljQm 〈〈 Xm,m+1 j \|Xr+1,r p 〉〉 − (m + 1) ∑ j tlj 〈 Xm,m+1 j 〉〈〈( Xm,m l − Xm+1,m+1 l ) \|Xr+1,r p 〉〉 . (14) The levels \|m〉, \|m+1〉 form a two-state subspace where we can introduce pseudospin formalism: Xm,m+1 l = S− l,m , Xm+1,m l = S+ l,m , 1 2 〈Xm+1,m+1 l − Xmm l 〉 = Sz l,m . (15) 541 I.V. Stasyuk, T.S. Mysakovych The mean value Qm = −2〈Sz m〉 is connected with the longitudinal component of the pseudospin, while the order parameter 〈bj〉 in the condensate phase is proportional to the transverse component 〈bj〉 ≈ √ m + 1〈S− m〉 = √ m + 1〈Sx m〉. The set of equations (12), (13), and (14) can be rewritten using pseudospin Green’s functions. The equations are analogous to those obtained in the RPA for the XXZ model. Their solution can be written in the form given in [14] 〈〈S+ m\|S− m〉〉 = 1 2π 〈σz m〉 Em(cos2 θm + 1) + 2ω cos θm − 4〈σz m〉Jq,m cos2 θm ω2 − (Em − 2〈σz m〉Jq,m)(Em − 2〈σz m〉Jq,m) cos2 θm . (16) In the mean field approximation the angle θm and parameter 〈σz m〉 are defined by the following equations: ∆m sin θm − 2(m + 1)t0〈Sx m〉 cos θm = 0, (17) 〈σz m〉 = 1 2 tanh βEm 2 , Em = √ ∆2 m + 4(m + 1)2t20〈Sx m〉2 , t0 ≡ tq=0 , Jq,m = (m + 1)tq , 〈Sx m〉 = −〈σz m〉 sin θm . (18) The solution sin θm = 0 (〈Sx m〉 = 0) corresponds to the normal (nonsuperfluid) phase; in the phase with condensate cos θm = − ∆m 2(m+1)t0〈σz m 〉 and 〈σz m〉 is given by the solution of the following equation 〈σz m〉 = 1 2 tanh (β(m + 1)t0〈σz m〉) . (19) The transition to the phase with condensate takes place at temperature T c m = ∆m ( ln 1 + ∆m J0,m 1 − ∆m J0,m )−1 . (20) This temperature reaches its maximum value T c max = J0,m/2 at ∆m = mU − µ = 0. The case of pair of the states \|0〉, \|1〉 corresponds to the hard core boson model; the peculiarities of the spectrum and density of states were investigated in [14]. The thermodynamics of the model was also studied in [14] in the case of the thermal activated hopping. 3. Results The phase transition from the Mott-insulator to superfluid phase is characterized by divergence of the Green’s function Gk=0(ω = 0) → ∞. This condition can be obtained using the following reasoning. Let us find the expression for the commutator [Xm,m+1 l , H ] in the Hubbard-I type approximation: [Xm,m+1 l , H ] ≈ ∆mXm,m+1 l − ∑ j tlj √ m + 1〈Xmm − Xm+1,m+1〉bj . (21) The mean value of this commutator should be equal to zero (because the average 〈Xm,m+1〉 does not depend on time in a stationary state). After decoulping we get: 〈bk〉[1 − tk ∑ m m + 1 ∆m Qm] = 0. (22) Thus, we can see that the nonzero solution 〈bk〉 6= 0 appears in the point of divergence of the Green’s function Gk(ω = 0) → ∞ (8). The condition (22) in the case t > 0 (considered in this paper) is fulfilled at the center of the Brillouin zone when k = 0 (an uniform phase case) and Gk=0(ω = 0) → ∞ in the phase transition point. In the case t < 0 this condition is satisfied at the Brillouin zone edge when k = (π, π, . . . ) (a modulation phase case). 542 Phase diagrams of the Bose-Hubbard model at finite temperature In our numerical calculations we used the diagonal Hamiltonian H0 to calculate the averages 〈Xmm〉. In this case the condition (22) coincides with that obtained in [13] in the framework of the diagrammatic hopping expansion at the calculation of temperature Green’s functions (the first order in the self-energy in the notations used in [13]). Let us consider the effect of thermal activation. We rewrite the parameter of the ion hopping in the following form: t = t0 exp(−β∆). Such an approximation was used in [8,14] to take into account the presence (in ionic conductors) of a static barrier over which a particle should hop. It can also be obtained as a result of renormalization due to the interaction with phonons. 0 0.1 0.2 0.3 0.4 -0.5 0 0.5 1 1.5 2 2.5 3 W /U µ/U (MI) (MI) (MI) (SF) (a) T=0.01U T=0.05U T=0.1U T=0.2U 0 0.1 0.2 0.3 0.4 0.5 0 0.5 1 1.5 2 2.5 3 W /U µ/U (MI) (MI) (MI) (SF) (b) ∆=0 ∆=0.01U Figure 1. The (W,µ) phase diagram at ∆ = 0 (a) and ∆ 6= 0 (b). In figure 1 the (W, µ) phase diagram is shown (W is a half width of the initial energy band, −W < t < W ). The system can pass from the MI to SF phase at the change of the chemical potential as well as at the change of the transfer parameter. The effect of thermal activation on phase transitions is illustrated in figure 1b. Thermal activation leads to the expansion of the region of the Mott-insulator phase. 0 0.1 0.2 0.3 0.4 0.5 0.6 0 0.5 1 1.5 2 2.5 T/ U µ/U (MI) (SF) (a) W=0.07U W=0.1U W=0.15U 0 0.2 0.4 0.6 0.8 0.8 1 1.2 T/ U µ/U (MI) (SF) ∆=0 ∆=0.01U ∆=0.025U (b) Figure 2. The (T, µ) phase diagram at ∆ = 0 (a) and ∆ 6= 0, W = 0.07U (b). 543 I.V. Stasyuk, T.S. Mysakovych The (T, µ) phase diagrams are given in figure 2. The effect of thermal activation is illustrated in figure 2b for the case of the one SF lobe (µ/U ≈ 1). At the increase of the activation energy this region narrows and then the superfluid phase disappears (there is a critical value ∆cr). It should be noted that at the fixed value of the chemical potential there are two critical temperatures Tc1 and Tc2 at which the phase transitions between the MI and superfluid phase occur (it takes place at ∆ 6= 0). One can think that the lower critical temperature Tc2 can be connected with the temperature of the transition to the superionic state in the case of superionic crystals. The simple estimate, based on the typical values of t and ∆ parameters for proton conductors (∆ ∼ 0.1 . . .0.4 eV; t ∼ 0.2 . . .0.5 eV) shows that Tc2 can be of the order of several hundreds of Kelvins, which corresponds to temperatures observed in the known family M3H(XO4)2 (M=NH4, Rb, Cs; X=S,Se) of hydrogen-bonded superionic crystals [15]. 0 4 8 12 16 0 4 8 12 U/W µ/ W (MI) (SF) (MI)(a) T=0.1W T=0.5W T=0.7W 0 4 8 12 16 0 4 8 12 U/W µ/ W (MI) (SF) (MI)(b) ∆=0 ∆=0.1W Figure 3. The (U, µ) phase diagram at ∆ = 0 (a) and ∆ 6= 0, T = 0.5W (b). From the (U, µ) phase diagrams shown in figure 3 we can see that the effect of thermal activation on the phase diagrams is similar to the effect of temperature change in the case ∆ = 0. In figure 4 the (T, µ) phase diagram calculated in the independent subband approximation using the equation (20) (where J0,m = (m + 1)t0 exp(−βm c ∆)) is shown (dashed line). Solid line is the line shown in figure 2a in the case W = 0.07U . We can see that the independent subband approximation is applicable at small values of m. In this region the expression (20) adequately describes the nearly linear growth of T m c with the increase of the number m. 4. Conclusions The phase transitions in the Bose-Hubbard model at finite temperature are investigated. Models of such type can be applied to the description of optical lattices, ionic conductivity, intercalation in crystals, and kinetics of ionic adsorption on the crystal surfaces. A single-particle Green’s function is calculated in the random phase approximation (which is an analogy to the Hubbard-I approximation in the case of the fermionic Hubbard model) and the formalism of the Hubbard operators is used. This Green’s function is obtained for the Mott- insulator and superfluid phases (in the last case we use the independent subband approximation, which is appropriate for the case of the noninteger (intermediate) mean number of the particle concentration per lattice site). The regions of the occurrence of the superfluid and Mott insulator phases are established and the phase diagrams in the plane (µ, t) (the chemical potential – transfer parameter) are built. The effect of temperature change on this transition is analysed and the phase diagram in the (T, µ) 544 Phase diagrams of the Bose-Hubbard model at finite temperature 0 0.05 0.1 0.15 0 0.5 1 1.5 2 2.5 T/ U µ/U (MI) (SF) Figure 4. Dashed line denotes the phase transition line obtained in the independent subband approximation using (20). Solid line denotes the phase transition line obtained using (22). W = 0.07U, ∆ = 0. plane is constructed. The role of thermal activation of the ion hopping is investigated by taking into account the temperature dependence of the transfer parameter (t = t0 exp(−β∆)). It is shown that thermal activation leads to the reconstruction of the Mott-insulator lobes (in particular, to the narrowing of the SF phase region). It is revealed that due to thermal activation of the ion hopping, there are two critical temperatures at which the MI-SF phase transition takes place (at the fixed value of the chemical potential). A simple estimation shows that the lower temperature can be quite high in proton conductors and can correspond to the temperature of the superionic phase transition in such compounds. References 1. Greiner M., Mandel O., Esslinger T., Hänsch T.W., Bloch I., Nature, 2002, 415, 39–44. 2. Sheshadri K., Krishnamurthy H.R., Pandit R., Ramakrishnan T.V., Europhys. Lett., 1993, 22, 257– 263. 3. Konabe S., Nikuni T., Nakamura M., Phys. Rev. A, 2006, 73, 033621–12. 4. Ohashi Y., Kitauri M., Matsumoto H., Phys. Rev. A, 2006, 73, 033617–5. 5. Freericks J.K., Monien H., Europhys. Lett., 1994, 26, 545–550. 6. Capogrosso-Sansone B., Söyler S.G., Prokof’ev N., Svistunov B., Phys. Rev. A, 2008, 77, 015602–4. 7. Byczuk K., Vollhardt D., Phys. Rev. B, 2008, 77, 235106–14. 8. Mahan G.D., Phys. Rev. B, 1976, 14, 780–793. 9. Tomoyose T., J. Phys. Soc. Jpn., 1997, 66, 2383–2385. 10. Stasyuk I.V., Dulepa I.R., Condens. Matter Phys., 2007, 10, 259–268. 11. Mysakovych T.S., Krasnov V.O., Stasyuk I.V., Condens. Matter Phys., 2008, 11, 663–667. 12. Reilly P.D., Harris R.A., Whaley K.B., J. Chem. Phys., 1991, 95, 8599–8615. 13. Ohliger M., Pelster A. Preprint arXiv: cond-mat.stat-mech/0810.4399, 2008. 14. Stasyuk I.V., Dulepa I.R., J. Phys. Stud., 2009, 13, No. 2, 2701. 15. Stasyuk I.V., Pavlenko N., J. Phys.: Condens. Matter, 1998, 10, 7079–7090. 545 I.V. Stasyuk, T.S. Mysakovych Фазовi дiаграми моделi Бозе-Хаббарда при скiнченнiй температурi I.В. Стасюк, Т.С. Мисакович Iнститут фiзики конденсованих систем НАН України, 79011 Львiв, вул. Свєнцiцького, 1 Отримано 8 липня 2009 р. Дослiджено фазовi переходи в моделi Бозе-Хаббарда. Одночастинкова функцiя Грiна розрахова- на в наближеннi хаотичних фаз iз використанням формалiзму операторiв Хаббарда. Встановлено областi iснування надплинної фази та фази моттiвського дiелектрика i побудовано фазовi дiагра- ми в площинi (µ, t). Проаналiзовано вплив змiни температури на фазовi переходи та побудовано фазовi дiаграми в площинi (T, µ). Дослiджено випадок термiчної активацiї перескоку частинок та проаналiзовано перебудову фазових дiаграм у цьому випадку. Ключовi слова: модель Бозе-Хаббарда, фазовi переходи PACS: 03.75.Lm, 05.30.Jp 546

Phase diagrams of the Bose-Hubbard model at finite temperature

Репозитарії

Схожі ресурси