Rydberg states of hydrogenlike ions in a braneworld model

Precise measurements of optical transition frequencies between Rydberg states of hydrogen-like ions could be used to obtain an improved value of the Rydberg constant, by avoiding the uncertainties about the proton radius. Motivated by this perspective, we investigate the influence of the gravitational interaction on the energy levels of Hydrogen-like ions in Rydberg states in a braneworld model. As it is known, in this scenario, the gravitational interaction is amplified in short distances. We show that, for Rydberg states, the main contribution for the gravitational potential energy does not come from the rest energy concentrated on the nucleus but from the energy of the electromagnetic field created by its electric charge. The reason is connected to the fact that, when the ion is in a Rydberg state with high angular momentum, the gravitational potential is not computable in zero-width brane approximation due to the gravitational influence of the electrovacuum in which the lepton is moving. Considering a thick brane scenario, we calculate the gravitational potential energy associated to the nucleus charge in terms of the confinement parameter of the electric field in the brane. We show that the gravitational effects on the energy levels of a Rydberg state can be amplified by hidden dimensions even when the compactification scale is shorter than the Bohr radius.


Introduction
In the past decades we have seen a renewed interest in extradimensional theories stimulated by the braneworld scenario, according to which our ordinary universe is a submanifold embedded in a higher-dimensional space [1][2][3][4]. The fundamental characteristic of these models is the assumption that matter and fields are confined in a three-dimensional space (the 3-brane), while gravity has access to the whole space, implying that the gravitational field could feel direct effects of the extra dimensions on a length scale much greater than the a e-mail: fdahia@fisica.ufpb.br scale where the standard model fields are directly affected. Since the gravitational interaction has been tested in short distance only recently, then braneworld scenario may allow the formulation of phenomenologically viable models in which the hidden dimensions can be much larger than the fourdimensional Planck length (10 −35 m), which is the characteristic size of the fifth direction postulated by the Kaluza-Klein model, the first modern higher-dimensional theory.
Direct laboratory tests of the inverse square law impose that the compactification radius (R) should be less than 44 µm [5,6], for instance. This is the most stringent constraint when the supplementary space has only one extra dimension. Regarding greater codimensions, experimental limits from astrophysics and colliders are tighter, but, according to current data, they are still greater than Planck length by 15 orders at least [7][8][9][10][11][12][13][14].
The original purpose of braneworld theories was to explain the hierarchy problem, i.e., why gravity is so much weaker in comparison to the other fundamental forces [1][2][3][4]. According to these models, the reason, roughly speaking, is that gravity is the only field that spreads in every direction, hence it seems weaker than the other interactions in the long distance scale, where the usual four-dimensional behavior is recovered.
On the other hand, in distances smaller than the size of the extra dimensions, these models predict that the gravitational force is amplified in comparison to the Newtonian three-dimensional force by a factor that depends on the number and on the size of the additional dimensions. This fact has motivated many authors to study the effects of the gravitational field in the atomic and molecular system as a way to obtain independent experimental bounds for parameters of the higher-dimensional theories [15][16][17][18][19][20][21][22].
Here we intend to study the influence of the extra dimensions on Rydberg states of hydrogen-like ions, motivated by recent developments in the spectroscopy area which suggest that precision measurements of optical transitions between Rydberg states can become an efficient method to test QED theory and also to determine the value of the Rydberg constant more accurately [23].
The current uncertainty of Rydberg constant is of the order of 10 −12 [24] and its value is extracted from the comparison of theoretical predictions and measurements of transition frequencies in hydrogen and deuteron, which includes transitions between S-states [23,24]. It happens that, regarding S-states, theoretical calculations are mainly limited by the uncertainty about the proton charge radius [23], whose value is also a source of a controversy known as the proton radius puzzle [27][28][29][30], which is a discrepancy between measurements of the proton size as inferred from the muonic hydrogen spectroscopy and the CODATA value [24], based on the proton-electron interaction.
In Rydberg states with high angular momentum, as the influence of the internal structure of the nucleus is negligible, the predictions are practically independent of the proton radius and then measurements would be free from those uncertainties. With recent calculations of QED corrections for higher-quantum number states, the accuracy of the predicted energy level could reach parts in 10 17 [23]. At the same time, advances on experimental techniques, such as optical frequency combs [25], promises to reduce the relative uncertainties to the order of 10 −19 in measurements of transition frequencies around the optical range [26]. Together, these theoretical and experimental advances may allow precise measurements of certain optical transitions between appropriate Rydberg states which, in principle, could lead to an improved value of the Rydberg constant [23].
Faced with such expectations, it seems relevant to investigate as how the extra dimensions could interfere in the energy level of the Rydberg states in the braneworld scenario. The effects of gravitational interaction on the S-states have already been considered before [22]. An interesting point to mention here is that when the space has more than two extra dimensions, the mean gravitational potential energy of a hydrogen-like atom in a S-state, U S , diverges if the brane structure is not taken into account. Indeed, this is connected to fact that the gravitational potential, ϕ, produced by the proton mass is not computable in the interior of the nucleus in the approximation of zero-width brane [9,22,31].
On the other hand, when the atom is in a Rydberg state the situation may be different. For a state with a high enough angular momentum, the leading term of the mean gravitational potential energy, U Ry , is finite in the thin brane limit. So the zero-width brane is a valid idealization here. However, under such circumstances, the amplification of the gravitational potential energy by extra dimensions is significant only if the compactification radius is greater than the atomic Bohr radius. Thus, considering the current constraints on R [14], at the first sight, we could be led to think that extra dimensions would have little influence on Rydberg states. Nevertheless, it is important to realize that the geometry around the nucleus should be similar to the Reissner-Nordstrom geometry, due to the nucleus electric charge. In a Reissner-Nordstrom spacetime, besides the gravitational potential, ϕ, produced by the rest mass, there is also the gravitational influence associated to the energy of the electromagnetic field created by the electric charge of the source, which, in the weak field regime, can be described by a certain potential χ .
Although, in the traditional three-dimensional picture, χ is smaller than the potential ϕ, this relation in the braneworld scenario can be different in a certain region outside the nucleus. As the electromagnetic field is spread in the threedimensional space, the potential χ diverges in every point of the brane, in the zero-thickness limit. Thus, in order to compute χ , the distribution of the electric energy inside the brane should be considered. Addressing this problem in the thick brane scenario, we consider the Green function for the gravitational potential associated to length scales smaller than R, where the effects of extra dimensions are stronger. With this approximated Green function, we determine the short-distance potential χ s . We find that it can be greater than the short-distance potential ϕ s outside the nucleus. This can lead to interesting consequences. Indeed, considering a Hydrogen-like ion in a Rydberg state, we show that the extra dimensions can amplify the gravitational potential energy of the ion even when the compactification scale R is smaller than the Bohr radius, due to the behavior of χ s .
Finally, we estimate the effects of higher-dimensional gravity on particular optical transitions of Hydrogen-like ions.

The gravitational field produced by the nucleus
In the original ADD-braneworld model [1], the spacetime has δ additional spacelike dimensions with the topology of a torus T δ . In the ground state of the model, the extra dimensions have a certain radius R, and it is assumed that the metric is flat. This means that the energy of the brane itself does not curve the bulk at long distances compared to the length scale of the brane, by the intervention of some mechanism [1,2]. It is also admitted that the gravitational field produced by the matter confined in the brane is governed by the Einstein-Hilbert action extended to higher dimensions: whereR is the scalar curvature of the bulk,ĝ is the determinant of the metric whose signature is assumed to be (−, +, . . . , +) and G D is the gravitational constant defined in the higher-dimensional space. The coordinates x and z refer to parallel and transversal directions with respect to the brane. Due to the topology of the supplementary space, the metric is periodic in the extra directions and can be expanded in a Fourier series with respect to the z-coordinates, giving rise to the so-called KK-modes. The zero-mode is supposed to reproduce the four-dimensional gravitational field at large distance in comparison to R. This condition requires that G D should be related to the Newtonian gravitational constant G according to the following formula [1,9]: Additionally, in order to get the correct Newtonian limit, it is also necessary that some mechanism ensures the stabilization of the volume of the supplementary space at large distance [35,36]. The extremization of the action (1) yields the higherdimensional version of the Einstein equations. In the weak field regime, which we assume to be valid in the atomic domain, the metric can be written as g AB = η AB + h AB . Here the capital Latin indices run from 0 to 3 + δ, η AB is the Minkowski metric and h AB is a small perturbation of the order of G D M. In a coordinate system in which the gauge condition is satisfied, the linearized equations reduce to the form where is the D'Alembertian operator associated to the Minkowski metric andT AB = T AB − (δ + 2) −1 η AB T C C is defined in terms of the energy-momentum tensor T AB of the source.
Considering the topology R 3 ×T δ , the solution of equation (4) for static sources is where X = ( x, z), K i = 2π R (0, 0, 0, k 1 , . . . , k δ ) and each k i is an integer number. If we consider T δ as a manifold embedded in R δ , then the vectors K i may be viewed as the localization of the mirror images of the source induced by the topology of the supplementary space on R δ (the covering space). The presence of the mirror images makes the solution periodic with respect to the z-coordinates as demanded by the topology. It can be shown that the higher-dimensional Green function recovers the four-dimensional behavior for long distances | x| >> R [33]. On the other hand, in short distance ( x − x < R ), the Green function is dominated by the first term of the expansion (5). In this paper, we are interested in studying the effects of the short-distance behavior of the gravitational interaction on the atomic energy levels.
Thus for the sake of simplicity, we are going to take only the first term ( K 0 = 0) of the series (5) as an approximation of the solution. Of course, this approximation gives a low estimate of the gravitational effect, since all terms of the series (5) we are neglecting have the same sign. It is also important to mention here that in order to be phenomenologically consistent the theory should provide some mechanism to stabilize the volume of the supplementary space. There are different mechanisms which can ensure an appropriate behavior of the radius of the extra dimension. In general, the effect of the stabilization is to generate some mass for the radion field [35][36][37][38]. This is necessary because long-range effects (r >> R) of a massless radion field would lead to an effective gravitational coupling with a different value from the Newtonian constant G N [1]. For a large class of stabilization mechanisms the radion mass m is smaller than the scale of the compactification 1/R [35][36][37][38]. This means that gravitational effects of the radion can extend to a length scale greater than that of KK-modes of the graviton and, in principle, could be tested by Cavendish-type experiments [32,36,39]. However in the short range r << R, the effects of the radion mass is negligible. Indeed, in the short-distance domain, the supplementary directions could be treated as a space with the non-compact topology R δ and, therefore, the green function for gravitational field reduces to the first term of (5). This is the usual procedure adopted in the study of black-hole production in colliders [40][41][42]. In this paper, as we are studying the effects of the short-distance gravitational potential, then the mass of the radion can be ignored in the following calculations.
In the context of the braneworld, it is assumed that the energy-momentum tensor of the confined fields has the form [9]: where T μν (x) is the ordinary energy-momentum tensor of the four-dimensional fields that live in the brane and f (z) is a normalized distribution very concentrated around the brane, which is approximately given by a delta-like distribution in the thin brane limit. The function f (z) describes the confinement of a field in the brane and, in principle, it could be a different function for each type of field.
In our system, the atomic nucleus is the source of the gravitational field. Due to its electric charge, it is reasonable to expect that the spacetime geometry around the nucleus should be similar to the Reissner-Nordstrom geometry. Based on this, it is convenient to decompose the energy-momentum tensor as a sum of two terms: μν describes the rest energy concentrated inside the nucleus and T

(E M) μν
represents the stress-energy tensor of the electromagnetic field created by the charge, which is spread in the space.
We deal with each term separately. The first tensor can be written as T where ρ is the mass density of the source. On its turn, in SI units, T where F μν is the electromagnetic tensor and 0 is the electric permittivity of the free space. In a first approach, let us ignore the contribution of the magnetic field produced by the proton's magnetic dipole. In this approximation, only the components of the electric field, F 0i = E i /c, are non-null (in the zero order of G D ). Take into account all this consideration, we can show that the metric can be written as: where x 0 = ct, the coordinates (r, θ, φ) are the usual spherical coordinates associated to the "almost Cartesian" coordi- . The function ϕ s plays the role of a Newtonian potential produced by the nuclear mass in R 3+δ : The sub-index s emphasizes that ϕ s is the potential given by the shortdistance Green function, which corresponds to first term of the series (5). Analogously, χ s is the short-distance gravitational potential produced by the energy of the electromagnetic field u = ( 0 E 2 /2) created by the electric charge: The spatial components, , of the electromagnetic stress-energy tensor, give rise to the functions λ 2,s and λ 1,s , which are defined by: As we shall see in the next section, in the gravitational sector of the atomic Hamiltonian (H G ), the leading term comes from the component g 00 , which depends on the potential ϕ s and χ s . To determine explicitly the functions ϕ s and χ s , it is important to look at the internal structure of the brane to see how the fields are localized inside. In a field-theoretic framework, topological defects are possible realizations of a brane, since these structures are capable of localizing fermions in their cores, as illustrated in Ref. [34], where it is shown that Dirac fields can be trapped inside a domain wall by means of a Yukawa-like interaction. In this context, usually known as thick brane scenario, a delta-like localization is replaced by a non-singular confinement where the fermion's states are described by a regular wave-function with a tiny width σ in the transversal directions.
Following these ideas, in Ref. [22], we calculated the gravitational potential, ϕ s , produced by the proton mass M p in the thick brane scenario. By admitting that the proton wavefunction in the transversal directions has a Gaussian profile, we have estimated the influence of the gravitational interaction on the energy of S-states of a Hydrogen atom. The leading contribution is proportional toĜ D m M p /a 3 0 σ δ−2 , where m is the electron mass and a 0 is the Bohr radius. It is clear, from this expression, that, the calculation diverges, even at the tree level, in the thin brane limit. Therefore, the mass distribution of the nucleus inside the brane cannot be neglected when the atom is in S-states. We should highlight that the mentioned term corresponds to the mean value of the gravitational potential energy integrated in the interior of the nucleus, which we shall denote as H G in . On its turn, outside the nucleus, the dominant contribution is proportional toĜ D Mm/a 1+δ 0 , if the compactification scale is greater than the Bohr radius. As, in realistic scenarios, σ << a 0 then the interior contribution H G in is much greater than the exterior contribution H G out for the energy of S-states.
However, if the atom is in a Rydberg states with large angular momentum, the inverse happens. Indeed, in states with angular momentum l, the internal contribution is reduced by a factor of the order of (r N /a 0 ) 2l , where r N is the radius of the nucleus. Therefore, for some angular momentum higher than the exterior contribution H G out become greater than H G in . Therefore, as the leading term does not depend on the brane thickness, we may say that a zero-width brane is a valid approximation for such states with high quantum numbers. Replacing f m (z) in Eq. (8) by a delta-Dirac distribution, we find that, in the exterior region, the short-distance gravitational potential produced by the nucleus mass is With respect to the complete potential given by the series (5), the term (13) is the leading contribution in the region r N << r << R. By using (2), we can conclude that the contribution of (13) to the gravitational potential energy of the atom in a Rydberg state will be of the order of (G Mm/a 0 ) (R/a 0 ) δ and, therefore, extra dimensions would amplify significantly the gravitational energy of the atom in a Rydberg state only if R >> a 0 . This conclusion is not necessarily valid when we consider the potential χ s produced by the electromagnetic energy.
If f e (z) is approximated by a delta-like distribution in Eq. (9), the potential χ s diverges everywhere, not only inside the nucleus, as happens with the potential ϕ s , since the electromagnetic field is spread in the space. Therefore, due to the behavior of potential χ s , the zero-width brane is not a valid idealization even when we are dealing with ions in Rydberg states.
Thus, in order to compute χ s , we have to consider the distribution of the electric energy inside the thick brane. With the purpose of obtaining some estimates, we are going to admit that the electric energy is uniformly distributed inside a compact region of the brane with a size ε, which may have the same order of the brane thickness. Thus, if V δ (ε) denotes the volume of a ball with a radius ε in the supplementary space, f e (z) can be defined as the step function: Taking this function in the expression (9) and integrating it with respect to the angular coordinates, the potential χ s evaluated in the brane ( z = 0) can be written as Here S δ is the hyper-area of a spherical hypersurface of δ dimensions with a unit radius. In (15), the integration interval of the transversal variable is 0 < z < ε. The divergent term in the thin brane limit comes from the function χ − as the integrating variable r passes through r . In the context of a thick brane scenario, this "dangerous" term can be isolated and calculated explicitly. As we are interested in studying Rydberg states with high angular momentum, we may restrict our analysis to points far from the nucleus (r >> r N ). In the calculation of χ s , it is convenient to separate the integration domain in two parts: the region inside the interval r − r ≤ R, where the short-distance behavior of the Green function is dominant, and the external region . Noticing that the zero-thick brane approximation is valid in , we can write: To proceed further we need now to specify the electromagnetic energy distribution in the three-dimensional space by means of the function u. To be consistent, the model should reproduce the ordinary behavior of electromagnetic field in a length scale greater than the brane thickness (r >> ε). It is important to remark at this point that the compactification scale (R) is admitted to be greater than ε in our approach, so the gravitational field is the only field directly affected by extra dimensions in the region r >> ε. Based on this considerations, we will assume that in the three-dimensional space the energy density of the electric field has the usual form given by: which is the well-known energy of electromagnetic field produced by a charge Q uniformly distributed inside a ball of radius r N in three-dimensional space. With this energy distribution, we find: where β δ can be written in terms of the gamma function as: and the corrections terms are of the following order: At this point, it is important to remark that the exact profile of the electromagnetic energy distribution inside the brane, given by the function f e (z), should depend on the localization mechanism of the gauge field. However, it is interesting to mention that for any normalized function f e (z), the potential χ s , in the leading order, has exactly the same form given by (18), provided that the parameter ε be defined as Therefore, from the phenomenological point of view, the gravitational potential χ s , in leading order, does not depend on many details of the localized fields, but only on a specific moment 1/z δ−2 f of the distribution f e , i.e., on the average value of 1/z δ−2 with respect to the distribution f e ( z). As f e (z) should be concentrated inside the brane, it is reasonable to expect that ε is related to the brane thickness. Now, we can compare the strength of the potentials χ s and ϕ s . We can see that χ s is greater than ϕ s if r > 4 β δ For realistic values of ε, the above condition is satisfied outside the nucleus. So, in this scenario, the gravitational potential produced by the electromagnetic field is greater than the short-distance term of the gravitational potential of the nuclear mass, in the exterior region. In three-dimensional space, this does not happen. Indeed, considering the same distribution (17), we find in the exterior region: where E is proportional to the energy of the electromagnetic field. The attractive term of χ (3) depends on the inverse of the distance r , therefore, it can be incorporated in the three-dimensional gravitational potential ϕ (3) by absorbing the electromagnetic energy as part of the rest energy of the system. It follows then that the potential, ϕ (3) = −G M/r , with M = M 0 + E/c 2 , is greater than the repulsive part of the potential χ (3) outside the nucleus.
At this point, let us mention that the mass distribution of the nucleus could be described by a continuous density ρ m , instead of a compact distribution, without changing our conclusion. In this new configuration, the potential ϕ s would diverge everywhere in the thin brane limit, but, in the thick brane scenario, the "dangerous" term would be proportional toĜ D ρ m /σ δ−2 . Thus, it will be smaller than χ s for r >> r N , provided that ρ m is a fast decreasing function, like an exponential.
In this section, we have calculated the gravitational field produced by the nucleus. In order to obtain finite results, it was necessary to take into account the distribution of the nuclear mass and the electromagnetic field inside the brane. In the next sections, we intend to discuss the complementary aspect of the picture, namely, the influence of gravity on the dynamics of the electromagnetic and the Dirac fields in this scenario.

The electrostatic potential on the brane
As we have mentioned before, in the ADD-model, it is assumed that the effects of extra dimensions appears first in the gravitational field, since it can spread into the extraspace, which has a size R that is much greater the brane width ε, where the other fields are confined. So it is reasonable to expect that direct effects of extra dimension on the other fields would be smaller than the gravitational effect by a factor of the order of (ε/R) δ . Thus, in a first approximation picture, we may assume, as usual, that the standard model fields can be treated as traditional four-dimensional fields confined in a 3-brane.
Based on this considerations, it follows that, in length scales above the brane thickness, the fields do not couple to the bulk geometry directly, but their dynamic in the spacetime is influenced by gravity through the brane geometry. According to the present approximation scheme, the geometry of the 3-brane is isometric to the geometry of thick brane center (i.e., the submanifold z = 0), since any deviations in the metric components inside the brane are negligible given that they are of the order of ε, at least. Thus, in the given coordinate system, the induced metric on the submanifold z = 0 can be directly obtained from (7). After a transformation to isotropic coordinates, the induced metric can be written in the following form: where In this context, let us determine the modification on the electrostatic interaction between the electron and nucleus caused by the geometry of the brane. The Maxwell equation in the curved space can be written as where g is the determinant of the metric tensor (22), J ν is the four-current and μ 0 is the vacuum magnetic permeability. If U μ = dx μ /dτ is the four-velocity field of the source, then J μ = −ρ e U μ , where ρ e is the proper charge density and τ is the proper time. For the sake of simplicity, let us admit that the nucleus is at rest in the given frame. Thus, the normalized four-velocity has the following components: U μ = cw −1 , 0 , in the given coordinate system. In the static regime, the four-potential reduces to A μ = (−φ/c, 0). Considering that F μν = ∂ μ A ν − ∂ ν A μ , we can derive from (25) the field equation for the electric potential φ: where ρ e,0 = v 3 ρ e is the charge density with the flat measure d 3 x. This means that in a hypersurface orthogonal to U μ , it satisfies the condition: As the gravitational field is weak, we expect that the gravitationally modified electric potential can be written as φ = φ 0 + φ G , where φ 0 is the solution in the flat space sourced by ρ e,0 and φ G is a small correction of the order of G D due to the spacetime curvature. In the first order of G D , the field Eq. (26) reduces to the form where 0 ∇ 2 is the Laplacian in the flat three-dimensional space and plays the role of a polarization charge density which works here as a source of the gravitational correction of the electrostatic potential. The solution is which is proportional toĜ D M Q 2 /c 4 in the leading order.

Dirac equation on the brane
The localization of spinors fields in the core of topological defects is a very interesting and involving problem [34,[43][44][45]. In Ref. [34], this localization mechanism is illustrated by a procedure that is capable of trapping a five-dimensional massless spinor in the center of a domain wall. Due to a Yukawa-type interaction between the spinor and scalar field that describes the wall, there is a chiral zero-mode Weyl spinor which is normalizable and localized in the center of brane. It is assumed that this localized spinor, which acquires some mass by means of another mechanism, describes the matter in the low-energy regime. The fundamental theory in high-energy regime could be very complex and possibly contains additional degrees of freedom (the KK-modes, for instance). However, whatever the theory in the high-energy regime, the localization mechanism must ensure that the long-range degrees of freedom should reproduce the known four-dimensional behavior [31,46].
Based on this phenomenological perspective, we are going to assume here that at a length scale greater than the brane thickness, the matter is described by four-dimensional Dirac fields confined to the center of the brane. Thus, to take into account the gravitational interaction between fermions, we shall admit that the confined particles do not feel the geometry of the bulk, but instead they interact with the metric of the (3 + 1)-spacetime given in (22).
Associated with that metric g μν , we find a set of four orthonormal vector fields, whose components, e Assuming a minimal coupling with the gravitational field, the Dirac equation (in a curved space) that describes the state of a fermion with mass m and electric charge q subject to gauge field A μ is given by: where the operator D μ = ∇ μ −iq A μ and ∇ μ is the covariant derivative of the spinor, which can be written in terms of the spinorial connection μ (x) as On its turn, μ (x) is defined as where σ ab = i 2 γ a , γ b is a representation of the Lorentz Lie Algebra in the spinor space. As usual, the symbol [, ] is the commutator operator and e bν;μ (x) is the covariant derivative of the vierbein fields which depend on the Levi-Civita connection λ νμ (the affine connection compatible with the metric) according to For the diagonal metric (22), a possible choice for the vierbein fields are: and all the other components equal to zero. By a direct calculation, the Dirac equation can be written in the form ih ∂ψ ∂t = H ψ. The operator H is the atomic Hamiltonian, which, in a convenient representation, assumes the following form in the first order of G D [47]: Here p is the usual three-dimensional momentum operator in flat spacetime, α i = γ 0 γ i and β = γ 0 . For the sake of simplicity, we are neglecting the effects of the potential vector A. It is also important to mention that, in the chosen representation, the Hamiltonian (36) is Hermitian in the Hilbert space of square-integrable functions endowed with the usual inner product calculated with the flat measure d 3 x.
Following the Foldy-Wouthuysen procedure, the nonrelativistic limit of the Hamiltonian (36) can be obtained, by admitting that the mean value of each term of H is much smaller than the rest energy of the test particle, mc 2 , in Rydberg states. In the gravitational sector of the Hamiltonian, i.e., in the part of H constituted by terms proportional to G D , we find that the leading terms are proportional to the test particle mass. Thus, in this approximation, we can write: As χ s is greater than ϕ s in the exterior region, then we can consider a further approximation: in order to find corrections on the energy levels of Rydberg states, coming from the short-distance behavior of the higherdimensional gravitational field produced by the energy of the electrovacuum surrounding the test particle in a Hydrogenlike ion.

Results and discussion
Treating H G as a small term of the total Hamiltonian, we can use the perturbation method to estimate the gravitational potential energy of a Hydrogen-like ion in a Rydberg state.
For an ion with an atomic number Z found in a state whose principal number is n and l is the angular momentum, the mean gravitational potential energy is approximately given by where e is the fundamental electric charge and Actually, it is important to stress, the energy (39) is provided by the short-distance behavior of the Green function for the gravitational potential. Without extra-dimensions, the gravitational potential energy of the ion in the same state is where A is the mass number. Here for the sake of simplicity, we admit that the mass of the nucleus is approximately equal to AM p . Comparing (39) and (41), we can find conditions under which the amount of gravitational potential energy coming from the short-distance behavior of the higher-dimensional potential surpasses the three-dimensional value. In this domain, extra dimensions may amplify significantly the gravitational effects on the ion. The condition H G > H (3) G can be expressed in terms of the size of the extra dimensions R in the following form: Notice that, in order to the gravitational potential energy be amplified by the extra-dimension, R should be some orders of magnitude greater than the brane thickness, but interestingly the compactification scale can be smaller than the Bohr radius. As an example, consider the Neon ion ( 20 Ne +9 ) in a six-dimensional space. In a state with n = 15 and l = 14, the extra dimensions would provide an amplification if R > 10 −14 m, when ε = 10 −20 m, for instance. Now let us consider the transition between the states (n, l = n − 1) and (n − 1, l = n − 2). The principal part of the energy gap between these levels is From (39), we can calculate the difference of the gravitational potential energy between the mentioned levels ( E G ). Relative to the principal part of the energy gap, we find If we consider muonic hydrogen-like ions, then the Bohr radius of the ground-state of the Hydrogen, a 0 , is replaced by the Bohr radius of the muonic Hydrogen, a 0,μ , and the electron mass m is substituted by the muon mass m μ . Due to this, the magnitude of expression (44) is multiplied by a factor a 0 /a 0,μ 3 m μ /m ∼ 10 9 .
According to Ref. [23], if the transition frequency lies in the optical band, the relative experimental precision could reach the order of 10 −19 [26]. Given this fantastic precision, it is interesting to investigate the effects of the hidden dimensions on the transitions between Rydberg states.
In a previous paper [48], it was suggested that the proton radius puzzle can be explained in the braneworld scenario. Briefly speaking, the puzzle is a discrepancy between the proton charge radius recommended by CODATA (r (C D) p ) − which is obtained from data of hydrogen and deuteron spectroscopy combined with data of electron-proton elastic scattering − and the radius inferred from the measurement of the Lamb shift of the muonic Hydrogen [27][28][29][30]. The point is that the experimental value of the energy difference between 2S and 2P states, in the muonic hydrogen, is approximately 0.3 meV greater than theoretical prediction calculated with r C D p . This excess of energy, which has no explanation within the standard model of physics, may be an indication of a missing interaction between the proton and muon.
In Ref.
[48], it was argued that the muon-proton gravitational interaction, modified by extra-dimensions, can account for the unexpected energy excess measured in the muonic Hydrogen Lamb shift. Considering a thick brane scenario, as described here, the mean gravitational potential of the muonic hydrogen in states 2S and 2P can be calculated. Then, to solve the proton radius puzzle, it is required that the difference of the gravitational potential energy be equal to  Fig. 1 In transitions between Rydberg states of muonic Hydrogenlike ions below the dashed lines, the influence of hidden dimensions, through the classical higher-dimensional gravity, can be greater than promised experimental uncertainty in measurements of transition frequencies around the optical band (100-1000 THz). Here we are assuming that G D /σ δ−2 has the appropriate value to solve the proton radius puzzle [48] energy excess of 0.3 meV. This leads to the condition: where γ δ is a factor whose values depend on the number of dimensions [48]. The Eq. (45) imposes a condition for the ratio,Ĝ D /σ δ−2 , between the fundamental gravitational constant G D (or equivalently the fundamental Planck mass M D ) and the parameter σ , which here can be seen as the confinement parameter of the nucleus in the brane.
Taking this value ofĜ D /σ δ−2 as reference in Eq. (44) and admitting that the matter and the electromagnetic fields have the same confinement parameter (ε σ ), we can make some estimates. We show in Fig. 1 what are the optical transitions between Rydberg states in which the effect of the higher-dimensional gravity exceeds the experimental precision promised by the optical metrology [26].
It is important to remark that these predictions are based on the classical behavior of the gravitational field. On the other hand, as pointed out in Ref. [9], quantum-gravity effects would become relevant at a length scale around the fundamental Planck length ( D ) defined in the higher dimensions or even in greater distances. If this is the case, then quantum-gravity effects could modify the classical result unpredictably. Nevertheless, as a definite quantum-gravity theory is not yet known, only experiments can answer this question.

Final remarks
The prospect of measuring optical transitions between Rydberg states with a relative uncertainty near to the impressive order of 10 −19 motivated us to investigate the influence of hidden dimensions on the energy levels of the Rydberg states of hydrogen-like ions.
In the braneworld scenario with large extra dimensions, gravity is the first interaction to be affected by the supplementary space. These models predict that, in comparison to the three-dimensional Newtonian interaction, the gravitational field can be hugely amplified in short distance. In this paper, we have investigated the effects of extra dimensions in the atomic domain by studying the influence of the higherdimensional gravitational field produced by the nucleus on the energy levels of the ion.
Due to the electric charge of the source, the spacetime around the nucleus is a brane-version of the Reissner-Nordstrom geometry. In this kind of spacetime, the metric, in the weak field regime, depends on the gravitational potential produced by the rest energy concentrated on the nucleus, ϕ, and also on the gravitational potential generated by the energy of the electric field created by the nucleus charge, χ .
In the zero-thickness approximation, the short-distance potential χ s diverges everywhere in the brane. Due to this, the potential gravitational energy is not computable in the thin brane limit even when the ion is in Rydberg states with high angular momentum.
Therefore, in order to estimate the influence of the potential χ s on the energy of the ion, we had to address this problem in the thick brane scenario. Considering a regular distribution of the electromagnetic energy inside a brane strip, we calculate the gravitational potential energy in terms of the confinement parameter of the electric field. We find that the extra dimensions are capable of amplifying the gravitational potential energy of the ion in a Rydberg state even when the compactification radius is smaller than the Bohr radius. This a consequence of the behavior of the potential χ s .
It is interesting that the short-distance behavior of the Green function associated to the gravitational potential plays an important role in the present situation. At the first sight, we could be led to think that the short-distance behavior is irrelevant when the ion is in the Rydberg state, since the probability to find the test particle close to the nucleus is very small in this state. Indeed this is the reason why the potential ϕ s is weak here. However, the source of the potential χ is spread in the space. Therefore, around any external point x, there is an interval smaller than R, in the electrovacuum, where the short-distance Green function is dominant. It is through this term that the gravitational influence of the electrovacuum surrounding the test particle can increase significantly the gravitational potential energy of the ion.