Landau Levels in a Gravitational Field: The Schwarzschild Spacetime Case

: We investigate the gravitational effect on Landau levels. We show that the familiar inﬁnite Landau degeneracy of the energy levels of a quantum particle moving inside a uniform and constant magnetic ﬁeld is removed by the interaction of the particle with a gravitational ﬁeld. Two independent approaches are used to solve the relevant Schrödinger equation within the Newtonian approximation. It is found that both approaches yield qualitatively similar results within their respective approximations. With the goal of clarifying some results found in the literature concerning the use of a third independent approach for extracting the quantization condition based on a similar differential equation, we show that such an approach cannot yield a general and yet consistent result. We point out to the more accurate, but impractical, way to use such an approach; a way which does in principle yield a consistent quantization condition. We discuss how our results could be used to contribute in a novel way to the existing methods for testing gravity at the tabletop experiments level as well as at the astrophysical observational level by deriving the corrections brought by Yukawa-like and power-law deviations from the inverse-square law. The full relativistic regime is also examined in detail.


Introduction
Bringing to light the quantum properties of gravity has recently been the subject of intense investigation. Many of the investigations conducted in the past were focused on the more modest attempt to observe the effect of gravity on a quantum particle by studying, more noticeably, the behavior of cold neutrons inside a gravitational field [1][2][3][4][5][6].
In such experiments the gravitational effect manifests itself on the quantum particle through the discrete gravitational energy spectrum the particle acquires while moving inside the gravitational field of Earth.
On the other hand, more recent proposals [7][8][9] have rather been more ambitious in the sense that they consist of experimental setups designed to make manifest the quantum character of gravity itself rather than the way a particle behaves under the influence of gravity. The quantum character investigated is the quantum superposition of the gravitational field. Notwithstanding the higher importance of these more recent proposals towards unraveling the quantum nature of gravity, we believe that any additional effect based on the quantum interaction between the classical gravitational field and a quantum particle would, not only take us closer towards understanding gravity at the quantum level as well, but might even give us new ways of testing classical gravity itself.
In more recent papers, new proposals have been put forward to bring to light this interaction of a classical gravitational field with a quantum particle. In Ref. [10], the possibility of probing the effects of gravity on a quantum particle in a superposition along different paths has been investigated. In Ref. [11] the effect of gravity on a macroscopic quantum state-the superfluid made of superconducting electrons-has been investigated through the Josephson effect. In both works, the key idea is based on quantum interference induced by gravity (see also Refs. [12][13][14]). In Ref. [15], the effect on the quantum states of cold neutrons of a gravitational field created by an oscillating massive object has been studied in detail. In Ref. [16], the possibility of stimulating gravitons with ultra cold neutrons has been suggested. In both of these latter works the investigation consisted of using the quantum behavior of neutrons in a known pure gravitational field-that of Earth. Now it would be interesting to investigate the gravitational effect on a particle, not through what is expected from its quantum behavior inside a given potential-here gravitational-but by searching for a novel effect the gravitational field could have on an otherwise usual quantum phenomena involving the particle that does not originally include gravity. Such an investigation could indeed serve two independent purposes; one which would be purely theoretical and one which might be oriented more towards the experimental side of the field. Indeed, our first purpose here is to make a new contribution to the existing proposals for making manifest the interaction of gravity with quantum particles. Our second purpose, which will be based on the outcomes of the first, is to put forward a new way for testing gravity itself using quantum particles.
We propose here to investigate the effect of the gravitational field on the well-known Landau quantized energy levels of a charged particle moving inside a uniform and constant magnetic field (see, for example, Ref. [17]). More specifically, we show that the gravitational field splits the Landau levels and removes their usual degeneracy (See Ref. [18] for a recent investigation on the effect of a linear electric field on Landau levels. In this regard, it should be remarked here that one could conduct a similar investigation using gravity rather than an electric field. In addition, one could then use a large massive plane with a circular hole in order to create the needed linear gravitational field that would lead to the simple harmonic oscillations of the charge in the vertical direction to the plane. The fundamental issue behind such an approach remains, of course, the radiation of the charge during its simple harmonic motion in the vertical direction). It is known that this simple quantum effect helped explain the very important quantum Hall effect behind which many important technological applications have sprung (see, for example, Refs. [19,20]).
It should be noted here, that, unlike the indifference of the quantum Hall effect to a difference in location within a gravitational field, as shown in Ref. [21], gravity does have an effect on the Landau levels of a charged particle as we are going to show. Our investigation aligns thus with what was pointed out in Ref. [22], based on general grounds, that Landau levels are split under the effect of an external potential which is monotonic in the radial distance along the plane perpendicular to the magnetic field. This does not, however, contradict what was found in Ref. [21], as the investigation in the latter reference was about Hall conductivity, which is indeed topological, and hence metric-independent. As shown in Ref. [23], the gravitational influence on the quantum Hall effect is due to the splitting of Landau levels caused by the potential gradient of gravity in the case of Earth's field and to the quadratic gravitational potential in the case of two weakly separated massive hemispheres.
The rest of this paper is organized as follows. In Section 2, we derive the differential equation governing the motion of a charged particle inside a uniform magnetic field from the Klein-Gordon equation of a charged particle minimally coupled to the electromagnetic field in the curved Schwarzschild spacetime around a spherical mass. We first show how one recovers the familiar Landau levels by setting the mass to zero in the differential equation, and then we extract an approximate equation for the case of a weak non-zero Newtonian gravitational field. In Section 3, we expose four different methods for using the equation to extract the quantization condition. We show that only the first two methods yield a consistent and general result. The third approach introduces an extra condition on the source of the gravitational field, whereas the last one is consistent but cannot be useful in practice. In Section 4, we rely on the first approach to show how to use the splitting of Landau levels induced by the gravitational field to test any departure from the square-law for gravity. We examine the widely investigated Yukawa-like and power-law potentials as concrete examples. In Section 5, we examine the full relativistic regime of a charged particle moving within the curved Schwarzschild background under the influence of the uniform magnetic field. We conclude this paper with a brief summary and conclusion section. Appendix A is devoted to a detailed derivation of the various mathematical formulas needed in the text.

A Particle inside a Magnetic Field in the Schwarzschild Spacetime
Landau quantization of the energy levels of a charged particle in Minkowski spacetime arises in the presence of a uniform and constant magnetic field B which is perpendicular to the plane in which the charged particle is moving. In the usual textbook treatment of the problem (see Ref. [17]) one writes down the Schrödinger equation for a particle, minimally interacting with the field, by substituting the partial spatial derivatives −ih∂ i by the operator −ih∂ i − qA i , where q is the charge of the particle and A i is the vector potential causing the magnetic field B.
In order to take into account the interaction of the particle with the gravitational field, however, we are going to use, instead of the Schrödinger equation, the Klein-Gordon equation, + m 2 c 2 ϕ = 0. We ignore in this paper the spin of the test particle. The equation is written in curved spacetime, where m is the mass of the particle, c is the speed of light, and = g µν ∂ µ ∂ ν is the d'Alembertian operator. Therefore, our particle is going to be governed by the following second-order differential equation (see Ref. [24]), Here, g µν is the inverse of the spacetime metric and g its determinant. We have used here the minimal prescription D µ = −ih∂ µ − eA µ to couple the particle to the four dimensional vector potential A µ of the electromagnetic field. In addition, in view of the possibility of using heavy ions as test particles to minimize the contribution of the intrinsic spin, we choose here to take the electron charge e as the charge of our particle.
We are interested in this paper in finding how the gravitational field of a spherical mass affects the Landau levels of a charged particle moving along the equatorial plane of the mass. For that purpose, one needs first to consider the curved spacetime created by the spherical mass M. Such a spacetime is represented by the Schwarzschild metric. As is well know in general relativity, however, one should a priori take into account the effect on the curvature of spacetime due to the uniform and constant magnetic field B as well. The combined effect of the mass M and such a B-field-chosen to lie along the z direction-gives rise to the so-called Schwarzschild-Melvin spacetime [25]. It is a special case given by Ernst in Ref. [26] of the full solution of the Einstein-Maxwell equations obtained in Ref. [27]. In the coordinates (t, r, θ, φ), such a metric reads [25], We have used here, after restoring the standard units, the convenient notation of Ref. [28] in which, ∆ = r 2 − 2GMr/c 2 and Λ = 1 + 2πG 0 c −2 B 2 r 2 sin 2 θ, with 0 being the vacuum permittivity constant. On the other hand, the vector potential corresponding to such a constant magnetic field within the metric (2) takes the form [28] (Note that we resorted back here to the covariant form of the potential vector as the tetrad form Aφ = 1 2Λ Br sin θ we displayed in our previous version of the manuscript leads to much confusion. In fact, the advantage of the tetrad expression displayed in the previous version of the manuscript is that it has the right dimensions for a potential vector and it allows one to straightforwardly extract the magnetic field from the spatial components using the usual curl formula, B = ∇ × A), Now, the classical motion of both neutral and charged test particles are affected by such a spacetime and have been studied analytically in detail in Ref. [28]. It was found that the magnetic field influences neutral particles geometrically and creates a potential barrier that prevents both neutral and charged particles from escaping away from the axis of symmetry to infinity. In addition, a numerical study of the trajectories of such classical particles has also been carried out in Ref. [29] with the conclusion that the particles' trajectories in such a spacetime are integrable in very special cases only.
In fact, unlike the Schwarzschild spacetime, the metric (2) does not reduce to the flat Minkowski metric even in the absence of the central mass M. In such a case, the metric (2) represents the so-called Melvin's magnetic universe [30,31]. The investigation of the behavior of a quantum particle in Melvin's universe has been carried out in Ref. [32]. The energy spectrum of the test particle was found to display indeed a Landau-like form which is corrected by a term arising from the geometric influence of the B-field on the particle. In addition, as argued in Ref. [32], such a correction term becomes significant only for a magnetic field intensity as high as 10 19 G and for a quantum number n of the order 10 30 .
However, our focus in this paper is, rather, the effect of the gravitational field created by a spherical mass on the usual Landau levels, not the geometric, i.e., the "gravitational", effect of the magnetic field itself. For that purpose, our setup here will be that of a spherical mass immersed inside a magnetic field the intensity of which brings a negligible geometric contribution to the background spacetime. Such an assumption is indeed amply sufficient for tabletop experiments on Earth as well as at the level of astrophysical objects for which the approximation 2πG 0 c −2 B 2 r 2 1 holds. For young neutron stars and magnetars the intensity of the magnetic field can indeed be low and range around 10 8∼12 G (see Ref. [33]). Under those conditions we may thus safely adopt the metric (2) with the approximation Λ ∼ 1, in which case the metric reduces to the pure Schwazschild spacetime. Therefore, in this paper we take the metric in Equation (1) to be that of the pure Schwarz-schild spacetime around a spherical mass M [34]: Here, dΩ 2 = dθ 2 + sin 2 θdφ 2 is the 2-metric on the unit sphere. Furthermore, we are going to use the vector potential A µ of expression (3) which reduces, for Λ ∼ 1, to its usual form in the axially symmetric (cylindrical) gauge (the equivalent of the Landau gauge) in spherical coordinates (t, r, θ, φ), in which it reads, A µ = (0, 0, 0, 1 2 Br 2 sin 2 θ). Thus, we assume, for definiteness, that the magnetic field B is indeed directed along the z-axis. Now, because of the symmetry of the planar motion of the particle around the z-axis, we expect the wavefunction of the particle to be of the form, ϕ(t, r, θ, φ) = e −i Et h e i φ R(r, sin θ), where E is the energy of the particle and is a positive integer-we assume, for definiteness, a specific projection of the orbital angular momentum of the particle along the z-axis, i.e., a positive projection. Therefore, the Klein-Gordon equation for the particle in this curved spacetime reads, explicitly, Furthermore, by having the particle move along the equatorial plane, along which θ = π 2 , the cylindrical symmetry of the system allows us to also expect the function R(r, θ) to depend only on the distance ρ = r sin θ of the particle from the z-axis which is perpendicular to the plane of motion (Note that we assume here the motion of the particle to be restricted to the equatorial plane for the sake of simplicity. However, any plane perpendicular to the magnetic field would be valid provided the latitude angle θ = π 2 of such a plane is taken into account. Around the poles of a magnetized astrophysical object, for example, one should consider θ 1. Indeed, we restrict here the particle's motion to a plane in view of the many applications such a configuration gives rise to. Such applications range from the two-dimensional motion of the electron gas in the quantum Hall effect [23] to the study of the equation of state of the dense matter confined to the thin envelop or crust of magnetized astrophysical objects, as will be discussed in Section 4.2. In fact, for a particle which is instead free to explore the three-dimensional space the crucial Landau levels do not arise, for then the particle's wavefunction would be described by spherical harmonics governed by a central gravitational potential inside a strong magnetic field, much like in the case of atoms inside large magnetic fields (see the discussion below Equation (16))). This would then make R a function of the form R(r, sin θ) = R(r sin θ) = R(ρ). Therefore, we have, ∂ θ R = r cos θ∂ ρ R. Thus, we also have, ∂ 2 θ R = −r sin θ∂ ρ R + r 2 cos 2 θ∂ 2 ρ R. Substituting these inside the previous equation, the latter takes the following simplified explicit form for θ = π 2 , This is the general equation that describes the planar motion of a charged particle inside a magnetic field within the spherically symmetric Schwarzschild spacetime. Before we examine how to extract the quantization condition from this equation for weak gravitational fields, we shall first set M = 0 and solve the equation for R(ρ) in order to make contact with what we already know about the dynamics of a charged particle in a uniform magnetic field inside the Minkowski spacetime. In fact, in addition to allowing us to check the correctness of Equation (6), this first step will provide us with the fundamental wavefunction to be used later in the more interesting case of M = 0. For M = 0, we know that we should recover the free particle of relativistic energy E = E + mc 2 , moving inside a uniform magnetic field of magnitude B and perpendicular to the plane θ = π 2 . Let us assume, for simplicity, a non-relativistic regime for the particle, i.e., E mc 2 . Then, Equation (6) becomes, after setting M = 0, as follows, In order to easily solve this equation, let us choose the following ansatz for the radial wavefunction (In the literature, the ansatz R(ρ) = ρ | | e − eB 4h ρ 2 v(ρ) is sometime chosen for the radial function R(ρ). This just makes the two possibilities of a left(right)-moving particle around the mass M), R(ρ) = ρ v(ρ) × exp(− eB 4h ρ 2 ), and denote by a prime a derivative with respect to ρ. The equation then becomes, Here, we have set, for convenience, α = 2mĒ h 2 and β = eB h . Next, we perform the following change of variable: z = 1 2 βρ 2 . This allows us, in turn, to rewrite the equation in the following canonical form [35], This equation is of the well-known form zv + (b − z)v − av = 0, called a confluent hypergeometric differential equation, the solution of which is a linear combination of two confluent hypergeometric functions 1 F 1 (a; b; z). These functions are also known as Kummer's functions [35]. The general solution to the canonical Equation (9) The constants A and A are the two constants of integration and, in our case, a = (1 − α/β)/2 and b = + 1. Combining this general solution with our ansatz for R(ρ), and using our definition of z, we finally get the particle's radial wavefunction as follows, We have discarded here the second solution that goes with the constant A as it would make R(ρ), and hence the wavefunction ϕ(r), diverge at the origin r = 0 for any azimuthal quantum number ≥ 0. Now, keeping only this first solution (10), the latter also diverges exponentially for ρ → ∞ unless we impose the following condition [35,36], for some positive integer n, in which case the confluent hypergeometric function in expression (10) becomes a finite-degree polynomial. The condition (11), in turn, implies that, By substituting the values of α and β we defined below Equation (8), we find the following more familiar quantization condition for the energy E of a charged particle inside a uniform magnetic field: These are the usual Landau quantized energy levels in which we clearly see the high degeneracy of the levels due to the freedom the particle has with the orbital quantum number . Let us now use these results to examine the case of the spherically symmetric curved spacetime in the Newtonian approximation. We are interested in finding the effect of the gravitational field on the Landau levels (13). Therefore, the case of a spherical mass M for which GM c 2 ρ will amply be sufficient for us here. In addition of simplifying greatly our calculations, this restriction is also greatly motivated by its practical side with regard to an eventual experimental setup. Indeed, since even for neutron stars, magnetars and magnetic white dwarfs, for which the mass could be of the order of a few solar masses (or a fraction thereof for white dwarfs) and for which the radius ranges from a few kilometers for neutron stars/magnetars to a few thousands of kilometers for magnetic white dwarfs, the approximation GM c 2 ρ is very realistic and, hence, serves well our main purpose in the present paper.
Therefore, we can now expand what is inside the square brackets in Equation (6) in powers of 2GM/(c 2 ρ) and keep only the leading order in such a ratio. Equation (6) then takes the following form: This equation is very general and applicable even at the astrophysical level, provided the restrictions on the mass and the radius of the astrophysical objects of interest, as mentioned above, are satisfied. On the other hand, for not too strong magnetic fields, like those used in tabletop experiments, we always haveheB/m mc 2 . Furthermore, for magnetic fields in the range 10 3∼9 G, like those on the surface of magnetic white dwarfs [37], the in-equalityheB/m mc 2 holds even for electrons, whereas for neutron stars/magnetars, for which the magnetic fields could reach the range 10 8∼12 G, such an inequality holds for protons and heavier ions. Therefore, by using E = E + mc 2 , the non-relativistic regime for whichheB/m mc 2 holds makes the approximation E mc 2 amply sufficient for us here. This then implies that E 2 ≈ 2mc 2 E + m 2 c 4 and Equation (14) becomes, after keeping only the leading first-order correction coming from the GM/(c 2 ρ) factor, as follows: It turns out, as discussed in the Introduction, that there are at least four ways of using this equation to extract the quantization condition of the energy of the particle. In what follows, we are going to present all four methods and show that only the first two are able to provide us with a consistent, practical, and general quantization condition of energy.

Four Methods Leading to Quantization
Differential equations of the form similar to Equation (15) arise in many areas of both physics and quantum chemistry, ranging from the study of the harmonium to the confinement potentials [38][39][40][41][42][43][44][45][46]. Various exact and approximate methods are known for solving such an equation with central potentials [47][48][49][50][51][52][53][54][55]. However, among the more familiar ones, only two yield the right quantization condition on the energy of the particle that is both general, consistent, and practical. We are going to examine first the ones that do yield a consistent quantization condition, but which, unfortunately, are the least used ones in the recent physics literature dealing with the quantization condition of a particle obeying a similar equation. The first one consists in using the time-independent perturbation theory, whereas the second one consists in approximating the system by a simple harmonic oscillator.

Using Perturbation Theory
To make contact with the time-independent perturbations, we should first rewrite Equation (15) in a Schrödinger-like form. To achieve that, let us set R(ρ) = ρ −1/2 ψ(ρ). Equation (15) then takes the form, This is just a Schrödinger equation with an effective potential V eff (ρ) made of two parts. The first part, consisting of the first three terms inside the parentheses, represents the potential of a particle inside a uniform magnetic field. This part of the effective potential is what would give rise to the familiar Landau levels (13) when solving Equation (16) without the last term inside the parentheses. The eigenfunctions ψ(ρ) of the corresponding Hamiltonian would then be found using expression (10). The last term inside the parentheses in Equation (16) represents the second part of the effective potential, and constitutes just a small perturbing potential V(ρ). In fact, recall that in our approximation of a weak gravitational field, we assume that ρ −1 GMm ehB/m, so that the Newtonian potential −GMm/ρ is indeed nothing but a small perturbation compared to the terms coming from the magnetic interaction. We are thus set for making use of the time-independent perturbation theory formalism (see Ref. [56]). Now, since Equation (16) describes a particle inside a constant magnetic field and a central Newtonian potential, the system we are dealing with is very reminiscent of the well-studied systems of atoms inside high magnetic fields (see Refs. [57,58] for a review). However, the works dealing with atoms inside a magnetic field are concerned with finding the energy-eigenstates and the electromagnetic transitions of an atom inside a magnetic field. As such, the problem dealt with in the literature on the subject is that of solving for the fully three-dimensional motion of the electron moving within the central potential of the nucleus. A specific approximation has to be used in that case (for strong magnetic fields it is called the "adiabatic approximation" [58]), which consists in expanding the wavefunction in terms of the Landau states weighed by longitudinal wavefunctions along the direction parallel to the magnetic field. As such, the result is that the usual Landau levels simply get augmented by the Coulomb bound state energies as displayed in Equation (8.7) of Ref. [58]. In our case, however, we shall use instead the time-independent perturbation theory to deal with the effect of the gravitational central potential on the planar motion of a particle moving around a massive object inside a magnetic field. For this reason our results will differ drastically from what is found in the case of an atom inside a strong magnetic field based on the adiabatic approximation. Our result will indeed display a product of the magnetic and gravitational contributions that has not been previously reported in the literature.
If we denote the unperturbed n th Landau energy level (13) by E (0) n then, as long as the spacing E m between these energy levels is greater than the Newtonian potential GMm/ρ, we are guaranteed that the perturbation expansion will be legitimate. The problem that might arise with this method is that the computation might become inaccessible as the Landau levels are infinitely degenerate. Fortunately, however, the fact that the Newtonian perturbation depends only on ρ, i.e., is rotational symmetric, means that there is no coupling between two different Landau orbitals within the same Landau level.
Since we already saw that only expression (10) converges for small ρ, we are going again to keep here only that expression and use it as the eigenfunction of the unperturbed Hamiltonian corresponding to Equation (16). Recalling then that ψ(ρ) = ρ 1/2 R(ρ), we have the following explicit expression for the unperturbed wavefunctions of the Schrödinger Equation (16), with n = 1 2 (α/β − 1). Therefore, our first task is to find the normalization constants A n which can be determined by imposing the completeness condition on the eigenfunctions, However, since for a more realistic setting the sphere of mass M has a finite nonzero radius ρ 0 , much larger than its Schwarzschild radius, the test particle's position would be limited to the interval ρ ∈ [ρ 0 , ∞). In other words, because of the finite size of the mass source we are considering here, the spherical mass imposes a physical restriction on the test particle so that the interior region of the former (located below ρ = ρ 0 ) remains physically inaccessible to the latter. Iimportantly, also, is the fact that our gravitational field is valid only for ρ ≥ ρ 0 , i.e., outside the spherical mass.
In the realistic case of a finite-radius spherical mass we should therefore distinguish two different regions when solving the Schrödinger equation. Region I, say, would repre-sent the outside of the spherical mass, for which ρ > ρ 0 , while region I I would represent the inside of the spherical mass, for which ρ < ρ 0 . As an ideal system, however, we assume the spherical mass to be completely reflective for the particle. This means the particle's wavefunction vanishes inside the sphere as the particle has no chance of penetrating inside the latter. Indeed, in this case our system consists effectively of a particle moving around a semi-infinite potential well, inside of which the potential is infinite and outside of which the potential is just that given in Equation (16). The wavefunction in the region I being then just expression (17), all we need to impose therefore is the continuity of the latter across the surface ρ = ρ 0 . This amounts to imposing the following requirement: Based on expression (17), this requirement, in turn, amounts to imposing the following condition involving the confluent hypergeometric function: The physical meaning of condition (19) can easily be understood as arising from the geometry of our system. In fact, solving condition (19) for the two unknown integers n and simply returns these as a function of β, i.e., the magnetic field, and the radius ρ 0 of the massive sphere. Indeed, by substituting the definition of the confluent hypergeometric function given below Equation (A11) into Equation (19), the latter yields a single algebraic equation involving both integers n and . By having the finite-radius mass sit at the center does indeed geometrically disturb the motion of the particle to which all the possible Landau levels n = 1, 2, 3, . . . and all the possible orbital numbers = 1, 2, 3, . . . would have otherwise been accessible. The existence of the forbidden region 0 < ρ ≤ ρ 0 implies that only certain values of n and are possible depending on the value of the product 1 2 βρ 2 0 . Thus, to make the particle's energies acquire the Landau levels, the magnetic field itself should be adjusted with the size of the massive sphere to allow for the condition (19) to be simultaneously satisfied. In the case of a point-like mass, i.e., for ρ 0 = 0, condition (18) is automatically satisfied and one does not need to impose (19) and, hence, no restriction is imposed on the quantum numbers n and either. Being interested here simply in the effect of the gravitational field on the Landau energy levels, however, we are going to assume in the remainder of this paper that a specific combination of the magnetic field and the size of the sphere is already such that the existence of Landau quantum levels and orbitals is guaranteed.
Going back to the normalization constants A n , the normalization condition that we should impose here to get these constants is then For this purpose, we make use of the integrals derived in Appendix A. We easily find that the constants A n are given by M −1/2 n , where M n is given by Equation (A10) after setting n = m there. Now, as noted above, although the Landau energy levels are infinitely degenerate, the fact that the gravitational interaction is rotational symmetric means that the perturbing potential V(ρ) does not couple between two different Landau orbitals of quantum numbers and . This implies that the matrix elements n, |V(ρ)|n, of the perturbation are diagonal. Therefore, the degenerate time-independent perturbation theory yields n is the unperturbed nth Landau level (13). This gives then the following more explicit first-order correction to the energy of the orbital belonging to the nth Landau level, To evaluate the improper integral, we first substitute expressions (17) for the wavefunctions and their normalization constants A n as given by Equation (A10). Then, using the result (A17) given in Appendix A we find, where P n is given by Equation (A17) after setting n = m. Although the product P n M n has a long and cumbersome expression, it is actually easy to conclude from such a product that the first-order correction to the energy levels is proportional to the square root of the magnetic field. In fact, in both infinite series (A10) and (A17) defining P n and M n , respectively, there appears the common constant factor (β/2) − . However, the M n -series has, in addition, the constant factor (β/2) − 1 2 , whereas the P n -series comes with the additional constant factor (β/2) − 1 2 . This implies that the product P n M −1 n gives rise to the constant factor (β/2) 1 2 . In light of this observation, a better formula for the general first-order correction to the Landau levels is then the following: We have introduced here the reduced seriesP n andM n which consist of expressions (A17) and (A10), respectively, without their respective factors (β/2) − − 1 2 and (β/2) − −1 . In order to see this more clearly, let us compute the explicit correction to the first Landau level by setting n = 1 in Equation (21). We find, The quantities M 1 and P 1 are given by Equations (A14) and (A18), respectively. This result shows that the splitting induced on the Landau levels by the gravitational field has actually a simple form at each level and for each orbital. We notice that the splitting is larger for stronger magnetic fields. The additional term on the right-hand side of expression (23) is still not fully transparent, however, as it involves inside the product P 1 M −1 1 incomplete gamma functions of the form Γ( + 1, β 2 ρ 2 0 ) and Γ( + 1 2 , β 2 ρ 2 0 ), which, in turn, involve infinite series made of and β 2 ρ 2 0 . For small values of , it is already obvious from the general definitions (A9) and (A16) of the terms M 1 and P 1 , respectively, that the product P 1 M −1 1 is of the order ρ −1 0 , so that the correction term in Equation (21) is of the order −GMm/ρ 0 , as expected. However, in order to extract a useful and a more transparent expression for the splitting of the n = 1 Landau level, we shall use the large-limits (A15) and (A19) of the exact expressions (A14) and (A18), respectively. Then, the splitting (23) takes on the elegant and more transparent form, Thus, for large-orbitals the correction to the first Landau level is due to the generic product of the gravitational field contribution and the square root of the magnetic field. The last step in Equation (24) comes from the asymptotic expansion for large arguments z of the gamma function [35], Therefore, the correction term in Equation (24) decreases like 1/ √ , and thus becomes gradually suppressed for large . In contrast, from Equations (A10) and (A17) giving M n and P n , respectively, we see that the first-order correction (22) does not get suppressed for large n. In addition, we notice that for large , the correction becomes insensitive to the radius ρ 0 of the massive sphere.
We would like to emphasize here the fact that, as explained in Appendix A, the largeapproximation (24) is valid for extremely large values of , for the limit was found by taking into account the already very large term β 2 ρ 2 0 inside the incomplete gamma functions in Equations (A14) and (A18). Therefore, contrary to what it might seem at first sight, the correction term on the right-hand side in Equation (24) is really small (as is required for a perturbation) in comparison to the first term even for masses M of the order of the solar mass. As derived in detail in Appendix A, this is in fact guaranteed provided that is bigger than βρ 2 0 2 . Therefore, the correction term in Equation (24) is indeed smaller than GMm/ρ 0 , which, as we already argued below Equation (16), is nothing but a small perturbation compared to the Landau energy represented by the first term. The same remark is also valid for the general splitting formula (23). In the latter, the smallness of the correction term compared to the first is less transparent but can, nevertheless, still be inferred from the less trivial expressions (A10) and (A17) of M n and P n , respectively. In fact, first of all, both series are exponentially suppressed by the term e − β 2 ρ 2 0 . More important, however, is that the series (A10) contains the factor Γ( + 1) in the numerator whereas the numerator of the series (A17) contains the factor Γ( + 1 2 ). It is now important to remark here that when setting B = 0, i.e., in the absence of the magnetic field, one does not find any quantization of energy coming from magnetism. Setting B = 0 in Equation (22), however, makes the energy vanish altogether. Actually, this is simply due to the fact that in this case both series M mn and P mn do not exist, for the integrals that gave rise to these series vanish since the Kummer's functions vanish for β = 0. A proper treatment of the case B = 0 consists indeed in solving the Schrödinger equation with the full central gravitational potential as the unique potential.
It is interesting to compute now the second-order correction to be able to fully appreciate the effect of the gravitational field. The second-order corrections to the energy levels are of the form, We have introduced here again the reduced seriesP kn andM n . The form of this correction is actually very familiar, for one should indeed recover at the second order in GM a form for the energy levels similar to that of a particle inside a central potential of the Coulombic 1/ρ-form. This is in analogy with the hydrogen atom for which the electron's energy levels are ∝ m(k e e 2 ) 2 /h 2 (see, for example, Ref. [56]). However, this expression cannot be used to find be energy levels in the case of a pure gravitational field, i.e., by setting B = 0 inside this formula either. In fact, a proper treatment in this case would be to set B = 0 instead in Equation (16), as the latter solves exactly just like for the hydrogen atom in terms of Laguerre polynomials [56].
Before we move on to the second approach, another short note is in order here. It is actually possible to start instead from the unperturbed Hamiltonian of a particle inside the gravitational potential V(ρ), together with the second term inside the parentheses in Equation (16), and consider the rest of the terms rising from the interaction of the particle with the magnetic field as being the perturbing potential. This approach would also easily work because the Laguerre polynomials-which constitute then the eigenfunctions of the unperturbed Hamiltonian-are also easy to integrate like the hypergeometric functions. The downside of this strategy for finding the effect of gravity on the Landau levels is that, experimentally, it does not make much sense to have a magnetic field so small that a mass of a few kilograms would overcome the force that such a magnetic field exerts on the charged particle. Earth's magnetic field would already affect the particle with a force that is greater than the gravitational force a few kilograms of iron would exert on the particle. Furthermore, although the reverse might be true for astrophysical processes, it is rather still the weak gravitational field relative to the magnetic field that matters most, as we shall discuss in Section 4.

Using a Harmonic Oscillator Approximation
This approach is based on finding the equilibrium distance of the particle from the spherical mass at which the potential energy of the particle is minimum [39]. In fact, the gravitational interaction of the particle with the spherical mass adds up to the interaction of the particle with the magnetic field to balance the centrifugal force due to the kinetic term and, hence, form the Landau bound states. This balance takes place at a specific radial distance ρ 0 from the center of the mass for each specific Landau level n and orbital . At this specific radial distance, the total potential of the particle can be approximated by that of a simple harmonic oscillator for which the quantized energy spectrum is well known.
Let us then start again from the Schrödinger Equation (16) and expand the effective potential V eff (ρ), contained inside the parentheses, in a Taylor series around the equilibrium position ρ * given by V eff (ρ * ) = 0. At the second order of the expansion, the effective potential then reads, where V 0 = V eff (ρ * ) and mω 2 = V eff (ρ * ). With such an approximate potential, Equation (16) takes the form of the usual Schrödinger equation of a simple harmonic oscillator for which the energy eigenvalues are given by, where n is again a positive integer. Our task now reduces then to solving for ρ * the condition V eff (ρ * ) = 0. That is, we need to solve the following quartic equation in ρ * , Notice that the case = 0 does not arise here as in such a case, Equation (29) does not admit any real solution. Now, before we solve this equation for the general case M = 0 and B = 0, it is instructive to examine first what would this approach give for the well-known cases of a particle inside a Coulombic potential, i.e., when B = 0, and then for a particle inside a magnetic field only, i.e., when M = 0. This will help us find out to what extent we can rely on this approach when tackling the general case of non-vanishing B and M.
Setting B = 0 in Equation (29) turns the latter into a first-degree equation which can easily be solved for ρ * . Substituting then the resulting expression of ρ * inside Equation (28) gives the following quantized energy levels of the particle inside the gravitational field, We recognize in this expression once again the main terms characteristic of the quantized energy levels of a particle inside a Coulombic potential [17,56]. This would give rise to the familiar proportionality ∝ m(GMm) 2 /(h 2 2 ) for large . The large-limit required for this approach to be accurate can be understood as enhancing our approximation of a very weak gravitational field, for then the Lorentz force becomes indeed much bigger than the Newtonian attraction. This actually agrees qualitatively with what we found in Equation (26) using the perturbation theory approach at the second order, provided of course one takes the large-limit there too.
Setting M = 0 in Equation (29), on the other hand, leaves the latter as a quartic equation, the solution of which is, however, very easily found. Substituting then the resulting expression of ρ * inside Equation (28) gives the following quantized energy levels of the particle inside a uniform magnetic field in Minkowski spacetime, We recognize in the first two terms of this formula the Landau quantized energy levels. However, the exact Formula (13) is recovered only in the large-limit again.
Let us now turn to the general case of a non-vanishing magnetic field in Equation (29). The four independent solutions to general quartic equations are well known, see Ref. [59]. However, instead of writing down the exact cumbersome expression of the physical solution for which ρ * is real and positive, we are going to content ourselves here by extracting only an approximation for it. In fact, since we are already in a weak gravitational field approximation, GMm 2 e 2 B 2 ρ 3 * , solving exactly Equation (29) is not necessary for us here. Thus, the approximate real and positive solution we find for ρ * is the following, Note that to get this approximate expression of ρ * , we have relied on the additional assumption that the dimensionless parameter x is very small. This assumption is amply satisfied only at laboratory-level experiments. Without such an approximation, the full expression of ρ * would become unwieldy and useless for neatly extracting the contribution of the gravitational field as we do in Equation (33) below. In fact, substituting expression (32) into Equation (28), the latter gives the sought-after quantized energy levels, We clearly see form this expression that we recover again the usual Landau levels plus the first-and second-order corrections we obtained using perturbation theory. Both firstand second-order corrections agree qualitatively with expressions (22) and (26), respectively. More importantly, however, is that for large the first-order correction does agree quantitatively as well with the expression obtained in Equation (24) using perturbation theory. In fact, in addition to displaying the generic product of the gravitational field contribution and the square root of the magnetic field contribution, the two corrections become actually identical for large . The numerical factors coincide and both imply a correction that decreases like 1/

√
. The second correction does not depend on the magnetic field. It is entirely due to the the gravitational field encoded in the Schwarzschild metric, and is quadratic in GMm as in Equation (26).

Using the Biconfluent Heun Equation: The Polynomial Approach
In contrast to the previous two methods for dealing with Equation (15), the polynomial approach does not consist in extracting the energy levels from the corresponding Schrödinger Equation (16) by approximately solving the latter. Instead, it is based on solving Equation (15) exactly and then requires that such a solution be physical by imposing a specific condition to be satisfied. The quantization condition on the particle's energy thus merely comes from imposing such a condition on the wavefunction.
In order to solve Equation (15) exactly, we begin, as we did in Section 2, by setting, R(ρ) = ρ v(ρ) exp(− eB 4h ρ 2 ), and denote by a prime a derivative with respect to ρ. The equation then takes the form, where, we have set, α = 2mĒ h 2 , β = eB h , and γ = 2m 2 GM h 2 . One way of dealing with this equation would be to expand the function v(ρ) in an infinite series in the variable ρ. Thus, one would set v(ρ) = ∑ k=+∞ k=−∞ a k ρ k and one would then plug this series inside Equation (34). The latter would be turned into the following identity involving an infinite series in ρ, Requiring that this identity holds for all ρ leads to the following three-term recursion relations: Unfortunately, there is no known general analytical solution and no simple convergence criterion for dealing with such a three-term recursion relation series. In addition, the truncation method one often uses to terminate the series when dealing with two-term series does not work here. A closely related but more accurate approach, that is indeed based on a truncation method, will be given in what follows and in the following subsection.
The right way to deal with Equation (34) is actually to solve exactly for the radial function v(ρ) as follows. First, we need to introduce the new variable z = β/2 ρ. Substituting this inside Equation (34) makes the latter take the following canonical form, known as the biconfluent Heun differential equation [60][61][62][63], with a = 2 , b = 0, c = 2α/β + 2 , and d = −2γ 2/β. The solution to this equation is given by a linear combination of two independent biconfluent Heun functions as follows, The constants C 1 and C 2 are the two constants of integration. Only the first term is convergent for z → 0, however. Although in our case the distance ρ, and hence z, does not go to zero because the mass used to create the Schwarzschild metric has a finite radius, we only keep the first term in order to guarantee the convergence of the series for arbitrary small, but non-zero, z.
The special function H(a, b, c, d; z) is called the biconfluent Heun function and it is given explicitly by the following infinite series [45,60], where (1 + a) k is the Polchhammer symbol (see Appendix A). The factors A k in the series satisfy the following three-term recursion relation, Two possibilities [45] are now available for extracting valid quantization conditions on the energy of the particle based on the biconfluent Heun function (39). The first will be exposed here and the second will be left for the next subsection.
The first possibility for finding a quantization of energy is to use the fact that the series (39) is highly divergent at infinity and, therefore, one necessarily needs to truncate the series to obtain an n-th order polynomial which would represent indeed a physical solution. For this purpose, one easily deduces that, according to the recursion relation (40), one has just to set [45], for some positive integer n. In fact, such a requirement guarantees that all the subsequent factors A n+k in the series (40) vanish for any positive integer k. This method has indeed been adopted by many authors interested in finding quickly a quantization condition involving Landau levels of various systems under a magnetic field. This approach has recently been used in Ref. [64] to study the energy levels of a galaxy moving in a Newtonian potential corrected by the cosmological constant. On the other hand, in Ref. [65] the effect on the Landau levels of a charged particle moving around a rotating cosmic string inside a magnetic field was also investigated by extracting a quantization rule for the energy of the charged particle based on this approach. Unfortunately, the downside of this approach is that it does not provide a consistent quantization rule for all values of the magnetic field-or of the mass-source of the gravitational field-in which the physical system is immersed. In fact, by substituting in the first condition of Equation (41) the expressions of a and c as given in terms of α, β, and , one just recovers a quantization condition similar to the usual Landau quantization; i.e., α = β(n + 1), from which one deduces that, This expression is indeed similar, but not exactly identical, to condition (13). In the latter the energy levels are proportional to half-integer multiples of the cyclotron frequency eB/2πm whereas according to condition (42) the quantized energy is any integer or halfinteger multiple of that frequency.
Furthermore, the second condition in Equation (41) constitutes, as already pointed out in Ref. [45], an additional constraint involving again the energy E of the particle as well as the magnetic field B. In fact, the coefficient A k+1 in the recursion relation (40) is given by the following matrix determinant [45], where, In our case, we have b = 0 so that D = 1 2 d. It is clear then from the matrix determinant (43) that requiring the term A n+1 to vanish will produce an (n + 1)-degree equation in the parameter d = −2γ 2/β = −2 eBh/(2m 2 GM). Such an equation involves, in addition, the term a = 2 as well as the term c = 2α/β + 2 = 4mE /(eBh) + 2 .
However, having already obtained the quantization condition (42) that relates the energy E to the magnetic field B, it is clear that the extra (n + 1)-degree equation A n+1 = 0 can only impose a specific constraint on the spherical mass M used to create the gravitational field. This means that this approach cannot be consistently applied for an arbitrary gravitational field but only for cases in which one manages to fine-tune the mass M in such a way to produce the quantization (42) itself. In addition, such a quantization would not therefore introduce any novelty as all it yields are quantized energies similar to, but not exactly the same as, the usual Landau levels in Minkowski spacetime. The only advantage of the method is thus to show that it might be possible to achieve a quantization of a particle's energy in a magnetic field even in the presence of a non-zero gravitational field by fine-tuning the mass-source of the gravitational field.

Using the Biconfluent Heun Equation: The Asymptotic Approach
The asymptotic approach is again based on the solution (39) to Equation (37) in terms of the biconfluent Heun function. However, in this approach one does not truncate the solution (39) by imposing the conditions (41) in order to recover a finite-degree polynomial to guarantee convergence. Instead, in this approach [45] one is rather concerned by the fact that for z → ∞, the asymptotic behavior of Heun's biconfluent function H(a, b, c, d; z) is given by [45,60,66], where N (a, b, c, d) is a constant. This asymptotic behavior renders indeed the function H(a, b, c, d; x) not square integrable, and hence unphysical as a wavefunction. Therefore, the only way to guarantee square integrability, and hence for the wavefunction to represent a physical system, is to impose N (a, b, c, d) = 0 [45]. This condition is what constitutes a real quantization condition free of any inconsistency as it does not involve two separate conditions involving the energy of the particle and the parameters of the system. The constant N (a, b, c, d) with b = 0, which is the case of interest to us here, is given by the following infinite series [45], where, λ = a 2 + c 2 + 1 and the factors A k (a, b, c, d) in the series are given by the determinant (43) in which the terms D k and A k are given this time by [45], It is obvious from these expressions that, while the requirement N (a, 0, c, d) = 0 would indeed give a genuine quantization condition, as it involves a single equation in all the parameters of the system, it is not at all useful in practice as it requires one to find the zeros of the infinite series (46) in order to be able to extract the quantization condition.

With a Yukawa-Like Deviation
In the previous two sections we were concerned only with bringing to light the effect of a known spherically symmetric gravitational field on the familiar Landau energy levels of a charged particle inside a magnetic field. The results we obtained in those sections may actually be put to use the other way around. In other words, we may now use what we have learned on how to deal with the simultaneous presence of a gravitational field and a magnetic field to unravel, or at least to test, a given unknown spherically symmetric gravitational field based on this specific splitting of the Landau energy levels of charged particles.
The field of testing gravity using elementary particles is also rich in intense recent investigations aimed at testing the inverse-square law for the gravitational attraction, either at short distances based on tabletop experiments [67][68][69], or at the observational astrophysical level [70]. The most widely investigated form of departure from the inverse-square law for gravity has the following Yukawa-like gravitational potential (see, for example, Refs. [71,72]), where δ is a dimensionless parameter that quantifies the relative strength of the additional Yukawa-like potential compared to the Newtonian potential, and λ represents the distance range of such an additional potential.
Our task now is to use again the time-independent perturbation theory based on this new gravitational potential. Having found the contribution E (1) n of the 1/ρ -perturbing term in Equation (22), all we need now is to evaluate the extra contribution of the perturbation coming from the second term inside the parentheses of expression (48). More specifically, to the first-order, the correction to the nth Landau level will read, (49) where, E (0)L n is the usual unperturbed Landau nth level, and E (1)N n is given by the second term on the right-hand side in Equation (22), and is due to the perturbation coming from the Newtonian potential. The presence of the exponential in this integral makes its evaluation less straightforward than those needed in Section 3. The detailed steps leading to the exact evaluation of the integral are exposed in Appendix A. Substituting Equation (A34) into the third term of the right-hand side of Equation (49), leads to the following first-order correction to the nth Landau level due to the modified gravity potential (48), where the first-order Yukawa perturbation E (1)Y n is given by the following expression: We have introduced here again the reduced seriesȲ n andM n which consist of expressions (A22) and (A10), respectively, but without their respective factors (β/2) − − 1 2 and (β/2) − −1 . The splitting of the first Landau level, n = 1, is found by substituting expressions (A23) and (A14) of Y 1 and M 1 , respectively, into (51).
While the result of (51) has the advantage of being exact, the fact that it involves, as displayed in Equation (A22), an infinite series in 1/λ (λ being a very small parameter) makes it of limited practical use for estimating the order of magnitude of such a correction in a given physical setting. For this reason, we shall, instead, look for an approximate estimate of the correction by replacing the exponential function e −ρ/λ in the integral (49) by the dominant value e −ρ 0 /λ , it takes within the whole interval of integration. In fact, the integrand in that integral, being already exponentially suppressed for large values of ρ by the factor e − β 2 ρ 2 as we saw in Section 3.1 for the Newtonian correction, the error we introduce in our estimate based on such a substitution will be at most of the order of λe −ρ 0 /λ times the exact correction. Therefore, by performing the integral in Equation (49) after replacing e −ρ/λ by e −ρ 0 /λ , we find the first-order correction to the nth Landau level to be given by Equation (50), where the Yukawa contribution to the correction is approximately given by, We have again used here the results of Appendix A where the integral in Equation (49) without the term e −ρ/λ is found in terms of the quantities M 1 and P 1 as given by Equations (A14) and (A18), respectively. We introduced, as done in Section 3.1, the barred quantitiesM 1 andP 1 to be able to extract the factor √ eB/2h.

With a Power-Law Deviation
A widely investigated power-law departure from the inverse-square law for gravity has the following form (see, e.g., Ref. [73]), where the parameter L has the dimensions of a length and quantifies thereby both the strength of the "fifth" force and the range of the latter. The power s is an integer which, in the extra spatial dimensions scenarios, is usually taken to be s ≤ 6. Following the same steps as those we took to arrive at the result (20) for the inverse-square law case, we find the following more explicit first-order correction to the energy of the orbital belonging to the nth Landau level due to a power-law deviation, In the second line we have used the formulas given in Appendix A to evaluate the integral. Just as we did in Section 3.1, we shall examine here, for concreteness, the correction brought to the first Landau level n = 1 and in the large-limit. Using the expressions (A15), (A19), and (A27) for M 1( 1) , P 1( 1) , and Q 1( 1) , respectively, from Appendix A, we find, From this result we clearly see the additional contribution of the power-law deviation from the inverse-square law. In contrast to the correction brought by a Yukawa-like deviation, the correction introduced by a power-law deviation is not exponentially suppressed. Furthermore, the correction is more significant, as one might have expected, when the power s is smaller. This last result can be effectively used as a means to measure the possible range of values for the parameter L for a given integer power s. Now, these two results can actually be used to probe any deviation from the squarelaw not only by using tabletop experiments, but even by using astrophysical observations. In fact, there is no better setup for a strong magnetic field combined with a relatively weak gravitational field required here-but still strong compared to Earth's standards-than within the environment of an astrophysical object, like a neutron star, a magnetar, or even a magnetic white dwarf. The quantized Landau states of the electrons in the magnetic field of a white dwarf undergo a change in their equation of state which influences the pressure and energy density of the white dwarf. The combined pressure and energy densities of matter and of the magnetic field determine the mass-radius relation of strongly magnetized white dwarfs [74][75][76][77]. On the other hand, it is known [78] that the strong magnetic field of a neutron star softens the equation of state of matter inside the star also due to Landau quantization. In addition, it was recently shown in Refs. [79,80] that as a result of Landau quantization of electrons' motion under the strong magnetic field of a magnetar, the neutron-drip transition in the crust of the star is shifted to either higher or lower densities depending on the magnetic field strength. The quantization of the motion of the electrons also makes the star's crust almost incompressible, leading to a density that remains almost unchanged over a wide range of pressures. With a splitting of the Landau levels due to the gravitational field of the star, however, we should observe a variation in any of these predicted patterns as a function of the star's mass. The extent of such a variation should betray any eventual deviation from the inverse-square law thanks to the factor e −ρ 0 /λ in Equation (52) or the factor (eBL 2 /2h ) s 2 in Equation (55). Our main approximation, ρ −1 GMm ehB/m, and the non-relativistic regime E mc 2 , on which our investigations in this paper are based, easily find common ground within the environment of young neutron stars and/or magnetars and magnetized white dwarfs for which the intensity of the magnetic field is high enough compared to the gravitational interaction of the particle and, yet, is low enough for the spacetime effect of the magnetic field itself to be negligible. We thus expect our results to very well find a good application in those extreme environments of these stellar objects where both gravity and the magnetic field strengths are beyond what one could hope to achieve at the laboratory level.
In fact, the gravitational effect simply adds up to the already existing thermal broadening of the Landau levels [81]. In this case, however, the coherence lengths in such environments are, of course, drastically reduced. Nevertheless, as discussed in the footnote below Equation (5), taking the horizontal plane to be around the poles of the star, where the longitudinal angle is θ 1, still allows us to rely on our approach to obtain an estimate of the gravitational effect by avoiding the need for coherence extending all along the equator. Indeed, around the poles of the star the cylindrical coordinate ρ is such that ρ = r sin θ r. The results obtained in Section 3.1 can then easily be adapted to such a case (as done in Ref. [23] for the inside of a sphere) by expanding here, instead, the non-relativistic Newtonian potential around the poles of the star into that of a simple harmonic oscillator. Denoting now the radius of the star by R, we have at the poles, r = R 2 + ρ 2 , with ρ R. Therefore, we have, up to a constant, −GM/r ∼ 1 2 GMρ 2 /R 3 . With such a potential, Equation (16) can easily be solved to yield the following counterpart of Equation (24) for the splitting of the Landau levels (see Equation (12) of Ref. [23], where the role of the constant K there is here played by GM/R 3 ): In contrast to the gravitational splitting given by Equation (24) for the case of motion along the equator, the splitting of the Landau levels for the case of motion around the poles is proportional to the orbital quantum number for 1, and it is inversely proportional to the magnitude B of the magnetic field. The reason being that the gravitational harmonic oscillator is stronger, the larger the circles made by the particle around the poles are. On the other hand, it is clear that the stronger the magnetic field, the more overwhelmed the gravitational harmonic oscillator becomes. Remarkably, this simple non-relativistic approximation is valid at the poles regardless of the smallness of the coherence length of the particles around the poles. Note, in addition, that, when it comes to astrophysical objects, the uniform-magnetic field approximation we adopted throughout this paper is certainly much more adequate for the case of motion around the magnetic poles of a star than elsewhere on the surface of the latter. Now, for the sake of generality and for wider applications, a relativistic extension of the analysis we have done so far is needed and is going to be the subject of the next section.

The Full Relativistic Regime
By working within the full relativistic regime, the full relativistic energy E should be kept and, as a consequence, terms involving products of E with terms containing the radial distance ρ in the general Equation (6) cannot be eliminated. Therefore, the equation cannot be brought to a Schrödinger-like form using the variable ρ. For this reason, we are going to proceed in three steps. First, we set R(ρ) = f (ρ)/ρ in Equation (6) so that the latter becomes, Next, in order to isolate the energy term E 2 and, at the same time, keep the two derivative terms without any factors containing GM/c 2 ρ, we use the Eddington-Finkelstein radial coordinate ρ * defined by [34]: With such a change of the radial coordinate, Equation (57) becomes, Finally, in order to cast this equation into a Schrödinger-like form, we set f [ρ(ρ * )] = ρ − 2GM/c 2 χ(ρ * ). The equation takes then the following exact Schrödingerlike form (To the best of our knowledge, this equation has never been reported in the literature before):

−¯h
This equation is indeed of the Schrödinger type and, hence, the quantization of the eigenvalue E 2 /2mc 2 and the corresponding perturbative corrections can easily be found using our results from Section 3. In order to be able to use perturbation theory, we are still going to work within the approximation 2GM c 2 ρ(ρ * ). Let us then first find the unperturbed relativistic energy E by setting M = 0 in Equation (60). In this case, we have ρ * = ρ and the equation takes exactly the same form as Equation (6) with the only difference that the term E on the right-hand side of Equation (6) is here replaced by (E 2 − m 2 c 4 )/2mc 2 . Therefore, the unperturbed wavefunction χ (0) n, (ρ) is given by expression (17) and the unperturbed relativistic energies are quantized and given by, This gives back the non-relativistic result (13) within the approximation E mc 2 . Now that we have the unperturbed relativistic energy, we are going to consider again Equation (60) in which we take all terms containing the ratio GM/c 2 ρ(ρ * ) as mere perturbing terms. Recalling, however, that our unperturbed wavefunction χ (0) n, (ρ) depends on the variable ρ * through its dependence on ρ, we should re-express Equation (60) in terms of ρ * as well. For that purpose, we need to express ρ in terms of ρ * . This could actually be done using Equation (58). However, as the variable ρ can only be extracted from Equation (58) as an infinite series in ρ * , and given that ρ appears in the denominator of many terms of Equation (60), we would have to truncate the infinite series by keeping only the zeroth order in GM/c 2 . According to Equation (58) this amounts to simply approximating ρ * by ρ. This approximation is in fact, according to Equation (58) The additional terms inside the square brackets will all supply corrections to the quantized relativistic energy (61). Using these terms, the first order correction to the energy reads, In the second line we have used again the formulas given in Appendix A to evaluate the integral. Just as we did in Section 3.1 for the non-relativistic case, we shall examine here, for concreteness, the correction brought to the first Landau level n = 1 and in the large-limit. In fact, in contrast to the general result (63), setting n = 1 and letting 1 has the advantage of showing explicitly the different -dependence of both the Newtonian and the general relativistic contributions. Using the expressions (A15), (A19), and (A31) for M 1( 1) , P 1( 1) , and R 1( 1) , respectively, from Appendix A, we finally find, This result gives the orbital splitting of the relativistic energy corresponding to the first Landau level for large orbitals . Comparing this expression with the non-relativistic result (24), we easily distinguish between the different contributions. In fact, we see from the second step leading to Equation (64) a distinct signature of both the special relativistic effects in the term mc 2 due to the fast motion of the particle as well as the general relativistic effects in the term GM due to the curvature of the spacetime background. The latter gives rise to the extra product GMeB which is absent in the Newtonian approximation of the Schwarzschild metric. This simultaneous contribution of the magnetic field and the curvature of spacetime in a single term arises thanks to the angular motion of the particle within a curved background. Furthermore, we see from the last line that the larger is the mass of the charged particle, the weaker is such a general relativistic correction. This is due to the fact that the mass-energy of the particle becomes then too large for the coupled magnetic and curved spacetime contribution to make much difference. The other important fact is the different -dependence between the Newtonian correction and the general relativistic correction. The former scales like 1/ √ whereas the latter scales like √ . This difference is due to the fact that the general relativistic setting becomes more important as the particle moves faster around the central mass. This is actually very reminiscent of the astrophysical observation that the general relativistic correction to the motion of planet Mercury is more important than those of the other planets that move slower as they are farther away from the Sun.
This purely relativistic result (64) can actually find a more immediate application for astrophysical situations like in the motion of electrons at the surface of magnetic white dwarfs and even neutron stars and/or magnetars. In fact, for a magnetic white dwarf of radius of about 0.015 R and a mass of about 0.4 M [76], the ratio GM/c 2 ρ is of the order of 10 −4 on the surface of the star. Our approximations done in this section are therefore suitable to fully apply our results (63) and (64) to study the behavior of the charged electron gas making the surface of such stars and confront them with observations regarding the mass-radius predictions [76]. Furthermore, as in the case of the Newtonian approximation done in Section 3.1, any significant deviation from the predicted corrections (63) and (64) within the range of our approximations would signal a possible deviation from general relativity. Of course, our present study does not take into account the additional corrections that would come from the spin of the charged particles but, nevertheless, the general relativistic corrections we found here would not deviate much by including the contributions of spin.

Summary and Conclusions
We have studied the effect of the gravitational field of a spherical mass on a charged particle inside a uniform and constant magnetic field. We started from the full Schwarzschild metric and derived Equation (14) that describes the dynamics of the particle for both laboratory-based experiments and astrophysical level observations. Although we subsequently simplified further the equation by keeping only the leading contribution from the Schwarzschild term GM/(c 2 ρ), going through Equation (14) was actually a crucial step for filtering out each individual contribution to the interaction potential. Furthermore, such a step allows us to clearly see how to adapt our present approach to other spacetime metrics as we show in Ref. [82] and how to extend it to include relativistic regimes as we did in Section 5.
We found that the infinite Landau degeneracy is removed as the Landau orbitals of the same Landau level split in energy. We used two independent methods to achieve our goal for the case of the Newtonian approximation of the Schwarzschild metric and both gave similar results up to the approximating scheme adopted under each approach. We pointed out that a third method, which seems at first sight to be able to yield the right quantization condition, does not actually yield a complete and consistent quantization of energy. We saw in detail where the shortcoming of such an approach occurs. We also indicated in detail why a fourth method, which is capable, in principle, of leading us to the desired quantized energies, is not really practical for finding the right quantization condition.
In addition to bringing to light the effect of the gravitational field on the Landau levels, we showed that our results could also serve another purpose which is to test any departure from the inverse-square law for gravity using quantum mechanical particles. In this paper, we restricted ourselves to the Schwarzschild spherically symmetric metric for which a Yukawa-like potential and a power-law potential are most easily incorporated. Moreover, such a spherically symmetric setup is easily achievable at the tabletop exper-iments level with, say, a two-dimensional gas of electrons or protons moving around a massive sphere immersed inside a strong magnetic field. The splitting of the Landau levels could be detected through the induced quantum Hall effect. In fact, the stationary states worked out in this paper find the most immediate laboratory application in experiments relying, for example, on the quantum Hall effect [23]. However, the applications of the stationary states explored in this paper are not limited to the laboratory-scale experiments, for the equation of state of matter making magnetized stars is studied using the Landau levels obtained using stationary states. Nevertheless, as already indicated at the end of Section 4.2, it should be stressed here that the turbulent environment around neutron stars drastically reduces the possibility of ever observing charged particles occupying stable Landau levels around the stars. The Schwarzschild vacuum geometry we used in this paper constitutes only a first idealized step towards a more rigorous study of, for example, the dynamics of the particles deep inside the stars, which require, instead, the use of the interior Schwarzschild solution [83]. In fact, the other very interesting investigation of this gravitational splitting of Landau levels would be to consider the superfluid phases of baryonic matter inside neutron stars (see, Ref. [84]). Such a topic should be the subject of our future work that will, indeed, heavily rely on the mathematical results we obtained, and methods we used, in this paper. Nonetheless, the idealized Schwarzschild geometry we started with in this paper clearly showed how gravity splits Landau levels, and it helped us extract the various approximations that might be relevant to the matter on the surface of neutron stars.
At the laboratory level, however, the Newtonian correction being of the order of GMm/ρ 0 , a spherical mass of the order of a thousand kilograms with a radius ρ 0 of one meter would only bring a splitting in the Landau levels of the order of 10 −16 eV for a proton. It should, nevertheless, be mentioned that such tiny corrections can easily be amplified if one uses instead of protons heavier quantum particles. In fact, such a possibility of using larger and heavier quantum objects has already been contemplated in relation to quantum superposition of nano-and micron-sized quantum objects. Quantum objects consisting of molecules of masses up to 10 4 amu have been successfully used in quantum interference experiments [85] and much more massive quantum objects have also been investigated in the literature [86][87][88][89][90][91][92][93][94][95], reaching masses of up to 10 13 amu [96]. With a charged molecule of such a mass, the correction to the first Landau level reaches up to about 10 −3 eV, which is easily measurable provided that one reaches down to sufficiently low temperatures. In fact, using low temperatures allows one not only to guarantee the longest coherence lengths possible, but to also prevent the gravitational splitting of the Landau levels we derived here from being overwhelmed by the inevitable temperature broadening of the Landau levels. Thus, according to our formula (24), the threshold temperature allowed at the laboratory level is given by, T 0 ∼ GMm eB/h K 2 B . For test particle masses of the order of 10 13 amu, the required temperature is of the order of 10 −4 K. On the other hand, for a Yukawa range λ of the order of the micrometer [67,71,73], the splitting caused by a Yukawa-like deviation from the Newtonian potential would be drastically reduced below the 10 −3 eV by a factor of the order of e −10 6 . For a power-law deviation, however, the correction would reach 10 −6 for an integer power s = 1 and for a parameter L of the order of a micrometer. Therefore, a significant deviation from what our results predict would certainly signal a possible deviation of gravity from the inverse-square law and would even suggest the possible form of the deviation.
We have also investigated in Section 5 the full relativistic regime that allowed us to extract the energy corrections brought to a relativistic particle moving inside the curved Schwarzschild background. We saw that the different corrections neatly separate between the purely Newtonian one and the general relativistic one. We found that the corrections brought by the latter are more important for particles with large orbital quantum numbers . Just as important is the fact that such corrections are not constrained anymore to slowly moving particles and, hence, can be confronted with observations at astrophysical levels. Indeed, as briefly discussed in Section 5, our results are valid for the case of charged particles moving on the surface of typical magnetic white dwarfs. It should be noted, however, that despite the neat results we obtained throughout this paper for the level n = 1 by relying on the large-limit, it will be of great interest to fully exploit our various results for general n and by relying, instead, on numerical methods. Now, in addition to exploiting numerical methods, our next tasks-the first of which is already presented in Ref. [82] for the non-relativistic case-are to (i) first consider the effect on Landau levels of more involved spacetime metrics, and then (ii) take into account the relativistic regime into consideration and investigate closer the implications for the highly magnetized astrophysical objects. The first allow us to use this splitting to probe in novel ways the modifications brought by GR to Newtonian gravity. The second promises to make it possible to probe any deviation from the inverse-square law when no restriction on the magnetic strengths of the astrophysical objects is imposed, for only within such general and hostile environments does the Yukawa correction have greater chances of being detected based on this splitting of the Landau levels. A closer and more specialized investigation on how to adapt our approach to the case of particles with spin and how to combine the results of Section 5 for the relativistic splitting of the Landau levels with astrophysical observations will be attempted elsewhere.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A. Evaluating Integrals Involving Kummer's Functions Using Laguerre Polynomials
In this appendix we show in detail how to evaluate the various improper integrals involving two Kummer' functions used in Sections 3 and 4. Note that an integral involving two Kummer's functions might be found solely in terms of a ratio between products of gamma functions [35]. However, for our purposes here a better strategy is to make a detour and use instead integrals involving Laguerre polynomials. In fact, given that the gamma function Γ(x) has simple poles at x = 0, −1, −2, . . . and that our Kummer's functions come with the negative integers −n, plugging our completeness condition directly into the integral formulas given in Ref. [35] would lead to ambiguous ratios between simple poles of the gamma functions.
Let us therefore start from the following relation between Kummer's hypergeometric function 1 F 1 (−n, b + 1; x) and the generalized Laguerre polynomial L n (x) (see Refs. [97,98]), On the other hand, an improper integral, consisting of the Laplace-like transform involving two generalized Laguerre polynomials, reads, [97,98], Therefore, substituting the expression (A1) of the generalized Laguerre polynomial in terms of Kummer's function, also gives, This result shows that our wavefunctions form a complete normalized set of basis for our eigenfunctions. On the other hand, we also have the following more general integral involving two generalized Laguerre polynomials [97,98]: Here, 3 F 2 (a, b, c; d, e; 1) is the generalized hypergeometric function, and it is defined by the following infinite series [35,99], where (a) k is the so-called Pochhammer symbol, defined by the product (a) k = a(a + 1) . . . (a + k − 1), and such that (a) 0 = 1. Therefore, combining Equations (A1) and (A4), we have also, Furthermore, we have the following useful addition theorem for Kummer's functions [35] that we are going to use: Therefore, by starting from result (A6) and shifting the integration boundary from 0 to x 0 and then using theorem (A7), we arrive, after a little bit of extra work, at the following important result, In the second step we have used the generalized binomial theorem to expand the term (x + x 0 ) c−1 in powers of x and x 0 for an arbitrary real exponent c, and we then used theorem (A7). In the last step we have used again result (A6).
Appendix A.1. Integrals Needed in "Using Perturbation Theory" Now, from this very general result (A8) we can extract all the useful integrals needed in this paper. In fact, performing the change of variable, x = ρ 2 , and setting b = , z = β/2, c = + 1, and x 0 = ρ 2 0 , the above result yields, where, (A10) Notice that the series in k in this last expression terminates at because in this case the exponent in (x + x 0 ) c−1 , to which we apply the binomial theorem, is the integer . The result (A10) is what is needed to find the normalization constants A n of the wavefunctions ψ n (ρ) used in Sections 3 and 4.
From Equation (A10), one can also find the needed expression M 1 by setting n = m = 1. Unfortunately, the presence of the three sums in expression (A10) makes the task very tedious. Therefore, a better strategy to get an explicit expression for M 1 is to just set n = m = 1 directly in the integral (A9) giving M mn and then evaluate the integral. The task becomes then very straightforward indeed and yields the following: In the second step, we used the definition, 1 F 1 (a, b; z) = ∑ ∞ k=0 (a) k (b) k z k k! , of Kummer's function [35]. To evaluate this integral (A11), recall the definition of the incomplete gamma function [35], ∞ t 0 t s−1 e −t dt = Γ(s, t 0 ), from which the variable re-definitions, t = β 2 ρ 2 and t 0 = β 2 ρ 2 0 , yield, In the second line we have used the explicit expression of the incomplete gamma function in terms of a sum of the usual gamma function and an infinite series [35]. Expanding the square brackets of the integral (A11) and then applying the result (A13) to each term of that sum individually, we easily find the sought-after expression, As explained in the text below Equation (23), we are interested in the large-limit of M 1 . From expression (A14) it is actually easy to estimate M 1 in the large-limit using the explicit expression of the incomplete gamma function as given by the square brackets in the second line of Equation (A13) as well as the recursion relation Γ(x + 1) = xΓ(x). In fact, we have, Note that to get to the middle step in Equation (A15), we had to discard the term β 2 ρ 2 0 from the incomplete gamma functions in Equation (A14). This becomes possible only for values of bigger than the dimensionless factor β 2 ρ 2 0 , as can be seen by examining the infinite defining series of the incomplete gamma function on the right-hand side of Equation (A13). In fact, for the ratio 1 s βρ 2 0 /2 s of the infinite series inside the square brackets to be negligible compared to the first term Γ(s) ∼ s!, one needs to have s > βρ 2 0 e/2, where e-not to be confused with the electron's charge-is the base of the natural logarithm. This condition is arrived at by evaluating the limit of the ratio 1 s! βρ 2 0 /2 s using Stirling's approximation formula for large n: ln n! ∼ n ln n − n. Thus, needs to be bigger than βρ 2 0 e/2 for the approximation (A15) to hold. The last term in expression (A15) is used in Section 3 to estimate the splitting of the energy levels for large orbitals .
Another useful integral can be extracted from the general result (A8) by performing again the change of variable x = ρ 2 and by setting this time b = , z = β/2, c = + 1 2 and x 0 = ρ 2 0 . We easily find then the following result: where, This result is what allows us to find in Section 3 the correction to the quantized energy levels E n due to the splitting of the orbitals. We would like to find from this result the expression of P 1 by setting n = m = 1. However, like for the expression M 1 , it is also here much easier instead to proceed by substituting n = m = 1 directly into the integral (A16) defining P mn and then evaluate the former. In fact, following similar steps as those that led us to the result (A14), we easily find the following expression for P 1 , From this expression it is also easy to estimate P 1 in the large-limit, as we did for the term M 1 above, using the explicit expression of the incomplete gamma function as given by the square brackets in the second line of Equation (A13) as well as the recursion relation Γ(x + 1) = xΓ(x). In fact, we have, The same remark, as the one made below Equation (A15), about how large should be for the middle step in the approximation (A19) to be justified, holds here also. The right-hand side of the last equality is used in Section 5 to estimate the splitting of the relativistic energy levels for large orbitals .