Spin-0 bosons near rotating stars

In this work we study the effects of rotating stars on the behavior of bosons close to their surfaces. For this task, metrics determined by the rotation of these stars will be taken into account. We will consider the Klein–Gordon equation in the Hartle–Thorne metric and in the one proposed by Berti et al. by considering some kinds of stars. Pions and Higgs bosons will be investigated as examples. By solving the equations, the main effect that may be observed is a rotational phase, δR\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta _R$$\end{document}, that depends on the angular momentum of the star, J, in the wave function of the boson. Corrections of higher orders in J are also investigated.


I. INTRODUCTION
The observation of the effects of the general theory of relativity in physical systems has always been a challenging task.Since the first observation of light deflection by the Sun or the effect of frame dragging due to the rotation of the earth [1], [2], the Lense-Thirring effect, which, after decades of experiments, was finally observed by the Gravity Probe B [3], or even the recent detection of gravitational waves [4], we have witnessed great experimental efforts in order to obtain reliable results.
In this work, we will study the effect of the dragging of the spacetime due to the rotation of stars in quantum systems, and for this purpose, we will use spin-0 bosons as test particles.
A description of the spacetime of "slowly" rotating relativistic stars was constructed during the 1960s by James B. Hartle and Kip S. Thorne [23], [24], for both the exterior and interior regions of a star.The limits of the approximation, or in other words, the "slowness" of the stars, are in fact not so narrow when considering physical stars, such that this criterion is not valid solely in the slow-rotation realm, yet it seems as if the solution is accurate for many systems, even for some that would not have been considered slow then.Therefore, in this paper, we will use the external Hartle-Thorne metric as given by [24] and also the corrections proposed in [25], and then write a Klein-Gordon equation.We will solve the equations and then show numerical results for some typical systems.This paper is divided into five sections beyond the introduction.In section II we will describe the line elements we will be using in our calculations and in Sec.III we will sketch the process by which we take the line element and expand it up to first order in the star's angular momentum  .Then, we will write the Klein-Gordon equation and solve it in various ways, in terms of several different special functions.In section IV, we will show an expansion that contains terms up to second order of the angular momentum  and the resulting power series solution to the Klein-Gordon equation obtained in this spacetime.
The resulting solution corresponds to the wave function for a free particle plus additional mass and angular momentum corrections.We will also compare this solution with the one obtained in section III, thereby showing that they have intrinsic similarities, despite having been obtained in conceptually different ways.In Sec.V the solution will then be studied for some physical systems with different values of the mass and angular momentum and in Sec.VI we will draw the conclusions of this work.

II. THE TWO METRICS
In this section, we will study the spacetime structure by means of metrics that take into account the rotation of stars.For this purpose we will consider the Hartle-Thorne metric [23], [24] and the metric proposed by Berti et al. [25].Since we are only interested in studying observable effects of bosons near rotating stars, we will only consider the solutions of the equations in an external region.
The Hartle-Thorne metric is a well-known spacetime metric that describes the geometry determined by slowly rotating stars.It was first obtained in a sequence of papers in the late 1960s and it has been used ever since in a variety of important works.The metric was initially constructed out of a very general canvas, that of an axially symmetric, stationary system that rotates uniformly and "slowly".This "slow" rotation refers to the fluid's angular velocity Ω, which must be slow enough for the changes in pressure, energy density, and gravitational field due to rotation to be small, much less than unity.These considerations lead to the criterion presented by the authors which determines a scale factor for a star of mass  and radius  that also means Ω ≪ , non relativistic particles inside the star.We note, however, that this criterion is actually wide-reaching, and, as opposed to what may be thought at first glance, it does not necessarily correspond to a slow rotation when compared with the values of real stars [26].
A general line element may then be proposed in order to study a rotating object by considering it as an axially symmetric stationary system in terms of spherical coordinates where  , , , and  are functions of  and  .Rotational effects are introduced as perturbations in these functions, and they become [24] where  is the star's mass quadrupole moment,  its angular momentum, and are associated Legendre functions of the second kind.In [27] it was pointed out that this metric is not consistently truncated, and for this reason, it presents some errors.A corrected metric was then proposed, but Berti et al. [25] studied that new metric and found some minor sign errors.They proposed the solution that will be the second metric we will be using in this paper, which may be considered a corrected version of the original Hartle-Thorne metric.From [25], the corresponding line element is given by where (10) and  = Ω/Ω * with Ω * = √︁ / 3 .Observing these results, it seems reasonable to consider both metrics in the study of quantum systems, and since the corrections proposed in [25] are very small, we will check whether they are important in the behavior of spin-0 bosons near some kinds of stars.

III. KLEIN-GORDON EQUATION IN THE HARTLE-THORNE METRIC
In this section we will study the Klein-Gordon equation in the metrics presented in Sec.
II up to first order in the angular momentum ratio  1 =  / 2 .We must remark that it is a suitable quantity for the analysis that we intend to do in this work since it presents small numerical values for the stars.For the sun, taking  =  for example,  1 = 9.7 × 10 −13 ; for a typical white dwarf,  1 ∼ 10 −10 ; and for a typical neutron star,  1 ∼ 10 −5 .As we can see, in a first approximation we may neglect terms of the order  2 1 since they become numerically very small in the region outside of the star.In this section, we will work out the solutions up to first order in  1 , and in the next section, we will study the corrections up to  2 1 .In [24] the line element (2) was expanded in powers of / , which, for the external region of the Sun, for example, is smaller than / ∼ 2 × 10 −6 , and when using numerical data relative to the Sun, an approximate form of this line element with an accuracy of the order of 1 part in 10 15 has been proposed We will study stars with nearly symmetric mass distributions by taking  = 0, and since  2 / 6 ≪ / for the external line element, we may write which is the metric that will be considered in the calculations.
The Klein-Gordon equation for an arbitrary spacetime may be obtained by replacing the regular derivatives in the usual Klein-Gordon equation with covariant derivatives [28], [29].
In so doing, we arrive at the following: Because the components of our metric are not functions of  or , we can apply the method of separation of variables in the equation and obtain where  = 0, ±1, ±2..., but we still need to solve the equations for ( ) and Θ( ).The whose solution is simply the set of associated Legendre functions    ( ) with  ′ =  ( + 1), where  is an integer.The radial equation is where ( ), making the change of variables  =  2 , and rearranging fractions, we arrive at the following equation where This equation is the normal form of the confluent Heun equation, the canonical form of which is given by where  () =  (, , , , , ) are the so-called confluent Heun functions with By means of an integrating factor, it is possible to relate the equations ( 17) and (19), and thus each of their coefficients, Having obtained these expressions, we have thus determined the exact analytical solution of the Klein-Gordon equation in the Hartle-Thorne metric given by eq. ( 12).It is possible to find a solution directly for eq.( 19) as a power series.We will use the asymptotic expansion for large  found in [30], which was based, among others, on [31], a book specialized in Heun's differential equations.In the vicinity of the irregular singularity at infinity, the two solutions of equation ( 19) can be written as 2 (, , , , , ) ∼ valid in the domain given by the following three-term recurrence relation: where assuming  0 = 1, and The negative power series given by ( 22) do not necessarily converge, and so we are required to truncate them at an appropriate spot.To choose said spot, we assume  +1 = 0, for some , and then take  0 = 0, which guarantees that, starting from a certain  ′ ,  +2 will be 0 for any  >  ′ .Doing this procedure allows us to discover a relation that places a constraint on the coefficients of the confluent Heun equation (19).
By calling to mind how these coefficients relate to the physical quantities defined in the Hartle-Thorne metric (12) and in the wave function ansatz (14), it is clear that eq. ( 30) defines an energy spectrum,  (), that may be solved numerically.
We may notice that some terms in eq. ( 16) are still very small, which makes it possible to obtain a simpler solution with the purpose of acquiring a better understanding of this system's physics by writing the solution in terms of more common special functions.With this purpose in mind, we can write eq. ( 17) as where By observing that  ′ / 2 and  ′ / are, numerically, several orders of magnitude smaller than  ′ /( −2) 2 or  ′ /( −2), for fixed parameters and  values (outside the stars), a reasonable approximation to be made is to neglect these terms in eq. ( 31).If we also perform a change of variables  = 2 ( − 2), the equation becomes where  =  ′ /(2) e  = √︁ 1/4 −  ′ and it is still relative to the line element (12), but with the approximations shown previously.The equation above is the so-called Whittaker equation, a modified version of the confluent hypergeometric equation, and it can be written in the form of Kummer's equation through a transformation of the function.This equation has two solutions, the Whittaker functions  , () and  , (), which can be written, respectively, in terms of a Kummer function  (, , ) and a Tricomi function  (, , ), which are themselves solutions to the confluent hypergeometric equation [32], Since the term proportional to  ′ is still much smaller than the remaining terms of eq. ( 31), and since in the region exterior to the stars  ≫ 2, we have  ≃  − 2, whereby the equation becomes, with good accuracy, Supposing a solution with the form and a change of variables  = 2 , we obtain the equation which has the solution Note that, for the physical problem we are studying,  2 = 0.The Kummer function may be defined as which asymptotically behaves like a Bessel function.Thus, for large values of the variable , which is precisely the region about which we are concerned, the solution may be written in the form or which explicitly shows the phase originated by the effect of the rotation of the star in the lowest order,   , which may be interpreted as a rotational phase.In the next section, we will investigate the solution for higher orders of  1 .

IV. EXPANSION OF THE HARTLE-THORNE METRIC TO SECOND ORDER IN J
Become some mistakes were discovered in the Hartle-Thorne metric given by eq. ( 2), since we are interested in studying terms of the order of  2 1 and higher in the solution we have to use the Berti et al. metric (9) in our calculations.It is interesting to note that, if we keep the terms up to first order in  1 and  = 0, eq. ( 11) is recovered, which means that the results presented in the last section are correct and the differences occur for higher terms in  1 .So, in this section, we will derive an equation taking into account the terms that appear in the corrected metric and write down a solution based upon an infinite series ansatz.In order to observe the behavior of the corrections in the solution an easier way to do that is to solve the equations for fixed values of the angle  .Taking these assumptions into account, the resulting radial equation becomes where After some manipulations, the equation may be written in the form where and This equation may be solved by the method presented in [33], where the solution has the form ( ) = e ± ( ) (48) and the exponent  ( ) is a sum of functions and a power series of / where the   coefficients must eventually be determined.It is possible to find the first few coefficients by simply substituting this ansatz back into equation ( 45), however, as  increases the expressions for the coefficients become increasingly larger.The first five coefficients are given by For large values of  we have which is the same result obtained in Sec.III, in a very different way, shown in eq.(41).

V. ANALYZING THE SOLUTIONS
In this section, we will investigate the solution (48) for several physical configurations and observe how it changes depending on the mass, angular velocity, and energy of the star or of the particle whose motion the Klein-Gordon equation purports to describe.We will also study high-energy and low-energy test particles.We will start by listing all physical configurations we chose in order to verify the size of the mass and angular momentum corrections.We selected three celestial bodies: the Sun [34], the neutron star PSR B1257+12 [35], [36], and the white dwarf PG 2131+066 [37].All numerical values used in this section are given in Planck units, where  =  = ℏ = 1.For the sun we have for the neutron star, and for the white dwarf, the corresponding values are We also chose two test particles with very different masses, the pion with mass   =0.139 GeV [34], and the Higgs boson with mass   = 125.1 GeV, which in Planck units are given by We separate eq.( 49) into four different parts: A is the real part of  containing only mass terms; B is the imaginary part of  containing only mass terms; Γ is the real part of  containing the angular momentum corrections, i.e., terms that are proportional to the angular momentum  or to its square  2 ; Δ is the imaginary part of  containing the angular momentum contributions.
These functions have the same structure as  ( ), given by eq. ( 49), and may be calculated by the sum of five terms determined by different dependencies on the variable  .Together, as they should, they completely make up (49).Then, by separating each type of contribution (mass and angular momentum) and also the real and imaginary parts of (49), we can see that these newly defined functions are written as and the coefficients inside of it can be found by separating the coefficients (50) in their real and imaginary parts and highlighting the ones with an angular momentum dependence, appearing in  2 ,  3 , and  3 .
We will show the results for the six combinations of stars and particles given by eq. ( 52)-(55) for particles with energies  = 2,  =  (particles near the surface),  = 1, and  = /4.These results are displayed in Tab.I-VI.
In the tables, each column is related to one kind of dependence on the variable  .Each line in each table below corresponds to one of the four functions previously defined: A(r), B(r), Γ( ) e Δ( ).Each column corresponds to the respective function's dependence on  , such that the first column corresponds to the contributions of  ( ) which are proportional to  , then to ln  , and so on.For example, if we look at Tab. I, the number 8.25031 × 10 23 found in the very first cell of the table, corresponds exactly to the term  .The tables above show the expected result that contributions diminish as the series goes and that the angular momentum contributions (corresponding to the functions Γ( ) and Δ( )) are present, albeit very small.As we can see, the star for which the corrections are more significant is the neutron star.However, even the largest angular momentum corrections are small if compared with the other terms.For this reason, we will continue this analysis considering neutron stars.By varying the physical quantities in question (such as the energy of the particle, and the angular momentum of the star) we will further explore our solution to the radial equation (45).We begin by defining three functions related to the solution which we will plot and which might give us some information on its status as a solution.These functions are The function  1 is precisely the exponent in the solution (48).The function  2 is eq.( 48), along with the integrating factor defined in (46); the end result is that  2 represents the radial solution of the wave function.The function  3 is defined as the difference between the exponent (49) and the mass contributions as defined in (56), divided by the absolute value of  .Since all that is left in the numerator of  3 is the sum of angular momentum contributions Γ( ) + Δ( ) as defined in (56), this function tells us the approximate order of magnitude of the angular momentum contributions when compared to the mass contributions.
Fig. 1 and Fig. 2 represent the function  1 in two regions, in Fig. 2 we show the results to the higher values of the energy.As we can see, the major difference between figures 1 and 2 is that the value of the function stabilizes as the energy increases and even big differences in energy account for little changes in  1 .For lower energies, even with small changes, the value of  1 changes considerably.In Fig. 3 and 4 we represent the function  2 , the radial wave function, in two regions according to the values of the energies of the pions.As we can see in Fig. 3, for large distances, far away from the surface of the star, the curves seem to coalesce, and this effect is especially visible for higher energies.
The last two graphs, Fig. 5 and 6, represent the function  3 , which shows the size of the angular momentum contribution as a variation in its value is done while maintaining the star's original mass and radius fixed.We observe that the largest angular momentum selected,  6 = 10, 000  , where   is the angular momentum of the neutron star as given in (53), most likely cannot be reached for a real star as it exceeds the maximum possible value given in the literature [26] for the angular velocity (and therefore angular momentum) of a neutron star.We used it regardless in order to show what happens to the solution for extreme rotations.and fixed angular momentum  =   .For small energies, the effect seems to be quite small, but it is possible that more terms in (49) may be needed to discover the full size of the low-energy contributions.
As we can see in figure 7, the angular momentum contribution is larger for higher energies.
Indeed, our solution seems to take the angular momentum into account more significantly when the star is rotating quickly (though still under the "slow" rotation regime as dictated by the original paper) and the particle is a high-energy zero-spin boson.We must remark that in the results presented in this section, the velocity of the particle   ≫   , where   is the rotation velocity of the star for  = .If we have small energies,   ∼   , more terms will be needed in eq. ( 49).
final state interactions of decaying particles or in the study of interacting particles.We must also remark that even if this kind of effect is small, an experimental verification would be proof of the frame-dragging by rotating objects and, consequently, of the theory of general relativity's influence on quantum systems.

FIG. 1 .FIG. 2 .FIG. 4 .
FIG.1.Absolute value of the real part of  1 for some values of the energy.For even small changes to the energy, the value of  1 changes substantially.

FIG. 5 .
FIG. 5. Imaginary part of  3 as a function of  for the actual angular momentum   of the neutron star and for variations of  .As previously stated,  3 gives an estimate of the size of the angular momentum corrections, and this graph shows the shape of the curves is independent of the value of the angular momentum, but its value changes as expected.

Fig. 5 FIG. 7 .
Fig.5 and 6show that the shape of the curves does not change considerably when varying the angular momentum, however, the absolute value of the function does.This is the expected result: the larger the angular momentum of the neutron star, the larger will be its contribution to the solution.We can also investigate what happens to variations of the energy  in the case where the angular momentum used was the neutron star's own.The results are shown in Fig.7.