Pith. sign in

REVIEW 1 major objections 5 minor 62 references

Exact QED vacuum birefringence shifts magnetar resonance by up to 10x

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · glm-5.2

2026-07-08 06:06 UTC pith:3D5WVMZC

load-bearing objection Solid quantitative study of finite-field QED corrections to magnetar polarization observables; deserves a serious referee. the 1 major comments →

arxiv 2607.06422 v1 pith:3D5WVMZC submitted 2026-07-07 astro-ph.HE hep-ph

Finite-Field QED Corrections to Vacuum Birefringence and Magnetar Polarization Transport

S. Abbassi , F. A. Chishtie , S. R. Valluri This is my paper
classification astro-ph.HE hep-ph PACS 12.20.Ds95.30.Gv97.60.Jd
keywords finite-fieldone-loopbirefringencecorrectionsfactorsmagnetarmagneticmaximum
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper argues that retaining the exact one-loop Heisenberg-Euler finite-field expressions for vacuum birefringence (keeping the refractive-index normalization gamma_s unexpanded) produces quantitatively significant corrections to magnetar vacuum-resonance observables. The standard weak-field Cotton-Mouton approximation, which treats birefringence as quadratic in B/B_cr, is shown to overestimate the birefringent phase near the surface of strongly magnetized neutron stars by up to a factor of 2.9 at 10^15 G. The polarization-limiting radius, however, is unchanged to better than one part in 10^12 because mode decoupling occurs at ~100 stellar radii where the field has fallen far below the critical field B_cr. At the plasma-vacuum resonance, the paper shows that all atmosphere-dependent factors cancel in the ratio of exact to weak-field birefringence, leaving a single atmosphere-independent correction factor R = Delta_n_weak / Delta_n_exact evaluated at the local field strength. For SGR 1806-20, the resonance density is reduced by a factor of 9.7 and the adiabatic conversion energy is raised by a factor of 2.13. The paper also identifies a broad maximum in the resummed one-loop parallel-mode photon magnetic response near B ~ 17 B_cr (~7.5 x 10^14 G), a feature that does not appear in the strictly truncated O(alpha) response and whose existence is robust but whose sub-percent profile is not order-controlled against omitted two-loop corrections.

Core claim

The central mechanism is the reduction of all finite-field corrections to the plasma-vacuum resonance to a single atmosphere-independent ratio R = Delta_n_weak / Delta_n_exact, made possible by the identity q_hat + m_hat = 2*Delta_n + O(alpha^2) which holds at the working order. This ratio, evaluated at the local field strength of the resonance layer, directly gives the shifts in resonance density (factor R^{-1}) and adiabatic conversion energy (factor R^{1/3}) without requiring knowledge of the atmosphere model. The paper shows this ratio departs significantly from unity for surface fields xi_s > 10, making vacuum-resonance observables in high-field magnetars the most sensitive channel forQ

What carries the argument

The Heisenberg-Euler effective Lagrangian in the one-loop, constant-field approximation, evaluated with the refractive-index normalization gamma_s kept exact (unexpanded). The two photon eigenmode refractive indices are built from the second derivatives of the effective Lagrangian (gamma_GG for the parallel mode, gamma_FF for the perpendicular mode), each normalized by gamma_s = 1 - gamma_F. The ratio R = Delta_n_weak / Delta_n_exact, where Delta_n = n_parallel - n_perpendicular, carries the finite-field correction into resonance observables. The extremum condition (1 + kappa_p) * kappa_p'' = (3/2) * (kappa_p')^2 locates the peak of the resummed response.

Load-bearing premise

The identity q_hat + m_hat = 2*Delta_n + O(alpha^2) reduces all finite-field corrections at the plasma-vacuum resonance to a single ratio R. This identity holds at leading order in the vacuum corrections (O(alpha)), which is the working order throughout the paper. If the O(alpha^2) corrections to this identity are non-negligible at the field strengths of interest (xi ~ 10-45), the atmosphere-independent ratios would acquire atmosphere-dependent corrections, weakening therob

What would settle it

If the O(alpha^2) correction to the identity q_hat + m_hat = 2*Delta_n is not negligible at magnetar field strengths (xi ~ 10-45), then the atmosphere-independent ratios for resonance density and conversion energy would acquire atmosphere-dependent corrections, and the clean source-by-source predictions of Table 2 would no longer hold at the stated precision.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Polarization-transport models for magnetars with surface fields above ~10 B_cr (including 1RXS J1708-4009, a leading near-term IXPE/eXTP target) require the exact finite-field birefringence rather than the Cotton-Mouton approximation; weak-field inputs misplace the resonance density by a factor of 2.6 and the conversion energy by 37% for this source.
  • The highest-field magnetars (SGR-class, B_s ~ 10^15 G or above) offer maximal contrast between exact and weak-field predictions, with resonance density shifts approaching an order of magnitude, making them the sharpest available test of finite-field QED vacuum birefringence.
  • The broadband polarization degree of magnetars, set at the polarization-limiting radius far from the star, is insensitive to finite-field corrections and tests the Heyl-Shaviv enhancement mechanism rather than finite-field QED; this cleanly separates two different physics channels.
  • The broad maximum in the photon magnetic response near 17 B_cr falls within the observed range of magnetar surface fields, suggesting that the known magnetar population straddles a theoretically identified transition scale in the QED vacuum response.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

1 major / 5 minor

Summary. This paper studies low-energy photon propagation in strong magnetic fields using the one-loop Heisenberg-Euler effective Lagrangian, retaining the refractive-index normalization gamma_s without expansion in B/B_cr. The authors derive finite-field refractive indices for both polarization modes, characterize the photon magnetic response over 0 <= xi <= 30, and propagate the resulting birefringence into magnetar polarization transport. The key results are: (1) the polarization-limiting radius r_pl is unchanged to better than 10^-12 because mode decoupling occurs where B << B_cr; (2) the weak-field Cotton-Mouton formula overestimates near-surface birefringent phase by up to a factor 2.9 at 10^15 G; (3) at the plasma-vacuum resonance, finite-field corrections reduce the resonance density by up to a factor 9.7 and raise the adiabatic conversion energy by up to a factor 2.13 for SGR 1806-20. The paper also identifies a broad maximum in the resummed parallel-mode magnetic response near xi ~ 17 B_cr and provides an analytic extremum condition.

Significance. The paper provides a careful, internally consistent finite-field treatment of vacuum birefringence within the one-loop Heisenberg-Euler framework. Several strengths deserve explicit credit: (1) The weak-field limit correctly reproduces the standard coefficients (14/45, 8/45, 2/15), fixing the normalization convention. (2) The internal consistency check on r_pl (agreement to 10^-12) is a convincing null result. (3) The atmosphere-independent ratios R = Delta_n_wf / Delta_n_exact at the plasma-vacuum resonance (Eq. 76) are a clean, falsifiable prediction. (4) The extremum condition (Eq. 66/88) is a general structural result, and its closed-form estimate (xi_peak ~ 16.73, within 1.4% of the exact 16.963) is a useful internal consistency check. (5) The paper is commendably honest about the perturbative status of the maximum: it clearly states that the non-monotonicity arises from the resummed normalization at the ~2% level, comparable to omitted two-loop corrections, and that only the existence and approximate location -- not the sub-percent profile -- are controlled statements. (6) Reproducible companion code is provided. The source-by-source quantification against specific IXPE/eXTP-

major comments (1)
  1. Sec. 5.5, Eqs. (75)-(76): The load-bearing identity q_hat + m_hat = 2*Delta_n + O(alpha^2) reduces all finite-field corrections to the single ratio R = Delta_n_wf / Delta_n_exact. The paper states this identity holds at O(alpha) because q_hat and m_hat are defined from the same one-loop mode-index corrections whose difference is Delta_n. This is plausible, but the identity is not explicitly derived or verified numerically in the manuscript. Given that the atmosphere-independence of the ratios in Eq. (76) -- a central robustness claim -- depends entirely on this identity, an explicit verification (even a brief numerical check at a few representative xi values, or a two-line derivation showing the cancellation explicitly) would substantially strengthen the paper. The paper acknowledges the O(alpha^2) caveat but does not demonstrate the O(alpha) identity itself.
minor comments (5)
  1. Sec. 5.4: The centered-dipole, radial-propagation, theta=pi/2 geometry is a simplification. The paper notes this but does not discuss how multipolar surface fields or non-radial propagation would affect the phase-accumulation results (as opposed to r_pl, which is argued to be robust). A brief sentence estimating the sensitivity would strengthen the near-surface phase-accumulation results.
  2. Eq. (28) and Sec. 5.5: The strong-field asymptote Delta_n ~ (alpha/6*pi)*xi is used for SGR 1806-20 (xi_s = 45.3), which lies beyond the numerically validated interval xi <= 30. The paper notes this and shows the values dashed in Fig. 5(a). The subleading corrections are O(C_0/xi) ~ 2% at xi = 45, which is small, but the reader must consult Eq. (28) to infer this; a brief quantitative statement in Sec. 5.5 would help.
  3. Fig. 2, panel (b): The caption states the response 'decreases by less than 1% toward xi = 30' but the y-axis range makes this difficult to verify visually. An inset zooming in on the region xi = 15-30 would improve clarity.
  4. Sec. 3.1, Eq. (31): The expression for N_perp(h) is given in terms of h = 1/(2*xi), while most of the paper uses xi. A brief note that N_perp(h) = N_perp(1/(2*xi)) would aid readability.
  5. The paper is lengthy; Sec. 5.3 mixes laboratory (PVLAS, ATLAS) and astrophysical context somewhat diffusely. The ATLAS light-by-light scattering result, while interesting, is in a different kinematic regime and could be condensed.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for a careful and constructive report, and for explicitly crediting the internal consistency checks and falsifiable predictions. The referee's single major comment is well-taken and will be addressed in revision.

read point-by-point responses

  1. Referee: Sec. 5.5, Eqs. (75)-(76): The load-bearing identity q_hat + m_hat = 2*Delta_n + O(alpha^2) is not explicitly derived or verified numerically. An explicit verification (numerical check or two-line derivation) would substantially strengthen the paper.

    Authors: We agree with the referee that this identity is load-bearing for the atmosphere-independence claim and that it should be demonstrated explicitly rather than merely asserted. The identity follows directly from the definitions of the vacuum polarizability coefficients q̂ and m̂ in terms of the O(α) mode-index corrections. In the standard convention (Lai & Ho 2002, 2003; van Adelsberg & Lai 2006), the vacuum contributions to the dielectric tensor at θ=π/2 are expressed through the same second derivatives of the Heisenberg–Euler Lagrangian that appear in our refractive indices: the parallel mode is controlled by γ_GG and the perpendicular mode by γ_FF. The quantities q̂ and m̂ are linear combinations of these same one-loop derivatives, and their sum q̂ + m̂ reproduces twice the birefringent splitting 2Δn = 2(n_∥ − n_⊥) at O(α), because the Maxwell normalization γ_s cancels identically in the difference. In the weak-field limit this reduces to the trivial check q̂ + m̂ = 3δ_V = 2Δn_wf already noted in the manuscript. At finite field, the cancellation persists because the same γ_s factor divides both mode indices. We will add a brief derivation (approximately half a page) to Sec. 5.5 showing this cancellation explicitly at the level of the definitions, together with a numerical verification table at several representative ξ values (e.g., ξ = 1, 5, 10, 20, 30) confirming that q̂ + m̂ − 2Δn_exact vanishes at the O(α²) level (i.e., at the ~10⁻⁴ relative level set by α/π). This directly addresses the referee's concern and makes the atmosphere-independence of the ratios in Eq. (76) transparent. revision: yes

Circularity Check

0 steps flagged

No significant circularity found; derivation is self-contained against external benchmarks

full rationale

The paper's derivation chain is largely self-contained and does not exhibit circularity. The one-loop Heisenberg-Euler Lagrangian (Sec. 2) is sourced from standard external references (Schwinger 1951 [10], Heisenberg & Euler 1936 [8]). The refractive indices (Sec. 3) are derived from field derivatives of this Lagrangian, with weak-field limits checked against known QED coefficients (Eqs. 24-26) and strong-field asymptotes validated numerically. The photon magnetic response (Sec. 4) is defined as μ_γ = -∂ω/∂B, a derivative of the refractive index — a definition the paper is transparent about, not a prediction. The central observable results in Sec. 5.5 rest on the identity q̂ + m̂ = 2Δn + O(α²), which the paper states holds because q̂ and m̂ are defined from the same O(α) dielectric tensor eigenvalues whose difference is Δn. This is a mathematical identity at the working order, not a fit or a definition of q̂+m̂ in terms of Δn; it reduces the finite-field correction to the ratio R = Δn_wf/Δn_exact, which is computed from the exact one-loop special-function expressions (Eqs. 12-14), not fitted to data. The atmosphere-independent ratios (Eq. 76) follow because both Δn_wf and Δn_exact carry the same sin²θ factor, which cancels — a mathematical fact, not circularity. The polarization-limiting radius null result (Sec. 5.4) is geometric: freeze-out occurs at ~10² R_NS where B << B_cr, so exact and weak-field Δn coincide by construction of the weak-field limit. No parameters are fitted to data and then 'predicted.' The magnetar field values come from the external catalog Olausen & Kaspi (2014) [51]. The vacuum-resonance formalism follows Lai & Ho (2002, 2003) [59, 60], external references. The self-citations to Valluri et al. ([20, 21, 22, 33]) are for historical context (higher harmonics, paramagnetic properties) and are not load-bearing for the central derivation. The maximum at ξ ≈ 17 (Sec. 5.2) is acknowledged as a structural feature of the resummed one-loop model, not order-controlled in α, and the paper is explicit about this limitation. Score 1 reflects the presence of non-load-bearing self-citations in a derivation that is otherwise self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 5 axioms · 0 invented entities

The paper introduces no new particles, forces, dimensions, or postulated entities. The 'photon magnetic response' is explicitly stated (Sec. 4.3) to be an effective property of the photon dispersion relation induced by vacuum polarization, not an intrinsic magnetic moment of a free photon. No free parameters are fitted to data; all numerical results follow from the exact one-loop special-function expressions with physical constants (alpha, m_e, e, B_cr) as inputs. The magnetar field strengths are taken from the McGill Online Magnetar Catalog (Olausen & Kaspi 2014).

axioms (5)
  • standard math One-loop Heisenberg-Euler effective Lagrangian in constant external electromagnetic fields (Eq. 1) provides the correct effective-field-theory description for low-energy photon propagation in slowly varying fields.
    Invoked in Sec. 2.1; standard since Heisenberg & Euler (1936) and Schwinger (1951). The low-energy condition |k|/m << 1 is stated in Eq. (59).
  • domain assumption The background magnetic field is constant or slowly varying on the electron Compton wavelength scale.
    Stated in Sec. 2.1 and Sec. 4.3 (Eq. 59). This is the standard constant-field approximation of the Heisenberg-Euler framework. For magnetar magnetospheres, field variation scales are much larger than the Compton wavelength.
  • domain assumption The magnetar magnetic field is well-approximated by a centered dipole with radial propagation at theta=pi/2 for the polarization-limiting radius calculation.
    Stated in Sec. 5.4. The paper acknowledges this is a simplification and notes that multipolar surface fields would move r_cr but not the freeze-out region.
  • domain assumption The identity q_hat + m_hat = 2*Delta_n + O(alpha^2) holds at the working order O(alpha) for the vacuum-resonance analysis.
    Invoked in Sec. 5.5. The paper states this holds identically in the weak-field limit and at finite field because q_hat and m_hat are defined from the same O(alpha) mode-index corrections whose difference is Delta_n. The O(alpha^2) correction is not quantified.
  • domain assumption The dielectric tensor for the polarization-limiting radius calculation is purely the vacuum one (no plasma contribution).
    Stated in Sec. 5.4 caveats. The paper notes this means vacuum-resonance mode conversion is not modeled in that subsection but is addressed separately in Sec. 5.5.

pith-pipeline@v1.1.0-glm · 30266 in / 3549 out tokens · 498895 ms · 2026-07-08T06:06:29.410022+00:00 · methodology

0 comments
read the original abstract

{We study low-energy photon propagation in a constant magnetic field within the one-loop Heisenberg--Euler theory, retaining the refractive-index normalization $\gamma_s$ without expansion. Here ``finite-field'' denotes exact dependence on $B/B_{\rm cr}$ within the one-loop, constant-field approximation. The resulting birefringence is propagated into magnetar polarization transport. In a centered-dipole model, the polarization-limiting radius is unchanged to better than $10^{-12}$ because mode decoupling occurs at $\sim10^2R_{\rm NS}$, where $B\ll B_{\rm cr}$. Near the surface, however, the weak-field Cotton--Mouton expression overestimates the accumulated birefringent phase by up to a factor $2.9$ at $10^{15}$~G. At the plasma--vacuum resonance, finite-field corrections reduce the resonance density by $32\%$ and raise the adiabatic conversion energy by $14\%$ for 1E~1547.0$-$5408; the corresponding changes are factors $2.6$ and $1.37$ for 1RXS~J1708$-$4009, and factors $9.7$ and $2.13$ for SGR~1806$-$20, the latter controlled by the strong-field asymptote. The resummed one-loop parallel-mode magnetic response remains positive and develops a broad maximum near $17B_{\rm cr}$. The strictly truncated $\mathcal O(\alpha)$ response is monotonic; therefore the maximum is a structural prediction of the resummed one-loop constitutive model, while its detailed profile and precise location require higher-loop validation. These results identify vacuum-resonance observables as the most sensitive channel for testing finite-field QED in magnetars.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

62 extracted references · 62 canonical work pages · 4 internal anchors

  1. [1]

    Born and L

    M. Born and L. Infeld, Proc. R. Soc. Lond. A144, 425 (1934)

    work page 1934

  2. [2]

    G. V. Dunne, Int. J. Mod. Phys. A27, 1260004 (2012). 23

    work page 2012

  3. [3]

    Baier and P

    R. Baier and P. Breitenlohner, Nuovo Cimento B47, 117 (1967)

    work page 1967

  4. [4]

    Dittrich and H

    W. Dittrich and H. Gies,Probing the Quan- tum Vacuum(Springer, Berlin, 2000)

    work page 2000

  5. [5]

    A. E. Shabad and V. V. Usov, Phys. Rev. D 83, 105006 (2011)

    work page 2011

  6. [6]

    Villalba-Ch´ avez and A

    S. Villalba-Ch´ avez and A. E. Shabad, Phys. Rev. D86, 105040 (2012)

    work page 2012

  7. [7]

    Altschul, Astropart

    B. Altschul, Astropart. Phys.29, 290 (2008)

    work page 2008

  8. [8]

    Heisenberg and H

    W. Heisenberg and H. Euler, Z. Phys.98, 714 (1936)

    work page 1936

  9. [9]

    McKeon, Can

    G. McKeon, Can. J. Phys.57, 615 (1979)

    work page 1979

  10. [10]

    Schwinger, Phys

    J. Schwinger, Phys. Rev.82, 664 (1951)

    work page 1951

  11. [11]

    G. V. Dunne, arXiv:hep-th/0406216 (2004)

    work page internal anchor Pith review Pith/arXiv arXiv 2004

  12. [12]

    Bialynicka-Birula and I

    Z. Bialynicka-Birula and I. Bialynicki-Birula, Phys. Rev. D2, 2341 (1970)

    work page 1970

  13. [13]

    S. L. Adler, Ann. Phys. (N.Y.)67, 599 (1971)

    work page 1971

  14. [14]

    V. I. Ritus, Zh. Eksp. Teor. Fiz.69, 1517 (1975) [Sov. Phys. JETP42, 774 (1975)]

    work page 1975

  15. [15]

    Gies and F

    H. Gies and F. Karbstein, J. High Energy Phys.03, 108 (2017)

    work page 2017

  16. [16]

    Karbstein, Phys

    F. Karbstein, Phys. Rev. Lett.122, 211602 (2019)

    work page 2019

  17. [17]

    J. S. Heyl and L. Hernquist, J. Phys. A30, 6485 (1997)

    work page 1997

  18. [18]

    J. S. Heyl and L. Hernquist, Phys. Rev. D55, 2449 (1997)

    work page 1997

  19. [19]

    J. S. Heyl and L. Hernquist, J. Phys. A30, 6475 (1997)

    work page 1997

  20. [20]

    Bhartia and S

    P. Bhartia and S. R. Valluri, Can. J. Phys. 56, 1122 (1978)

    work page 1978

  21. [21]

    S. R. Valluri and P. Bhartia, Can. J. Phys. 58, 116 (1980)

    work page 1980

  22. [22]

    W. J. Mielniczuk, D. R. Lamm, and S. R. Val- luri, Can. J. Phys.66, 692 (1988)

    work page 1988

  23. [23]

    Zavattiniet al., Phys

    E. Zavattiniet al., Phys. Rev. D77, 032006 (2008)

    work page 2008

  24. [24]

    Bregantet al.(PVLAS Collaboration), Phys

    M. Bregantet al.(PVLAS Collaboration), Phys. Rev. D78, 032006 (2008)

    work page 2008

  25. [25]

    Cantatore (PVLAS Collaboration), Lect

    G. Cantatore (PVLAS Collaboration), Lect. Notes Phys.741, 157 (2008)

    work page 2008

  26. [26]

    Della Valleet al., Eur

    F. Della Valleet al., Eur. Phys. J. C76, 24 (2016)

    work page 2016

  27. [27]

    Ejlliet al., Phys

    A. Ejlliet al., Phys. Rep.871, 1 (2020)

    work page 2020

  28. [28]

    Ahmadiniazet al., High Power Laser Sci

    N. Ahmadiniazet al., High Power Laser Sci. Eng.13, e7 (2025)

    work page 2025

  29. [29]

    Aaboudet al., Nat

    ATLAS Collaboration, M. Aaboudet al., Nat. Phys.13, 852 (2017)

    work page 2017

  30. [30]

    Aadet al., Phys

    ATLAS Collaboration, G. Aadet al., Phys. Rev. Lett.123, 052001 (2019)

    work page 2019

  31. [31]

    A. I. Ibrahimet al., Astrophys. J. Lett.574, L51 (2002)

    work page 2002

  32. [32]

    R. P. Mignaniet al., Mon. Not. R. Astron. Soc.465, 492 (2017)

    work page 2017

  33. [33]

    S. R. Valluriet al., Mon. Not. R. Astron. Soc. 472, 2398 (2017)

    work page 2017

  34. [34]

    Tavernaet al.(IXPE Collaboration), Sci- ence378, 646 (2022)

    R. Tavernaet al.(IXPE Collaboration), Sci- ence378, 646 (2022)

    work page 2022

  35. [35]

    Lai, Proc

    D. Lai, Proc. Natl. Acad. Sci. USA120, e2216534120 (2023)

    work page 2023

  36. [36]

    R. E. Stewartet al., arXiv:2509.19446 (2025)

    work page arXiv 2025

  37. [37]

    H. P. Rojas and E. R. Querts, inProc. MG10 Meeting, edited by M. Novelloet al.(World Scientific, 2006), p. 2241

    work page 2006

  38. [38]

    H. P. Rojas and E. R. Querts, Int. J. Mod. Phys. D16, 165 (2007)

    work page 2007

  39. [39]

    P´ erez Rojas and E

    H. P´ erez Rojas and E. Rodr´ ıguez Querts, Eur. Phys. J. C74, 2899 (2014). 24

    work page 2014

  40. [40]

    Karbstein and R

    F. Karbstein and R. Shaisultanov, Phys. Rev. D91, 085027 (2015)

    work page 2015

  41. [41]

    Lundin, Europhys

    J. Lundin, Europhys. Lett.87, 31001 (2009)

    work page 2009

  42. [42]

    Lundin, Ph.D

    J. Lundin, Ph.D. thesis, Ume˚ a Universitet (2010)

    work page 2010

  43. [43]

    V. S. Adamchik, Comput. Phys. Commun. 157, 181 (2004)

    work page 2004

  44. [44]

    Dittrich, Phys

    W. Dittrich, Phys. Rev. D19, 2385 (1979)

    work page 1979

  45. [45]

    F. W. Olver, D. W. Lozier, R. F. Boisvert, and C. W. Clark,NIST Handbook of Math- ematical Functions(Cambridge Univ. Press, 2010)

    work page 2010

  46. [46]

    Tsai and T

    W.-y. Tsai and T. Erber, Phys. Rev. D12, 1132 (1975)

    work page 1975

  47. [47]

    C. M. Kim and S. P. Kim, arXiv:2202.05477 (2022)

    work page internal anchor Pith review Pith/arXiv arXiv 2022

  48. [48]

    Zaneet al., Astrophys

    S. Zaneet al., Astrophys. J. Lett.944, L27 (2023)

    work page 2023

  49. [49]

    The long quest for vacuum birefringence in magnetars: 1E 1547.0-5408 and the elusive smoking gun

    R. Tavernaet al., Astrophys. J. (2026), arXiv:2601.15452

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  50. [50]

    Detectability of Magnetar-Induced Vacuum Birefringence with IXPE and eXTP

    F. Abu-Ajamieh, arXiv:2605.02541 (2026)

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  51. [51]

    S. A. Olausen and V. M. Kaspi, Astrophys. J. Suppl.212, 6 (2014)

    work page 2014

  52. [52]

    Berceau, R

    P. Berceau, R. Battesti, M. Fouch´ e, and C. Rizzo, Can. J. Phys.89, 153 (2010)

    work page 2010

  53. [53]

    Cad` ene, P

    A. Cad` ene, P. Berceau, M. Fouch´ e, R. Bat- testi, and C. Rizzo, Eur. Phys. J. D68, 16 (2014)

    work page 2014

  54. [54]

    Eng.12, e26 (2024)

    HIBEF Collaboration, High Power Laser Sci. Eng.12, e26 (2024)

    work page 2024

  55. [55]

    J. S. Heyl and N. J. Shaviv, Mon. Not. R. Astron. Soc.311, 555 (2000)

    work page 2000

  56. [56]

    J. S. Heyl and N. J. Shaviv, Phys. Rev. D66, 023002 (2002)

    work page 2002

  57. [57]

    J. S. Heyl, N. J. Shaviv, and D. Lloyd, Mon. Not. R. Astron. Soc.342, 134 (2003)

    work page 2003

  58. [58]

    J. S. Heyl and I. Caiazzo, Galaxies6, 76 (2018)

    work page 2018

  59. [59]

    Lai and W

    D. Lai and W. C. G. Ho, Astrophys. J.566, 373 (2002)

    work page 2002

  60. [60]

    Lai and W

    D. Lai and W. C. G. Ho, Astrophys. J.588, 962 (2003)

    work page 2003

  61. [61]

    van Adelsberg and D

    M. van Adelsberg and D. Lai, Mon. Not. R. Astron. Soc.373, 1495 (2006)

    work page 2006

  62. [62]

    A. Y. Potekhin, D. Lai, G. Chabrier, and W. C. G. Ho, Astrophys. J.612, 1034 (2004). 25

    work page 2004