arXiv:1901.04971v1 [hep-ph] 15 Jan 2019 pcrma al si 98[ 1928 hydrogen in the as of early structure as fine field, spectrum Coulomb correct electrostatic the predicted an has in immersed electron an ec fDrcwv ucin ti fe er htsince that heard diver- of often integral near-the-origin is spatial It this the function? wave behind Dirac of physics gence the is What corresponding the origin, the at are hydrogen regular of always functions wave Schrödinger the Whereas as scales equation Dirac (rel- of triumph greatest of mechanics. one quantum as ativistic) viewed achievement, often historic is this which of importance the derestimate of solution Coulomb The physics. equa- equation, of Dirac celebrated time the top all the in among spin- tions ranked the often for framework electron, mechanical quantum ativistic Introduction. a eairo ia aefnto rnspronounced brings small function extremely singu- in wave only probabilistic this Dirac effect violate Moreover, of to behavior function. not lar wave appear of origin interpretation the at gence egn tteoii [ origin the at di- vergent (logarithmically) weakly become functions wave Dirac at9 er.Teeeitsm prdcwr htat- that work sporadic some exist There the years. in 90 addressed satisfactorily past been never has divergence origin. weak the this near of underlying behavior myth is divergent the it unravel angle, to theoretical curiosity the great practi- from any However, cause disaster. to cal not nucleus. seems the divergence of this radius Therefore, charge the than smaller magnitude erteoii iegneo ia aefntoso hydr of functions wave Dirac of divergence Near-the-origin j si ieykon h olm ouin fthe of solutions Coulomb the known, widely is As upiigy h uzerltdt h near-the-origin the to related puzzle the Surprisingly, + 2 1 2 − Z 2 nmlu cln eairehbtdb h ia aefunc wave Dirac the by exhibited behavior scaling anomalous stepyisgvrigtewal iegn near-the-ori divergent weakly the the of governing physics the is ru qain n a eu h edn oaihst l o all W to logarithms NREFT. leading the the resum in ca can The one term equation, fields. (zitterbewegung) group nucleus Darwin the and and electron correction the of (OPE) expansion the by accounted be olo hswr st eosrt ht ihntenonrelat di the logarithmic within universal that, the demonstrate Kle atoms, Coulombic to and for is Schrödinger tailored work this atomic of the goal of behavior near-the-origin P e pi h eaiitcQD eaeteeyefre to enforced thereby are We when meaningful QED. be relativistic to the cease must in function up set logarithm OPE universal these that counterintuitive somewhat α ia qain[ equation Dirac hr saln-tnigpzl ocrigteCuobsolu Coulomb the concerning puzzle long-standing a is There 2 nttt fHg nryPyis hns cdm fScienc of Academy Chinese Physics, Energy High of Institute [ 6 colo hsc,Uiest fCieeAaeyo Science of Academy Chinese of University Physics, of School , nS 7 e.g. .Acroscnnrmte arises. then conundrum curious A ]. 1 Ψ / 2 6 † eaiitcwv qainfor equation wave relativistic a , .Hwti atr a arise? can pattern this How ]. Ψ r yrgn sasqe forpeeigwr htamt demyst to aim that work preceding our of sequel a As hydrogen? ν − tl ean nt,ti diver- this finite, remains still 1 r 3 erteoii,with origin, the near 4 culymn orders-of- many actually , rvdsascesu rel- successful a provides ] , etraieycalculable perturbatively 5 .I sdffiutt un- to difficult is It ]. ighn Huang, Yingsheng Dtd aur 6 2019) 16, January (Dated: nS ν expansion 2 1 r isnceceteegn rmteoeao product operator the from emerging coefficient Wilson / = ∗ 2 ssotrta h lcrnsCmtnwavelength. Compton electron’s the than shorter is 1 2 uJia, Yu rae ehia details. technical greater [ (OPE) expansion the product by operator accounted hydrogen be calculable the can batively of functions wave divergences (NREFT), Dirac logarithmic theory field universal effective the frame- nonrelativistic the within of that, field-theoretical show work We modern physics. some particle in to tools be appealing can by equation wave resolved relativistic in puzzle long-standing real- not have nature. and perturbative nuisance, its a ized as divergence this regard odnwv ucin,fre rlg.Tecur- The Klein-Gordon the resembles trilogy. [ aspects case a many in formed work rent Klein- functions, and wave Schrödinger atomic Gordon the of co- electron-nucleus behavior universal alescence the for account to OPE and [ work recent together two function, our wave with Dirac of behavior short-distance Comp- meaningless. electron become the may than function wavelength shorter wave ton scale length Dirac a the at that probed is clarification conceptual h ia qaini olm edo xdnucleus fixed a of field Coulomb charge a with in for equation Dirac functions The wave Dirac of Divergence tvsi olm on tt omls fra incom- rel- [ an of see (for context list, formalism plete the state in bound myth Coulomb this ativistic understand to tempted h olo hsLte st eosrt htthis that demonstrate to is Letter this of goal The hsLte,wiham tdcpeigteuniversal the deciphering at aims which Letter, This † cdvrecscnntb cone ythe by accounted be not can divergences ic i eairo h ia aefunctions wave Dirac the of behavior gin egneo h ia aefntoscan functions wave Dirac the of vergence 2 n u Yu Rui and ,t hc erfrteitrse edr for readers interested the refer we which to ], infrthe for tion vsi ffciefil hoy(NREFT) theory field effective ivistic nGro aefntos[ functions wave in-Gordon dr in rders in fteDrcequation, Dirac the of tions − s sdet h eaiitckinetic relativistic the to due is use i Ze s ejn 009 hn and China 100049, Beijing es, ~ c 8 t h i frenormalization of aid the ith ,Biig104,China 100049, Beijing s, α ocueta h ia wave Dirac the that conclude – reads 12 Zα · nS ∇ isncecetalae ihthe with affiliated coefficient Wilson ‡ ) notntl l fte edto tend them of all Unfortunately ]). gnadoeao product operator and ogen n eoe the recover and 1 + / 2 βmc yrgn tappears It hydrogen. figteuniversal the ifying 1 , 2 2 ,wihas pl NREFT apply also which ], − Ze 4 πr 1 i.e. 2 r , − 13 2 Z what , = Ψ ,the ], 2 .A important An ]. α nS 2 / 2 1 E / 2 Ψ , hydrogen. pertur- (1) 2
2 2 α ∇ e where ,β are the Dirac matrices, m is the electron mass. Here H0 = 2m 4πr is the nonrelativistic Coulomb In the nonrelativistic limit (Zα 1), the energy spec- Hamiltonian− (for brevity,− we adopt the ~ = 1 unit from ≪ trum for the hydrogen-like atoms are approximately now on). The kinetic p4 correction, Darwin term (re- sponsible for the zitterbewegung motion) and the spin- 2 2 4 4 2 Z α Z α n 3 orbit terms comprise the (1/c2) relativistic corrections Enj = mc 1 + , − 2n2 − 2n4 j + 1 − 4 ··· O 2 to H0: (2) ∇4 1 Zα H = , H = 2 . 2 kin 3 2 Dar 2 2 where α e signifies the fine structure constant. −8m c −8m c ∇ r ≡ 4π~c Zα From now on we will concentrate on the Dirac wave Hs.o. = σ L. (8) 1 4m2c2r3 · functions for j = 2 , positive parity (nS1/2) hydrogen. The normalized wave functions read [7]: The fine structure of hydrogen level in (2) is readily recovered by incorporating (8) as small perturbations. 1 Here we explore the first-order perturbative correction Fn(r) 4π ξm Ψn 1 m(r)= , (3) to the Schrödinger wave function at the origin for the 2 G (r) q3 σ ˆr ξ n 4π · m 1S state, again treating (8) as the perturbation. We q wonder whether the correction is finite or not. Following 1 0 the standard first-order perturbation theory in quantum where the Pauli spinors ξ 1 = and ξ 1 = . 2 0 − 2 1 mechanics, we have Near the origin, Fn(r) roughly scales as [6] ∞ (1) dk Sch k0 ∆H 10 Z2α2 ∆R10 (0) Rk0 (0)h | | i. (9) − 2 ∼ ˆ0 2π E10 Ek0 Sch 2r − Fn(r) Rn (0) , (4) ≈ na0 To our purpose, we need only consider the contribution from the continuum states when summing over the in- ~ Sch where a0 = /(mcZα) is the Bohr radius, and Rn (0) termediate states. Since ∆H is rotation scalar, those represents the radial Schrödinger wave function at the surviving intermediates states must have l =0. origin for the nS hydrogen state. We have also taken the It turns out that the integration in (9) are dominated nonrelativistic approximation √1 Z2α2 1 Z2α2/2 by the the high-frequency intermediate states. To reg- in the exponent. The singularity− has noticeable≈ − effect 2 ularize the occurring UV divergence, we impose a hard only when r / na0 exp 2/Z2α2 na0 10−16300/Z [14], momentum cutoff mZα Λ m in (9). Accordingly, 2 − ∼ 2 which is even many orders shorter than the length scale each piece of ∆H renders≪ the following≪ divergent correc- related to the QED Landau pole! tion to the 1S wave function at the origin: A more careful expansion of the Dirac radial wave func- tion near the origin in first few powers of Zα reveals that (1) Sch ZαΛ 2 2 ∆R10 (0) = R10 (0) + Z α ln Λ , (10a) kin πm Z2α2 2 2 Sch (1) Sch ZαΛ Z α lim Fn(r)= Rn (0) (1 r/a0) 1 ln r + . ∆R (0) = R (0) ln Λ , (10b) r→0 − − 2 ··· 10 10 Dar − 2πm − 2 (5) and Hs.o. yields a null correction for S-wave hydrogen. Assuming Λ mZα, one then verifies that each of ∆H The wave function is manifestly logarithmically divergent ∼ 2 2 2 at the origin. does generate a correction of order v Z α , compati- ble with the 1/c2 scaling in (8). Summing∼ each entry in First-order quantum-mechanical perturbation theory. Let (10), we find us first try to build some insight via a useful exercise 2 2 in the nonrelativistic quantum mechanics. The electron (1) Sch ZαΛ Z α ∆R10 (0) = R10 (0) + ln Λ . (11) in the hydrogen-like atom is undoubtedly nonrelativistic, ∆H 2πm 2 with v/c Zα 1). As can be found in many texts in 2 2 ∼ ∼≪ Z α quantum mechanics [7, 15], it is convenient to conduct a Intriguingly, the 2 ln Λ term coincides with what ap- nonrelativistic reduction (v/c expansion) to Dirac equa- pears in the expanded Dirac wave function (5), provided tion, and end up with an effective Schrdinger¨ equation: Λ is replaced with 1/r. This might be a strong indica- tor that the kinetic and Darwin terms are responsible for Heff ψ(r)= Eψ(r), (6) generating the mild logarithmic divergence of the Dirac wave function near the origin. The linear power diver- where ψ is a Pauli spinor. The effective Hamiltonian is gence in (11) is clearly absent in the original Dirac wave made of function (5). In our opinion, these UV power divergences arise from our specific use of the sharp momentum cutoff Heff = H0 + ∆H = H0 + Hkin + HDar + Hs.o.. (7) as the UV regulator, which may be viewed as an artifact. 3
OPE in relativistic QED fails. To reformulate the rel- r 0 q r 0 ativistic Dirac equation in an external Coulomb poten- q tial in a field-theoretical fashion, one may naturally at- = − − UVCT k tempt to couple relativistic QED with a external elec- p k p p k trostatic field. Analogous to our previous work [1, 2], rather than introduce the classical electromagnetic field, we introduce the heavy nucleus effective theory (HNET) + + UVCT (very similar to the famous heavy quark effective the- p k ory [16, 17]), and assign a HNET field N to take into account the quantum fluctuation around a positively- charged and static nucleus. The relativistic UV theory Figure 1. The next-to-leading Wilson expansion for the prod- tailored to describe atoms is then obtained by combining uct of the Ψ and N fields. The original diagram is reexpressed as the sum of two terms, by subtracting and adding the renor- the QED and HNET: malized Green function inserted with a local ΨN operator. 1 = Ψ(¯ iD/ m)Ψ + N †iD N F F µν , (12) LUV − 0 − 4 µν reads where the covariant derivatives acting on electron and nucleus are Dµ = ∂µ + ieAµ and Dµ = ∂µ iZeAµ, re- − Zα spectively. (1) only concerns the MN limit, which ln r + , Feynman is insensitive to the nucleus’s detailed→ properties ∞ such as 2π ··· c(r)= (15) magnetic moment and charge radius. We thereby retain Zα ln r + , Coulomb only the leading HNET lagrangian in (12), and treat N π ··· as a bosonic field because the spin degree of freedom de- couples. in Feynman and Coulomb gauge, respectively. Following [1, 2], we tentatively express the Dirac wave function of a hydrogen atom in term of the following Notwithstanding its gauge dependence, the Wilson co- vacuum-to-atom matrix element: efficient c(r) in (15) does contain the ln r term, which un- fortunately arises prematurely at order Zα, rather than Ψ (r) 0 Ψ(r)N(0) njm; Π , (13) the desired Z2α2 ln r in (5). njmΠ ≡ h | | i where Π= 1 signifies the parity. Note both field oper- Therefore, we have to admit that the Wilson expan- ators are defined± in an equal time t =0. This spatially- sion formulated in our UV theory for atoms, (12), fails nonlocal matrix element clearly is not gauge invariant. to account for the near-the-origin behavior of the hydro- The universal near-origin behavior of the Dirac wave gen Dirac wave function. This failure may cast some function may be transparently understood from (13), doubt on the legitimacy of interpreting the matrix ele- once the Wilson expansion of the operator product [13] ment (13) as the Dirac wave function. The very short distance (r 1/m) behavior of the Dirac wave func- for φ(x)N(0) is conducted. In the r 0 limit, the an- ≪ ticipated OPE relation reads → tion cannot be correctly accounted by the OPE based on (12), which may hint that one should not associate phys- Ψ (r)N (0) = (1+ c(r)) [ΨN] (0) + , (14) ical significance with the Dirac wave function literally as R R R ··· r 0. → where r r , the notion “[...] ” implies the renor- ≡ | | R malized operator, and the ellipsis indicates the higher- NREFT lagrangian and the desired OPE. Inspired by the dimensional operators carrying more derivatives. Were success of quantum mechanical perturbation theory en- 2 2 the Wilson coefficient c(r) exhibit the desired Z α ln r coded in (11), one might like to reexamine the OPE for- structure as in (5), we would successfully accomplish our mulated in the nonrelativstic EFT rather than the UV goal. theory (12), but including a certain class of relativistic The renormalization of ΨN resembles very much with corrections. that of the heavy-light quark current q¯Γhv in HQET [18], where the logarithmic UV divergence is first encountered A natural option is to replace the QED lagrangian at order Zα. This immediately indicates the Wilson co- in (12) by the nonrelativistic QED (NRQED) [19]. To efficient should exhibit Zα ln r structure. Following the mimic the dynamics underpinning the Coulomb solution instruction expounded in [1], by analyzing the asymptotic of the Dirac equation, we deliberately drop all the occur- behavior of the four-point Green function in the r 0 rences of the A field. In the lowest order of 1/MN and limit as guided by Fig. 1, one finds that c(r) at order-→Zα through the order of v2/c2, the nonrelativistic EFT for 4 atoms [20] then reduces to relation turns out to read
lim ψ(r)N(0) = (r)[ψN]R(0) + , (18) ∇2 ∇4 1 † r→ m C ··· NREFT = ψ iD0 + + 3 L ( 2m 8m Z2α2 (r)=1 mZαr (ln µr + const) + (Z3α3). ∇2A0 σ ∇A0 ×∇ C − − 2 O c e + ic e · + ψ (16) D 2 S 2 This formula constitutes the key result of this work. The − 8m 4m ··· ) (Zα) Wilson coefficient was first derived in the nonrel- † c4 † † 1 2 O + N iD0N + ψ ψN N + A0 + . ativistic Schrödinger EFT [1], which is invoked to explain m2 2 ∇ ··· the universal electron-nucleus coalescence behavior of the The positron sector in NRQED has been neglected. At atomic Schrodinger wave function [21, 22]. 2 2 tree level, the short-distance coefficients associated with The (Z α ) coefficient in (18) is highly desirable, which reproducesO the logarithmic scaling of the Dirac the Darwin and spin-orbit terms are simply cD = cS =1. Note that an intrinsic UV cutoff Λ m must be imposed wave function at small r in (5). Moreover, it also de- to make (16) meaningful. ≤ pends on an artificial scale mZα < µ m, which is the renormalization scale associated with≤ the local operator Eq. (16) is no longer gauge invariant. However, this [ψN] . does not necessarily imply a nuisance. In fact, Coulomb R In [2], we have shown that the kinetic p4 term plays the gauge A =0 is a convenient choice for tackling non- vital role in generating the Z2α2 ln r coefficient for spin-0 relativistic∇ · bound state problem, where the effects of A electron. For spin- 1 electron, the Darwin term is as im- and A0 can be cleanly separated. Unlike the physics un- 2 portant as the kinetic term, both of which represent the derpinning the Lamb shift [20], Dirac equation (1) per se (1/c2) relativistic corrections in (16) (Spin-orbit term does not involve the effect of the physical photon excited yieldsO a vanishing contribution). The extra contribution by the A field. The A0 field is the very agent to mediate from Darwin term in the Dirac case explains why the the instantaneous Coulomb interaction, which provides Z2α2 ln r coefficient in (18) differs by the factor of two the crucial binding mechanism all atoms. We will stay from its counterpart in the Klein-Gordon case. exclusively with the Coulomb gauge. The OPE relation in the momentum space can be ob- One can descend from (12) to (16) via the perturba- tained by Fourier transforming (18). In the q m limit, tive matching procedure. There may potentially arise one has → a contact interaction in NREFT due to integration out + − of the shortly-lived e e intermediate states. Neverthe- ψ(q)N(0) d3r e−iq·rϕ(r)N(0) less, through an explicit one-loop matching calculation, ≡ ˆ we find no room for the contact interaction to contribute. e (q)[ψN]R(0) + , (19) Thus we conclude c = (Z3α3). → C ··· 4 O 8πmZα + (Z3α3) π2Z2α2 + (Z4α4) (q)= e O O , Due to the single-pole structure of the nonrelativistic C q4 − q 3 propagators and the instantaneous Coulomb interaction, | | self-energy diagrams for the ϕ and N fields, as well as wheree we have suppressed the free-theory Wilson coef- crossed ladder diagrams simply vanish in arbitrarily high ficient (2π)3δ(3)(q). In contrast to the renormalizable order. As a gratifying simplification, we only need con- field theory, here the coexistence of both q−4 and q −3 | | sider the uncrossed Coulomb ladder diagrams in NREFT scalings in the Wilson coefficient is peculiar in NREFT, sector. which is dubbed the double-layer form of OPE [2]. Analogous to (13), we tentatively define the large com- Renormalization of local Operator in NREFT. To better 2 2 ponent of the Dirac wave function for the nS1/2 as the understand why the Z α ln r arises in the Wilson coef- following nonlocal matrix element: ficient, let us first examine the renormalized composite operator ψN, which appears in the right hand side of 1 OPE in (18). Fn(r) ξm = 0 ψ(r)N(0) nS1/2,m , (17) √4π h | | i It is known that the logarithmic UV singularity arises first at two-loop order in NRQED (NRQCD) vertex di- with the Dirac field Ψ in (13) replaced by the Pauli spinor agrams, first discovered by matching the vector current field ψ. from QCD onto NRQCD [23, 24]. Later an explicit two- In [2], we have successfully accounted for the near-the- loop calculation is also conducted in the nonrelativistic origin behavior of the Klein-Gordon hydrogen wave func- EFT of QCD [25]. It is a certain class of relativistic ef- tion with the aid of OPE formulated in the NREFT. fects that play a vital role in generating logarithmic UV Hearteningly, the universal near-the-origin behavior of divergence. This knowledge hints that the local operator the Dirac wave function can also be understood by ap- ψN may well be subject to a multiplicative renormaliza- plying the OPE to (17). In the small r limit, the OPE tion at order Z2α2. 5
are automatically absent in this scheme. At order Z2α2, only two diagrams in Fig. 2 generate logarithmic UV di- vergence:
Z2α2 1 = + ln µ2 + finite , (20a) G2a) 2 ǫ (a) (b) (c) (d) Z2α2 1 = + ln µ2 + finite , (20b) G2c) − 4 ǫ The first term comes from Fig. 2a) which involves a p4 insertion in the bottom loop [1]. The second term is new, which arising from Fig. 2a) containing a Darwin vertex. The divergence of order-v0 is specific to NRQED, because 2 (e) (f) (g) (h) the Coulomb enhancement factor (1/v) compensates the relativistic correction v2. ∝ Figure 2. Representative diagrams for local operator renor- The total UV divergence∝ at two loop order is then is malization. The cap represents the insertion of the operator 2 2 1 4 Z α 4ǫ . Define operator renormalization constant via ψN, cross stands for the p relativistic correction, solid square [ψN] = Z ψN, whose value in the MS scheme is then for the Darwin vertex, while empty square for spin-orbital ver- R S tex, which making vanishing contribution in this case. The Z2α2 empty circle represents the contact interaction. The last two ZS =1 + . (21) diagrams are beyond the prescribed accuracy of O(Z2α2). − 4ǫ ··· The anomalous dimension of the operator ψN then reads
2 2 We use dimensional regularization (spacetime dimen- d ln ZS Z α γS = . (22) sion d = 4 ǫ) to regularize UV divergence, which re- ≡ d ln µ 2 spects the all− the cherished symmetry. Power divergences
r 0 r 0 q q = − − UVCT + + UVCT p k p k p k p k r 0 r 0 r 0 = − − UVCT + − + + UVCT r 0 r 0 r 0 = − − UVCT + − + + UVCT
Figure 3. Illustration of the OPE structure of the four-point Green functions through order Z2α2. The first line is for the Wilson coefficient C(1)(r), the two bottom lines for the Wilson coefficient C(2)(r). The thick line indicates the corresponding loop momentum to be “hard” (∼ q).
Sketch of derivation of OPE relation (18). To deduce the property demanded by (18), order by order in Zα. In aforementioned OPE relation, it is convenient to con- [2], we have presented a detailed derivation of the OPE sider the four-point Green function formed by the vac- relation through order Z2α2 for Klein-Gordon case, with uum matrix element of ψ(r)N(0) multiplied with the the technicalities expounded. The calculation in [2] can Fourier-transformed fields ψ†(p) and N †(k). Assuming be directly carried over here, supplemented merely with p, k mZα 1/r m are soft momenta, our task is to the extra complication of including the Darwin term. ∼ ≪ ∼ verify this Green function hase the desirede factorization In Fig. 3, we pictorially show the procedure to deduce 6 the intended coefficient (r) through order Z2α2, by sub- Summary. The Coulomb solution to Dirac equation tracting and adding backC the Green function inserted marks an milestone in the history of quantum mechanics. with the local operator ψN, and expanding the integrand Until very recently, the physics underlying the divergent in the hard sub-loop region (our method is analogous to near-the-origin behavior of Dirac wave function for nS1/2 the Wilson expansion for the product of two φ fields in hydrogen is rather obscure. the λφ4 theory [26]). Here we just present the final re- Together with our recent work [2], we employ the mod- sults in MS scheme: ern field-theoretical machinery such as EFT and OPE to 1 solve the long-standing conundrums related to hydrogen = Z2α2 ln µr + ln 4πeγE−1 , (23) KG and Dirac wave functions. The key is to construct C3b) − 2 the nonrelativistic EFT for atoms by combining NRQED Z2α2 1 3 = ln µr + γ + 2 ln 2π ψ , and HQET. By incorporating the relativistic kinetic cor- C3c) 2 4 E − 2 rection and Darwin term, we show that the universal log- arithmic near-the-origin divergence can be simply inter- where γ is the Euler constant, ψ is the di-Γ function. E preted as the Wilson coefficients arising from the product The first term has already been given in [2], the second of the electron and nucleus fields. RGE further allows us is new. Summing up these two terms, we then obtain the 2 2 to recover the anomalous scaling behavior observed in Wilson coefficient Z α ln µr recorded in (18). − 2 the KG and Dirac wave functions at small r. By our RGE for Wilson coefficient (r). It seems advantageous approach, we can readily trace the different scaling be- C to invoke the powerful renormalization group equation havior of these two types of relativistic wave functions is (RGE) to resum large logarithms appearing in the Wilson due to Darwin term. coefficient in (18) to all orders in Zα, effectively explain A virtue of our field-theoretical approach is that, since the anomalous scaling in (4) for Dirac wave function. OPE holds at operator level, our main result (18) can Demanding the left hand side of (18) is independent of apply to any multi-electron atom as well, not necessarily µ, we then write down the RGE for (r): confined for hydrogen-like atoms. C The folklore says that [14], as the solutions to the rela- ∂ (r, µ) µ C + γS (r, µ)=0, (24) tivistic wave equations, the KG and Dirac wave functions ∂µ C in Coulomb field can be probed in arbitrarily small dis- with the anomalous dimension given in (22). Since µ and tance. Our OPE analysis based on relativistic QED ut- r are glued together inside the logarithms, dimensional terly fails to reproduce the true short-distance divergence analysis leads to (µ ∂ r ∂ ) = 0. We can transform of the KG and Dirac wave functions. In our opinion, this ∂µ − ∂r C the RGE (24) into the evolution equation that controls symptom may indicate that the KG and Dirac equations the scaling of in small r reads of hydrogen are actually secretely nonrelativistic bound- C state equations. It does not make sense to talk about ∂ the near-the-origin behavior of the KG and Dirac wave r C + γS =0. (25) ∂r C functions when r 1 . When probing the short-range ≪ m The solution of (25) becomes behavior of these wave functions, the minimum value of r should be frozen around 1 , since the electron Comp- Z2α2 m r − 2 ton length of constitutes the smallest length scale where (r, µ)= (r0,µ) . (26) nonrelativistic EFT may still apply. C C r 0 Acknowledgment. One of the authors (Y. J.) thanks As anticipated, after resumming the leading logarithms, Xiangdong Ji for spurring his interest in investigating the RGE-improved Wilson coefficient (r) fully agrees the near-the-origin divergence of hydrogen Dirac wave C with the near-the-origin behavior of Dirac wave function function during his visit at University of Maryland in in (4). 2011, and especially for the emphasis that this phe- It is enlightening to pursue the closer contact be- nomenon should be linked with renormalization. This tween our field-theoretical machinery and (4). Let us work is supported in part by the National Natural Sci- na0 define a reference length as r0 = 2 for the nS hy- ence Foundation of China under Grants No. 11875263, drogen state, and µ0 = 1/r0. The boundary condition No. 11621131001 (CRC110 by DFG and NSFC). is (r = r0; µ = µ0)=1. In accord with (17)), pro- videdC that we identify the following local operator matrix element with the nonrelativistic Schrodinger-Pauli wave function at the origin: ∗ [email protected] † 1 Sch [email protected] 0 [ψN]R(0; µ0) nS1/2,m R (0) ξm, (27) ‡ h | | i≈ √4π n0 [email protected] [1] Y. Huang, Y. Jia and R. Yu, arXiv:1809.09023 [hep-ph]. we are then able to fully recover (4). [2] Y. Huang, Y. Jia and R. Yu, arXiv:1812.11957 [hep-ph]. 7
[3] P. A. M. Dirac, Proc. Roy. Soc. Lond. A 117, 610 (1928). Addison-Wesley Publishing Company, 1992. [4] C. G. Darwin. The wave equations of the electron. Pro- [16] E. Eichten and B. R. Hill, Phys. Lett. B 234, 511 (1990). ceedings of the Royal Society A: Mathematical, Physical [17] H. Georgi, Phys. Lett. B 240, 447 (1990). and Engineering Sciences, 118(780):654–680, apr 1928. [18] A. V. Manohar and M. B. Wise, Heavy quark physics, [5] W. Gordon. Die energieniveaus des wasserstoffatoms Camb. Monogr. Part. Phys. Nucl. Phys. Cosmol. 10, 1 nach der diracschen quantentheorie des elektrons. (2000). Zeitschrift für Physik, 48(1-2):11–14, jan 1928. [19] W. E. Caswell and G. P. Lepage, Phys. Lett. 167B, 437 [6] H. A. Bethe and E. E. Salpeter, Quantum Mechanics of (1986). One- and Two-Electron Atoms, Springer, 1977. [20] A. Pineda and J. Soto, Phys. Lett. B 420, 391 (1998) [7] V. B. Berestetskii, E. .M. Lifshitz and L. P. Pitaevskii, [hep-ph/9711292]. Quantum electrodynamics, (translated by J. .B. Sykes [21] P.-O. Löwdin, Phys. Rev. 94, 1600 (1954). and J. S. Bell), Pergamon press, 1982. [22] T. Kato, Communications on Pure and Applied Mathe- [8] B. Durand and L. Durand, Phys. Rev. D 28, 396 (1983) matics 10 (2): 151-177 (1957). Erratum: [Phys. Rev. D 50, 6662 (1994)]. [23] A. Czarnecki and K. Melnikov, Phys. Rev. Lett. 80, 2531 [9] P. Castorina, P. Cea, G. Nardulli and G. Paiano, Phys. (1998) [hep-ph/9712222]. Rev. D 29, 2660 (1984). [24] M. Beneke, A. Signer and V. A. Smirnov, Phys. Rev. [10] J. Malenfant, Phys. Rev. D 36, 863 (1987). Lett. 80, 2535 (1998) doi:10.1103/PhysRevLett.80.2535 [11] J. Malenfant, Phys. Rev. D 38, 3295 (1988). [hep-ph/9712302]. [12] R. Chen, Int. J. Theor. Phys. 47 (2008), 881-890. [25] M. E. Luke, A. V. Manohar and I. Z. Rothstein, Phys. [13] K. G. Wilson, Phys. Rev. 179, 1499 (1969). Rev. D 61, 074025 (2000) [hep-ph/9910209]. [14] C. Itzykson and J. B. Zuber, New York, Usa: Mcgraw-hill [26] J. C. Collins, Renormalization : An Introduction to (1980) 705 P.(International Series In Pure and Applied Renormalization, The Renormalization Group, and the Physics) Operator Product Expansion, Cambridge, University [15] B. R. Holstein, Topics in Advanced Quantum Mechanics, Press (1984).