arxiv: v1 [physics.atom-ph] 22 Jan 2010

Transcription

1 Self Energy Correction to the Bound Electron g Factor of P States Ulrich D. Jentschura Departent of Physics, Missouri University of Science and Technology, Rolla, Missouri , USA arxiv: v1 [physics.ato-ph] 22 Jan 21 The radiative self-energy correction to the bound-electron g factor of 2P 1/2 and 2P 3/2 states in one-electron ions is evaluated to order α(zα) 2. The contribution of high-energy virtual photons is treated by eans of an effective Dirac equation, and the result is verified by an approach based on long-wavelength quantu electrodynaics. The contribution of low-energy virtual photons is calculated both in the velocity and in the length gauge and gauge invariance is verified explicitly. The results copare favorably to recently available nuerical data for hydrogenlike systes with low nuclear charge nubers. PACS nubers: 12.2.Ds, 31.3.js, p, 6.2.Jr I. INTRODUCTION When a bound electron interacts with an external, unifor and tie-independent agnetic field (Zeean effect), the energetic degeneracy of the atoic energy levels with respect to the agnetic projection quantu nuber is broken, and the different agnetic sublevels split according to the forula E = g j µ B Bµ (1) where g j is the bound-electron (Landé) g factor, µ B = e/(2) is the Bohr agneton, and B is the agnetic field whichisassuedtobe orientedparalleltothe quantization axis. Finally, µ is the agnetic projection quantu nuber ofthe electron; i.e. the projectionofits total angular oentu (divided by ) onto the quantization axis. In leading order, the bound-electron g factor is deterined by nonrelativistic quantu theory and is equal to a rational nuber for all bound states in a hydrogenlike ion. Both relativistic atoic theory as well as quantu electrodynaics (QED) predict deviations fro the nonrelativistic result. The relativistic effects follow fro Dirac theory and can be expressed in ters of a power seriesin the paraeterzα, wherez is the nuclearcharge nuber and α is the fine-structure constant. The QED effects are caused ainly by the anoalous agnetic oent of the electron, which is turn in caused by the exchange of high-energy virtual photons before and after the interaction with the external agnetic field. Here, by high-energy we refer to a virtual photon with an energy of the order of the electron rest ass. A second source for QED effects are exchanges of virtual photons with an energy coensurate with the atoic binding energy scale, which is saller than the electron rest ass energy by a factor (Zα) 2. Here, the electron eits and absorbs a virtual photon before and after the interaction with the external agnetic field, undergoing a virtual transition to a excited atoic state in the iddle. For P states, the latter effects lead to a correction to the bound-electron g factor of order α(zα) 2. The coplete result for the correction of order α(zα) 2 is obtained after adding the anoalous agnetic oent correction (high-energy part) and the low-energy photon contribution of the sae order. Previous studies of the bound-electron g factor for P states in hydrogenlike systes include Refs. [1 5]. Quite recently, the proble has received renewed interest [6, 7]. For few-electron ions, the bound-electron g factor has been investigated in Refs. [6 1]. For the 2 3 P states of heliu, there is still an unresolved discrepancy of theoretical and experiental results (see Refs. [1 12]). The expansion of the quantu electrodynaic radiative correction to the electron g factor, which is an expansion in powers of α for a free electron, is intertwined with an expansion in powers of Zα for a bound electron (this fact has been stressed in Ref. [13]). For an np j state in a hydrogenlike syste, we can write down the following intertwined expansion in powers of α and Zα, δg(np j ) = g +(Zα) 2 g 2 n 2 +O(Zα)4 + α π { b +(Zα) 2 b } 2 n 2 +O(Zα)4. (2) The coefficients g and g 2 characterize the relativistic effects, whereas b and b 2 are obtained fro the oneloop radiative correction. The nonrelativistic result for the Landé g factor reads g (np 1/2 ) = 2 3, g (np 3/2 ) = 4 3. (3) The relativistic correction follows fro Breit theory and the Dirac equation in an external agnetic field [3, 14], g 2 (np 1/2 ) = 2 3, g 2(nP 3/2 ) = (4) The leading correction due to the anoalous agnetic reads as (see Refs. [1, 3]), b (np 1/2 ) = 1 3, b (np 3/2 ) = 1 3. (5) We are concerned here with the evaluation of the b 2 coefficient of np j states, which is deterined exclusively by self-energy type corrections(vacuu polarization does not contribute).

2 2 We adopt the following outline for this paper. In Sec. II, we reexaine the contribution of high-energy virtual photons (see also Ref. [3]). Two alternative derivations are presented, which are based on an effective Dirac equation (Sec. II A) and on an effective low-energy longwavelength quantu electrodynaic theory (Sec. II B) which is obtained fro the fully relativistic theory by a cobined Foldy Wouthuysen and Power Zienau transforation [1]. The low-energy part is also treated in two alternative ways. The velocity-gauge calculation in Sec. III A is contrasted with the length-gauge derivation in Sec. III B. Conclusions are reserved for Sec. IV. Natural units ( = c = ǫ = 1) are used throughout the paper. II. HIGH ENERGY PART A. Effective Dirac Equation In Ref. [3], the contribution to b 2 due to high-energy virtual photons was obtained on the basis of the twobody Breit Hailtonian. Here, we perfor the calculation using a siple approach, based on an effective Dirac Hailtonian (see Ch. 7 of Ref. [15]). For an electron interacting with external electric and agnetic fields, this equation reads [ H rad = α p ef 1 ( 2 ) A ] +β+f 1 ( 2 )V +F 2 ( 2 e ( ) i γ 2 E βσ B ). (6) We here take into account the Dirac for factor F 1 and the Pauli for factor F 2. The atrices α = γ γ and β = γ are the standard Dirac atrices in the Dirac representation [15], is the electron ass, and e = e is the electron charge. Up to the order relevant for the current calculation, we ay approxiate both for factors in the liit of vanishing oentu transfer as F 1 ( 2 ) F 1 () = 1, F 2 ( 2 ) F 2 () κ α 2π. (7) The vector potential A corresponds ( to ) a unifor external agnetic field, i.e. A = 1 B 2 r, and the electric field E is that of the Coulob potential (ee = V). Finally, V = Zα/r is the binding potential. So, H rad α p+β+v e 2 α ( B r ) iκ 2 γ V e 2 κβ Σ B. (8) Dirac eigenstates fulfill ( α p+β+v)ψ = E D ψ, where E D is the Dirac energy. The first few ters in the perturbative expansion of E D in a agnetic field read E = ψ e α ( ) B ψ r 2 e 2 κ ψ + eκ 2 βσ B ψ ( ψ i γ V ) Q [ α ( B 2 r)] ψ. (9) where Q = 1 2 (1 γ ) is a projectoronto virtual negativeenergy states, and ψ is the relativistic wave function. An evaluation of the first ter on the right-hand side of (9) with Dirac wave functions confirs the results for g and g 2 given in Eqs. (3) and (4). The second ter on the right-hand side of (9) yields the result for b as given in Eq. (5). When the Dirac wave functions are properly expanded in powersof Zα, the second and third ters on the right-hand side of (9) yield the following high-energy contribution b (H) 2 to the b 2 coefficient defined in Eq. (2), b (H) 2 (np 1/2) = 1 2, b(h) 2 (np 3/2) = 1 1. (1) Theseresultsareinagreeentwith thosegivenineq.(5) of Ref. [3]. B. Long Wavelength Quantu Electrodynaics It is instructive to copare the fully relativistic approach outlined above to an effective nonrelativistic theory. In Ref. [16], a systeatic procedure has been described in order to perfor a nonrelativistic expansion of the interaction Hailtonian for a light atoic syste with slowly varying external electric and agnetic fields. This procedure involves two steps, (i) a Foldy Wouthuysen transforation of an interaction of the type (6), suitably generalized for any-electron systes, and (ii) a Power Zienau transforation to express the vector potentials in ters of physically observable field strengths. The result is an interaction, given in Eq. (3) of Ref. [16], which describes a nonrelativistic expansion of the ato-field interaction in powers of Zα and can be used in order to identify ters which contribute at a specified order. If we are interested in evaluating the corrections to the g factor up to order α(zα) 2, i.e. all corrections listed in Eq. (2), the relevant effective interactions for a one-

3 3 electron syste are H ag = H M + 3 H i, H M = µ B ( L+ σ) B, (11a) i=1 H 1 = µ B 2 2 p2 ( L+ σ) B, H 2 = µ B(1+2κ) Zα 4 r 3 ( r σ) ( r B), H 3 = µ Bκ 22( p σ)( p B), (11d) (11b) (11c) where µ B = e/(2) is the Bohr agneton. We denote the Schrödinger Pauli two-coponent wave function by φ in order to distinguish it fro the fully relativistic wave function ψ. Specifically, φ reads as φ( r) = R(r)χ µ κ(ˆr) in the coordinate representation, where R(r) is the nonrelativistic radial coponent of the wave function and χ µ κ (ˆr) is the standard two-coponent spin-angular function [17]. An evaluation of the perturbation E = φ H ag φ (12) confirs the results of Eqs. (3), (4), (5) and (1) for the high-energy part. No second-order effects need to be considered in this foralis up to the order in the Zα-expansion relevant for the current study. III. LOW ENERGY PART A. Velocity Gauge The ost econoical approach to the calculation of the low-energycontributionoforderα(zα) 2 totheg factorof P states consists in a calculation of the orbital g l factor, with a conversion of the orbital g l factor to the Landé g j factor in a second step of the calculation. In the order α(zα) 2, one ay indeed convert the spin-independent correction to the orbital g l factor to a spin-dependent correction to the g j factors of the 2P 1/2 and 2P 3/2, as described in Appendix A. However, a ore systeatic approach to the proble, which is also applicable to higherorder (in Zα) corrections, is based on a perturbation of the nonrelativistic self-energy of the bound electron by the agnetic interaction Hailtonian (11). This is the approach outlined below. We thus investigate the perturbation of the nonrelativistic bound-electron self-energy [18] due to the agnetic interaction Hailtonian H M = e 2 ( L+ σ ) B (13) given in Eq. (11a). In the velocity gauge, the interaction of the electron with the vector potential A of the quantized electroagnetic field is given by the ter p A/, where p is the electron oentu. In an external agnetic field, it is the physical oentu p e 2 ( B r), not the canonical oentu p, which couples to the quantized electroagnetic field. This aount to a correction δ J to the electron s transition current given by δj i = e ( ) i B e r = 2 2 ǫijk B j r k, (14) Because of the syetry of the proble (Wigner Eckhart theore), we ay fix the axis of the B field to be along the quantization axis(z axis) and the projection of the reference state to be µ = 1 2. This procedure allows one to siplify the angularalgebra. It is inspired by the separation of nuclear and electronic tensors that are responsible for the hyperfine interaction. Such a separation has been used in Eq. (1) of Ref. [19] and in Eqs. (1) and (11) of Ref. [2]. In the case of the g factor, the agnetic field of the nucleus is replaced by the external, hoogeneous agnetic field of the Zeean effect. Wedivide outaprefactor e/(2)froboth theagnetic Hailtonian H M and fro the current δj i and obtain the perturbative Hailtonian h M, and the scaled current δj i. In the spherical basis, this procedure leads to the operators h M, = L +σ, δj i = ǫi3k r k, (15) which are aligned along the quantization axis of the external agnetic field. The index zero of the operators, in the spherical basis, denotes the z coponent in the Cartesian basis (see Ref. [17]). The following shorthand notation for the atoic states with agnetic projection µ = 1 2 proves useful, j 1 2 np j(µ = 1 2 ), h M, j 1 2 h M, j 1 2. (16) Finally, we can proceed to the calculation of the perturbed self-energy. The nonrelativistic (Schrödinger) Hailtonian of the ato is H NR = p2 +V, (17) 2 and the nonrelativistic self-energy reads δe = 2α 3π 2 ǫ j 1 2 p 1 H NR E NR +ω p j 1 2. The wave-function correction to the self-energy reads ǫ (18) δe ψ = 4α 3π 2 (19) ( ) j 1 2 p 1 H NR E NR +ω p 1 h M, E NR H NR j 1 2. Here, [1/(E NR H NR )] is the reduced Green function, with the reference state being excluded fro the su over

4 4 virtual states. The contribution of virtual np j states (with j being equal to that of the reference state and n 2) vanishes because of the orthogonality of the nonrelativistic radial wave functions. The interaction h M, couples 2P 1/2 and 2P 3/2 states, but the contribution of virtual states with different j as copared to the reference state vanishes after angular algebra [17] because the self-energyinteractionoperator{ p[1/(h NR E NR +ω)] p} is diagonal in the total angular oentu. Virtual states with different orbital angular oentu than the reference state are not coupled at all to the reference state by the action of the perturbative Hailtonian h M,. Because all contributions vanish individually, we can thus conclude that δe ψ =. Hence, we have to evaluate first-order corrections to the Hailtonian and to the energy corresponding to the replaceents H NR H NR + h M, and E NR E NR + h M, in Eq. (18). Furtherore, we have a correction to the current corresponding to p/ p/+δ j. The energy correction reads as δe E = 2α 3π 2 h M, (2) ǫ j 1 2 ( ) 2 p 1 p H NR E NR +ω j 1 2, where ǫ is an upper cutoff for the photon energy [18] (the scale-separation paraeter ǫ is diensionless). It corresponds to the following g factor correction, δg E (np j ) = g (np j ) δe E h M,. (21) A nuerical evaluation of this correction according to established techniques [21] yields δg E (2P 1/2 ) = 2α ( 3π (Zα)2 1 ( ) ) ǫ 6 ln (Zα) , (22a) δg E (2P 3/2 ) = 4α ( 3π (Zα)2 1 ( ) ) ǫ 6 ln (Zα) For the correction to the Hailtonian, we get δe H = 2α 3π 2 j 1 2 ǫ ( ) 2 p 1 h M, p H NR E NR +ω j1 2 (22b), (23) where we have used the relation [H NR,h M, ] =. This translates into the following correction for the g factor, δg H (np j ) = g (np j ) δe H h M,. (24) A nuerical evaluation leads to the following results, δg H (2P 1/2 ) = 2α ( ( ) ) 1 ǫ 3π (Zα)2 6 ln (Zα) , (25a) δg H (2P 3/2 ) = 4α ( ( ) ) 1 ǫ 3π (Zα)2 6 ln (Zα) The correction to the current is given by δe C = 4α 3π 2 j 1 2 ǫ (25b) 1 pi H NR E NR +ω δji j 1 2, (26) where we take into account the ultiplicity factor due to the current acting on both sides of the propagator. The corresponding correction to the g factor is δg C (np j ) = g (np j ) We obtain the following nuerical results, δg C (2P 1/2 ) = 2α 3π (Zα)2.2467, δg C (2P 3/2 ) = 4α 3π (Zα) δe j h M,. (27) (28a) (28b) Suing all low-energy corrections, the spurious logarithic ters cancel, and we obtain b (L) 2 (2P 1/2) =.98428, (29a) b (L) 2 (2P 3/2) =.49216, (29b) as the spin-dependent low-energy contribution to the bound-electron g factor. The result for 2P 3/2 is equal to half the correction for 2P 1/2 (this fact is independently proven also Appendix A). We denote the low-energy contribution to the b 2 coefficient defined in Eq. (2) by b (L) 2. B. Length Gauge In the length gauge, the interaction with the quantized electroagnetic field is given by the dipole interaction e x E, where E is the electric-field operator [22]. The gauge-invariant[22] nonrelativistic self-energy in the length-gauge reads δe = 2α 3π ǫ j x 1 H NR E NR +ω x j 1 2. (3)

5 5 In the length gauge, the contribution of the wave-function correction vanishes because of the sae reasons as for the velocity gauge. Also, there is no correction to the transition current, because the canonical oentu does not enter the interaction Hailtonian in the length gauge. We only have corrections to the Hailtonian and to the energy. We start with the energy perturbation, δe E = 2α 3π h M, (31) ǫ ( ) 2 j x 1 x H NR E NR +ω j 1 2. The subscript E instead of E serves to differentiate the length-gauge as opposed to the velocity-gauge for of the correction. Indeed, the nuerical results for the corresponding correction to the bound-electron g factor are different fro those given in Eq. (22) and read δg E (2P 1/2 ) = 2α ( 3π (Zα)2 1 ( ) ) ǫ 2 ln (Zα) , (32a) δg E (2P 3/2 ) = 4α ( 3π (Zα)2 1 ( ) ) ǫ 2 ln (Zα) For the correction to the Hailtonian, we get δe H = 2α 3π j 1 2 ǫ 3 (32b) ( ) 2 x 1 h M, x H NR E NR +ω j 1 2. (33) A nuerical evaluations leads to δg H (2P 1/2 ) = 2α ( 1 3π (Zα)2 2 ln ( ǫ (Zα) 2 ) ) , (34a) δg H (2P 3/2 ) = 4α ( ( ) ) 1 ǫ 3π (Zα)2 2 ln (Zα) (34b) Adding the length-gauge corrections, the logarithic ters cancel, and it is straightforward to nuerically verify the gauge-invariance relation and with the help of a soewhat lengthy calculation, it is possible to show analytically that the velocity-gauge and the length-gauge fors of the low-energy contributions are equal. The calculation follows ideas outlined in detail in Ref. [23] where the ore coplicated case of a relativistic correction to a transition atrix eleent was considered. Here, we are interested ainly in the nuerical value of the correction, for which the gauge invariance provides a highly nontrivial check. Note that the atrix eleents governing the transitions of the reference to the virtual states are copletely different in the length and in the velocity gauges, and the final results are obtained after suing over the discrete and continuous parts of the spectru of virtual states. IV. CONCLUSIONS In our approach to the calculation of the boundelectron g factor of P states, the contribution due to high-energy virtual photons can be obtained using two alternative approaches, based either on an effective Dirac equation or on a low-energy effective Hailtonian. The contribution due to low-energy photons is treated as a perturbation of the bound-electron self-energy [18] due to the interaction with the external unifor agnetic field. Corrections to the Hailtonian, to the bound-state energy and (in the velocity gauge) to the transition current have to be considered. The final results for the lowenergy parts in the velocity- and length-gauges agree although the individual contributions differ (including the coefficients of spurious logarithic ters). Adding the high-energy contribution to the g factor correction given in Eq. (1) and the low-energy effect given in Eq. (29), we obtain the following results for the self-energy correction of order α(zα) 2 to the boundelectron g factor of 2P states (b 2 = b (H) 2 +b (L) 2 ), b 2 (2P 1/2 ) =.48429, b 2 (2P 3/2 ) =.59214, (37a) (37b) wheretheb 2 coefficienthasbeendefinedineq.(2). Both above results copare favorably with recently obtained nuerical data for low-z hydrogenlike ions [24] (see also Appendix B). An obvious generalization of the foralis outlined here to the 3P and 4P states yields the results b 2 (3P 1/2 ) =.45, b 2 (3P 3/2 ) =.5525, (37c) (37d) δg E +δg H +δg C = δg E +δg H (35) b 2 (4P 1/2 ) =.31331, (37e) and thus, the nuerical results already given in Eq. (29). Let us finally discuss the analytic proof of the gauge invariance. Using the coutator relation p i = i[h NR E NR +ω, x i ], (36) b 2 (4P 3/2 ) = (37f) The two ain results of the current investigation can be suarized as follows. First, in Sec. III we forulate a generalizable procedure for the calculation of lowenergycorrectionsto the Landé g j factors in one-electron

6 6 ions, applicable to P states and states with higher angular oenta. This procedure is based on choosing a specific reference axis for the external agnetic field. In the future, it ight be applied to include higher-order ters fro the Hailtonian (11) which couple the orbital and spin degrees of freedo. Second, we resolve the discrepancy reported in Ref. [24] regarding the low Z liit of the α(zα) 2 correction to the g j factor with previous results reported in Ref. [5] for this correction [see Eq. (37) and Appendix B]. In Appendix A, it is shown that the discrepancy to the results of Ref. [5] can be traced to the final evaluation of the logarithic sus over virtual states, while the angular oentu algebra is in agreeent. That is a further reasonwhy the cross-checkofour calculation in the length and velocity gauges appeared to be useful. Regarding the experiental usefulness of the obtained results, we can say that recent proposals [25] concerning easureents of the bound-electron g factor of low Z hydrogenlike ions are based on double-resonance schees that also involve transitions to 2P states in the presence of the strong agnetic fields of Penning traps. In order to fine-tune the double-resonance setup, the results obtained here ight be useful. Also, the results reported here serve as a general verification for the analytic foralis used in the theoretical treatent of α 3 corrections to the g factor of 2 3 P states in heliu, for which an interesting discrepancy of experiental and theoretical results persists (see Ref. [12] and Sec. V of Ref. [16]). We conclude with the following reark. Our calculation concerns the g factor of P states, and we confir that in the order α(zα) 2, low-energy virtual photons yield an iportant contribution. Fro the treatent in Sec. III, we can understand physically why there is no such low-energy effect of order α(zα) 2 for S states. Naely, a Schrödinger Pauli S 1/2 state happens to be an eigenstate of the Hailtonian h M,. We have h M, ns 1/2 (µ = ± 1 2 ) = ± ns 1/2(µ = ± 1 2 ). This property holds because an S state carries no orbital angular oentu, and therefore is an eigenstate of the third coponent of the spin operator σ, and b 2 therefore vanishes for S states [26]. For P states and states with higher orbital angular oenta, the situation is different: thesestatesarenoteigenstateofh M, irrespectiveoftheir angular oentu projection, even though h M, coutes with the nonrelativistic Hailtonian H NR. Therefore, there is a residual effect of order α(zα) 2 due to low-energy virtual photons for states with nonvanishing angular oenta. Institutes of Standards and Technology (NIST). TABLE I: Higher-order reainder functions i 1/2 (Zα) and i 3/2 (Zα) for the self-energy correction to the g factor of 2P 1/2 and 2P 3/2 states, respectively, as obtained recently in Ref. [24]. The value of α eployed in the calculation is α 1 = , and the nuerical uncertainty of the all-order (in Zα) calculation due to the finite nuber of integration points in the nuerical calculation is indicated in brackets. Z g 1/2 (Zα) g 3/2 (Zα) (21).14814(21) (22) (24) (19) (24) (14) (22) (1) (18) (7) (14) (5).1535(11) (4).15933(7) (2) (4) (2) (4) FIG. 1: (Color online) The higher-order reainder function i 1/2 (Zα) is shown as a function of Z. Nuerical values for i 1/2 (Zα) are given in Table I. The point at Z = is given by the coefficient 1 4 b2(2p 1/2) and is approached soothly. Acknowledgents Valuable and insightful discussions with V. A. Yerokhin and K. Pachucki are gratefully acknowledged. This project was supported by the National Science Foundation (Grant PHY ) and by a precision easureent grant fro the National FIG. 2: (Color online) Sae as Fig. 1, but for the 2P 3/2 state. The higher-order reainder function i 3/2 (Zα) is plotted as a function of Z, with nuerical values for i 3/2 (Zα) given in Table I. The point at Z = is li Zα i 3/2 (Zα) = 1 4 b2(2p 3/2). The liit is approached soothly.

7 7 Appendix A: Rearks on the Low Energy Part First of all, let us reark that our results reported in Sec. III can be expressed as logarithic sus over the spectru of atoic hydrogen. After perforing the angular algebra [17], we can write the total low-energy correction δg (L) (np j ) = δg E (np j )+δg H (np j ) as follows, δg (L) (np 1/2 ) = 8α 9π (E S E np ) 2 S r np 2 ( ) 2 ES E np ln (Zα) 2 8α (E D E np ) 2 9π and ( ) np r nd 2 2 ED E np ln (Zα) 2 δg (L) (np 3/2 ) = 4α 9π (A1) (E S E np ) 2 S r np 2 ( ) 2 ES E np ln (Zα) 2 4α (E D E np ) 2 9π ( ) np r nd 2 2 ED E np ln (Zα) 2. (A2) Here, S r np and D r np are reduced atrix eleents in the notation of Ref. [27]. Because (E S E np ) (Zα) 2 and S r np (Zα) 1, the above correctionstotheg factorareanifestlyoforderα(zα) 2. The sus over extend over both the discrete as well as the continuous part of the hydrogen spectru and can conveniently be evaluated using basis-set techniques [28]. Because it ay not be copletely evident fro the presentation in Ref. [5], we reephasize here that the authors of the cited article evaluate a correction δg l to the orbital g factor according to the definition H eff = g s µ B S B +gl µ B L B (A3) for the effective interaction of a bound electron with an external agnetic field ( S = 1 2 σ easures the electron spin, and we have g S 2, g l 1). Here, an interesting analogy to the description hyperfine splitting can be drawn, because the interaction with the nuclear agnetic field can also be separated into distinct coponents, naely the orbital, spin-dipole, and Feri contact ters given, e.g., in Eqs. (24) (33) of Ref. [29]. Returning to the discussion of the g factor, we see that as a cobination of g s and g l, the Landé g factor is obtained as g j = g l j(j +1)+l(l+1) s(s+1) 2j(j +1) +g s j(j +1)+s(s+1) l(l+1) 2j(j +1) 3j(j +1) l(l+1)+s(s+1) 2j(j +1). (A4) In our notation, the analytic result given in Eq. (14) of Ref. [5] for the correction to the orbital g factor reads δg l (np) = α (E S E np ) 2 S r np 2 3π ( ) 2 ES E np ln (Zα) 2 α (E D E np ) 2 3π ( ) np r nd 2 2 ED E np ln (Zα) 2. (A5) The prefactors ultiplying g l in Eq. (A4) read 4 3 for P 1/2 and 2 3 for P 3/2 states, and therefore the analytic forula obtained in Eq. (14) of Ref. [5] is in agreeent with our approach for both P 1/2 and P 3/2 states, after the correction to g l is converted into the corresponding odification of g j. However, their nuerical result δg l =.24α 3 disagrees with our result both in sign and in agnitude. Indeed, expressed in ters of our b (L) 2 coefficient, the results indicated in Ref. [5] would iply that b (L) 2 (2P 1/2) = 4.2 and b (L) 2 (2P 3/2) = 8.4. Finally, let us note that the calculation of the orbital correction to the g j factor is only applicable to order α(zα) 2, not α(zα) 4, because the higher-order ters in the agnetic interaction (11) couple the orbital and spin degrees of freedo. The foralis outlined in Sec. III generalizes easily to the calculation of higher-order corrections that couple spin and orbital angular oentu, which ight be needed in the future, whereas the separation into spin and orbital g factors only holds up to order α(zα) 2. Appendix B: Coparison to Nuerical Data We paraeterize the one-loop self-energy correction δ (1) g to the g factor of P states as δ (1) g(2p j ) = α π { b +(Zα) 2 i j (Zα) }, (B1) where i j (Zα) is the nonperturbative (in Zα) reainder function. The reainder functions i 1/2 (Zα) and i 3/2 (Zα) for 2P 1/2 and 2P 3/2, respectively, have recently been evaluated in Ref. [24]. Nuerical values of i j (Zα) for Z = 1,...,1 are given in Table I. Note that these

8 8 data iply a spin-dependence of the higher-order correction ter b 2 beyond the spin-dependence of the highenergy part given by Eq. (1). The liit as Zα of the reainder function i j (Zα) is As shown in Figs. 1 and 2, this liit is being consistently approached by the nuerical data. li Zα i j(zα) = 1 4 b 2(2P j ). (B2) [1] S. J. Brodsky and R. G. Parsons, Phys. Rev. 163, 134 (1967). [2] R. A. Hegstro, Phys. Rev. A 7, 451 (1973). [3] H. Grotch and R. Kashuba, Phys. Rev. A 7, 78 (1973). [4] H. Grotch and R. A. Hegstro, Phys. Rev. A 8, 2771 (1973). [5] J. Calet, H. Grotch, and D. A. Owen, Phys. Rev. A 17, 1218 (1978). [6] J. M. Anthony and K. J. Sebastian, Phys. Rev. A 48, 3792 (1993). [7] I. Gonzalo and E. Santos, Phys. Rev. A 56, 3576 (1997). [8] Z. C. Yan and G. W. F. Drake, Phys. Rev. A 5, R198 (1994). [9] Z. C. Yan, Phys. Rev. A 66, 2252 (22). [1] K. Pachucki, Phys. Rev. A 69, 5252 (24). [11] M. L. Lewis and V. W. Hughes, Phys. Rev. A 8, 2845 (1973). [12] C. Lhuillier, J. P. Faroux, and N. Billy, J. Phys. (Paris) 37, 335 (1976). [13] U. D. Jentschura and J. Evers, Can. J. Phys. 83, 375 (25). [14] G. Breit, Nature (London) 122, 649 (1928). [15] C. Itzykson and J. B. Zuber, Quantu Field Theory (McGraw-Hill, New York, 198). [16] K. Pachucki, Phys. Rev. A 71, 1253 (25). [17] D. A. Varshalovich, A. N. Moskalev, and V. K. Khersonskii, Quantu Theory of Angular Moentu (World Scientific, Singapore, 1988). [18] H. A. Bethe, Phys. Rev. 72, 339 (1947). [19] J. Bieroń, F. A. Parpia, C. Froese Fischer, and P. Jönsson, Phys. Rev. A 51, 463 (1995). [2] V. A. Yerokhin, A. N. Arteyev, V. M. Shabaev, and G. Plunien, Phys. Rev. A 72, 5251 (25). [21] U. Jentschura and K. Pachucki, Phys. Rev. A 54, 1853 (1996). [22] U. D. Jentschura and C. H. Keitel, Ann. Phys. (N.Y.) 31, 1 (24). [23] B. J. Wundt and U. D. Jentschura, Phys. Rev. A 8, 2255 (29). [24] V. A. Yerokhin and U. D. Jentschura, Phys. Rev. A 81, 1252 (21). [25] W. Quint, B. Nikoobakht, and U. D. Jentschura, Pis a v ZhETF 87, 36 (28), [JETP Lett. 87, 3 (28)]. [26] K. Pachucki, U. D. Jentschura, and V. A. Yerokhin, Phys. Rev. Lett. 93, 1541 (24), [Erratu Phys. Rev. Lett. 94, (25)]. [27] A. R. Edonds, Angular Moentu in Quantu Mechanics (Princeton University Press, Princeton, New Jersey, 1957). [28] S. Saloonson and P. Öster, Phys. Rev. A 4, 5559 (1989). [29] V. A. Yerokhin, Phys. Rev. A 78, (28).