Nonlinear Klein-Gordon equation and the Bose-Einstein condensation

The interest in the Klein-Gordon equation with different potentials has increased in recent years due to its possible applications in Cosmology, Hadron Physics and High-Energy Physics. In this work we investigate the solutions of the Klein-Gordon equation for bosons under the influence of an external potential by using the Feshbach-Villars method. We present detailed results for two cases: the Coulombic potential and the harmonic potential. For the latter case, we studied the effects of self-interacting particles by adopting a mean-field approach. We show that our results converge smoothly to the solution of the Schr\"odinger equation for the same systems as the relativistic effects diminish.

Introduction: The Klein-Gordon (KG) equation has been investigated by various methods due to its application in different fields. The case of a scalar field has an increasing interest because of its connections with the axion, a hypothesized particle that may have a relevant role in Astrophysics and in studies of Dark-Matter [1,2,3,4].
In Astrophysics, the investigations of a possible bosonic phase in massive stars has triggered the study of the solutions of the relativistic equation with an external potential. The possibility of a condensate regime inside the stars has been considered. Several methods have been used to solve the equation for different potentials [5,6]. In the case of self-interacting particles, the problem becomes particularly difficult, and in Ref. [7] the authors claim that it is not possible to find a solution for this case. Nevertheless, they investigated the axion star and the possibility of the formation of a condensate.
Systems formed by interacting bosons are also of interest in High-Energy Physics and in Hadron Physics, where the possibility of a gluonic field in the condensate regime has been considered [8,9,10]. With the advent of colliders such as the EIC, FAIR and NICA, with the capability to investigate the hadron structure in details, we can expect important advances in our understanding of the complex hadronic structure. The effects of the interaction of the components of the system with an external potential has been extensively studied in connection with the Bose-Einstein condensate [11,12]. The interaction can favour or disfavour the formation of the condensate. The correlations and interactions among the particles do not prevent the formation of the condensate [13], and can lead to non-extensive forms of the condensate described in terms of the Tsallis statistics [14,15]. The characteristics of the non extensive condensate has been described in Ref. [16,17].
In this work, we investigate the solutions of the KG equation for a system of bosons under the influence of an external potential, in particular for the Coulombic potential and for the harmonic potential, and then we include the self-interaction through a mean-field approach. We investigate the ground-state solution and obtain the chemical potential for the case of the external harmonic potential, and with the self-interaction regulated by the scattering length. We use the Feshbach-Villars formalism, which allows us to evaluate the transition from the relativistic to the non-relativistic regimes. This is particularly useful in the case of self-interacting particles in a harmonic potential, for which the Schrödinger equation is known [18,19]. We show that our results with the KG equation smoothly approach those for the Schorödinger equation as the systems become non-relativistic.
The Feshbach-Villars formalism: The KG equation for free particles is This equation can be written in a convenient form by using the Feshbach-Villars formalism [20] (see also Ref. [21]). In this case, we define the spinor where ϕ(t, r) and χ(t, r) are wave-functions related to ψ(t, r) by In the Feshbach-Villars formalism, the KG equation is written for the spinor as a Schrödinger-like equation, i.e. where In the equation above σ i (i = 1, 2, 3) are the Pauli matrices, andˆ p = −ih ∇. The conserved charge and current are given by respectively. The spinor representation is equivalent to the coupled differential equations The Feshbach-Villars formalism becomes especially convenient to study the non-relativistic regime, in which the relativistic wave equation reduces to the solution of the Schrödinger equation. We will study below some solutions of the KG equation for different potentials. The time-independent form of the KG equation is and the corresponding solution of the time-dependent equation, Eq. (4), is Introduction of the chemical potential: In the present work we are interested in the investigation of the Bose-Einstein condensate formed in a relativistic system. We include the chemical potential, by assuming that the ground-state energy coincides with the chemical potential. The stationary-state solution is where ± corresponds to the positive(negative) chemical potential solution. Then, the original KG equation, Eq. (1), leads to a time-independent KG equation which writes whose solution in the spinor representation is where the spinor  24) and (33)) We display the wave functions for the first three levels: n = 1, 2, 3. Middle panel: Reduced wave functionū(r) and the spinor components (ū ϕ (r),ū χ (r)) for n = 1. In both panels we have considered λ = 0.1 and we have normalized the wave functions such that u(r) r→0r β + (cf. Eq. (26)). Right panel: −δ M n as a function of λ for the first three levels (cf. Eq. (35)).

has the amplitudes
Notice Introduction of an external potential: The Bose-Einstein condensation can happen even in the presence of an external potential.
To study the effects of an external potential, let us solve the time-independent KG equation in different cases. This equation has been obtained and studied in the literature in several contexts, see e.g. Refs. [22,23,24,25,26,27,28,29,30]. It can be derived by applying the correspondence principle in quantum mechanics: E → ih ∂ ∂t , p → −ih ∇; to the relativistic expression for the energy E = p 2 c 2 + m 2 c 4 +V ( r). Using the Feshbach-Villars formalism we obtain the relation with ψ = ϕ + χ. Then, we find from Eq. (16) that the spinor components ϕ and χ can be computed from the wave function ψ as follows In the massless case, one has to consider m = 0 in Eq. (15), and replace m → N in Eqs. (16) and (17), where N is an arbitrary real nonzero parameter [31]. The Schrödinger equation can be obtained from the KG equation (15) by considering the non-relativistic limit, i.e. c → ∞. To see that, let us define the energy with the rest-mass energy subtracted,μ, as Then, we can write the KG equation (15) in a form similar to the Schrödinger equation, but with an explicit relativistic correction. This is which is completely equivalent to Eq. (15). Then the Schrödinger equation is obtained from the KG equation in the limit |V ( r) − µ| mc 2 . Notice that in this limit one has from Eq. (17) that ϕ( r) c→∞ ψ( r) and χ( r) c→∞ 0.
Coulombic potential: Let us study the case of a Coulombic potential, We consider Eq. (15) in spherical coordinates, and assume = 0 (zero angular momentum), obtaining the equation This equation can be written in a simpler form by introducing the reduced radial wave function Then one gets The solution of Eq. (23) with the potential of Eq. (20) can be obtained analytically as a linear combination of two independent solutions, i.e.
where L a n (z) and U(a, b, z) are the generalized Laguerre polynomial and the confluent hypergeometric function, respectively, while we have defined ∆ ≡ m 2 c 4 − µ 2 and The solution at short distances behaves as with +(−) standing for the dominant behavior of the L a n (z) (U(a, b, z)) solution. The physical (regular) solution corresponds to r β + , so that the coefficient C − should vanish and then we should consider C 2 = 0 in Eq. (24). From the behavior of the L-solution at large distances, one finds that the square-integrable condition for the wave function leads to the exact result for the eigenvalue problem of Eq. (23) which we will denote by µ KG n , corresponding to the solutions of the equation leading to the chemical potential Notice that the form of Eq. (28) implies that |µ KG n | < mc 2 , so that the eigenfunction associated to the eigenvalue µ KG n behaves as u(r) ∼ r→∞ e −κr with κ ≡ ∆ n /(hc) > 0, thus corresponding to a bound state. To consider the non-relativistic limit, we make c to be large enough. Then we have We can compare our results in the non-relativistic limit with the results obtained with the Schrödinger equation. In this latter case, the eigenvalues areμ NR n = − mg 2 2h 2 n 2 , n = 1, 2, 3, · · · .
Let us define δ µ n as the difference between the result for the eigenvalue obtained with the KG equation (after subtracting the rest-mass energy), and the non-relativistic result of Eq. (30), i.e.
µ KG n = mc 2 +μ NR n + δ µ n . In the limit of c 1, the expression of Eq. (28) has the following expansion so that one can easily identify the correction δ µ n at order O(c −2 ) from a comparison with Eq. (31). Before displaying the results, it is convenient to define the dimensionless variables Then, the dimensionless eigenvalues are functions only of the parameter λ . In particular, the decomposition of Eq. (31) writes in dimensionless units as The non-relativistic limit corresponds to λ → 0, or equivalently c → ∞. We display in Fig. 1 (left and middle panels) the dimensionless reduced wave functionū(r), as well as the spinor componentsū ϕ (r) ≡ hc mc 2 · rϕ(r) andū χ (r) ≡ hc mc 2 · rχ(r). We display also in Fig. 1 (right)  showing that the result of the Schrödinger equation is recovered in the non-relativistic limit.
Finally, let us note that the eigenvalue µ KG n is real only for real values of β + , and this happens for i.e. for λ < 1/2. One can see that in the regime g >hc/2, the asymptotic behaviors of the reduced wave function are so that while it is a square-integrable function, it is not a regular solution at r → 0.
System of N interacting particles: Another interesting case that we will study is a system of N interacting particles under the influence of an external potential. In the non-relativistic limit, such a system is usually described in terms of the Gross-Pitaevskii nonlinear Schrödinger equation, which is given by In Ref. [18] (see also [19]) this system was studied with an harmonic external potential and a mean-field potential for the interaction term given in terms of a scattering-length, a, that is, Here we study this case in the relativistic regime, and evaluate its non-relativistic limit comparing our results with the results from Ref. [18]. In the relativistic case, from Eq. (23) and using the potential of Eq. (39), one has Let us define the dimensionless variables Using Eq. (18) one finds β = 1/α +β . Then, Eq. (40) can be written in the equivalent form The non-relativistic limit corresponds to α → 0, or equivalently c → ∞. In the non-relativistic case, we define the normalization as N = d 3 r|ψ(r)| 2 = 4π ∞ 0 dr r 2 |u(r)| 2 where N is the number of atoms or partons. The corresponding normalization for the dimensionless wave function is The length s is the characteristic length of the system. In the relativistic case, it is better to define the normalization in terms of the conserved charge, i.e.
which for the dimensionless spinor components of the wave function Φ f (x) ≡ 8π|a| · r f (r) , f = ϕ, χ, leads to Obviously Q −→ c→∞ N and q −→ c→∞ n. We show in Fig. 2 the results ofβ and q by solving the Gross-Pitaevskii nonlinear equation in the relativistic case, Eq. (42), and compared with the non-relativistic result (α = 0) of Ref. [18]. Even though the curves of Fig. 2 contain all the relevant information, it would be interesting to show explicitly the dependence ofβ with the scattering length a, by using thatβ (q) =β (Q|a|/ s ). To do so, we display in Fig. 3 the results ofβ as a function of Q for several values of |a|/ s (left panel), and as a function of |a|/ s for several values of Q (middle panel). One can see that for fixedβ , the larger the value of |a|, the smaller the value of Q. This can be easily understood by noticing that the relationβ =β (q) given by Fig. 2 (middle) can be inverted to give the one-valued function q = q(β ), so that using Eq. (45) one finds Q = s q(β )/|a|. Notice that the normalization ofβ given by Eq. (41) doesn't contain any dependence in a, so that the dependence ofμ with a is proportional to the dependence ofβ with a, i.e.μ(a) =hω ·β (Q|a|/ s ). Finally, we display in Fig. 3 (right) the dependence of N/Q with the parameter α for Φ (0) = 2.5, where N is computed in the relativistic case by using Eq. (43) with Φ(x) = Φ ϕ (x) + Φ χ (x) the solution of Eq. (42). Notice that in the non-relativistic limit (α → 0) both quantities are equal, but in the relativistic regime one finds N < Q.
Discussion and conclusions: Concluding, in this work we analysed the time independent solutions of the non-linear Klein-Gordon equation. We used the Feshbach-Villars approach with potential, and analysed the behaviour of the system under different external potentials, namely, the Coulombic potential and the harmonic potential.
With the introduction of the chemical potential as the eigenvalue of the Hamiltonian operator for the ground-state, we obtained the relevant conditions for the formation of the Bose-Einstein condensate in a relativistic system. The use of the Feshbach-Villars formalism allowed us to easily investigate the transition of the system from the relativistic to the non-relativistic regimes. We compare the results obtained for the relativistic system with a mean-field potential and an harmonic external potential, with those for the non-relativistic system under the same conditions. We observe a smooth transition of several parameters as the relativistic effects diminish, for both the Coulombic potential and the harmonic potential.
The effects of the interaction between the particles of the system, described by the scattering length in a mean-field approach, is analysed, and in particular, we study the variation of the chemical potential with the scattering length.
The study of the nonlinear Klein-Gordon equation to investigate condensation of relativistic particles invites us to pursue extensions of this work. One possibility is to use the same formalism of the present work to investigate condensation in the relativistic hadronic matter, and another option is to study a modified form of the Gross-Pitaevskii equation with fractional derivatives.