arXiv:2107.01612v1 [hep-th] 4 Jul 2021 oaihial iegn nerl nlo igasaerpae ihco with replaced are t diagrams with loop increases in eigenvalues integrals eigenvalu momentum divergent momentum logarithmically the discrete the between over separation summation the with replaced be will eoi h hscllmtt eosrt hscancellation. this demonstrate to out. final limit cancel physical a cross-section the in total in exist the zero may in divergences that infrared photons the (low-energy) soft for accounting eainaaoost htfrteaglrmmnu components, momentum the angular that the Such for that [8]. to commute analogous not coupled straig relation is do most torsion components where The four-momentum theory, v Einstein–Cartan the finite [9]. the is torsion moment a torsion angular requires with to total mechanics the leads quantum for and law relativistic conservation electrodynamics the of quantum consistency in integrals divergent sdvret u ti o measurable. not char is redefinition mass, Such it (bare) but experiment. original divergent, the kno in the is measured a redefining physically to is by according that absorbed value infinity regular to be dimensional tends can Hooft–Veltman loop form metho four-momentum ’t common a the most in two and integral The [4] finite. regularization become Villars to integrals the modifying mt eodpoo,adte ebob h rt[] thigher- At ver [3]. a first is the diagrams. third fundamental reabsorbs The three then these photon. and from photon, virtual va second a is a reabsorbs anoth first emits and to The amplitud emits annihilates which diagrams. electron probability pair of an electron-positron kinds the virtual three a exist creates of there level, Feynma calculation one-loop by At represented a corrections perturbation electrodynamics, higher-order quantum In Introduction. 1. hs ag-nain at r oaihial iegn 1 2]. [1, divergent logarithmically are parts over gauge-invariant integrated whose are propagators These propagators. Feynman † ∗ [email protected] [email protected] nnncmuaiemmnu pc,gnrtdb oso,inte torsion, by generated space, momentum non-commutative In provides [7], torsion spacetime on based regularization, Torsional h nerl ersnigtesl-nryadtevre loha also vertex the and self-energy the representing integrals The ntefu-oetmsae olwn h rnil fsuperpositio of principle the following space, four-momentum the In . hsnnhsclrsl,cle h lrvoe iegne streate is divergence, ultraviolet the called result, nonphysical This h bevd enraie masses. re-normalized observed, the nrrddvrec.I h n,w aclt h nt bare finite n the small calculate a we photon 0 end, a re the lepton: gives torsional In torsion apply divergence. that thi show infrared we that also show article, We We this divergence. electrodynamics. quantum In in lepton charged Feynm sums. in convergent integrals into divergent magn w the regularizes in gravity, with replacement increases of this eigenvalues theory eigenva the Einstein–Cartan momentum between the separation the In over momentum. summation theory, the field quantum in osoa euaiaino efEeg n aeMs fEl of Mass Bare and Self-Energy of Regularization Torsional ntepeec fsaeietrin h oetmcomponen momentum the torsion, spacetime of presence the In nvriyo e ae,30Bso otRa,Ws ae,CT Haven, West Road, Post Boston 300 Haven, New of University . 39MV 90 MeV, 4329 1 eateto hmclEgneig alaeaCleeof College Tagliatela Engineering, Chemical of Department 2 eateto ahmtc n hsc,Cleeo rsadS and Arts of College Physics, and Mathematics of Department ihe e Grosso Del Michael . 5MV n 53MV epciey hs ausconstitute values These respectively. MeV, 1543 and MeV, 95 1 ∗ n ioe Pop lawski Nikodem and none nryadmmnu,rsligi integrals in resulting momentum, and energy unbounded igaswt lsdloso ita atce 1 2]. [1, particles virtual of loops closed with diagrams n ml asi ie otepoo n hntknto taken then and photon the to given is mass small A mo reDrcpril ncre pctm with spacetime curved in particle Dirac free a of um rpoo.Tescn ssl-nry hc swhen is which self-energy, is second The photon. er ndarm ihlosb unn them turning by loops with diagrams an e rwv ucin edn oafiie(dressed) finite a to leading function, wave or ge, rto vrtemmnu nFymndiagrams Feynman in momentum the over gration tt 1 ] hnadn otbaigradiation, braking soft adding When 2]. [1, state re ees l ena igasaecomposed are diagrams Feynman all levels, order oriaefaecnb hsn[11]. chosen be can frame coordinate a ultraviolet- eliminates that mechanism physical a osi 1] ntepeec ftetrintensor, torsion the of presence the In [10]. spin to ea nrrddvrec,wihaie rmnot from arises which divergence, infrared an ve td ftemmnu.Consequently, momentum. the of itude emgiueo h oetm Consequently, momentum. the of magnitude he asso h lcrn un n tau and muon, electron, the of masses umplrzto cag cenn) photon a screening): (charge polarization cuum yrglrzto:amteaia ehdof method mathematical a regularization: by d scle enraiain[,6.Abr quantity bare A 6]. [1, re-normalization called is e,wihi hna lcrneisaphoton, a emits electron an when is which tex, oetmcmoet aif commutation a satisfy components momentum naypoi omi h hscllmt This limit. physical the in form asymptotic wn rcdr lmntsteultraviolet the eliminates procedure s towr n otntrlter fgravity of theory natural most and htforward uaiaint h efeeg fa of self-energy the to gularization ,teecretosivleteitgainof integration the involve corrections these n, st lmnt hsdvrec r h Pauli– the are divergence this eliminate to ds ihtrini ope osi,the spin, to coupled is torsion hich s olwn h isenCra equations, Einstein–Cartan the Following es. neoms,wihrglrzsthe regularizes which mass, onzero vretsm.R-omlzto eoe a becomes Re-normalization sums. nvergent zto 5.We euao aihs the vanishes, regulator a When [5]. ization uswl elc nerto over integration replace will lues leo h aeeeti hre[] The [8]. charge electric bare the of alue sd o omt;therefore, commute; not do ts o hsclpoeswl involve will process physical a for e niern and Engineering 2 † 61,USA 06516, ciences, bu 5 of 85% about ectron 2

finite procedure and gives finite forms of the gauge-invariant vacuum polarization tensor, the running coupling, and the bare electric charge of the electron (and other charged particles) [8]. In this article, we apply torsion regularization to the self-energy of a charged lepton. We review this procedure [8] in Section 2. Next, in Section 3, we review the self-energy of the electron in quantum electrodynamics [2]. We show how torsion eliminates the ultraviolet divergence of the self-energy in Section 4, which constitutes the main part of this work. Then, in Section 5, we show how torsion gives the photon a small mass and thus eliminates an infrared divergence. Lastly, we determine the bare masses of the charged leptons in Section 6. Torsion may therefore physically eliminate the infinity problem in quantum field theory. This elimination would be a remarkable feature of the torsion tensor because torsion may also solve the singularity problem in general relativity [12]. Also, torsion may provide fermions with effective spatial extension [13] and explain the observed dynamics of the early Universe [14].

2. Torsion and Non-commutative Momentum.

This section reviews torsional regularization, following [8]. In the presence of torsion, the components pi of the four-dimensional momentum operator satisfy a commutation relation ~ k [pi,pj ]=2i S ij pk, k where S ij is the torsion tensor [8]. For the Dirac fields, the spin tensor and thus the torsion tensor are completely antisymmetric [10]. In a coordinate frame, in which the torsion pseudovector has only the temporal component S0 = Q/(2~), this commutation relation becomes − [px,py]= iQpz, [py,pz]= iQpx, [pz,px]= iQpy. We use c = 1. According to the Cartan equations, Q depends on the four-momentum and is proportional to p3, where p2 = (p0)2 p2, where p0 is the energy and p is the momentum: − Q = Up3. The constant U, which can be taken positive without loss of generality, is on the order of the squared inverse of the Planck mass MP. The cyclic commutation relations for the momentum can be written as

[nx,ny]= inz, [ny,nz]= inx, [nz,nx]= iny (1) for the spatial components ni of a vector p n = . (2) Q The relations (1) are analogous to those for the angular momentum. Consequently, the eigenvalues of n are given by n = n = j(j + 1) and nz = m, where an integer j 0 is the orbital quantum number and an integer m [ j, j] is the magnetic| | quantum number [11]. ≥ ∈ − Torsional regularizationp [8] replaces integration over the momentum in Feynman diagrams with summation over the momentum eigenvalues related to the eigenvalues of n: j ∞ ∞ ∞ dn dn dn f(n) 4π f(n) n =4π f(n) m =4π f(n) j(j +1)=4π f(n) n2, (3) x y z → | z| | | eigenstates j=1 m= j j=1 j=1 ZZZ X X X− X X where f(n) is an arbitrary scalar function of n. According to this prescription, such a function is multiplied by the absolute value of the commutator of two integration variables, [n ,n ] = n , and integration over continuous | x y | | z| variables nx, ny, nz is replaced with summation over the eigenvalues of n. For j = 0, m = 0 and thus nz = 0, which does not contribute to the sum in (3).1

1 The prescription proposed in [8] is based on the correspondence between the classical and quantum partition functions: 1 dq dp ↔ , dq dp f(H(q,p)) ↔ 2π f(E) |[q,p]|, π~ ZZ ZZ 2 eigenstatesX eigenstatesX

where f(x) = e−x/kT . The quantum commutation relation between the integration variables q and p leads to a discrete spectrum of energy eigenstates. Therefore, integration over a continuous phase space in a classical partition function is replaced with summation over energy eigenstates in the corresponding quantum partition function. 3

A Feynman diagram with a loop involves calculating an integral in the four-momentum space that has a form d4l/(l2 ∆+ iǫ)s, where ∆ > 0 does not depend on the four-momentum li, s is an integer, and ǫ 0+ [1, 2]. − 0 → 0 Applying the Wick rotation, in which the temporal component of the four-momentum l is replaced with ilE and R 2 0 2 2 thus lE = (lE) + lE, turns the integration in spacetime with the Lorentzian metric signature to the integration in the s 4 2 s s 0 2 s four-dimensional Euclidean space. The integral becomes i( 1) d lE/(lE + ∆) = i( 1) dlEdl/(lE + ∆) , where −3 − dl = dlxdlydlz. Omitting the subscript E and using Q = Ul with U > 0 and (2) gives R R l2 = (l0)2 + l2 = (l0)2 + U 2n2l6. (4)

The integration over l0 can be replaced with the integration over l:

dl0 1 3U 2n2l4 dl0 = dl = dl − . dl (1 U 2n2l4)1/2 − The integration over l can be replaced with the integration over n = l/Q and then changed to the summation over j in the non-commutative momentum space:

∞ dl dl dl f(l2) dn dn dn Jf(Q2n2) 4π Jf(Q2n2) n2, (5) x y z → x y z → j=1 ZZZ ZZZ X where J = ∂(lx,ly,lz)/∂(nx,ny,nz) is the Jacobian of the transformation from l to n. Differentiating (4) with respect 2 2 2 2 2 2 5 2 6 to nx, using n = nx + ny + nz, gives 2l(∂l/∂nx)=6U n l (∂l/∂nx)+2U l nx, which is equivalent to

∂l U 2l5n = x . ∂n 1 3U 2n2l4 x − The transformation derivatives are thus ∂l ∂(Qn ) ∂l Q x = x = Q +3Un l2 = [1 3U 2l4(n2 + n2)], ∂n ∂n x ∂n 1 3U 2n2l4 − y z x x x − ∂l ∂(Qn ) ∂l Q x = x =3Un l2 = (3U 2l4n n ), ∂n ∂n x ∂n 1 3U 2n2l4 x y y y y − and similarly for other components. They give the Jacobian:

∂lx/∂nx ∂lx/∂ny ∂lx/∂nz Q3 J = det ∂ly/∂nx ∂ly/∂ny ∂ly/∂nz = .   1 3U 2n2l4 ∂lz/∂nx ∂lz/∂ny ∂lz/∂nz −   Consequently, integrating over the Euclidean four-momentum li and using the prescription (5) gives

dl0 Q3 √1/(Un) U 3l9 dl0dl = dl J dn = dl dn =2 dl dn dl (1 U 2n2l4)1/2 (1 U 2n2l4)1/2 Z Z Z − Z0 − 1/(Un) ∞ √ U 3l9 8π dl n2, → (1 U 2n2l4)1/2 j=1 0 X Z − with n = j(j + 1). The integral dl0dl/(l2 + ∆)s in the presence of the Einstein–Cartan torsion is therefore [8] p R 1/(Un) 1 1 ∞ √ U 3l9 1 ∞ U 3ξ9n2(Un) 5 dl0dl 8π dl n2 =8π dξ − (l2 + ∆)s → (1 U 2n2l4)1/2 (l2 + ∆)s (1 ξ4)1/2[ξ2/(Un) + ∆]s j=1 0 j=1 0 Z X Z − X Z − 1 1 π/2 ∞ 8πU 2ξ9n 3(Un)s ∞ 4πU s 2ζ4ns 3 ∞ U s 2 sin4 φns 3 = dξ − − = dζ − − =4π dφ − − (1 ξ4)1/2(ξ2 + U∆n)s (1 ζ2)1/2(ζ + U∆n)s (sin φ + U∆n)s j=1 0 j=1 0 j=1 0 X Z − X Z − X Z π/2 ∞ U s 2 sin4 φ (j(j + 1))(s 3)/2 =4π dφ − − , (6) [sin φ + U∆(j(j + 1))1/2]s j=1 0 X Z 4 using a series of substitutions: Unl2 = ξ2 = ζ = sin φ. For U < 0, U in the sum-integral (6) is replaced with U . 3 | | At large values of j, the sum-integral (6) behaves as j∞=1 j− for any s, so it converges, regularizing the original integral. ∼ In the limit of continuous momentum space, U 0,P the separation between adjacent values of j has no effect on the integral, so the summation over j can be replaced→ with the integration over y, where y = j(j +1). For s = 3, this limit is equal to the original, finite value of the integral. For s = 2, the integral in this limit diverges as ln(U), equivalently to the Pauli–Villars regularization. p ∼

3. Electron Self-energy.

This section reviews the self-energy of the electron in quantum electrodynamics, following [2]. The Feynman propagator for an electron with bare mass m0 and four-momentum pµ is given by i( p + m ) = 6 0 , G (p2 m2 + iǫ) − 0 where p = γµp and γµ are the Dirac matrices. We use ~ = 1. The propagator in the second-order in e is given by 6 µ i( p + m ) i( p + m ) Σ = 6 0 [ iΣ (p)] 6 0 , G G (p2 m2 + iǫ) − 2 (p2 m2 + iǫ) − 0 − 0 where Σ2 represents the electron self-energy at one-loop level: d4k i( k + m ) i iΣ (p) = ( ie)2 γµ 6 0 γ − . − 2 − (2π)4 k2 m2 + iǫ µ (p k)2 µ2 + iǫ Z − 0 − − Infrared divergence is regularized by adding a small photon mass µ, which will be determined in Section 5. This propagator describes an electron that emits a photon and reabsorbs it. Introducing a Feynman parameter x combines two denominators into one: 1 1 1 1 = dx k2 m2 + iǫ (p k)2 µ2 + iǫ [(k2 2xk p + xp2 xµ2 (1 x)m2 + iǫ)]2 − 0 − − Z0 − · − − − 0 Consequently, one obtains 1 d4l 2x p +4m iΣ (p)= e2 dx − 6 0 , 2 − (2π)4 [l2 ∆+ iǫ]2 Z0 Z − where ∆= x(1 x)p2 + xµ2 + (1 x)m2. − − − 0 Applying the Wick rotation gives an integral in the four-dimensional Euclidean space:

1 d4l 2x p +4m Σ (p)= e2 dx E − 6 0 . 2 (2π)4 [l2 + ∆]2 Z0 Z E The full propagator of an electron is equal to the sum of terms that involve powers of one-particle-irreducible diagrams: i + Σ + Σ Σ + = . G G G G G G ··· p m Σ( p) 6 − 0 − 6 The physical (observed) mass m is located at a pole of the above denominator that is a solution of

[ p m0 Σ( p)] p=m =0. 6 − − 6 |6 The full propagator has a simple pole shifted from m0 by the self-energy Σ( p). Near the pole, the denominator is 6 1 equal to ( p m)(1 dΣ/d p) p=m, which gives the re-normalization constant of the electron wave function: Z2− = 6 − − 6 |6 (1 dΣ/d p) p=m. Therefore, the physical mass is shifted from the bare mass in the second order in e by − 6 |6 δm = m m =Σ ( p = m) Σ ( p = m ). − 0 2 6 ≈ 2 6 0 5

Consequently, one obtains

1 d4l 2 x δm =2m e2 dx E − , (7) 0 (2π)4 [l2 + ∆]2 Z0 Z E where

∆=(1 x)2m2 + xµ2. (8) − 0

4. Torsional Regularization of Ultraviolet Divergence.

Now, we proceed to the main calculation of this work. In the presence of torsion, integration over the mo- mentum is replaced with summation over the momentum eigenvalues according to (6):

π/2 4 ∞ 4 (s 3)/2 d lE s 2 sin φ [l(l + 1)] − 4πU − dφ . 2 s s (lE + ∆) → 0 [sin φ + U∆ l(l + 1)] Z Xl=1 Z For s = 2, this replacement gives p

π/2 d4l ∞ sin4 φ [l(l + 1)] 1/2 E 4π dφ − . (9) 2 2 2 (lE + ∆) → 0 [sin φ + U∆ l(l + 1)] Z Xl=1 Z With the prescription (9), the mass shift (7) becomes p

2 π/2 1 4 1/2 8πm e ∞ sin φ(2 x)[l(l + 1)]− δm = 0 dφ dx − . (10) 4 2 (2π) 0 0 [sin φ + U∆ l(l + 1)] Xl=1 Z Z 2 2 p Since U MP− , U∆ is on the order of (m0/MP) 1. In this limit, the separation between adjacent values of l does not affect∼ significantly the integrand and thus the≪ summation over l in (10) can be replaced with the integration over n, where n = l(l + 1). Consequently, we obtain

p 2 π/2 1 4 1 8πm e ∞ sin φ(2 x)n− δm = 0 dφ dx dn − . (2π)4 [sin φ + U∆n]2 Z0 Z0 Z√2 Substituting from n to a new variable w = U∆n gives

2 π/2 1 8πm e ∞ dw δm = 0 dφ sin4 φ dx(2 x) (2π)4 − w[sin φ + w]2 Z0 Z0 Z√2U∆ 1 π/2 3 2m0α sin φ 2 sin φ = 2 dx(2 x) dφ + sin φ ln +1 , (11) π 0 − 0 −sin φ + √2U∆ √2U∆ Z Z h  i where α = e2/4π is the fine structure constant. This integral is convergent as long as U = 0, that is, in the presence of torsion, and as long as ∆ = 0. 6 6

5. Torsional Regularization of Infrared Divergence.

According to the Einstein–Cartan theory of gravity, a fermion field is a source of torsion. When a photon is coupled to a lepton, it is also minimally coupled to the torsion tensor:

F 2 (∂A + SA)2, L∼ ∼ where is the Lagrangian density, F symbolizes the electromagnetic field tensor, A symbolizes the electromagnetic potential,L and S symbolizes the torsion tensor. Therefore, the term S2A2 arising from this coupling generates a mass 6 term µ2A2, so the mass of a photon is on the order of magnitude of the torsion tensor, which is on the order of Up3. Using p = m0, as in (7), gives

3 1 µ m0U,U 2 . (12) ∼ ∼ MP For the entire range of x [0, 1], ∆ in (8) is thus different from zero. Consequently, torsion removes infrared divergence of a lepton self-energy. For∈ a Feynman diagram that converges in the limit µ 0, the small nonzero mass µ does not significantly affect its amplitude. →

6. Calculation and Discussion.

Since U∆ 1, the integral (11) can be approximated as ≪ 1 π/2 3 2m0α sin φ 2 sin φ δm = 2 dx(2 x) dφ + sin φ ln π 0 − 0 − sin φ √2U∆ Z Z h  i 2m α 1 π/2 ∆ 0 2 2 √ 2 2 = 2 dx(2 x) dφ sin φ sin φ ln 2Um0 2 + sin φ ln(sin φ) π 0 − 0 − − m0 Z Z h   i 2m α 1 π/2 ∆ π/2 0 2 √ 2 2 = 2 dx(2 x) dφ sin φ 1 ln( 2Um0) ln 2 + dφ sin φ ln(sin φ) , π 0 − " 0 − − − m0 0 # Z Z h  i Z where ∆ is given by (8). After integrating over φ, we obtain

2m α 1 π ∆ 0 √ 2 δm = 2 dx(2 x) ( 1 ln( 2Um0) ln 2 + N , π 0 − " 4 − − − m0 # Z h  i where π/2 N = sin2 φ ln(sin φ) dφ 0.151697. ≈− Z0 Consequently, the mass shift is equal to

2m α 3 π 2m α 1 π ∆ 0 √ 2 0 δm = 2 1 ln( 2Um0) + N 2 dx(2 x) ln 2 π 2 4 − − − π 0 − 4 m0 h i Z   3m α 4N
m α 1 µ2 0 √ 2 0 2 = 1 ln( 2Um0) dx(2 x) ln (1 x) + x 2 . (13) 4π π − − − 2π 0 − − m0   Z   With the substitution of U and µ (12), we use (13) to determine the observed mass m = m0 +δm as a function of the bare mass m0. Conversely, we obtain the bare mass from the observed mass. For the electron, m =0.510999 MeV gives m0 =0.432928 MeV. For the muon, m = 105.658 MeV gives m0 = 90.9514 MeV. For the tau lepton, m = 1776.86 MeV gives m0 = 1542.63 MeV. Our results depend on the orders of magnitude of the constants U and µ estimated in (12), but they are not significantly affected by changing their exact values. For the electron, increasing U by a factor of 10 gives m0 = 0.434408MeV and decreasing it by the same factor gives m0 =0.431457MeV. Increasing or decreasing µ by a factor of 10 gives m0 = 0.432928 MeV. Accordingly, the value of µ is less significant than the value of U in the shift from the bare mass to the re-normalized mass. We conclude that the bare masses of charged leptons constitute about 85% of their observed, re-normalized masses. To compare, the bare electric charge is about 22% higher than the observed, re-normalized charge [8]. These numbers might change slightly when the contributions from the W and Z bosons are included. This work was funded by the University Research Scholar program at the University of New Haven.

[1] R. P. Feynman, Phys. Rev. 76, 769 (1949); J. D. Bjorken and S. D. Drell, Relativistic Quantum Fields (McGraw-Hill, 1965); C. Nash, Relativistic Quantum Fields (Academic Press, 1978); F. Mandl and G. Shaw, Quantum Field Theory (Wiley, 7

1993); W. Greiner and J. Reinhardt, Quantum Electrodynamics (Springer, 1994); M. Maggiore, A Modern Introduction to Quantum Field Theory (Oxford University Press, 2005); C. Itzykson and J.-B. Zuber, Quantum Field Theory (Dover, 2006); K. Huang, Quantum Field Theory: From Operators to Path Integrals (Wiley, 2010). [2] M. E. Peskin and D. V. Schroeder, An Introduction to Quantum Field Theory (Addison-Wesley, 1995). [3] J. Schwinger, Phys. Rev. 75, 651 (1949); Phys. Rev. 82, 664 (1951). [4] W. Pauli and F. Villars, Rev. Mod. Phys. 21, 434 (1949). [5] G. ’t Hooft and M. Veltman, Nucl. Phys. B 44, 189 (1972). [6] F. J. Dyson, Phys. Rev. 75, 486 (1949); M. Gell-Mann and F. E. Low, Phys. Rev. 95, 1300 (1954). [7] E. Schr¨odinger, Space-time Structure (Cambridge University Press, 1954). [8] N. Poplawski, Found. Phys. 50, 900 (2020); arXiv:1807.07068. [9] F. W. Hehl, P. von der Heyde, G. D. Kerlick, and J. M. Nester, Rev. Mod. Phys. 48, 393 (1976). [10] T. W. B. Kibble, J. Math. Phys. 2, 212 (1961); D. W. Sciama, in Recent Developments in General Relativity, p. 415 (Pergamon, 1962); Rev. Mod. Phys. 36, 463 (1964); Rev. Mod. Phys. 36, 1103 (1964); F. W. Hehl and B. K. Datta, J. Math. Phys. 12, 1334 (1971); E. A. Lord, Tensors, Relativity and Cosmology (McGraw-Hill, 1976); V. de Sabbata and M. Gasperini, Introduction to Gravitation (World Scientific, 1985); K. Nomura, T. Shirafuji, and K. Hayashi, Prog. Theor. Phys. 86, 1239 (1991); V. de Sabbata and C. Sivaram, Spin and Torsion in Gravitation (World Scientific, 1994); N. Poplawski, Classical Physics: Spacetime and Fields, arXiv:0911.0334. [11] P. A. M. Dirac, The Principles of Quantum Mechanics (Oxford University Press, 1930); L. D. Landau and E. M. Lifshitz, Quantum Mechanics: Non-Relativistic Theory (Pergamon, 1977); J. J. Sakurai, Modern Quantum Mechanics (Addison- Wesley, 1994). [12] F. W. Hehl, P. von der Heyde, and G. D. Kerlick, Phys. Rev. D 10, 1066 (1974); B. Kuchowicz, Gen. Relativ. Gravit. 9, 511 (1978); M. Gasperini, Phys. Rev. Lett. 56, 2873 (1986); N. J. Poplawski, Phys. Lett. B 694, 181 (2010); Phys. Lett. B 701, 672 (2011); Gen. Relativ. Gravit. 44, 1007 (2012); N. Poplawski, Phys. Rev. D 85, 107502 (2012); G. Unger and N. Poplawski, Astrophys. J. 870, 78 (2019); J. L. Cubero and N. J. Poplawski, Class. Quantum Grav. 37, 025011 (2020); N. Poplawski, Zh. Eksp. Teor. Fiz. 159, 448 (2021); J. Exp. Theor. Phys. 132, 374 (2021). [13] N. J. Poplawski, Phys. Lett. B 690, 73 (2010); Phys. Lett. B 727, 575 (2013). [14] N. Poplawski, Astrophys. J. 832, 96 (2016); S. Desai and N. J. Poplawski, Phys. Lett. B 755, 183 (2016).