The Scattering of Dirac Spinors in Rotating Spheroids

There are many stars that are rotating spheroids in the Universe, and studying them is of very important significance. Since the times of Newton, many astronomers and physicists have researched gravitational properties of stars by considering the moment equations derived from Eulerian hydrodynamic equations. In this paper we study the scattering of spinors of the Dirac equation, and in particular investigate the scattering issue in the limit case of rotating Maclaurin spheroids. Firstly we give the metric of a rotating ellipsoid star, then write the Dirac equation under this metric, and finally derive the scattering solution to the Dirac equation and establish a relation between differential scattering cross-section, $\sigma$, and stellar matter density, $\mu$. It is found that the sensitivity of $\sigma$ to the change in $\mu$ is proportional to the density $\mu$. Because of weak gravitational field and constant mass density, our results are reasonable. The results can be applied to white dwarfs, main sequence stars, red giants, supergiant stars and so on, as long as their gravitational fields are so weak that they can be treated in the Newtonan approximations, and the fluid is assumed to be incompressible. Notice that we take the star's matter density to be its average density and the star is not taken to be compact. Obviously our results cannot be used to study neutron stars and black holes. In particular, our results are suitable for white dwarfs, which have average densities of about $10^{5}-10^{6}$\,g~cm$^{-3}$, corresponding to a range of mass of about $0.21-0.61 M_{\bigodot}$ and a range of radius of about $6000-10000$\,km.


Introduction
As we know, a physical model usually contains a set of most basic physical quantities, some of which can be used to describe the structure and evolution of matter, while others can be used to describe the motion state of an object. A mathematical model of the motion of matter is usually composed of fundamental equations, fundamental physical quantities and definite solution conditions such as initial conditions, boundary conditions and joint conditions. The research method of physics is to abstract the objective world into a physical model, and then convert it into a mathematical model. By solving the mathematical model, we can obtain the functional relationship between physical quantities and space-time coordinates, then shall examine the physical model by observation or experiment, and improve the physical model based on comparisons.
Assuming that three semi-axes of an elliposoid are a, b and c, respectively, we define a spheroid as an elliposoid with a = b > c. As the most general physical model, a rotating spheroid is composed of fluids (magnetohydrodynamic, turbulence and perfect fluid) with the pressure, density and temperature in a gravitational field, particle groups in the non-thermal equilibrium (distribution function), and boson and fermion fields. Especially for scalar fields, the fundamental equations are the Klein-Gordon equations in flat space-time and curved space-time, respectively [1], [2], [3], [4]. A theory of gravity will be tested by the gravitational observation effect in different celestial bodies. For the same theory of gravity, different coupling coefficients give different solutions to the basic equations. There may exist significant differences between different theories of gravity, e.g., the theory of gravity with torsion, superstring theory and supergravity [5] [6] [7], [8].
In the nearly 80 years since the Dirac spinors were defined, their roles in describing fermion fields have been well established and widely accepted. For convenience, let us review the previous research work conducted by other authors related to spinors: Villalba and de Fisica [9] solved the massless Dirac equation in a non-stationary rotating causal Godel-type cosmological universe using the method of separation of variables; Dolan et al. [10] studied the scattering of massive spin-half waves by a Schwarzschild black hole using analytical and numerical methods; Deleo & Rotelli [11] developed the potential scattering of a spinor within the context of perturbation field theory; In 2011, Stefano et al. [12] studied Dirac spinors in Bianchi type-I cosmological models within the framework of torsional f (R)-gravity, while Bordé et al. [13] presented a secondquantized field theory of massive particles with a spin of a half or antiparticle in the presence of a weak gravitational field treated as spin two external field in a flat Minkowski background. In 2012, Poplawski [14] made use of Einstein-Cartan-Sciama-Kibble (ECSK) theory of gravity to discuss non-singular, big-bounce cosmology from spinor-torsion coupling, and Daude and Kamran [15]considered massive Dirac fields evolving in the exterior region of a 5-D Myers-Perry black hole and studied their propagation properties; In 2014, Brihaye et al. [16] studied Dirac equation for spherically symmetric space-time and application to a boson star in EGB gravity, and Ambrus & Winstanley [17] discussed Dirac fermions on an anti-desitter background; In 2015, Bini et al. [18] discussed massless Dirac particles in the vacuum C-metric; In 2017, Röken [19] showed the separability of the massive Dirac equation in the nonextreme Kerr geometry in the horizon-penetrating advanced Eddington-Finkelstein coordinates; Dzhunushaliev & Folomeev [20] investigated Dirac star in the presence of Maxwell and Proca fields, and Oliveira [21] investigated the influence of non-inertial and spin effects on the 2−D Dirac oscillator interacting with a uniform magnetic field and with Aharonov-Bohm effects in the cosmic string space-time.
As we know, there are two kinds of methods to study Dirac equation describing fermions in the vicinity of different types black holes. One is the Newman-Penrose (NP) formalism, and another is direct methods. In 2016, Batic et al. [22] [25] used the direct methods to obtain the Dirac equation in the background of the Garfinkle-Horowitz-Strominger black hole. In 2020, Ahmad et al. [26] took advantage of the NP formalism to study the Dirac equation around a regular Bardeen black hole surrounded by quintessence. The above-mentioned articles are related to the Dirac equation in curved space-time.
The remainder of this paper is organized as follows. In Sec. 2, we present the calculation method. In Section 3, we define connections and give the derivations of related quantities. In Section 4, we derive the Dirac equation in the weak gravitational field approximation. In Sec. 5, we give the scattering solutions to the Dirac equation and discuss the scattering issue in a limit case of rotating Maclaurin spheroids, and summarize the whole paper and look forward to the future work in Sec. 6. In an appendix we give concrete expressions of the Dirac equation in the weak gravitational potential and write the matrix M and the scattering solutions of the Dirac equation.

Calculation method
Let us emphasize that there are two approaches to study gravitational properties: (1) moment method, and (2)Dirac field scattering, particle geodesic motion and scalar field scattering. We will use the second approach, because the advantage of this approach is that the gravitational properties are closely related to the observations, and a physical model will be constrained by observations. In this paper, we will discuss the physical effect of a rotating spheroid: the scattering of Dirac spinors, which has never been done before. A specific method is first to give the metric of a rotating spheroid, then write out the Dirac equation under the metric,and find the scattering solution, and finally give a relationship between the scattering cross-section and the stellar density. According to the observed scattering amplitude, the density is expected to be determined to study the gravitational characteristics of spheroids, and the sensitivity of the scattering cross-section to the change of the density will be discussed.
In this paper, we will discuss the physical effect of a rotation spheroid: the scattering of Dirac spinors, which has never been done before. The specific method is to first give the metric of the spheroid, then write the Dirac equation under the metric, find the scattering solution, and finally give the relation between the scattering cross-section and the star density. According to the observed scattering amplitude, the gravitational properties of a rotating spheroid will be studied by determining the density, and the sensitivity of the scattering cross-section to the density change will be discussed.
Here, employing the quantum field theory of curved space-time, we will study the scattering of spinors of the Dirac equation by taking the following three steps: (1)In the first step, we get the metric (2)In the second step, we will study the Dirac equation in a curved space-time (weak gravitational field) where a = 0, 1, 2, 3 [27][28] [29]. The quantity D µ Ψ (a) is given as , which is dependent on the Pauli matrix σ k . Also we get the second term on the right side of Eq. (4), namely where the summation convention is used. At the same time, the term of γ (c) e µ (c) will be calculated, then we obtain the term of γ (c) e µ (c) D µ Ψ (b) . Therefore, we can investigate the Dirac equation in curved space-time.
Rearranging the above Dirac equation leads to For the physical meanings of quantities in Eq (6), to see the following sections. (3)In the last step, we will obtain the scattering solution of the Dirac equation in a curved space-time i.e., where the Green function , with the first-order Hamiltonian Using the first-order weak gravitational field approximation, the Ψ (x) on the right side of Eq. (4) is replaced by the solution of free particle equation, and the scattering solution is obtained.

Connections and related quantities
The Newtonian limit corresponds to the following line element [30] From Eq. (9), it is obtained that According to the metric above, we get nonvanishing connections of Γ λ µν (µ, ν, λ = 0, 1, 2, 3, and there is no summation over repeated indices) Based on Eq. (2), the quantities ω (a) µ(b) can be calculated as follows. In the following cases of (1-4), we set a, b, c = 1, 2, 3, and there is no summation over the indices (1)Case of ω

The Dirac equation in a weak gravitational field approximation
As shown in Section 2, the Dirac equation in curved space-time in a weak gravitational field approximation can be described by Eq. (3) and Eq. (4), with the wave function 2 ) T , where a = 0, 1, 2, 3. There are also mounting concerns that mΨ (a) = m φ (a) χ (a) and For convenience, by making the following substitution it is easy to obtain γ (c) e 0 (c) = Notice that we reconsider Eq. (17), and the second term on its right side is denoted as Π are written in a more compact way as follows (Here, we don't adopt Einstein convention in order to avoid confusion.) In order to study the Dirac equation in curved spacetime, we also calculate the quantity of 1 In quantum mechanism the three components of the Pauli matrix under the 2 × 2 spin representation are where b = 0, 1, 2, 3 and the Einstein convention is used. From the analysis above, the concrete expressions of the Dirac equation in a weak gravitational potential are derived, as shown in Appendix A. For convenience, Eq.(A1)-(A5) are expressed in a uniform and compact form: where σ k is the Pauli matrix and P is a [8×8] matrix, which is related to gravitational potential, the matrix elements are given as following: 1. For the first row, 2. For the second row, 3. For the third row, ∂Φ ∂x i σ j ), P 33 = im, P 34 = P 35 = P 37 = 0; (24) 4. For the fourth row, 5. For the fifth row, , P 55 = im, P 56 = P 57 = P 53 = 0; (26) 6. For the sixth row, 7. For the seventh row, 8. For the eighth row, where ε kij is the Levi-Civita symbol. By defining the matrix of P * ≡ P − P 0 that is a small quantity, we get ∂ 0 Π + HΠ = 0, and In order to get the scattering solution to the above Dirac equation using a perturbation theory, for convenience, we first set Φ = e iwt , (H 0 − E) = −P * and make the substitutions of → |Ψ >, P * → V and iw → E, then we obtain (H 0 − E)|Ψ >= −V |Ψ > and H 0 |ϕ >= E (0) |ϕ >.

The scattering solution to the Dirac equation and rotating Maclaurin spheroids
The main research objects of quantum mechanics are divided into two types: bound states and scattering states. Theoretically, the scattering state is a non-bound state, which involves the continuous region of the energy spectrum of a system. One can freely control the energy of incident particles, which is different from dealing with particles in the bound state. The bound state theory mainly involves the eigenvalues and eigenstates of discrete, quantized energies of the system [31], [32], [33].. The scattering theory mainly deals with the redistribution of scattering particles and their properties (such as polarization, correlation, etc.) in the scattering process. By analyzing the scattering results, one can find the structures inside particles, which promotes the development of basic theories. From trapped atoms to liberated quarks, a better understanding of the structure of matter depends largely on the study of scattering [34].

The scattering solution to the Dirac equation in a weak gravitational field of a rotating spheroid
In this subsection, at first, we focus on a solution to the equation of Dirac spinors of free particles, |ϕ >= |ϕ 0 > e −ik·r , which satisfies the following equation From the above expression, it is obvious that the secular equation is written as There are four different solutions to Eq. (32): E (0) = im ± ik and E (0) = −im ± ik. As the space is limited, we only choose the solution of E (0) = im + ik. It should be noted that the choice of E 0 depends on the observation of the scattering cross section (see the end of this section). Supposing k x = k y = 0 and k = k z and solving the above eigenvalues, we obtain According to the scattering formula in quantum mechanics, we have Based on the expressions of , we first set ε → 0, and define the matrices of H, H 01 , and H 02 as follows: Then we have Accordingly, their inverse matrices are written as Thus, the Green's function is given as and the scattering formula of Eq. (35) becomes where ϕ(x, t) = e iwt ϕ 0 (x) is the eigenfunction mentioned above. Letting S ≡ 1 0 , we obtain Next, we're going to derive V (x)ϕ 0 (x) ≡ P * (x)ϕ 0 (x). It should be noticed that P * is a [8×8] matrix, while ϕ 0 (x) is a [8 × 1] matrix. If both of S and σ are assumed to be numbers, the matrix elements of [P * (x)ϕ 0 (x)] would be given as follows.
1. For the first row, 2. For the second row,

For the fifth row,
6. For the sixth row 7. For the seventh row, 8. For the eighth row, We write the scattering solution of the Dirac equation , which is a [8 × 1] matrix, the eight matrix entries are given in Appendix B.
Let us continue exploring the matrix entries M 1 ∼ M 8 . Firstly, the relations between Cartesian coordinates x, y, and z and oblate elliptic coordinates ξ, η and ϕ are given as follows where ρ 0 is the focal length of a rotating spheroid, and the azimuth ϕ ∈ (0, 2π). By making a coordinate transformation (x, y, z) → (ξ, η, ϕ), we get d 3 x ≡ dxdydz = ρ 3 0 (ξ 2 + η 2 )dξdηdϕ. When ξ > ξ 0 , the gravitational potential Φ of a rotating spheroid with constant density µ is given by where M is the mass of an ellipsoid star, which can be estimated as M = 4/3πa 2 cµ = 4/3 × πρ 0 3 µ(1 + ξ 0 2 )ξ 0 ,(a = b = ρ 0 1 + ξ 2 0 is the semi-axis length in the equatorial plane, and c = ρ 0 ξ 0 is the semi-axis length in the axis of rotation). When 0 ≤ ξ ≤ ξ 0 , the gravitational potential Φ becomes where V 0 = A(ξ 0 ), Ω = 2B(ξ 0 ) and For details about the gravitational potential Φ, see Chandrasekhar [32]. From Eq. (B6), we get the matrix entry For convenience, we introduce the notation T to denote the following integral, Thus, the matrix entry of M 1 becomes In the same way, we get the other seven matrix entries T, (59) and Since the system we are studying is a rotating ellipsoid star, its gravitational potential Φ does not involve azimuth ϕ. Seeing from Eq. (55)-(62), the matrix entries of M and the symbol T are functions of p x , p y and p z . In the following, the scattering solution to the Dirac equation will be reexpressed with M or T . Then we have where the expression of Ψ (x, t)= (φ (0) , χ (0) , φ (1) , χ (1) , φ (2) , χ (2) , φ (3) , χ (3) ) is used.

A special case: Maclaurin spheroids
In this part, we will study the scattering of Dirac spinors in a special limiting case of Maclaurin spheroids. To understand the properties of Maclaurin spheroids, at first, let us discuss the motion of scalar particles in a rotating Maclaurin spheroid belonging to a special class of boson stars. When the distance between the center of the spheroid and a scalar particle, r, is larger than the stellar radius R, the exterior potential is given by Φ ex = −M/r, corresponding to its partial derivative (−∇Φ ext ) r = M/r 2 and the acceleration of gravity, a r = −M/r 2 [30]. From the Kepler's law, we get the angular velocity of the particle ω 2 = M/r 3 (or M/r 2 = rω 2 from the Newton Second Law of Motion). Similarly, when (r ≤ R, the interior potential is given by Φ in = −2πµ(R 2 − r 2 /3), corresponding to (∇Φ in ) r = 4/3 × πµr and ω 2 = 4/3 × πµ.
According to statistical mechanics, the scattering cross section, also known as the collision section, is a physical quantity describing the scattering probability of microscopic particles. The dimension of the scattering cross-section is the same as that of the area. One core issue of scattering theory is to solve the scattering amplitude by studying the probability of the particles being scattered to the unit solid angle in the direction of (p, θ, ϕ). This probability can be expressed by the scattering differential cross-section σ(p, θ, ϕ), determined by the amplitude of the spherical scattered wave f (p, θ, ϕ), i.e. σ(p, θ, ϕ) = |f (p, θ, ϕ)| 2 .
In order to evaluate T in a rotating Maclaurin spheroid, we define another new quantity as p 0 = −p x e x −p y e y − (p z − k z )e z here e x , e y and e z are the unit vectors in the x, y and z directions, respectively, then we have Since the scattering field is anisotropic, we can give the values of the scattering amplitude, f i (p, θ, ϕ) (i = 0, 1, 2, 3, 4, 5, 6, 7). Here we consider a simple case of i = 0, i.e. ϕ (0) . The scattering amplitude f 0 (p, θ, ϕ) is determined by the following expression: We also obtain the average value of f 0 (p, θ, ϕ), The probability of scattering particle in the interval of p ∼ p + dp and θ ∼ θ + dθ is determined by Taking the long-wavelength approximation k z → 0, and paying attention to p 0 = −p x e x − p y e y − (p z − k z )e z , we get the first component off 0 (p, θ), and its corresponding scattering cross-section, In the same way, the second component of the scattering amplitude is given as corresponding to the scattering cross-section Substituting Eq. (64) into the scattering solution, we find that the scattering cross-section is proportional to µ 2 , and depends on the radius of the rotating star. In addition, we can determine the constant density µ from observations,σ ± 0 (p, θ). From the above expression, it is obvious that the higher the density µ, the larger the sensitivity of scattering amplitudeσ ± 0 (p, θ) with regard to µ. Notice that above we have chosen the longwavelength approximation E 0 = −im ± ik ≈ −im, and the star is not too compact. The properties of white dwarfs (WDs) in Einstein-∧ gravity were investigated by Liu and Lü [35]. Considering the temperature effects [36], the massive WD, with a mass of about 0.61 M and a radius of about 6000 km, has an average density of about 10 6 g cm −3 , while the lower-mass WD, with a mass of about 0.21 M and a radius of about 10,000 km, has an average density of about 10 5 g cm −3 , where M is the mass of the sun. Thus, our results are reasonable for WDs.
From the above, we find that scattering cross sections σ ± 0 (E = im ± ik ≈ im) are independent of the mass of particles, m, while other scattering cross sections σ ± 0 * (E = −im ± ik ≈ −im) depend on m. Therefore the observations of scattering cross sections can determine which form of energy (E = ±im ± ik ≈ ±im) should be adopted in our physical systems.

Summary and outlook
In this work, we have studied the scattering of spinors in the Dirac equation in detail, and discussed the scattering issue in the limit case of rotating Maclaurin spheroids, and established a relationship between the scattering cross section σ and the density of matter µ. We also found: the higher the density µ, the higher the sensitivity of the scattering cross section σ to the change of µ. The results can be applied to all stars that can be treated with the Newtonian approximation approximately. In the future we'll determine the constant density of the star from the observed values of σ + 0 (p, θ, ϕ). Here we provide the prospect of the follow-up work, and study the direction continuing being advanced. First, we'll find out the other seven scattering cross sections σ ± i (p, θ)(E 0 = im ± ik ≈ im), andσ ± i * (p, θ)(E 0 = −im ± ik ≈ −im), (i = 1, 2, 3, 4, 5, 6, 7), thus the average density of the star will be determined; Secondly, we'll consider a more complex model on rotating spheroids with electromagnetic field. The scattering cross sections with and without electromagnetic(EM) fields maybe different. Thirdly, we will investigate the motion of particles inside an accretion disk around a rotating spheroid. Due to the gravity of a rotating spheroid, the disk might be possible. Finally, it is expected that we study the physical effects of a compact rotating spheroid, especially to solve the scattering solutions to the geodesic equation,scalar-field equation and spinor-field equation, respectively.For the expansion of the Newtonian analysis, we will be concerned with a metric perturbation h νυ away from flat space-time, defined as g µν = η µν + h µν , and explore the gravitational properties of rotating spheroids, as well as the scattering solution under supergravity.