Asymptotic Generalized Extended Uncertainty Principle

We present a formalism which allows for the perturbative derivation of the Extended Uncertainty Principle (EUP) for arbitrary spatial curvature models and observers. Entering the realm of small position uncertainties, we derive a general asymptotic EUP. The leading 2nd order curvature induced correction is proportional to the Ricci scalar, while the 4th order correction features the 0th order Cartan invariant Psi^2 (a scalar quadratic in curvature tensors) and the curved space Laplacian of the Ricci scalar all of which are evaluated at the expectation value of the position operator, i.e. the expected position when performing a measurement. This result is first verified for previously derived homogeneous space models and then applied to other non-trivial curvature related effects such as inhomogeneities, rotation and an anisotropic stress fluid leading to black hole"hair". Our main achievement combines the method we introduce with the Generalized Uncertainty Principle (GUP) by virtue of deformed commutators to formulate a generic form of what we call the Asymptotic Generalized Extended Uncertainty Principle (AGEUP).


I. INTRODUCTION
The standard uncertainty principle of quantum mechanics in its fundamental form does not take into account effects which are expected to arise from an underlying theory of quantum gravity. Taking inspiration from string theory [1,2], the Heisenberg uncertainty principle was generalized by the inclusion of the gravitational photon-electron interaction (acceleration) leading to the Generalized Uncertainty Principle (GUP) [3][4][5][6].
As usually assumed, the GUP takes into account the gravitational uncertainty of position related to the minimum fundamental length scale in physics. However, it may not be the only gravitationally induced change. In fact, the curvature of space-time does exert an influence over quantum mechanical uncertainty relations. This is the regime of the Extended Uncertainty Principle (EUP) [7][8][9][10] which takes into account the uncertainty related to the background space-time.
Both components can be formulated in terms of the standard deviations of position x and momentum p and combined to yield the most general Generalised Extended Uncertainty Principle (GEUP) [11,12] σ x σ p ≥ 2 1 + α 0 l 2 where l p denotes the Planck length, the Planck constant, r c some curvature scale related to the background space-time and α 0 and β 0 the GUP-and EUPparameters, respectively. * Mariusz.Dabrowski@usz.edu.pl † fabian.wagner@usz.edu.pl The GEUP as given by (1) has its heuristic Newtonian analogue [12] which represents the acceleration of an electron induced by both the gravitational interaction of a photon of energy E and the cosmological Hubble horizon r H = c/H = (Λ/3) 1/2 of de Sitter spacë with the cosmological constant Λ, the photon-electron distance r, the Hubble parameter H, the gravitational constant G and the speed of light c (the latter two will be set equal to 1 throughout the remainder of this paper). A derivation of the EUP based on the notion of geodesic balls was performed in Refs. [13,14]. It reflects the influence of spatial curvature on quantum-mechanics on 3-dimensional spacelike hypersurfaces of space-time. The method was applied to homogeneous and isotropic geometries of constant curvature K and the corresponding EUP was calculated.
In our recent paper [15] we applied this method to calculate the EUP for Rindler and Friedmann horizons and obtained corrections to the Hawking temperature and Bekenstein entropy of black holes. In said derivation emphasis was put on the local observer in an accelerated frame or at the center of symmetry seeing the effects related to such a specific choice of frame. Yet, a coordinate-independent (covariant) relation for the EUP was the same as derived earlier in [14]. This was further commented on in Ref. [16].
The purpose of this paper is to derive a general expression for the EUP in the case of small position uncertainties (to be specified below) and thus consider non-trivial (e.g. non-homogeneous) effects of curvature on the EUP.
The paper is organized as follows. Section II outlines a method of deriving the EUP as it was applied to homogeneous (constant curvature) spaces in Ref. [14]. In Section III we present our perturbative approach to apply it further to small position uncertainties in Section arXiv:2006.02188v2 [gr-qc] 9 Jun 2020 IV. Moreover, in Section V several types of perturbations are discussed to provide exemplary applications of the asymptotic EUP to various geometries of non-trivial curvature. Finally, Section VII is intended to summarize our results.

II. BACKGROUND GEOMETRY DETERMINED EUP
The reasoning behind Refs. [13,14] restricts our perspective to spacelike hypersurfaces of the underlying space-time (it is impossible to even define a mathematically sound uncertainty relation in standard quantum mechanics otherwise). Thus, as other effects which combine quantum mechanics and gravity, this method is observer dependent.
Consider a free wave function ψ living on said Riemannian manifold but confined to a geodesic ball (B ρ ) of radius ρ. Due to their invariance under diffeomorphisms, geodesic balls are the natural generalization of the Heisenberg slit to curved manifolds where ρ corresponds to the slit width thus providing a measure of position uncertainty. Confinement to this domain is ensured by imposing Dirichlet boundary conditions on the wave function.
As the Laplace-Beltrami-operator ∆ (representing the squared momentum operatorp 2 = − 2 ∆) is hermitian with respect to the measure dµ of the Hilbert space in question H = L 2 (B ρ ⊆ IR 3 ,dµ) and thus its eigenvalues provide an orthonormal base of the latter, deriving the uncertainty relation basically boils down to solving the eigenvalue problem Choosing the wave function to be real (the eigenvalue problem is the same for the real and the imaginary part), the Dirichlet boundary conditions ensure that p = 0. Hence, the momentum uncertainty reads with the 1 st eigenvalue of the eigenvalue problem (3). Multiplying by ρ, the uncertainty relation is obtained. In Ref. [14] it was found that the uncertainty relation for Riemannian 3-manifolds of constant curvature K reads Note that the uncertainty relation derived this way is not of the same kind as the one described by (1) because it features the characteristic length of confinement ρ, a generalisation of Heisenberg's slit width. Thus, ρ should rather be interpreted as uncertainty and does not represent the standard deviation of position.

III. WEAK CURVATURE EUP FORMALISM
Consider a generic weak curvature effect on spacelike 3D-hypersurfaces, i. e. for a spatial metric splitting where g (0) ij denotes the flat space Riemannian metric in some set of coordinates x i , i = 1, 2, 3. The perturbation will be treated at 2 nd order throughout this section. Correspondingly, the inverse metric can be approximated as where the last line was just given for notational reasons, i. e. to show that the unperturbed metric may be used to rise and lower indices. Note that the sub-and superscripts (0), (1) and (2) denote the perturbation order and should not be understood as covariant or contravariant indices.
In order to obtain an EUP for the metric (7), a general construction of geodesic balls in the weak curvature regime has to be given. Furthermore, the corresponding eigenvalue problem has to be solved in these domains.
This will be discussed in the following subsections.

A. Geodesic balls
In order to solve the given quantum mechanics problem, the domain to which the corresponding wave function ψ has to be confined needs to be found. As it is 3-diffeomorphism invariant, the most natural generalisa-tion of Heisenberg's slit width is a geodesic ball. This object is defined by its boundary on which every point p ∈ ∂B ρ has a constant geodesic distance σ(p 0 , p) from its center p 0 . Physically, the point p 0 describes the expectation value of the position operator.
Fixing p 0 , the biscalar σ(p 0 , p) becomes a scalar field σ(p) (related to Synge's world function [17]) which measures the distance along the shortest geodesic connecting the points p 0 and p. Hence, it has to satisfy the differential equation [18] g ij ∂ i σ∂ j σ = 1 (9) with the boundary condition σ(p 0 ) = 0.
As it conveniently describes surface normals of the spheres limiting the geodesic balls, we will interpret this scalar as coordinate. These geodesically spherical boundaries can then be conveniently defined by the relation σ = ρ. Note that the shape of geodesic balls depends on the coordinates and the background geometry which are chosen and does not need to even closely resemble ordinary balls in flat space and Cartesian coordinates. An example of their change is given in figure 1. In this illustration the background consists of constant time slices of the Schwarzschild static patch described by Schwarzschild coordinates. In order to illuminate the distortion, three geodesic balls of equal geodesic radius ρ but different coordinate distances r 0 between their center (p 0 ) and the center of symmetry of the space (origin) are compared. The corresponding calculations were done numerically thereby not invoking the small ball approximation made below. It can be seen quite clearly that the distortion of the geodesic balls increases with increasing curvature i. e. decreasing coordinate distance from the center of symmetry.
The coordinates which are normal to σ can be understood as surface parameters of geodesic spheres. In this paper we denote these as χ and γ to construct the geodesic coordinate system σ i = (σ, χ, γ). Imposing orthogonality as motivated above, the "angular" coordinates have to satisfy From a geometric point of view these coordinates by definition describe the background in such a way that the geodesic balls are not distorted. This might be quite unnatural for arbitrary metrics but simplifies a perturbative treatment in quantum mechanics enormously. Linearising (9) and (10) as indicated above, we can write with the geodesic coordinates to n th order σ i n . Similarly, we perturbatively construct derivatives with respect to the employed coordinates as which works up to 2 nd order because As a starting point for the perturbative quantum mechanical treatment we have to express the given metric in terms of geodesic coordinates. Then (9) and (10) become simply while the remaining four components of the metric contain all the information about the underlying space.
If the background is flat (as will be assumed later), this can be interpreted as a quasi-spherical (i. e. spherical up to corrections of order ) coordinate system. Then the unperturbed coordinates σ i 0 = (σ 0 , χ 0 , γ 0 ) in (11) furnish a spherical coordinate system constructed around p 0 .
For the remainder of this section we assume the metric (7) to be given in terms of geodesic coordinates.

B. Weak field approximation and gravitational waves
Our geodesic coordinates (11) are reminiscent of the weak field approximation applied in the case of gravitational waves [19]. The coordinate systems chosen in both cases are similar. In fact, the coordinate describing the geodesic direction of propagation of the gravitational waves is analogous to the geodesic distance coordinate employed in this paper. The other two coordinates are normal to the geodesic motion i. e. they are extended along hypersurfaces of constant geodesic distance and can thus be considered analogous to our "angular coordinates".
However, the representation which is usually used to describe gravitational waves is given in the so-called transverse-traceless (TT) gauge which provides four conditions in the 4-dimensional Lorentzian system while in our 3-dimensional approach there are only three conditions. Therefore the metric is not traceless in the geodesic coordinates defined above -it is merely expressed in a transverse gauge.

C. Eigenvalue problem
As in the original approach by Schürmann [14], the eigenvalue problem to be solved reads where λ n and ψ n (n stands for nlm in 3D) describe the eigenvalues and -functions (there are countably infinitely many) and ∆ denotes the Laplace-Beltrami operator which has to be evaluated on the curved background on which it reads Subjecting the background manifold to the metric splitting (7), the Christoffel symbols become (up to 2 nd order) with the unperturbed symbol (0) Γ k ij and the higher order expressions Correspondingly, the Laplace-Beltrami operator can be split as with the unperturbed operator ∆ (0) (on a flat background this is just the Laplacian in spherical coordinates) and the higher order corrections If coordinate sets are employed which naturally describe the underlying metric, the Dirichlet boundary conditions most probably become -dependent as well (due to the -dependence of ∂B ρ ) which makes a perturbative treatment unnecessarily complicated. This problem can be circumvented employing geodesic coordinates (11). As described above, then dσ denotes the geodesic radial form with respect to the origin p 0 and therefore describes the normals of surfaces of constant σ while dχ and dγ are the corresponding tangents.
On a flat background the 0 th order approximation to the metric is just flat space in spherical coordinates as expected and the boundary condition (16) has become homogeneous thus simplifying the calculation enormously. Now the problem can be separated order by order expanding the eigenvalues and -functions as There is a complication though: Perturbing the metric goes hand in hand with perturbing its determinant which is part of the measure dµ with respect to which we define the scalar product of the Hilbert space of our problem H.
In fact, dµ = √ gd 3 x but the square root of the determinant of the metric reads This translates to a perturbative measure dµ = dµ (0) + dµ (1) + O( 2 ) which in turn leads to a perturbation of the Hilbert space scalar product , = , 0 + , 1 +O( 2 ). The unperturbed operator ∆ (0) is self-adjoint with respect to dµ 0 i. e. it is only almost (meaning up to higher order corrections) self-adjoint with respect to dµ. Thus, the solutions of the unperturbed eigenvalue problem only furnish an almost orthonormal base of H. How to deal with this unusual kind of perturbation is explained in appendix A. There this method is applied to an abstract operatorÔ which in this case corresponds toÔ = −∆.
In the main text we will content ourselves with giving the perturbative corrections to the eigenvalues. To 1 st order we obtain which resembles the usual result. Note that the expectation value is taken with respect to the unperturbed measure. There is a significant change, though, to 2 nd order. The corresponding correction reads Then, restoring all quantum numbers nlm the momentum standard deviation of the state ψ nlm becomes to 2 nd order in In principle this approach can be used to get perturbative corrections to any already known solution. Nevertheless, we will concentrate on flat backgrounds.

D. Flat background perturbations
In this paper we are particularly interested in perturbing around a flat background i. e. flat space in spherical coordinates constructed around p 0 . Then, the 0 th order problem has general solutions characterized by the integers n, l, m : with the spherical harmonics Y l m , the spherical Bessel function of 1 st kind j l (x) and the n th zero of the spherical Bessel function of 1 st kind j l,n .
It is well-known that the lowest eigenvalue of the unperturbed operator equals λ 100 = π 2 /ρ 2 . Thus, the uncertainty acquires the lower limit The state in question is then described by the wave function which is independent of χ and γ. Using this and the characteristic of the geodesic coordinates (14), the relevant terms of the Laplace-Beltrami operator (i. e. the ones corresponding to derivatives with respect to σ denoted as ∆ (n) σ ) become order by order Thus, the problem is much more tractable than it appears at first glance. This is how far we can get without additional assumptions about the problem. These will be made in the following section.

IV. GENERAL SOLUTION FOR SMALL POSITION UNCERTAINTIES
The only assumption that this calculation is based on demands that the geodesic balls be sufficiently small i. e. that we stay in a coordinate neighbourhood of p 0 which is small enough such that we can approximate any metric in terms of Riemann normal coordinates x a as [20] g ij e a i e b j δ ab − with the dreibein e a i (i, j, . . . denote spherical coordinates σ 0 , χ 0 and γ 0 constructed around the point p 0 ) and the Riemann tensor R acbd . Thus, to 0 th approximation space is flat and there is no first order correction.
Correspondingly, the inverse metric reads where the flat space metric δ ab is used to rise and lower indices.
As the 1 st order correction to the metric vanishes, the 2 nd order can be treated similarly to a 1 st order correction.
Firstly, it will be shown that Riemann normal coordinates are "geodesic Cartesian coordinates" in the sense that spherical coordinates σ i 0 constructed around the center of the balls p 0 already equal geodesic coordinates σ i (c. f. (11)). This simplifies the problem enormously thus opening up the possibility to solve the general problem analytically.
In fact, it suffices to show that (14) is satisfied. In general, we have Note that we can assign x a = σ 0 l a with the unit radial vector l a = (sin χ 0 cos γ 0 , sin χ 0 sin γ 0 , cos χ 0 ) and similarly, e a σ0 = l a . Plugging this in, we obtain which by the symmetries of the Riemann tensor and the orthonormality of spherical coordinates becomes In particular, the same applies for higher order corrections showing that Riemann normal coordinates are in fact geodesic Cartesian coordinates. This makes an additional change of coordinates unnecessary and we can simply to continue with the formalism developed above. The radial part of the 2 nd order Laplace-Beltramioperator reads in this case where the 2 nd equality is valid because the metric at p 0 is in fact the identity. This translates into the 2 nd order correction to the eigenvalue with the area element of the geodesic two-sphere dΩ. The two integrals can be evaluated independently to yield Hence, we obtain the 2 nd order correction which leads to the uncertainty relation A higher order treatment of the problem is possible, yet tedious. An account to 4 th order is given in appendix B. In short, the 3 rd order correction vanishes and plugging in (B33) the final formula for the uncertainty up to 4th order reads with the 0 th order Cartan invariant Ψ 2 which satisfies in three dimensions [21] and the numerical constant Hence, the uncertainty relation is also altered on manifolds with vanishing Ricci scalar such as hypersurfaces of constant time in the static Schwarzschild patch (c. f. Section V A). Yet, the influence is significantly weakened.

V. NON-UNIFORM APPLICATIONS
In the preceding sections we introduced a formalism that makes an analytical treatment of the EUP possible without assuming further symmetries. Now the result will be applied to several examples. We will examine spherically symmetric metrics (constant curvature, spatial sections of Schwarzschild, Schwarzschild-(anti-)de Sitter, "hairy" and extremal Reissner-Nordström black holes) and spacelike hypersurfaces of the Gödel universe.

A. Spherical symmetry
Written in terms of the coordinates r i = (r, θ, φ) in Lemaître-Tolman-Bondi form [22], every spherically symmetric line element can be expressed as The form of the function E(r) indicates that we deal with spatially dependent curvature. This ansatz was used in many cosmological attempts to model dark energy due to a spherically symmetric void of matter in the universe [23] leading to a dipolar distribution of objects in the sky [24]. We denote the distance of the geodesic ball from the center of symmetry by r 0 . Due to the symmetry of the background space-time this completely characterizes the position of the center of the ball p 0 because the curvature function and its derivatives satisfy E (n) 0 ≡ E (n) | p0 = E (n) (r 0 ) where the subscript in this case denotes the n th derivative with respect to r.
First, say R| p0 = 0. Then, the leading order contribution reads Plugging in e. g. E(r) = (−1/2)Kr 2 , that is assuming constant curvature K, this equation recovers the result that was already obtained non-perturbatively in Ref. [14] and stated in (6). Additionally, the 4 th order contribution vanishes thereby validating the formalism -a perturbative expansion of a polynomial trivially equals itself confirming the result cited above. Secondly and conversely, the relation R| p0 = 0 is only satisfied by constant-time sections of the static Schwarzschild patch. The corresponding result describes the uncertainty relation in the presence of massive bodies (e. g. on the surface of planets or outside the horizon of black holes). It becomes to 4 th order Combining the results, the uncertainty relation for the spatial sections of the static patch of the Schwarzschild-(anti-)de Sitter-space-time can be obtained. The corresponding curvature function reads where K = Λ/3 with the cosmological constant Λ which can be expressed in terms of curvature scales as Λ/3 = R 2 H in dS (where it is interpreted as the Hubble radius) and Λ/3 = −R 2 H in AdS. Correspondingly, the uncertainty relation can be expressed as This equation also indicates the expected position at which both effects are equally strong. This happens when with R S = 2M. If R H R S , in the de Sitter-case these quantities still describe the positions of the horizons accurately. Thus, the relative strength of the effects depends on the position uncertainty.
Practically, the constant-time slices of our universe are not or at least just very slightly curved. Say the radius of curvature is of the order of the Hubble radius of our universe and an object is situated in the vicinity of the earth for which R S ∼ 10mm. Additionally assume (as for Hawking radiation) that ρ ∼ R S . Then the earth's influence is stronger than the constant curvature influence up to a distance of r 0 ∼ 10 4 km from the earth's center i. e. in the atmosphere. Hence, for laboratory experiments on the surface of the earth both effects need to be taken into account.
For Sagittarius A*, the black hole in the center of the Milky Way, this distance has already grown to r 0 ∼ 10 4 AU ∼ 1ly which extends far out to the orbits of the closest stars.
This could lead to interesting results in neutron-star and near black-hole physics. The Chandrasekhar mass would be one of such applications (see [25]).
Another possible extension are the spatial sections of a simple model of a hairy black hole based on the energy function E(r) = −M/r + Q 2ω /r 2ω [26] with the constant "hairy charge" Q and the particle number per unit area ansatz parameter ω. This geometry is a black hole solution of Einstein's equations in 4-dimensions with an anisotropic fluid stress energy tensor. For ω = 1 it describes the Reissner-Norström metric, while it develops short hair for ω > 1.
From (59) and (60) we obtain where the O(ρ 4 ) contribution for the hairy charge part was omitted for simplicity. This can be done for Q M and ω ≥ 1.
On the other hand, the result for the extremal Reissner-Nordström black hole (Q = M and ω = 1) can be given in terms of the 2 nd order contribution as Thus, in contrast to the Schwarzschild black hole the contribution is positive, thereby increasing the uncertainty. This aligns well with the fact that adding charge adds negative curvature which in turn reflects the relative sign between the charge and mass terms in the energy function.

B. Gödel universe
Describing the Gödel universe on spacelike hypersurfaces of the slicing given in [27,28], the 3-metric reads with the rotation radius a which is related to the angular velocity as Ω G = 1/( √ 2a). Due to the symmetry of the problem the position of the center of the ball is characterized by the polar distance from the rotation axis r 0 .
The corresponding leading order uncertainty relation reads This example concludes this section on non-uniform applications of the asymptotic EUP. However, we have not yet shown how to reconcile this relation with the well-known GUP from the quantum gravity side. This will be done in the following section.

VI. ASYMPTOTIC GENERALIZED EXTENDED UNCERTAINTY PRINCIPLE -AGEUP
In this section we will shortly comment on how to incorporate quantum gravity based effects into the problem. The (linear and quadratic) generalised uncertainty principle is often implemented in 3-D quantum mechanics by virtue of a deformed commutator (c. f. [5,29,30]) where the parametersα and α are dimensionful. More general deformed commutators which involve the coordinates are also possible [31], but we will not consider them here. We will only take into account the quadratic GUP here thus setting α = 0 and denote suggestivelyα = αl 2 p / 2 with the dimensionless GUP-parameter α and the Planck length l p .
Thus, the deformed commutator becomes With this choice of parameters noncommutativity of space coordinates appears only at order l 4 pp 4 / 4 and can thus be neglected.
We can write the momentum operator in terms of an auxiliary quantityp i which satisfies the usual commutation relations [x i ,p j ] = i i. e. it represents the normal quantum mechanical momentum operator. Following this ansatz, the thus called deformed momentum op-eratorP i readsP All the results obtained in this paper were written in terms of eigenvalues ofp 2 . Asp i andP i commute, we can simply translate the results usingP i terms of the results with respect top i . This works as follows: As the wave-functions solving the eigenvalue problem (3) furnish an orthonormal basis of the Hilbert space in question, we can write .
In every other case it leads to a GEUP in the given space once the EUP is known.

VII. SUMMARY
We have presented a formalism which allows to derive the Extended Uncertainty Principle (EUP) asymptotically on 3-dimensional Riemannian manifolds of arbitrary curvature as given by our formula (54). The result represents the uncertainty an observer (defining the foliation of spacelike 3-hypersurfaces) measures on a curved background. The formalism has first been applied to already derived homogeneous and isotropic space models, to be subsequently broadened in scope by inclusion of models which allow for other curvature related effects such as inhomogeneities, rotation and an anisotropic fluid leading to black hole "hair". The corresponding results influence neutron star and near black hole physics as well as cosmology, opening up an avenue for further investigation by observations. Interestingly, the leading (2 nd order) curvature induced correction is proportional to the Ricci scalar, while the 4 th order correction is proportional to the 0 th order Cartan invariant Ψ 2 and the curved space Laplacian of the Ricci scalar all evaluated at the expectation value of the position operator. We have shown that in the case of vanishing Ricci curvature (e. g. for the Schwarzschild metric), the lowest order perturbation is of 4 th order, while for the rotating Gödel metric as well as for anisotropic stress "hairy" black holes it is already of 2 nd order, though higher-order terms contribute as well.
Additionally, we have compared the influence which is exerted upon the uncertainty relation by cosmological horizons to the ones caused by local massive objects (planets, stars, black holes and so on) and proved the necessity to take into account the latter in many sensible applications on small (Earth) and large (Sagittarius A*) scales.
Finally, the formalism has been extended phenomenologically to combine the result with the Generalized Uncertainty Principle (GUP) by virtue of deformed commutators. Thus, we have presented an asymptotic form of the Generalized Extended Uncertainty Principle in the low energy limit given by the formula (86) which is our main achievement.
Appendix A: Measure changing non-singular perturbation theory to 2 nd order Consider an operatorÔ with domain D acting on the Hilbert space H = L 2 (D ⊆ IR 3 ,dµ) which is self-adjoint with respect to the measure dµ. Then it satisfies the eigenvalue equationÔ with its eigenvalues λ n and eigenvectors ψ n (n may stand for multiple quantum numbers) which are assumed to be normalized such that Further assume that the operator can be expanded perturbatively asÔ 1 and where all operatorsÔ (N ) areÔ (0)bounded. The latter is then called the unperturbed operator whose eigenvalues and -vectors are known. This is the starting point for perturbation theory in quantum mechanics.
A perturbative treatment of the eigenvalue problem (A1) can be complicated by the fact that the unperturbed operator might not be self-adjoint with respect to the measure µ. In fact, it may instead be self-adjoint with respect to a measure dµ 0 . In usual introductions to perturbation theory dµ = dµ 0 is assumed.
Yet, this is not generally the case. In principle, the measures may be related as with a positive non-vanishing function of the coordinates a 2 (x).
A simplifying assumption that has to be made here requires that the function a 2 (x) may be expanded perturbatively as It is questionable whether this problem is solvable or even mathematically well-defined without this additional assumption. Anyway, it certainly applies to the physical example studied in section III C. Expanding a 2 (x) as in (A6) is equivalent to introducing a perturbed scalar product as Additionally, we assume the spectrum ofÔ (0) to be discrete i. e. we consider non-singular perturbation theory. Of course, this could be generalised to the singular case.
By virtue of (A6) the unperturbed operator is almost self-adjoint meaning that the non-self-adjointness is of higher order i. e.
The first of these relations (A12) translates to which yields without loss of generality (discarding a global phase) At this point we see the first influence of the non-selfadjointness ofÔ (0) : The first order eigenvectors cannot be taken to be orthogonal to the unperturbed ones because they themselves are only orthogonal to 0 th order i. e.
Perturbing (A1) yields order by order These equations will now be contracted with ψ m | for some m. The case m = n is yields the corrections to the eigenvalues while the coefficients c (1) nm can be determined from the opposite case Note that said contraction shifts the order of perturbation of some terms not only due to the expansion of ψ m but also of the scalar product.
The 0 th order part then really corresponds to solving the unperturbed eigenvalue problem because takes place entirely within H 0 .
Furthermore, the 1 st order corrections to the eigenval-ues read (applying (A17)) which with (A12) becomes Hence,Ô (0) not being self-adjoint does not significantly change the 1 st order contribution to the eigenvalues. However, the expectation value is taken with respect to the measure dµ 0 .
Contracting the eigenvalue problem with ψ m | assuming n = m yields (after some algebra and applying (A18)) Again, this is equivalent to the case whereÔ (0) is selfadjoint and again the amplitude is calculated with respect to the unperturbed measure. Thus, to 1 st order the only clearly visible change is in the coefficient c (1) nn according to (A15). Analogously, the 2 nd order correction to the eigenvalues reads after application of (A13) and (A17) Finally, expanding ψ (1) n as in (A11) and applying (A15), (A16) and (A24), this can be expressed in terms of already known quantities as The first two terms of this equation appear in usual perturbation theory as well (which is indicated by the fact that they are evaluated with respect to the unperturbed measure) while the last two are entirely new a fact that can be inferred from the higher order scalar products they are featuring. Thus, the change in the measure has a strong effect on the 2 nd order eigenvalues and cannot be neglected.
Appendix B: Higher order corrections to the asymptotic uncertainty relation In this appendix we expand upon the perturbative derivation of the asymptotic uncertainty relation in the main text to get to 4 th order thereby explaining the origin of the mathematical constant ξ introduced in section IV. The required math is in principle the same as in the text although the calculations get rather lengthy and can in their entirety only be performed numerically. The Riemann normal coordinate expansion of the metric reads to 4 th order [20] g ij e a i e b j δ ab − Furthermore, the root of the determinant of the metric equals to 2 nd order which induces a perturbation in the measure with dµ 0 = σ 2 dσdΩ. As in section IV, the symmetries of the Riemann tensor imply that spherical coordinates constructed around p 0 are automatically geodesic coordinates. In other words (B1) already satisfies (14).
Accordingly, the higher order contributions to the ra-dial part of the Laplace-Beltrami-operator read where we defined for simplicity As there is no 1 st order contribution to the metric. the 3 rd order contribution can be treated as if it was of 1 st order while the 4 th order contribution can be treated as if it was of 2 nd order in the sense of appendix A. Hence, the 3 rd order contribution reads where the last equality stems from the fact that the angular integral vanishes. Thus, there is no 3 rd order contribution. However, the 4 th order eigenvalue is not as trivial. According to (A26), it reads (B10) In order to be able to compute the values of the terms which are not summed over, we need radial integral Furthermore, in three dimensions the Riemann tensor can be expressed in terms of the Ricci tensor and the metric (note that at p 0 we have g ab = δ ab ) as R abcd | p0 =δ ac R bd | p0 − δ ad R bc | p0 − δ bc R ad | p0 + δ bd R ac | p0 − R| p0 2 (δ ac δ bd − δ ad δ bc ) (B13) and by the contracted Bianchi identity ∇ c ∇ d R ab | p0 δ ab δ cd + δ ac δ bd + δ ad δ bc = 3∆R| p0 .
Combining these, we obtain (B17) Using this notation, we can express the 4 th order contri-bution to the eigenvalue as which is a scalar as expected.