PHYSICAL REVIEW D 98, 016006 (2018)

Vacuum stabilized by anomalous magnetic moment

Stefan Evans and Department of Physics, The , Tucson, Arizona 85721, USA

(Received 9 May 2018; published 10 July 2018)

An analytical result for Euler-Heisenberg effective action, valid for electron spin g−factor jgj < 2, was extended to the domain jgj > 2 via discovered periodicity of the effective action. This allows for a simplified computation of instability modified by the electron’s measured g ¼ 2.002319. We find a strong suppression of vacuum decay into electron positron pairs when magnetic fields are dominant. The result is reminiscent of mass catalysis by magnetic fields.

DOI: 10.1103/PhysRevD.98.016006

I. INTRODUCTION The here presented results are a step towards under- We explore the effect of anomalous electron spin g-factor standing of nonperturbative QED vacuum structure in E ð 2 − 1Þ on vacuum instability in Euler Heisenberg (EH) effective ultrastrong fields at the scale EH= g= . Here we action [1–3]. In the wake of theoretical development of EH accommodate the nonperturbative character of modifica- ≠ 2 action, studies of the nonlinear QED vacuum have almost tions arising from g explicitly and quantify how for ≠ 2 always relied on an electron spin g-factor of exactly 2. a given value of g modifications of the instability of Since both real and imaginary parts of the effective action the EH effective action arise in presence of external fields. g ≠ 2 In order to complete QED vacuum structure study in are modified by an anomalous -factor, the rate of “ ” α ≃ 1 137 vacuum decay in strong fields by pair production is two loop order: second order in = but to all affected. We study this effect here. orders in external (constant) fields, it is necessary to Employing an anomalous g ≠ 2-factor is an improvement consider (i) magnetic field dependent g, along with (ii) direct self-energy modifications of the particle mass on the EH effective theory in its current form, serving as an – introduction of higher order corrections into the theory. EH m by (constant) external fields [24 27]. Our results j j 2 – presented in this work allow to insert the field dependence action has been extended to g < [4 7]. A finite regular- → ð Þ → ð Þ ized expression for all g was obtained by exploiting perio- g gf a; b ;m mf a; b (subscript f reminds of field dicity of the action [8]. While a magnetic field added to an dependence) to obtain from the one loop EH effective electric field enhances the vacuum decay rate at g ¼ 2 [9,10], action the corresponding two loop result accounting in we will show that for strong magnetic fields a significant our approach for the nonperturbative behavior of EH suppression to particle production arises for g ≠ 2,aneffect effective action as a function of g. reminiscent of magnetic mass catalysis [11–15] in that it can be rather precisely described by modification of particle mass II. MODIFIED FORM OF EFFECTIVE ACTION induced by the magnetic field. A. jgj < 2 Strong magnetic fields are found on surfaces of magnet- B The Schwinger proper time method at g ¼ 2 was ized neutron stars (magnetars), suggested to possess fields j j 2 ∼102 E ¼ 2 extended to g < and produces effective action [3,7] above EH critical field ( EH m =e in Lorentz- Z ∞ Heaviside natural units, where m is electron mass) 1 ds 2 L ¼ −im s – EH 2 3 e [16 18]. Nonlinear QED phenomena in such scenarios have 8π 0 s been studied extensively, see [19,20] and references therein. eas cosh½g eas ebs cos½g ebs Even stronger B fields are produced for ultrashort time 2 2 − 1 ð Þ × ½ ½ ; 1 intervals in relativistic heavy ion collisions [19,21–23]. sinh eas sin ebs where a2 − b2 ¼ E2 − B2 ¼ 2S; ab ¼ E · B ¼ P ð2Þ Published by the American Physical Society under the terms of and α ¼ e2=4π (as used in [3]). We note that charge the Creative Commons Attribution 4.0 International license. renormalization subtraction which modifies in Eq. (1) Further distribution of this work must maintain attribution to 2 2 the author(s) and the published article’s title, journal citation, the term −1 → −1 þða − b Þ=3 is modified for g ≠ 2 and DOI. Funded by SCOAP3. to read −1 → −1 þða2 − b2Þðg2=8 − 1=6Þ.

2470-0010=2018=98(1)=016006(5) 016006-1 Published by the American Physical Society STEFAN EVANS and JOHANN RAFELSKI PHYS. REV. D 98, 016006 (2018) L ð Þ¼L ð Þ ð Þ For a pure electric field, meromorphic expansion and EH gk EH gk−1 ; 9 regularization produces temperature representation [7,28] for any real integer k where Z 2 ∞ X m T 2 2 πg − −2 þ 4 2 þ 4 ð Þ L ¼ ½ − ½1 þ i 2 E=T k

2 2 We can now write the effective action for all values of gk as m− ¼ m − iϵ ð5Þ Z ∞ 1 ds 2 L ¼ −im s EH 2 3 e to offset poles. These two representations of EH action are 8π 0 s only valid for jgj < 2: Eq. (1) diverges if jgj > 2, readily −4 −4 eas cosh½gk k eas ebs cos½gk k ebs shown when writing cosh and cos terms as exponentials, 2 2 − 1 × ½ ½ ; and meromorphic expansion in Eq. (3) fails for jgj > 2. sinh eas sin ebs ð12Þ B. jgj > 2 avoiding any divergences encountered if jg j > 2 appeared j j 2 k To extend to g > , we take advantage of EH action alone in the arguments of cos and cosh. In addition, being periodic in g. We repeat in abbreviated format the meromorphic expansion can now be applied to obtain B arguments and derivations seen in [8]. For a constant temperature representation for any g: Eq. (3) is also field pointing in zˆ the Landau orbit spectrum is obtained generalized to arbitrary gk. Whenever using in the EH using the so called Klein-Gordon-Pauli generalization of action an anomalous magnetic moment, of any magnitude, the Dirac equation and any B-dependence, the periodic reset is implied qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi according to 2 2 En ¼ m þ pz þ QjeBj½ð2n þ 1Þ ∓ g=2; ð6Þ g → g − 2 ð Þ 2 2 k; 13 where Q ¼1. Summing orbit quantum number n, we see that a shift so that the value of g is always in the principal domain −2 2 may never be used in g → g þ 4k ð7Þ Eq. (1) as the EH effective action is ill defined in that case corresponds to a shift in n, leaving the summed states [8]. Equation (13) explains how any value of g [including (Casimir energy, [2]) unchanged. Noting this periodicity gðBÞ] has to be periodically reset, and only the reset value is and using the usual method of discrete summation [1,2] to be used in all expressions involving the EH effective to obtain an integral representation, a periodic function action. introduces the Bernoulli polynomials [29,30]. With this periodic extension of the Bernoulli polynomials, a modified III. STABILIZATION OF THE VACUUM meromorphic expansion is used (Eq. (9) in [8]), agreeing in E limit g → 2 with [10,28,31]. A. Vacuum decay in a pure field At magnetic field strengths beyond We first show the standard case: vacuum decay in a pure electric field. The imaginary part of action (Eq. (14) in [7]) m2 jeBj < ∼ 862m2;g¼ g ð0Þ¼2.002319; ð8Þ becomes jg=2 − 1j f ∞   m2T2 Xð−1Þn g the energies in Eq. (6) become imaginary, that is the ℑ½L ¼ π −nm=T ð Þ EH 8π 2 cos 2n e : 14 associated eigenstates disappear from the spectrum and n¼1 n self-adjointness is lost. Since magnetic fields cannot do work, a compensating modification of mass may stabilize Comparing this result at g ¼ 2.002319 to that at g ¼ 2, the 1% 60E 1% the vacuum in presence of ultrastrong magnetic fields. We two are different by <: for a< EH, and < for 60 100E postpone this question to future study and limit the present

016006-2 VACUUM STABILIZED BY ANOMALOUS MAGNETIC MOMENT PHYS. REV. D 98, 016006 (2018)

B. Imaginary part of action: General form We compute the imaginary part of effective action for nonzero a, b. Using Eqs. (1) and (13), which now produce a nondivergent action for the electron anomalous g-factor, we plug in a normalized expansion of the csch functions, typically used in computing scalar QED action [10]: Z     ∞ −1 ds 2 g g L ¼ −im s EH 3 3 e cosh x cos y 4π 0 s 2 2 ∞ X ð−1Þn csch½nπy=x csch½nπx=y x3y − xy3 ; × 2 þ 2π2 2 − 2π2 n¼1 n x n y n ð15Þ plus renormalizing terms, and where FIG. 1. Ratio ℑ½L ða; b; g ¼ 2.002319Þ=ℑ½L ða; b; g ¼ 2Þ ¼ ¼ ð Þ EH EH x eas; y ebs: 16 is plotted against a, at fixed values of b in units EEH. We can neglect the second term as it does not contribute to the imaginary part of action, and write the first part as: enough b=a dominates the conventional g ¼ 2 vacuum Z decay enhancement by magnetic fields: the eþe− produc- −1 ∞ L ¼ dsð ½ 2 − ½ 2 Þ tion is maximized for b=a ∼ 102.Atb=a ∼ 103, all effects EH 3 3 cos m s isin m s 4π 0 s ∼1 E     cancel and we are left with a ratio of (pure field X∞ n ∼ 100E g g ð−1Þ 3 csch½nπb=a result). For fields of order a EH there is further ×cosh x cos y x y þ; g 2 2 2 þ 2π2 suppression, lowering the peak in Fig. 2 (cos½2 nπ con- ¼1 n x n n tribution discussed in Sec. III A). ð17Þ To better understand the results shown in Fig. 2, we look at the asymptotic limit b=a → ∞ in Eq. (18): where the second term gives after summing residues     X∞ ð Þð Þ X∞ ð−1Þn ðeaÞðebÞ ð−1Þn g ℑ½L ¼ ea eb g π ℑ½L ∼ cos nπ EH 2 cos n EH 8π2 2 8π n 2 ¼1 n n¼1 n g X∞ cosh½ nπb=a 2 ˜ 2 2 −nπm =ea ð Þ × 1 þ 2 e−2rnπb=a e−nπm =ea ð20Þ × ½ π e : 18 sinh n b=a r¼1 → 0 In the limit b we recover Eq. (14):   X∞ n   ðeaÞðebÞ ð−1Þ g 2 2 2 X∞ n −nπm˜ =ea e E ð−1Þ g 2 ∼ cos nπ e ; ð21Þ ℑ½L ¼ cos nπ e−nπm =eE 8π2 n 2 EH 8π2 2π 2 n¼1 n¼1 n ∞   m2T2 X ð−1Þn g ¼ π −nm=T ð Þ 8π 2 cos 2 n e ; 19 n¼1 n the result for pure E fields from [7].

C. Evaluation and asymptotic behavior In contrast to the result for a pure E field, Sec. III A, the result for arbitrary a and b obtained evaluating Eq. (18) exhibits vacuum stabilization by the anomalous magnetic moment, see Fig. 1, significant for large b=a. Pair pro- duction is completely suppressed as b=a → ∞. Figure 2 offers the same result as in Fig. 1, but with B field dependence of pair production made more visible by FIG. 2. ℑ½LEHða; b; gÞ=ℑ½LEHða; b ¼ 0;g¼ 2Þ is depicted as taking the ratio with imaginary action for a pure E field. We a function b=a. Dashed line shows rescaled mass function from see that the anomalous g-induced suppression at large Eq. (21) with g ¼ 2.002319.

016006-3 STEFAN EVANS and JOHANN RAFELSKI PHYS. REV. D 98, 016006 (2018) where rescaled mass We note that for field strengths b=a ∼ 103, the vacuum decay is suppressed by factor ∼ × 10−3. In limit b=a → ∞, j − 4 j 2 2 gk k pair production is completely removed. Our result is a step m˜ ¼ m þ − 1 eb; ð22Þ 2 toward accounting for extreme magnetic field physics effects with b=a ≫ 1 arising in heavy ion collisions and where we have written explicitly the transformation to g0, on magnetars. Eq. (13), to demonstrate the periodicity as a function of g While not the focus of this work, we recognize a growing of our result. The outer absolute value in Eq. (22) assures experimental interest in the high-intensity laser frontier that the mass always increases as is seen in Eq. (20). The [33–35]. Seen our results, the modification of vacuum inner absolute value is due to evenness of the cos and cosh polarization for anomalous g-factors [36–38] should be terms in Eq. (18). The result obtained using rescaled mass reconsidered: our results differ due to recognition of the according to Eq. (21) is shown as a dashed line in Fig. 2 and periodicity in g, potentially resulting in experimentally is valid for b=a > 1. discernible differences between the cited results, and a This form of mass rescaling shares similarities with mass periodic in g result for the effect of vacuum polarization. – catalysis computations [11 15], since the leading correc- Another concern we have is that in the prior work on the α tion to g is linear in fine-structure constant . self-energy in leading OðαÞ, B-dependent corrections were B E found to be negative within the range < EH, while being IV. CONCLUSION B ≫ E – positive for EH [24 27]. These results will need to be An extension of the anomalous moment modification to reinspected to account for the nonperturbative character of two-loop effective action is of interest [32]. We have made g-factor periodicity. here a step in this direction and have demonstrated that To conclude: we presented QED vacuum instability in due to the electron’s anomalous magnetic moment, vac- strong magnetic fields allowing for g ≠ 2. Our work is the uum instability by decay into eþe− pairs experiences large first step toward a complete characterization at two loop modification for large b=a. This effect can be described in order of strong B-field behavior, valid to all orders in compact form by rescaling mass, Eq. (21). We have shown magnetic field. We recognize new nonperturbative effects that the effect of the gyromagnetic anomaly compensates related to periodicity as function of g. Our results apply in and overwhelms the Schwinger g ¼ 2 result where B straightforward fashion to consideration of B-dependence → ð Þ → ð Þ enhances pair production. We thus believe that an evalu- in g gf a; b and m mf a; b . We have shown that ation of electron-positron pair production for nearly pair production suppression occurs irrespective of the sign constant fields on the scale of the order of electron of the g-anomaly. We believe that the periodicity in g Compton wave length λC, such that EH action is valid, exploited here will facilitate future computation of effective must include this effect which dominates for strong B. action including self-energy corrections.

[1] W. Heisenberg and H. Euler, Folgerungen aus der [8] J. Rafelski and L. Labun, A cusp in QED at g ¼ 2, Diracschen Theorie des Positrons, Z. Phys. 98, 714 arXiv:1205.1835. (1936). [9] L. Labun and J. Rafelski, Vacuum decay time in strong [2] V. Weisskopf, The electrodynamics of the vacuum based on external fields, Phys. Rev. D 79, 057901 (2009). the quantum theory of the electron, Kong. Dan. Vid. Sel. [10] Y. M. Cho and D. G. Pak, Effective Action: A Convergent Mat. Fys. Med. 14, 1 (1936). Series of QED, Phys. Rev. Lett. 86, 1947 (2001). [3] J. S. Schwinger, On gauge invariance and vacuum polari- [11] V. P. Gusynin, V. A. Miransky, and I. A. Shovkovy, Dynami- zation, Phys. Rev. 82, 664 (1951). cal chiral symmetry breaking by a magnetic field in QED, [4] R. F. O’Connell, Effect of the anomalous magnetic moment Phys. Rev. D 52, 4747 (1995). of the electron on the nonlinear Lagrangian of the electro- [12] V. P. Gusynin, Magnetic catalysis of chiral symmetry break- magnetic field, Phys. Rev. 176, 1433 (1968). ing in gauge theories, Ukr. J. Phys. 45, 603 (2000). [5] W. Dittrich, One loop effective potential with anomalous [13] E. J. Ferrer and V. de la Incera, Magnetic catalysis in the moment of the electron, J. Phys. A 11, 1191 (1978). presence of scalar fields, Phys. Lett. B 481, 287 (2000). [6] S. I. Kruglov, Pair production and vacuum polarization of [14] E. Elizalde, E. J. Ferrer, and V. de la Incera, Beyond constant arbitrary spin particles with EDM and AMM, Ann. Phys. mass approximation magnetic catalysis in the gauge Higgs- (N.Y.) 293, 228 (2001). Yukawa model, Phys. Rev. D 68, 096004 (2003). [7] L. Labun and J. Rafelski, Acceleration and vacuum temper- [15] I. A. Shovkovy, Magnetic catalysis: A review, Lect. Notes ature, Phys. Rev. D 86, 041701 (2012). Phys. 871, 13 (2013).

016006-4 VACUUM STABILIZED BY ANOMALOUS MAGNETIC MOMENT PHYS. REV. D 98, 016006 (2018)

[16] R. C. Duncan and C. Thompson, Formation of very strongly [28] B. Müller, W. Greiner, and J. Rafelski, Interpretation magnetized neutron stars—implications for gamma-ray of external fields as temperature, Phys. Lett. A 63, 181 bursts, Astrophys. J. 392, L9 (1992). (1977). [17] C. Kouveliotou et al., An X-ray pulsar with a superstrong [29] L. Euler, Methodus succincta summas serierum infinitarum magnetic field in the soft gamma-ray repeater SGR 1806- per formulas differentiales investigandi, M´emoires de l’A- 20., Nature (London) 393, 235 (1998). cad´emie Imp´eriale des Sciences de St.-P´etersbourg 5,45 [18] A. K. Harding and D. Lai, Physics of strongly magnetized (1815). neutron stars, Rep. Prog. Phys. 69, 2631 (2006). [30] T. M. Apostol, An elementary view of Euler’s summation [19] R. Ruffini, G. Vereshchagin, and S. S. Xue, Electron- formula, Am. Math. Mon. 106, 409 (1999). positron pairs in physics and astrophysics: From heavy [31] G. V. Dunne, Heisenberg-Euler effective Lagrangians: nuclei to black holes, Phys. Rep. 487, 1 (2010). Basics and extensions, in From Fields to Strings, edited [20] S. P. Kim, H. K. Lee, and Y. Yoon, Pair production density by M. Shifman et al. (World Scientific, Singapore, 2005), in strong QED, Int. J. Mod. Phys. Conf. Ser. 01, 303 (2011). Vol. 1, p. 445. [21] K. Rumrich, K. Momberger, G. Soff, W. Greiner, N. Grun, [32] M. Reuter, M. G. Schmidt, and C. Schubert, Constant and W. Scheid, Nonperturbative Character of Electron- external fields in gauge theory and the spin 0, 1=2, 1 path Positron Pair Production in Relativistic Heavy Ion Colli- integrals, Ann. Phys. (N.Y.) 259, 313 (1997). sions, Phys. Rev. Lett. 66, 2613 (1991). [33] A. Di Piazza, C. Müller, K. Z. Hatsagortsyan, and C. H. [22] G. Baur, K. Hencken, and D. Trautmann, Electron-positron Keitel, Extremely high-intensity laser interactions with pair production in relativistic heavy ion collisions, Phys. fundamental quantum systems, Rev. Mod. Phys. 84, 1177 Rep. 453, 1 (2007). (2012). [23] J. Rafelski, J. Kirsch, B. Müller, J. Reinhardt, and W. [34] B. M. Hegelich, G. Mourou, and J. Rafelski, Probing the Greiner, Probing QED vacuum with heavy ions, in New quantum vacuum with ultra intense laser pulses, Eur. Phys. Horizons in Fundamental Physics, edited by S. Schramm J. ST 223, 1093 (2014). and M. Schäfer, FIAS Interdisc. Sci. Ser. (Springer [35] B. M. Hegelich, L. Labun, and O. Z. Labun, Finding International Publishing, Cham, Switzerland, 2017), p. 211. quantum effects in strong classical potentials, J. Plasma [24] B. Jancovici, Radiative correction to the ground-state energy Phys. 83, 595830301 (2017). of an electron in an intense magnetic field, Phys. Rev. 187, [36] M. J. G. Veltman, Two component theory and electron 2275 (1969). magnetic moment, Acta Phys. Pol. B 29, 783 (1998). [25] R. G. Newton, Atoms in superstrong magnetic fields, Phys. [37] R. Angeles-Martinez and M. Napsuciale, Renormalization Rev. D 3, 626 (1971). of the QED of second order spin 1=2 fermions, Phys. Rev. D [26] D. H. Constantinescu, Electron self-energy in a magnetic 85, 076004 (2012). field, Nucl. Phys. B44, 288 (1972). [38] C. A. Vaquera-Araujo, M. Napsuciale, and R. Angeles- [27] W. Y. Tsai and A. Yildiz, Motion of an electron in a Martinez, Renormalization of the QED of self-interacting homogeneous magnetic field-modified propagation function second order spin 1=2 fermions, J. High Energy Phys. 01 and synchrotron radiation, Phys. Rev. D 8, 3446 (1973). (2013) 011.

016006-5