Strongly Magnetized Sources: QED and X-ray Polarization

Radiative corrections of quantum electrodynamics cause a vacuum threaded by magnetic field to be birefringent. This means that radiation of different polarizations travels at different speeds. Even in the strong magnetic fields of astrophysical sources the difference in speed is small; however, it has profound consequences for the extent of polarization expected from strongly magnetized sources. We demonstrate how the birefringence arises from first principles, show how birefringence affects the polarization state of radiation and present recent calculations for the expected polarization from magnetars and X-ray pulsars.


Introduction
The second prediction of quantum electrodynamics (QED) was the birefringence of the vacuum [1,2], in particular that light travels through a magnetic field at different speeds depending on its polarization. This came only eight years after the formulation of relativistic quantum mechanics [3] and five years after the first prediction of QED, the existence of the positron [4]. Of course, the positron was discovered soon after [5], but vacuum birefringence has not yet been observed definitively, although recently astronomers have discovered strong evidence for it [6]. During the next fifteen years, as QED was formalized and place on a rigorous foundation, the prediction of vacuum birefringence was also made more rigorous [7].
Why has it taken so long to observe directly the second prediction of QED? The primary reason is that, for magnetic field strengths achieved in the laboratory, the effect is incredibly small. At one Tesla, the difference in the index of refraction between the two modes is 4 × 10 −24 . Despite the incredible smallness of the effect, terrestrial experiments are reaching closer and closer to measuring it [8,9] in strong magnetic fields. The development of more powerful lasers promises to probe the closely related effect of photon-photon scattering [10,11]. Even in the magnetic fields of the most strongly magnetized objects in the Universe, the magnetars, the difference in the index of refraction is only a few percent; however, the vacuum birefringence can increase the observed extent of polarization in the X-rays by a factor of ten even for more weakly magnetized objects.
Although the difference in the index of refraction is small between the modes, this difference means that, for radiation of sufficiently high frequency, the polarization states will not mix as they travel through the magnetosphere of the source [12]. If the emission originates from a region over which the local magnetic field varies in direction, the decoupling of the polarization modes by vacuum polarization makes the final polarization fraction up to ten times larger than without it. This was first discovered by [13] and fueled a resurgence in the interest of using astrophysics to test QED.
Rather than provide a detailed review of all of the possible QED effects, we will focus on motivating why QED is important for X-ray polarization from magnetized objects. We will start from first principles, leading to some recent results for magnetars and X-ray pulsars that can be probed with instruments currently under development [14][15][16]. We will leave out other interesting phenomena, such as the role that the competition between vacuum and plasma birefringence plays in the formation of spectral lines [17]. In particular, we will derive the effective Lagrangian of QED from the classical Lagrangian of the electromagnetic field coupled to a Dirac field using modern functional techniques ( § 2.1), we will calculate the index of refraction ( § 2.2) and how it affects the propagation of polarization radiation ( § 2.3), and finally we will present recent results for polarized emission from magnetars and X-ray pulsars ( § 3).

Theoretical Treatment : Effective Action and Vacuum Birefringence
If an external field breaks the symmetry of the vacuum, we must develop new tools to account for its presence. The Lagrangian of QED is where the interaction between the electromagnetic fields and the external field is given by the This rule must be taken into account in all fermion propagators, including internal lines such as in the vacuum polarization and in photon splitting processes. The symbols γ µ are the Dirac matrices that span the spinorial components of the fermion fields. If the external field is sufficiently weak, the interaction with the field may be treated perturbatively as a series of discrete interactions. On the other hand, at a field strength of B QED = m 2 c 3 /(eh) = 4.4 × 10 13 G, the gyration energy of an electron or, equivalently, the potential energy drop across its Compton wavelength, is equal to its rest mass. For this reason, when the field exceeds a critical value of approximately B QED /2, this series fails to converge. Essentially, each term in the sum of diagrams is equally large in this limit.
In the next section, we derive the effective action of a general field configuration to one-loop order using the QED Lagrangian (Eq. 1) and techniques from statistical mechanics. The key results for X-ray polarization are the index of refraction ( § 2.2) and how polarization changes as radiation propagates through an inhomogeneous birefringent medium ( § 2.3).

Effective Action : Formal Derivation
The connections between the theory of quantum and statistical fields are manifold. Our derivation of the effective action and Lagrangian will exploit these connections. The final results that we present here are well known in the specialized literature for quantum field theory in strong fields [e.g 18,19] and older introductory texts [e.g. 20,21], but they are typically absent from recent introductory texts [e.g. 22,23]. We present a derivation of the effective action using functional techniques familiar from modern treatments of quantum field theory [e.g. 23] and statistical mechanics.
First, the partition function consists of the sum of the quantum phases (or statistical weights) of each possible state of the system (the analysis in this section draws on § §11.3-11.4 of [23]), where J µ and A µ are the electromagnetic current and vector potential respectively, η andη are the fermionic currents and ψ andψ are the fermionic fields. The variables η,η, ψ andψ are anti-commuting Grassman numbers. This means that ηψ = −ψη and η 2 = 0. Grassman numbers behave differently from commuting numbers under integration and differentiation as well. Brackets are used to denote functionals (integrals of the fields over the entire spacetime) while parentheses indicate functions. Each field configuration receives a phase proportional to the integral of the Lagrangian over spacetime, i.e. the action. The constant of proportionality is i/h; in statistical physics, the constant of proportionality is −1/(kT), and the states are weighted by energy and not action. Drawing the analogy further, the functional E[J µ ,η, η], would correspond to the Helmholtz free energy F(T, V, N) in statistical mechanics; it is the vacuum energy as a function of the external sources J µ ,η and η. The partial derivative F(T, V, N) with respect to the volume of the system yields the thermodynamic conjugate to V, the pressure. For the quantum field, the result proceeds similarly, where we use the symbol δ to denote a functional derivative and where |Ω denotes the vacuum state.
The functional derivative of E[ ] with respect to one of the currents yields the classical field, i.e. the vacuum expectation value of the corresponding field, which we denote as A 0 µ (x). Generally, when one considers the properties of the magnetized vacuum, it is the fields that are specified, not the currents. The effective action is related to E[ ] through a Legendre transformation, just as the Gibbs free energy G is related to F, i.e.
Analogously, the effective action is The functional derivative of the effective action, Γ[ ], with respect to one of the classical fields yields the distribution of the corresponding current. Using the analogy with thermodynamics, the effective action is the vacuum energy with the distribution of the fields fixed.

Functional integration
Computing the effective action begins with the expression for Z[ ], the partition function, specifically by expanding the classical action with currents about the values of the classical fields, (∆ψ(y)) + (∆A µ (x)) δ 2 L δA µ (x)δψ(y) (∆ψ(y)) + (∆ψ(x)) δ 2 L δψ(x)δA µ (y) (∆A µ (y)) + Higher Order Terms (8) where , the difference between the electromagnetic field including the quantum fluctuations and the classical electromagnetic field, and similarly for the other fields. Since the functional derivatives will be evaluated at ψ 0 (y) =ψ 0 (y) = 0, the last two terms vanish. Furthermore, second derivatives with respect to the same Grassman field also vanish because of the anti-commutative nature of the fields. Let's evaluate for an example the term The ← ∂ notation indicates that the resulting functional derivative is an operator that differentiates something to the left. Because the function derivative lives within a integral over all of spacetime, we can use integration by parts to simplify the result further if we assume that the boundary terms vanish: We see that the first-order derivatives vanish when the fields satisfy the field equations. Although we do not have an explicit relationship that connects the currents to the classical fields that they generate, we will impose that the currents J µ ,η and η along with the classical fields A 0 µ ,ψ 0 and ψ 0 satisfy the field equations and evaluate all of the functional derivatives at the values of the classical fields, so the first order terms in the expansion vanish and the vacuum energy E[ ] to lowest order is a Gaussian functional integral, We integrated over all of the possible field configurations by assuming that a particular field configuration is a point in an infinite dimensional space and by looking at the functional derivatives as infinite dimensional matrices. Let us look at the following term in detail to see how this works: where λ i are the eigenvalues of the matrix δ 2 L/(δA µ (x i )δA ν (y j )) , and we have chosen the eigenvectors as a basis for performing the integral over the field configurations, ∏ dA µ (x i ) . Now let us perform the integration over each of the dA µ (x i ) to yield The constant terms absorb the constant prefactor in front of the integral, −i/(2h), as well as some divergent terms. The symbol Det denotes the functional determinant over both the spacetime and the spin space in the case of the Dirac fields. One subtlety is the plus sign in front of the functional derivative involving the Grassman fieldsψ(x) and ψ(y). Simply, the integral of where the unexpected result comes from the anti-commuting nature of the fields. Performing the Legendre transformation yields an expression for the effective action to lowest order (one loop), where we have dropped the constant terms from the expression. We have not been concerned with renormalizing the effective action so far, but the functional determinants are probably divergent. Let us insist that the effective action vanishes as the classical field vanishes, so we have to subtract two terms corresponding to the functional determinants in the absence of an external field. This renormalizes the zero-point energy and yields, where The functional derivative of the Lagrangian with respect to the vector potential is same for all values of the classical vector potential as along as the fermionic classical field vanishes. The effective action contains the classical Maxwell action of electrodynamics and an additional term that quantifies the effects of the vacuum fluctuations of the Dirac (here electron-positron) fields. We can use the linear algebra result, ln DetA = Tr ln A, to simplify the expression for the effective action further. We use the convention that Det and Tr span both coordinate and spin space, while tr and det just cover the spinorial components. (the analysis is this subsection builds upon §4-3-3 and §4-3-4 of [20]), We would like to put the logarithm in a more manageable form.

Tr ln
The first equality holds since the charge conjugation matrix C satisfies Cγ µ C −1 = −γ T µ , so CΠ /C −1 = −Π / T (similarly for p /), and the trace of an operator is invariant under transposition. The second equality results from summing the first two expressions.

Effective Action : Proper-time Integration
We use the identity to expand the logarithm (29) and obtain the proper-time expression for the effective Lagrangian density [7], where U(s) is the time-evolution operator governed by the Hamiltonian, where is the analogous operator for vanishing external fields. Eq. 30 forms the basis of the worldline numerics techinque that facilitates the calculation of the effective action for arbitrary field configurations [24,25]

Results for a uniform field
We will select a particular frame and gauge to calculate the trace and obtain an expression for the effective Lagrangian from a uniform electromagnetic field. We begin with since σ µν F 0 µν commutes with Π µ Π µ for constant fields.
We will choose a frame such that E B, and a ≡ |E| and b ≡ |B| where E and B are the classical electric and magnetic fields respectively. In this frame, the eigenvalues of i 2 esσ µν F 0 µν are ±es(a ± ib). In a general frame, we can define which demonstrates that a and b are Lorentz invariants of the field. We obtain Evaluation of the first trace is more complicated. With no loss of generality we can assume that the magnetic and electric fields point in the z-direction and select a gauge with A 3 = −at and A 1 = −by; this yields, Using the commutation relation [x, p x ] = −i, we can define the shift operator which allows us to write To evalulate the trace itself we use the momemtum representation tr x| exp(isΠ µ Π µ )|x = dp 3 dp 1 (2π) 4 dp 0 dp 0 dp 2 dp 2 exp i( Let us examine the last of the two integrals in detail.
where H is the Hamiltonian of a harmonic oscillator with unit mass and spring constant k = e 2 b 2 .
Using the known eigenvalues of the system yields an expression for the integral, The result for the first integral is similar except here k = −e 2 a 2 , so the complete expression for tr x| exp(isΠ / 2 )|x = −i e 2 ab (2π) 2 coth(eas) cot(ebs).
Taking the limit of this expression as a and b vanish yields

Effective Lagrangian in a constant field
Substituting this result into Eq. 30 yields an expression for the effective Lagrangian, For small values of s, the integrand diverges as e 2 (a 2 − b 2 )/(3s). Since this is proportional to the classical Lagrangian, we can absorb this infinite term through a renormalization which yields the renormalized expression for the effective Lagrangian, We can get a more convenient expression by performing a Wick rotation, and substituting ζ = sm 2 where α = e 2 /(hc) and we can now take → 0.

Index of refraction for Low-Energy Photons
The index of refraction for low-energy photons is obtained most simply by defining the macroscopic fields as the generalized momenta conjugate to the fields [21], and linearizing these relations about the background field [26]. For an external magnetic field this yields [27] where ξ = B/B QED , n ⊥ and n are the index if refraction for the perpendicular and parallel modes respectively (see below), and ln A = 1 12 − ζ (1) (−1) ≈ 0.248754477. The functions X 0 (1/ξ) and X 1 (1/ξ) are related to the effective action and its derivative with respect to a in the limit of a → 0 [28].
Our naming convention is the following: if µναβ F µν F αβ = 0 or E · B = 0, the photon is in the parallel mode, otherwise it is in the perpendicular mode, where F µν is the sum of field tensors of the wave and external field. In the weak field limit, n − 1 ∝ ξ 2 ; while in the strong field limit n ⊥ − 1 ∝ ξ and n approaches a constant [27]. In particular, in the weak field we have and a birefringence of

Propagation through a Birefringent Medium
As polarized radiation propagates through a birefringent medium, the direction of polarization changes. If we use the Stokes parameters S 0 , S 1 , S 2 and S 3 (I, Q, U, V), we can define the normalized Stokes vector s = (S 1 , S 2 , S 3 )/S 0 and follow the evolution of the polarization using [29] and where λ measures the length of the photon path and ν is the frequency of the radiation. The value of ∆k is the difference in the wavenumber for the two polarization states, and the final equality holds in the weak-field limit of QED. The direction ofΩ points toward the polarization of the perpendicular mode on the Poincaré sphere of polarization states. This equation may seem more familiar if one considers the Faraday rotation of polarized light passing through a weakly magnetized plasma. In this case,Ω points toward the s 3 -direction, corresponding to the circular polarization, so the polarization direction of linearly polarized light will rotate. In general, if the direction ofΩ is constant, the vector s will circle the direction ofΩ. If the magnetic field is sufficiently strong so that |Ω| is sufficiently large, the vector s will circle the direction ofΩ even in the case in whichΩ changes direction and magnitude, if it does so sufficiently gradually. In particular, if the polarization state is initially parallel toΩ, that is, the initial polarization is parallel or perpendicular to the magnetic field, the polarization state will remain nearly parallel toΩ as long as [30] Ω d ln |Ω| dλ If this condition holds, the polarization states evolve adiabatically, and the polarization direction will follow the direction of the birefringence. If we specialize to the dipole field surrounding a neutron star where B ≈ µr −3 , where µ is the magnetic dipole moment of the star and r is the distance from the center of the star, we find the following condition where β is the angle between the dipole axis and the line of sight.
If we define the polarization-limiting radius (r PL ) to be the distance at which the equality holds, we find that the polarization will follow the direction of magnetic field out to  Fig. 1 illustrates the propagation of radiation away from the surface of the neutron star toward a distant observer. For X-ray photons coming from near the surface of a neutron star with a surface field of 10 12 G, the polarization-limiting radius is much larger than the star, according to Eq. 57, so the observed polarization of the photons will reflect the direction of the magnetic field at a Radiation leaving the surface of a neutron star follows geodesics so that the bundle of rays that reaches the distant observer is approximately cylindrical. The three-dimensional coordinates (x, y, z) are given in terms of the radius of the neutron star (R). If the polarization-limiting radius is small, the final polarization will reflect the magnetic field structure near to the star where the bundle covers a large fraction of hemisphere so the field structure varies a lot over the bundle at this point, and the final polarization will also vary a lot over the image. On the other hand, if the polarization-limiting radius is large, the field structure over the ray bundle is simpler, and the polarization direction will not vary much over the image. Since the magnetic field is assumed to be that of a dipole (aligned with the z-axis), it has axial symmetry and different images will be distinguishable by the observer's magnetic inclination angle i. Adapted from [31] large distance from the star and not at the surface. For a much more weakly magnetized star, the polarization-limiting radius will be comparable to the radius of the star, so the observed polarization will reflect the field structure close to the star. Fig. 2 depicts the final polarization states across the image of the neutron star surface assuming that the radiation is initially in the ordinary mode, that is, the electric field is parallel to the local magnetic field. The left panel shows the case where the vacuum birefringence is neglected, and the right panel show the case where the surface field is about 10 12 G and the frequency is 10 17 Hz or an energy of about 0.4 keV. This is appropriate for a thermally emitting neutron star such as one of the X-ray dim neutron stars (XDINS). The effect of the vacuum polarization is to comb the polarization direction to be aligned with the direction of the magnetic axis of the star and dramatically increase the observed total polarization from about 13% to about 70%. For more strongly magnetized neutron stars the effect is more dramatic.
The studies of the effects of QED on the net polarization from astrophysical sources are more mature for the magnetars and the XDINS. However, it may be important for other sources as well. Tab. 1 lists several classes of possible sources. For the magnetars and XDINS, the polarization-limiting radius is much larger than the radius of the star. Because we expect the emission in these sources to come from a larger region of the stellar surface or magnetosphere (in the case of the non-thermal emission from magnetars [32]), we expect a large increase in the observed polarization fraction due to QED. Although the ratio of the polarization-limiting radius to the stellar radius is also large for the X-ray pulsars (XRP), as we shall see in § 3, the effect for these objects is more subtle. The QED effects for more weakly magnetized stars such as millisecond XRPs (ms XRPs) and strongly magnetized For comparison, if one assumes that the entire surface is emitting fully polarized radiation, the net linear polarization on the left 13% while it is 70% on the right. Table 1. The expected polarization-limiting radii for various sources; the typical observing times are for eXTP at 2-8 keV for the magnetars (4U 0142+61) and XRPs (Her X-1) from the text and for RedSOX [34] at 0.2-0.8 keV for RX J1856.5-3754 to make a four-sigma detection. −3 ] r pl at 4 keV [cm] r pl /R t obs Magnetar 10 6 10 15 10 33 3.0 × 10 8 300 10 ks XDINS 10 6 10 13 10 31 4.7 × 10 7 50 * 1 ks XRP 10 6 10 12 10 30 1.9 × 10 7 20 100 ks ms XRP 10 6 10 9 10 27 1.2 × 10 6 1.2 AM Her 10 9 10 8 10 35 1.9 × 10 9 1.9 Black Hole 10 6+ ? N/A See [33] * XDINS have little emission at 4 keV, and therefore will be difficult to observe with eXTP and IXPE. The value of r pl /R at 0.4 keV is 32, so vacuum birefringence is important for observations with soft-X-ray polarimeters.

Radius [cm] Magnetic Field [G] µ [G cm
white dwarfs have not yet been explored. However, QED has been found to be generically important in the interpretation of X-ray polarization from accreting black holes [33, and Caiazzo & Heyl this volume].

Results
To demonstrate the importance of QED vacuum birefringence on the X-ray polarization from neutron stars, we focus on two particular objects: the magnetar 4U 0142+61 and the X-ray pulsar Hercules X-1. For the magnetar, the effect is especially dramatic. The magnetar 4U 0142+61 exhibits a X-ray spectrum from 2-8 keV (the range of sensitivity of IXPE and eXTP) that is dominated by thermal emission [35]. As argued by [e.g. [36][37][38][39][40], this thermal emission is expected to be nearly fully polarized as it is emitted. We have used the phase-resolved spectral fits of [35] to simulate the polarized spectra in XIMPOL [41] for a 10.0-ks observation with eXTP. Because we expect, when QED is included, that the final polarization will follow the magnetic field direction at a large distance from the star, the final polarization fraction with QED will be near 100% [36]. On the other hand without QED, the expected final polarization will be lower as apparent from the left panel of Fig. 2 and will depend on where the emission originates. Using the fits of [35] as a guide, we assume that the hot blackbody component originates from a hot-spot of ten degree radius and that the cool blackbody originates from the entire stellar surface. We assume that the hard-power-law component will not contribute in the 2-8 keV range. Even without QED birefringence, we expect the first component to remain highly polarized because it originates from a small region near the polar cap where the projected field is well aligned (see the left panel of Fig. 2), but the second component will not, because it comes from the rest of the surface. Fig. 3 presents the results of this calculation along with a short simulated observation with eXTP. Even a relatively short observation of this magnetar can detect the effect of QED vacuum birefringence dramatically. The emission from other magnetars such as AXP 1RXS J170849.0-400910 is dominated by non-thermal magnetospheric emission even in the 2-8 keV band of eXTP and IXPE so to understand the X-ray polarization from these objects requires a magnetospheric emission model [e.g. 32]. However, because the QED-induced polarization-limiting radius is much larger than the emission region even in this case, the results are similar: QED dramatically increases the expected polarization of these objects.
The second object that we consider is the X-ray pulsar Her X-1. We consider that the radiation comes from a slab at the polar cap with a radius of six degrees. In reality the X-ray emission comes from both the accretion disk and from near the star, and the contribution from the disk is expected to be more important at lower energies. Furthermore, even the structure of the emission near the star is uncertain: it could be a slab or an accretion column, and there are non-thermal contributions as well. For simplicity here we just treat the slab geometry and argue that with phase-resolved polarimetry,  x Figure 4. The polarization from Her X-1 as function of photon energy using the emission models of [12,42], averaged over the rotation of the pulsar. A positive value of the polarization degree indicates that the polarization direction is perpendicular to the projection of the rotation axis onto the plane of the sky. In the left panel, the dashed lines give the result without vacuum birefringence (we have reversed the polarization direction in this case for ease of comparison, see the text) and the solid lines include the QED effect. A positive value of Q/I indicates that the net polarization is perpendicular to the projected spin axis of the stars in the QED-on case, and parallel in the QED-off case. The crosses give the results of a 100-ks simulated observation with eXTP. The results for IXPE are similar with a 300-ks observation.
the geometric ambiguities and the location of various emission regions can be disentangled. We use the results of [42] to estimate the total flux from the region and the results of [29] to estimate the polarized fraction and the intensity as a function of direction from the slab. We choose a geometry that results in a large pulsed fraction. We again use XIMPOL [41] to perform the instrumental simulation including effects of QED birefringence by integrating Eq. 54 through the magnetosphere. We assume a range of possible stellar radii from 9 to 15 km. From Fig. 2, we see that when the emission is restricted to the region near the magnetic pole, the effect of QED is more subtle than for the magnetar. The projected magnetic field direction is well aligned near the pole even without vacuum birefringence in the left panel. The final polarization in the right panel is well aligned in general but substantial circular polarization can be generated near the pole [43]. Because the instrument only detects linear polarization, the circularly polarized radiation does not contribute to the observed polarization fraction, so the net effect of the vacuum birefringence is to reduce the polarization fraction integrated over the emission region and the rotational phase as shown by the left panel of Fig. 4. The solid lines trace the extent of linear polarization including QED, and the dashed lines neglect it. Across the 2-8-keV energy range, the trend in Q/I can be explained by looking at the opacities in the atmosphere. At low energies the opacity for photons in the parallel polarization is larger than for the perpendicular polarization, so the bulk of the radiation emerges in the perpendicular mode. However, as the photon energy approaches the cyclotron energy (here taken to be 38 keV), the opacity for the perpendicular mode increases and one gets more emission in the parallel mode. The right panel of Fig. 4 really highlights the role of QED. Below about 400 eV the polarization swings as well. All of the photons coming through the atmosphere pass through the vacuum resonance region in which the linear contribution to the birefringence from QED cancels that of the plasma [17]. In this region the value ofΩ swings from pointing along a particular direction in the s 1 − s 2 −plane up toward s 3 and back onto the s 1 − s 2 −plane in the opposite direction. As we know from Eq. 56 if this happens slowly enough the photon polarization will follow the direction ofΩ, and this is in fact what happens for photons with energies greater than about 350 eV. The polarization state is switched from perpendicular to parallel. Because this QED effect reverses the polarization direction for all of the photons above this energy, we have switched the sign of the value of Q/I for QED-off case for ease of comparison. Without an independent measurement of the projection of the spin axis of the star into the plane of the sky, measurements just in the 2-8-keV band cannot measure this polarization flip. The key effect of the QED birefringence in this harder band is to slightly reduce the polarization fraction.
We can apply Eq. 56 to determine how the conditions in the atmosphere, in particular the density scale height, determine the critical energy above which the polarization switches as the radiation passes through the vacuum resonance [17]. In the case of Her X-1 we estimate the temperature of the atmosphere to be 8 keV [42]. With the temperature fixed, the density scale height only depends on the composition of the atmosphere which is that of the donor star and the surface gravity; therefore, the photon energy at which the polarization direction flips depends on the surface gravity and can be used to measure the radius of the star. For the curves in Fig. 4 we have assumed a mass of 1.4 M for the neutron star.

Discussion
We have presented an ab initio derivation of the effective Lagrangian of quantum electrodynamics and the induced vacuum birefringence. We have found that QED vacuum birefringence dramatically increases the total polarization fraction of the the thermal emission from a magnetar. The bright, hot magnetar 4U 0142+61 is an excellent candidate to probe the QED vacuum with the current generation of X-ray polarimetry missions: IXPE and eXTP. We presented the first calculation of the expected polarization from an X-ray pulsar, Her X-1, to include vacuum birefringence in the magnetosphere. For these more weakly magnetized objects, the effect of QED is more subtle. In the IXPE and eXTP energy band, QED reduces the total polarization fraction slightly and the resulting extent of polarization has a modest dependence on the radius of the neutron star. If the emission at lower energies indeed comes from the stellar atmosphere, the direction of the polarization will swing by ninety degrees at about 400 keV. This occurs because higher energy photons can pass through the vacuum resonance in the stellar atmosphere adiabatically [17]. The transition energy between adiabatic and non-adiabatic evolution depends on the density scale height in the atmosphere, so a measurement of this energy could constrain the radius of the neutron star.

Materials and Methods
We integrate the equations of the polarization evolution through the neutron-star magnetosphere using an adaptive Runge-Kutta method as outlined in [13].