Accurate Thermodynamic Properties of Ideal Bosons in a Highly Anisotropic 2D Harmonic Potential

One can derive an analytic result for the issue of Bose–Einstein condensation (BEC) in anisotropic 2D harmonic traps. We find that the number of uncondensed bosons is represented by an analytic function, which includes a series expansion of q-digamma functions in mathematics. One can utilize this analytic result to evaluate various thermodynamic functions of ideal bosons in 2D anisotropic harmonic traps. The first major discovery is that the internal energy of a finite number of ideal bosons is a monotonically increasing function of anisotropy parameter p. The second major discovery is that, when p≥0.5, the changing with temperature of the heat capacity of a finite number of ideal bosons possesses the maximum value, which happens at critical temperature Tc. The third major discovery is that, when 0.1≤p<0.5, the changing with temperature of the heat capacity of a finite number of ideal bosons possesses an inflection point, but when p<0.1, the inflection point disappears. The fourth major discovery is that, in the thermodynamic limit, at Tc and when p≥0.5, the heat capacity at constant number reveals a cusp singularity, which resembles the λ-transition of liquid helium-4. The fifth major discovery is that, in comparison to 2D isotropic harmonic traps (p=1), the singular peak of the specific heat becomes very gentle when p is lowered.


Introduction
The Bose-Einstein condensation (BEC) of boson systems is always a hot-spot of research in physics.Three independent teams realized the BEC of trapped ultra-cold alkali atoms in 1995 [1][2][3].Quasi-particles of the Bose type in solids can also be in the BEC state.Recently, scientists have also observed the BEC of excitons [4,5], exciton-polaritons [6,7], and magnons [8,9] in several solid-state systems.Furthermore, the BEC of photons in low-dimensional optical microcavities has become a hot-spot of research in optics [10].Thermalization of a 2D photon gas has been realized in a dye-filled curved-mirror microcavity [11].The BEC of 2D photons in a dye-filled curved-mirror microcavity has been observed by Weitz and colleagues [12,13].In a previous work [14], we have investigated the accurate thermodynamic properties of ideal bosons in a 2D isotropic harmonic potential.In another previous work [15], we have proposed an exact analytic result for ideal boson gases in a highly anisotropic 2D harmonic potential.In this paper, we shall study the thermodynamic properties of ideal Bose atoms in a 2D anisotropic harmonic potential.The exact thermodynamic theory of the BEC state of such ideal Bose atoms will be expounded.
In 1995, the BEC of ultra-cold alkali atoms was realized in 3D harmonic potentials.This achievement kindled great interest in the BEC of ultra-cold bosons in 2D harmonic potentials.The 2D quantum system has been an attractive object since the creation of quantum theory [16].The 2D characters bring about unbelievably rich phenomena.Some accurate analytic results in 2D systems can be acquired using special methods.First, Salasnich analyzed the thermodynamics of ideal gases in a generic power-law potential [17].
Recently, Weitz and colleagues investigated both experimentally and theoretically the statistical mechanics of a gas of massive photons in a 2D harmonic potential [12,18].Further, Stein and Pelster investigated the thermodynamic properties of trapped ideal photons in dimensional crossover from 2D to 1D [19].However, these authors did not compute the internal energy, the entropy, the Helmholtz free energy, and the heat capacity of ideal photons in 2D anisotropic harmonic traps.The Berezinskii-Kosterlitz-Thouless (BKT) transition was realized in ultra-cold Bose atoms in a homogeneous quasi-2D optical trap [20].For a homogeneous 2D boson system with interactions, either a BEC or a BKT transition occurs depending upon whether the interaction strength is small or large, respectively.In fact, a crossover between BEC and BKT occurs for increasing interaction strength [21].
The object under study is an ideal boson gas in a 2D anisotropic harmonic potential.In reality, a non-interacting Bose-Einstein condensate can be produced by means of Feshbach resonance [22].In order to deal with an ideal boson gas in a 2D anisotropic harmonic potential, in all the references [23], researchers have utilized various approximate methods.In order to keep off these approximate methods, we propose an exact analytic result in an ideal boson gas in a 2D anisotropic harmonic potential.Using this accurate analytic solution, we elaborate upon the accurate thermodynamic behavior of an ideal boson gas in a 2D anisotropic harmonic potential.This thermodynamic theory is effective for arbitrary temperature, boson number, and anisotropy parameter.The first major discovery is that the internal energy of a finite number of ideal bosons is a monotonically increasing function of the anisotropy parameter p.The second major discovery is that, when p ≥ 0.5, the changing with temperature of the heat capacity of a finite number of ideal bosons possesses the maximum value, which happens at critical temperature T c .The third major discovery is that, when 0.1 ≤ p < 0.5, the changing with temperature of the heat capacity of a finite number of ideal bosons possesses an inflection point, but when p < 0.1, the inflection point disappears.The fourth major discovery is that, in the thermodynamic limit, at T c and when p ≥ 0.5, the heat capacity at constant number reveals a cusp singularity, which resembles the λ-transition of liquid helium-4.The fifth major discovery is that, in comparison to 2D isotropic harmonic traps (p = 1), the singular peak of the specific heat becomes very gentle when p is lowered.These accurate thermodynamic properties of ideal boson gases in a 2D anisotropic harmonic potential can be confirmed in current physics laboratories.
The remainder of this paper is arranged as follows.In Section 2, one derives the thermodynamic potential of ideal boson gases in a 2D anisotropic harmonic potential.In Section 3, we describe the exact thermodynamic behavior of ideal boson gases in a 2D anisotropic harmonic trap.Section 4 depicts the phase transitions of ideal boson gases in the thermodynamic limit.The general conclusions are given in Section 5.

Many-Particle State in the Particle-Number Representation
The present paper deals with an ideal boson gas with zero spin.The model system is comprised of N non-interacting bosons moving in a 2D anisotropic harmonic trap [24].In the beginning, one must bring out the two operators bn x n y and b † n x n y .bn x n y and b † n x n y denote the destruction and production operators of an atom occupying the oscillator level with n x and n y , respectively.bn x n y and b † n x n y obey the Bose commutation relations: Here, we can bring out the number operator Nn x n y of atoms occupying the oscillator level with n x and n y by Nn x n y = b † n x n y bn x n y . ( The eigenvectors of Nn x n y are |N n x n y and Nn x n y obeys the eigenvalue equation where N n x n y denotes an eigenvalue of Nn x n y .The eigenvector |N n x n y can be represented as |N n x n y = 1 Hence, the total number operator N of bosons is presented as According to the program of second quantization, we can immediately obtain the Hamiltonian of the boson system as where E n x n y is the level of a 2D anisotropic harmonic oscillator given by where ω x and ω y signify the angular frequencies of the trap along the x and y axes, respectively.Because the number operator N commutes with the Hamiltonian Ĥ, the common eigenstates of N and Ĥ are presented as The eigenstate vector in Equation ( 8) is symmetric under the interchange of any two production operators, consistently with the Bose-Einstein statistics.

Thermodynamic Potential
The eigenstate vector in Equation ( 8) represents a multimode number state of bosons, which is a pure state and hence far from thermal equilibrium.Nevertheless, the boson system in a 2D anisotropic harmonic potential goes into thermal equilibrium.This equilibrium is constructed by way of the successive collisions between bosons.The boson system in thermal equilibrium is represented by a definite temperature T. Since the number of bosons in the potential is conserved, the 2D boson system has a non-vanishing chemical potential µ.In order to represent the thermal equilibrium state of the boson system, one must devise a grand canonical ensemble of bosons.
In the beginning, one needs to derive the thermodynamic potential Ω. Ω is only a function of temperature and chemical potential.In the grand canonical ensemble at temperature T, the grand partition function Z is given by where k B denotes Boltzmann's constant.The basis states employed in the trace correspond to the eigenstates of the Hamiltonian Ĥ, which are presented by Equation (8).Here, Ĥ signifies the Hamiltonian of the boson system and is presented by Equation (6).The thermodynamic potential is linked to the Hamiltonian of the boson system via the grand partition function, All the macroscopical thermodynamic functions may be calculated from the thermodynamic potential.
Substituting Equations ( 5) and ( 6) into Equation ( 9), the grand partition function is acquired as If the trace in Equation ( 11) is written out at length with the complete set of eigenstates |N n x n y of number operator Nn x n y , one obtains Consequently, the grand partition function can be acquired as If we take the logarithm of Equation ( 13), we can obtain the thermodynamic potential as

Ushering of the Reduced Chemical Potential x a
For temperature T, N n x n y signifies the thermal average of the number of atoms occupying the oscillator level with n x and n y .The mean occupation number N n x n y of oscillator levels may be acquired from the thermodynamic potential as Substituting Equation (14) into Equation (15), one immediately finds that Equation ( 16) represents the famous Bose-Einstein distribution.The chemical potential µ is decided by the limitation that the total number of bosons in the system equals N: In the present paper, we assume that ω x ≥ ω y .It is necessary to define an anisotropy parameter p by the relation p = ω y /ω x .The anisotropy parameter varies in the range of 0 ≤ p ≤ 1.When p = 0, the trap is a 1D harmonic trap along the x axis.When p = 1, the trap is a 2D isotropic harmonic trap.
To decide µ, one must usher the fugacity z by the relation z = exp(µ * /k B T), where we have ushered an effective chemical potential µ * = µ − 1 2 h(ω x + ω y ).Further, we can usher the parameters q x = exp(−hω x /k B T) and x a = 1 − µ * /hω x .As the temperature T occurs in the definition z = exp(µ * /k B T), the fugacity z can not embody the chemical potential µ * by much and, hence, z is a bad physical quantity.The quantity x a embodies the chemical potential µ * a lot and, therefore, x a corresponds to a good physical quantity.Consequently, the quantity x a is known as the reduced chemical potential.In the same method, the quantity q x embodies the temperature T a lot and, hence, q x corresponds to a good physical quantity.By means of the good physical quantities x a and q x , Equation ( 17) is adapted into an equation of state: where the first term signifies the number of condensed bosons and H q x (x a ) represents the number of uncondensed bosons and is presented by where the sum upper limit ∞ is substituted for an upper cutoff k c and, in reality, one can let k c = 200.The numeral calculation displays that the upper cutoff k c = 200 is enough for a high-precision computation.ψ q (x) stands for the q-digamma function defined by ψ q (x) = d[ln Γ q (x)]/dx, where Γ q (x) stands for the q-gamma function defined by when |q| < 1 and x = 0, −1, −2, • • • .Jackson ushered the q-gamma function [25] and Krattenthaler and Srivastava ushered the q-digamma function [26].In past decades, the q-gamma function and the q-polygamma function have played an important role in science and technology [27].The reduced chemical potential x a can be calculated numerically from Equation (18).x a is a function of temperature T, boson number N, and anisotropy parameter p.To meet Equation ( 18), we require that x a ≥ 1.When x a = 1, a 2D atomic gas goes into the state of BEC.

Internal Energy
In the following, we shall investigate the thermodynamic properties of ideal bosons in a 2D anisotropic harmonic potential.A major thermodynamic function in the ideal boson system corresponds to the internal energy E, as defined by The internal energy E is decided collectively by Equations ( 16), (18), and (23).To calculate E, we first write out E n x n y = n x hω x + n y hω y + 1 2 h(ω x + ω y ).On substituting the last equation into Equation ( 16), Equation ( 16) can be adapted as follows: where β = 1/k B T. The system tends to go into the ground state with n x = n y = 0.Because we have shifted the zero-point energy into the effective chemical potential, the groundstate energy has been taken to be zero.Substituting Equation (24) into Equation ( 23), we thus acquire n y hω y ze −β(n x hω x +n y hω y ) where q y = exp(−βhω y ).Notice that z and β are two independent Lagrange multipliers.Consequently, we keep z constant for the partial derivative of the brace with respect to β. Equation ( 25) suggests that the condensed bosons have no donation to the internal energy E. Finishing the partial derivative with respect to β, one can simplify Equation ( 25) into the following form: where we have employed the relation z = exp(βµ * ) and x a = 1 − µ * /hω x .The infinite series in Equation ( 26) converges very rapidly and, therefore, it can be computed numerically.
We can do a numeral computation about E and, therefore, we let ω y /2π = 10.0 Hz, which is available to an actual experiment [28].For the sake of contrast, one can usher a 2D condensation temperature T 2D by the definition k B T 2D = hω y .Consequently, we acquire T 2D = 0.4799 nK.In the quantum statistical mechanics, the temperature of a 2D boson gas is meaningful only when the boson number is large enough.The internal energy E is a function of temperature T, boson number N, and anisotropy parameter p.When calculating E, we must combine Equation (18) with Equation ( 26).The reduced chemical potential x a can be calculated from Equation (18) to obtain the internal energy.According to Equation (26), in Figure 1 we show changing of the internal energy E with the boson number N for various T. We select the anisotropy parameter p = 0.5.One can find that, at T = 1 nK, E/hω x ≈ 7.580 for all N. Put differently, when p = 0.5 and T = 1 nK, a 2D boson gas is always in the state of BEC for arbitrary boson number N.An interesting discovery is that, when T ≤ 100 nK and N ≥ 3 × 10 4 , E/hω x is a constant dependent of temperature.Put differently, when T ≤ 100 nK and N ≥ 3 × 10 4 , a 2D boson gas is always in the state of BEC.These findings present an insight that, for a fixed T, there is a critical boson number N c , above which x a = 1.Based on Equation (26), in Figure 2 we reveal changing of the internal energy E with the temperature T for various N.We select the anisotropy parameter p = 0.5.When N is fixed, E is a monotonically increasing function of temperature T. When N is fixed, however, changing with temperature of the internal energy of a finite number of ideal bosons possesses an inflection point.The calculation shows that the inflection point corresponds to a critical temperature T c . Figure 2 reveals that the internal energy of a finite number of ideal bosons possesses a classic limit E = 2Nk B T as T T c .The computation also shows that, when p < 0.5 and for a fixed N, the inflection point on the curve of E versus T disappears.According to Equation (26), in Figure 3 we display the variation of the internal energy E with the anisotropy parameter p for various T. In Figure 3, we take N = 10 4 .An interesting discovery is that, for a fixed N and when T is finite, E = 0 at p → 0. When N and T are fixed, E is a monotonically increasing function of anisotropy parameter p.These findings present an insight that, when T and N are fixed, there is a critical anisotropy parameter p c , below which x a = 1.

Entropy and Helmholtz Free Energy
In statistical mechanics, entropy (conventional sign S) is commonly interpreted as a degree of disorder.The entropy S of the ideal boson system can be derived from the thermodynamic potential as Substituting Equation ( 14) into Equation ( 27), we immediately discover that where we have ushered an effective chemical potential µ * = µ − 1 2 h(ω x + ω y ) and, therefore, E n x n y = n x hω x + n y hω y .Utilizing Equations ( 16), (17), and (23), we can simplify Equation (28) as where q x = exp(−βhω x ), q y = exp(−βhω y ) and we have introduced a reduced chemical potential x a = 1 − µ * /hω x .E is the internal energy presented by Equation (26).Now, we can evaluate the entropy S, which is a function of temperature T, boson number N, and anisotropy parameter p.If one evaluates S, one must associate Equation (29) with Equations ( 18) and ( 26).The reduced chemical potential x a can be computed from Equation (18) to acquire the internal energy and the entropy.Based on Equation (29), changing of the scaled entropy S/k B with the boson number N is presented in Figure 4 for various T. We select the anisotropy parameter p = 0.2. Figure 4 demonstrates that, for a fixed temperature T, the entropy S is a slowly increasing function of the boson number N when N > N c .According to Equation (29), changing of the scaled entropy S/k B with the temperature T is presented in Figure 5 for various N.We select the anisotropy parameter p = 0.2. Figure 5 shows that, when T < 0.4 nK, S/k B is a constant dependent of the boson number N. At T = 4 nK, the curves corresponding to various N cross together.Moreover, Figure 5 reveals that, when T > 4 nK, the entropy S is a rapidly ascending function of temperature T for a fixed boson number N. Based on Equation ( 29), in Figure 6 we show changing of the scaled entropy S/k B with the anisotropy parameter p for various T. In Figure 6, we take N = 10 4 .An interesting discovery is that, for all N and T, S = 0 at p → 0. For a fixed N and T, S is a monotonically ascending function of anisotropy parameter p.In the next step, we make a Legendre transformation of the internal energy E to the Helmholtz free energy F. The Helmholtz free energy F is defined by At constant temperature, the Helmholtz free energy is minimized in thermal equilibrium.Putting Equation (29) into Equation (30), one immediately finds that  The Helmholtz free energy F given by Equation ( 31) is a function of temperature T, boson number N, and anisotropy parameter p.If one evaluates F, one must associate Equation (18) with Equation ( 31).The reduced chemical potential x a can be computed from Equation (18) to acquire the Helmholtz free energy.According to Equation (31), changing of the scaled free energy −F/hω x with the boson number N is presented in Figure 7 for various T. We select the anisotropy parameter p = 0.2. Figure 7 demonstrates that the free energy F is always negative and is a slowly descending function of the boson number N for a fixed temperature T. Based on Equation (31), changing of the scaled free energy −F/hω x with the temperature T is presented in Figure 8 for various N.We select the anisotropy parameter p = 0.2. Figure 8 reveals that the free energy F is a rapidly descending function of temperature T for a fixed boson number N. Notice that F = 0 at T = 0 K.At T = 5 nK, the curves corresponding to various N cross together.According to Equation (31), in Figure 9 we show changing of the scaled free energy −F/hω x with the anisotropy parameter p for various T. In Figure 9, we take N = 10 4 .An interesting discovery is that, for all N and T, F = 0 at p → 0. When N and T are fixed, F is a monotonically descending function of anisotropy parameter p.The plot is drawn for T = 1, 15, 50, 100 nK.We select the axial angular frequency ω y /2π = 10.0 Hz and the anisotropy parameter p = 0.2.

Heat Capacity at Constant Number
The heat capacity is the quantity of heat it takes to lift the temperature of a substance by one degree centigrade.This obviously depends on the conditions under which the heating happens.Since the internal energy E is a function of temperature T, boson number N and anisotropy parameter p, we study the heat capacity at constant number, which is specified by  The internal energy E given by Equation ( 26) is a composite function of intermediate variables q x , q y , and x a , where q x = exp(−hω x /k B T) and q y = exp(−hω y /k B T).The reduced chemical potential x a is an implicit function of temperature T, boson number N, and anisotropy parameter p and is determined by Equation (18).
Employing Equation (26), one can readily derive the two partial derivatives of E as follows: 1 where If we utilize Equation ( 18) and make some arrangements, then we can infer the partial derivative of x a as follows: where ψ q (x) = d[ψ q (x)]/dx denotes the q-trigamma function and The evaluation displays that the system possesses a critical boson number N c at which the heat capacity at constant number reaches the maximum value.An interesting discovery is that, in Equation ( 32), (∂E/∂T) x a is positive but (∂E/∂x a ) T is negative.(∂E/∂T) x a signifies the effect of temperature, which damages BEC, while (∂E/∂x a ) T stands for the effect of chemical potential, which keeps BEC.Consequently, there is a contest mechanism between the first and second terms in Equation (32).If the two effects equalize, the system possesses a critical boson number N c , below which the 2D system goes into the normal state, but above which the 2D system enters the BEC state.
Afterwards, the reduced chemical potential is employed in Equation (32) to acquire the heat capacity at constant number.C N is a function of temperature T, boson number N, and anisotropy parameter p.If one evaluates C N , one must associate Equation (18) with Equation (32).Based on Equation (32), changing of the scaled heat capacity C N /k B with the boson number N is presented in Figure 10 for various T. We select the anisotropy parameter p = 0.1.An interesting discovery is that, at T = 1 nK, C N /k B ≈ 26.0 for all N. Figure 10 demonstrates that, for a finite number of 2D bosons, the system possesses an approximate critical boson number N c .If T is fixed, the heat capacity C N is an ascending function of the boson number N when N is smaller than N c .As N ≥ N c , C N /k B corresponds to a constant dependent of T. According to Equation (32), changing of the scaled heat capacity C N /k B with the temperature T is presented in Figure 11 for various N.We select the anisotropy parameter p = 0.1.Obviously, C N = 0 at T = 0 K. Figure 11 demonstrates that, for p = 0.1 and at a fixed N, there is an inflection point on the curve of C N versus T, which corresponds to the critical temperature T c .If N is fixed, the heat capacity C N is a fast ascending function of temperature T when T is smaller than T c , but at T ≥ T c , it is a slowly increasing function of temperature T. In Figure 12, we also display changing of the scaled heat capacity C N /k B with the temperature T for various p at N = 1000.Figure 12 reveals that, for p ≥ 0.5, the system possesses a critical temperature T c at which C N (T) attains a maximal value.The specific-heat maximum is a monotonically increasing function of anisotropy parameter p.The circumstance of phase transitions grows clearer and clearer as the anisotropy parameter approaches unity.The specific-heat maximum represents an important marker in this transition.When T T c , C N /k B = 2N, which corresponds to the classic limit.According to Equation (32), in Figure 13 we show the variation of the scaled heat capacity C N /k B with the anisotropy parameter p for various N.The temperature is fixed at T = 50 nK.It is interesting to note that, for N = 10,000, C N possesses the minimum at p = 0.04 and the maximum at p = 0.51.For N ≤ 1000, C N is a monotonically descending function of anisotropy parameter p.A fine character of the accurate solutions in the above is that they are effective for arbitrary T, N, and p.

Critical Temperature and Reduced Chemical Potential
In the following, we intend to study the thermodynamic limit in the 2D anisotropic harmonic potential.At first, one must observe that the thermodynamic limit in the anisotropic potential corresponds to taking N → ∞ and √ ω x ω y → 0 with the product Nω x ω y kept constant [24].Then, we can adapt Equation (18) in the form, where N 0 signifies the boson number in the ground state n x = n y = 0. We all know that, as N → ∞, x a = 1.It can be verified that the above program of the thermodynamic limit in the anisotropic potential amounts to letting N 0 = 0 and x a = 1 in Equation ( 43).The operation leads to the following expression for the critical temperature T c , where q xc = exp(−hω x /k B T c ).In the trap thermodynamic limit, one can derive the solution of Equation (43) as Even though Equations ( 44) and (45) are acquired in the condition of N → ∞, they are effective for arbitrary T and large N (N ≥ 10 3 ).As 10 2 ≤ N < 10 3 , Equations ( 44) and (45) are acceptable too.In reality, one first uses Equation (44) to resolve the transition temperature T c and then employs Equation (45) to resolve the reduced chemical potential x a .The critical temperature T c given by Equation ( 44) is an implicit function of boson number N and anisotropy parameter p.As p is fixed, the critical temperature T c is a monotonically ascending function of boson number N. As N is fixed, the critical temperature T c is a monotonically descending function of anisotropy parameter p.When the number of harmonically trapped 2D bosons is finite, the properties of phase transitions at some critical temperature are very clear too.The above investigation obviously demonstrates that the properties of phase transitions in the anisotropic potential are more vivid than those in the isotropic potential.

Internal Energy, Entropy, and Helmholtz Free Energy
The statistical mechanics of phase transitions can portray phase transitions on the basis of thermodynamic functions.At critical temperatures, the Helmholtz free energy of the two phases should be successive.Nevertheless, phase transitions could be sorted into two kinds based on the properties of derivatives of the Helmholtz free energy.When the first derivatives of the Helmholtz free energy are not successive, the accompanied phase transitions are known as first-order phase transitions.When the second derivatives of the Helmholtz free energy are not successive, the accompanied phase transitions are known as second-order phase transitions.
In the anisotropic trap, the Helmholtz free energy F is a function of temperature T, boson number N, and anisotropy parameter p, that is, F = F(T, N, p).The entropy S is linked to the first derivative of F: Accordingly, when first-order phase transitions occur, the entropy of two phases a and b is discontinuous, S a = S b .When a transition between two phases a and b happens, the latent heat L is presented by L = T(S a − S b ).A first-order phase transition possesses two significant features: at the transition, the system has a latent heat, and the system makes a jump in the entropy.From gas to liquid phase, from liquid to solid phase, and from gas to solid phase, these transitions are all first-order phase transitions.The heat capacity C N at constant N is linked to the second derivative of F: Accordingly, when second-order phase transitions occur, the heat capacity of two phases a and b is not successive, C a = C b .A second-order phase transition possesses two significant features: at the transition, the system has no latent heat, and the system makes a jump in the heat capacity.When there is no magnetic field, the transition between superconductive and normal phases of a metal corresponds to a second-order phase transition.
At this point, we start to depict the transition behavior of ideal bosons in a 2D anisotropic harmonic potential.Based on Equations ( 26), (44), and (45), Figure 14 reveals the internal energy E versus the temperature T for N = 10 4 and p = 0.5.Owing to Equation (44), we discover that, at N = 10 4 and p = 0.5, T c = 51.90nK.In Figure 14, as the temperature T traverses T c = 51.90nK, the internal energy E is successive and flat.Computations also demonstrate that, at T c = 51.90nK, both the entropy S and the free energy F are successive and flat and do not reveal any clue of phase transitions.

Heat Capacity at Constant Number
In what follows, we can inspect the transition properties of the heat capacity at constant number.As T ≤ T c and since x a = 1, (∂x a /∂T) N = 0.As T > T c , (∂x a /∂T) N is resolved by the partial derivative of the expression H q x (x a ) = N.As a result, we obtain the heat capacity at constant N as follows: As T > T c , one can acquire the partial derivative of x a as follows: When a phase transition happens, the heat capacity of two phases a and b is not successive, C a = C b .Therefore, at transition temperature T c , the system makes a jump in the heat capacity and this jump is presented by Since, at T c , the heat capacity of the BEC state is larger than that of the normal state, this leap takes place.As a result, the transition between normal and BEC states corresponds to a second-order phase transition.
When calculating C N in the thermodynamic limit, we must combine Equation (48) with Equations ( 44) and (45).According to Equations (44), (45), and (48), Figure 15 demonstrates the heat capacity C N against the temperature T for N = 10 4 and 10 5 at p = 0.5.From Equation (44), we find that, at N = 10 4 , 10 5 and p = 0.5, T c = 51.90 and 166.07 nK, respectively.The cusp singularity of a second-order phase transition of BEC in 2D anisotropic traps is displayed in Figure 15.Such a critical property resembles the λ-transition of liquid helium-4.An up-to-date BEC experiment with 2D ideal photons has showed this cusp singularity of the specific heat [29].Consequently, our program for the thermodynamic limit is totally proper.Owing to Figure 15, we can find that the heat capacity of the BEC state is extremely larger than that of the normal state.At critical temperature T c , this leap in the heat capacity is immense.Further, Figure 16 reveals the heat capacity C N against the temperature T for N = 10 4 and p = 0.3.Owing to Equation (44), one can discover that, at N = 10 4 and p = 0.3, T c = 66.69 nK.Comparing Figure 15 with Figure 16, one can see that, for a fixed N and when p is lowered, the singular peak of the specific heat becomes very gentle.Many properties of BEC of 2D anisotropic traps are exhibited in the present paper.

Conclusions
This thesis investigates the thermodynamic behavior of ideal bosons in 2D anisotropic harmonic traps.In this study, we find that the issue of BEC in 2D anisotropic harmonic potentials can be figured out analytically.The analytical solution is connected with an analytical expression, which includes a series expansion of q-digamma functions.The q-digamma function was ushered in mathematics thirty years before and, at present, plays an important role in science and technology.For this probe, we usher a new thermodynamic function that is known as the reduced chemical potential to substitute the fugacity.In the issue of BEC of harmonically trapped 2D bosons, the fugacity z is a bad thermodynamic function, but the reduced chemical potential x a is a fine thermodynamic function.We construct a hypothetical model of quantum thermodynamics to evaluate various thermodynamic functions of ideal bosons in 2D anisotropic harmonic traps.These thermodynamic quantities are functions of temperature T, boson number N, and anisotropy parameter p.The graphs of their variation with T, N, and p are given.
An interesting observation is that, when T and N are fixed, the entropy of a finite number of ideal bosons is a monotonically ascending function of anisotropy parameter p but the free energy is a monotonically descending function of anisotropy parameter p.The first major discovery is that the internal energy of a finite number of ideal bosons is a monotonically ascending function of anisotropy parameter p.The second major discovery is that, when p ≥ 0.5, the changing with temperature of the heat capacity of a finite number of ideal bosons possesses the maximum value, which happens at critical temperature T c .The third major discovery is that, when 0.1 ≤ p < 0.5, the changing with temperature of the heat capacity of a finite number of ideal bosons possesses an inflection point, but when p < 0.1, the inflection point disappears.The fourth main finding is that, for a fixed T and when N < 10 4 , the heat capacity of a finite number of ideal bosons is a monotonically descending function of anisotropy parameter p, but when N ≥ 10 4 , C N possesses a minimum and a maximum.Therefore, a finite number of ideal bosons in 2D anisotropic harmonic traps possess larger heat capacity C N than those in 2D isotropic harmonic traps.The fifth major discovery is that, in the thermodynamic limit, at T c and when p ≥ 0.5, the heat capacity at constant number reveals a cusp singularity, which resembles the λ-transition of liquid helium-4.The sixth major discovery is that, in comparison to 2D isotropic harmonic traps (p = 1), the singular peak of the specific heat becomes very gentle when p is lowered.Since the observation of the BEC of ultracold dilute atomic gases in 1995, one has thought that the thermodynamic properties of 2D anisotropic harmonic traps can be evaluated exactly.This paper accomplishes this task successfully.
This thesis has studied phase transitions of ideal bosons in 2D anisotropic harmonic traps in the thermodynamic limit.Nevertheless, one can find that, when the boson number in 2D anisotropic harmonic traps is finite, the properties of phase transitions at some critical temperature are very clear as well.This thesis demonstrates that, in the trap thermodynamic limit, the system possesses an accurate critical temperature T c settled by a strict expression, which includes a series expansion of q-digamma functions.The calculation displays that, when p ≥ 0.5, the 2D critical temperature T c is identical to the specific-heat maximum of a finite number of harmonically trapped 2D bosons.Calculations also display that, at T c both the entropy S and the free energy F are successive and flat and do not reveal any clue of phase transitions.However, one can discover that, at T c , the heat capacity at constant number is not successive.Therefore, at T c , the system possesses a leap in the heat capacity at constant number.The leap in the heat capacity at constant number is enormous.Since, at T c , the heat capacity of the BEC state is much larger than that of the normal state, this leap takes place.The dissection of 2D anisotropic harmonic traps sheds more light on the behavior of BEC.
To sum up, the boson system in the BEC state is comprised of condensed bosons and uncondensed bosons.We ascertain that, in 2D anisotropic traps, the number of uncondensed bosons is represented by an analytical expression, which includes a series expansion of q-digamma functions in mathematics.One can utilize this analytic result to evaluate various thermodynamic functions of ideal bosons in 2D anisotropic harmonic traps.The first major discovery is that, in the thermodynamic limit, at T c and when p ≥ 0.5, the heat capacity at constant number reveals a cusp singularity, which resembles the λtransition of liquid helium-4.Another major discovery is that, in comparison to 2D isotropic harmonic traps (p = 1), the singular peak of the specific heat becomes very gentle when p is lowered.In this probe, we usher a new thermodynamic function that is known as the reduced chemical potential to substitute the fugacity.

Figure 1 .
Figure 1.A plot of the scaled internal energy E/hω x versus the boson number N in a 2D ideal boson gas.The plot is drawn for T = 1, 15, 50, 100 nK.We select the axial angular frequency ω y /2π = 10.0 Hz and the anisotropy parameter p = 0.5.

Figure 2 .
Figure 2. A plot of the scaled internal energy E/hω x versus the temperature T in a 2D ideal boson gas.The plot is drawn for N = 50, 200, 1000, 10,000.We select the axial angular frequency ω y /2π = 10.0 Hz and the anisotropy parameter p = 0.5.

Figure 3 .
Figure 3.A plot of the scaled internal energy E/hω x versus the anisotropy parameter p in a 2D ideal boson gas.The plot is drawn for T = 50, 55, 60, 65 nK and N = 10 4 .We select the axial angular frequency ω y /2π = 10.0 Hz.

Figure 4 .
Figure 4.A plot of the scaled entropy S/k B versus the boson number N in a 2D ideal boson gas.The plot is drawn for T = 1, 15, 50, 100 nK.We select the axial angular frequency ω y /2π = 10.0 Hz and the anisotropy parameter p = 0.2.

Figure 5 .
Figure 5.A plot of the scaled entropy S/k B versus the temperature T in a 2D ideal boson gas.The plot is drawn for N = 50, 200, 1000, 10,000.We select the axial angular frequency ω y /2π = 10.0 Hz and the anisotropy parameter p = 0.2.

Figure 6 .
Figure 6.A plot of the scaled entropy S/k B versus the anisotropy parameter p in a 2D ideal boson gas.The plot is drawn for T = 50, 55, 60, 65 nK and N = 10 4 .We select the axial angular frequency ω y /2π = 10.0 Hz.

Figure 7 .
Figure 7.A plot of the scaled free energy −F/hω x versus the boson number N in a 2D ideal boson gas.The plot is drawn for T = 1, 15, 50, 100 nK.We select the axial angular frequency ω y /2π = 10.0 Hz and the anisotropy parameter p = 0.2.

Figure 8 .
Figure 8.A plot of the scaled free energy −F/hω x versus the temperature T in a 2D ideal boson gas.The plot is drawn for N = 50, 200, 1000, 10,000.We select the axial angular frequency ω y /2π = 10.0 Hz and the anisotropy parameter p = 0.2.

Figure 9 .
Figure 9.A plot of the scaled free energy −F/hω x versus the anisotropy parameter p in a 2D ideal boson gas.The plot is drawn for T = 50, 55, 60, 65 nK and N = 10 4 .We select the axial angular frequency ω y /2π = 10.0 Hz.

Figure 10 .
Figure 10.A plot of the scaled heat capacity C N /k B versus the boson number N in a 2D ideal boson gas.The plot is drawn for T = 1, 15, 50, 100 nK.We select the axial angular frequency ω y /2π = 10.0 Hz and the anisotropy parameter p = 0.1.

Figure 11 .
Figure 11.A plot of the scaled heat capacity C N /k B versus the temperature T in a 2D ideal boson gas.The plot is drawn for N = 50, 200, 1000, 10,000.We select the axial angular frequency ω y /2π = 10.0 Hz and the anisotropy parameter p = 0.1.

Figure 12 .
Figure 12.A plot of the scaled heat capacity C N /k B versus the temperature T in a 2D ideal boson gas.The plot is drawn for p = 0.05, 0.1, 0.5, 1.We select the axial angular frequency ω y /2π = 10.0 Hz and the boson number N = 1000.

Figure 13 .
Figure 13.A plot of the scaled heat capacity C N /k B versus the anisotropy parameter p in a 2D ideal boson gas.The plot is drawn for N = 50, 200, 1000, 10,000.We select the axial angular frequency ω y /2π = 10.0 Hz and the temperature T = 50 nK.

5 Figure 14 .
Figure 14.In the thermodynamic limit, changing of the scaled internal energy E/hω x with the temperature T. We select the axial angular frequency ω y /2π = 10.0 Hz, the boson number N = 10,000, and the anisotropy parameter p = 0.5.

5 Figure 15 .
Figure 15.In the thermodynamic limit, changing of the scaled heat capacity C N /Nk B with the temperature T at N = 10 4 and 10 5 .We select the axial angular frequency ω y /2π = 10.0 Hz and the anisotropy parameter p = 0.5.

Figure 16 .
Figure 16.In the thermodynamic limit, changing of the scaled heat capacity C N /Nk B with the temperature T at N = 10 4 and p = 0.3.We select the axial angular frequency ω y /2π = 10.0 Hz.