Data availability
NICER observations (ObsIDs: 8020300101, 8020300102, 8020300103 and 8020300104) and IXPE observations (04003801) are readily accessible in the HEASARC data archive: https://heasarc.gsfc.nasa.gov/W3Browse (https://doi.org/10.25504/FAIRsharing.979d22). Radio observations made by Murriyang/Parkes (ObsIDs: r070819_124237, uwl_250326_173709, uwl_250330_171906 and uwl_250331_131207; doi10.4225/08/52292AE9B2D80, doi10.25919/rzdr-pw25 and doi10.25919/v5hn-4v34) are publicly available from the CSIRO Data Access Portal (https://data.csiro.au/) following an 18-month proprietary period starting on the observation date.
Code availability
Data reduction and analysis of X-ray products were performed using publicly available software HEAsoft v.6.35.0 (https://heasarc.gsfc.nasa.gov/docs/software/lheasoft/) from the High Energy Astrophysics Science Archive Research Center (HEASARC), particularly FTOOLs v.6.35.1, SAOImage DS9 v.8.4b1 and Xspec v.12.15.0. Generation and calibration of the NICER event lists was also performed by NICERDAS v.12 of HEASoft. The simulation and analysis framework ixpeobssim v.31.1.0 was used to generate high-level IXPE data products (https://ixpeobssim.readthedocs.io/en/latest/). Moreover, the software filterbackground.py was used for the treatment of the IXPE background, found at GitHub (https://github.com/aledimarco/IXPE-background). Timing analysis was performed using tempo2 (https://github.com/mattpitkin/tempo2), PINT (https://github.com/nanograv/PINT) and CRIMP (https://github.com/georgeyounes/CRIMP/tree/main). PyXspecCorner (https://github.com/garciafederico/pyXspecCorner) and corner.py (https://corner.readthedocs.io/en/latest/) were used to generate the X-ray spectro-polarimetric and radio RVM corner plots, respectively. Additional custom code for generating figures and performing analysis is available at GitHub (https://github.com/rae-stewart/Polarimetric-Analysis-of-1E-1547.0-5408). Custom code for the MAGTHOMSCATT Monte Carlo simulation is available upon reasonable request.
References
Kouveliotou, C. et al. An X-ray pulsar with a superstrong magnetic field in the soft γ-ray repeater SGR1806 - 20. Nature 393, 235–237 (1998).
Article ADS CAS Google Scholar
Lai, D. & Ho, W. C. G. Polarized X-ray emission from magnetized neutron stars: signature of strong-field vacuum polarization. Phys. Rev. Lett. 91, 071101 (2003).
Article ADS PubMed Google Scholar
van Adelsberg, M. & Perna, R. Soft X-ray polarization in thermal magnetar emission. Mon. Not. R. Astron. Soc. 399, 1523–1533 (2009).
Article ADS Google Scholar
Taverna, R., Turolla, R., Suleimanov, V., Potekhin, A. Y. & Zane, S. X-ray spectra and polarization from magnetar candidates. Mon. Not. R. Astron. Soc. 492, 5057–5074 (2020).
Article ADS CAS Google Scholar
Heisenberg, W. & Euler, H. Folgerungen aus der Diracschen Theorie des Positrons. Z. Phys. 98, 714–732 (1936).
Article ADS CAS Google Scholar
Schwinger, J. On gauge invariance and vacuum polarization. Phys. Rev. 82, 664–679 (1951).
Article ADS MathSciNet Google Scholar
Duncan, R. C. & Thompson, C. Formation of very strongly magnetized neutron stars - Implications for gamma-ray bursts. Astrophys. J. Lett. 392, L9–L13 (1992).
Article ADS CAS Google Scholar
Paczynski, B. GB 790305 as a very strongly magnetized neutron star. Acta Astronaut. 42, 145–153 (1992).
Google Scholar
Kaspi, V. M. & Beloborodov, A. M. Magnetars. Annu. Rev. Astron. Astrophys. 55, 261–301 (2017).
Article ADS CAS Google Scholar
Harding, A. K. & Lai, D. Physics of strongly magnetized neutron stars. Rep. Prog. Phys. 69, 2631–2708 (2006).
Article ADS CAS Google Scholar
Taverna, R. & Turolla, R. X-ray polarization from magnetar sources. Galaxies 12, 6 (2024).
Article ADS Google Scholar
Adler, S. L. Photon splitting and photon dispersion in a strong magnetic field. Ann. Phys. 67, 599–647 (1971).
Article ADS Google Scholar
Tsai, W.-Y. & Erber, T. Propagation of photons in homogeneous magnetic fields: index of refraction. Phys. Rev. D 12, 1132–1137 (1975).
Article ADS Google Scholar
Pavlov, G. G. & Gnedin, Y. N. Vacuum polarization by a magnetic field and its astrophysical manifestations. Astrophys. Space Phys. Res. 3, 197 (1984).
ADS Google Scholar
Adam, J. et al. Measurement of e+e− momentum and angular distributions from linearly polarized photon collisions. Phys. Rev. Lett. 127, 052302 (2021).
Article ADS PubMed CAS Google Scholar
Ho, W. C. G., Lai, D., Potekhin, A. Y. & Chabrier, G. Atmospheres and spectra of strongly magnetized neutron stars. III. Partially ionized hydrogen models. Astrophys. J. 599, 1293–1301 (2003).
Article ADS CAS Google Scholar
Mignani, R. P. et al. Evidence for vacuum birefringence from the first optical-polarimetry measurement of the isolated neutron star RX J1856.5-3754. Mon. Not. R. Astron. Soc. 465, 492–500 (2017).
Article ADS CAS Google Scholar
Taverna, R. et al. Polarized X-rays from a magnetar. Science 378, 646–650 (2022).
Article ADS PubMed CAS Google Scholar
Ho, W. C. G. & Lai, D. Atmospheres and spectra of strongly magnetized neutron stars - II. The effect of vacuum polarization. Mon. Not. R. Astron. Soc. 338, 233–252 (2003).
Article ADS CAS Google Scholar
Ejlli, A. et al. The PVLAS experiment: a 25 year effort to measure vacuum magnetic birefringence. Phys. Rep. 871, 1–74 (2020).
Article ADS Google Scholar
Lai, D. IXPE detection of polarized X-rays from magnetars and photon mode conversion at QED vacuum resonance. Proc. Natl Acad. Sci. USA 120, e2216534120 (2023).
Article MathSciNet PubMed PubMed Central CAS Google Scholar
Heyl, J. S., Shaviv, N. J. & Lloyd, D. The high-energy polarization-limiting radius of neutron star magnetospheres - I. Slowly rotating neutron stars. Mon. Not. R. Astron. Soc. 342, 134–144 (2003).
Article ADS Google Scholar
Heyl, J. S. & Shaviv, N. J. Polarization evolution in strong magnetic fields. Mon. Not. R. Astron. Soc. 311, 555–564 (2000).
Article ADS Google Scholar
Camilo, F., Ransom, S. M., Halpern, J. P. & Reynolds, J. 1E 1547.0–5408: a radio-emitting magnetar with a rotation period of 2 seconds. Astrophys. J. 666, L93–L96 (2007).
Article ADS Google Scholar
Coti Zelati, F. et al. The long-term enhanced brightness of the magnetar 1E 1547.0–5408. Astron. Astrophys. 633, A31 (2020).
Article Google Scholar
Lower, M. E. et al. The 2022 high-energy outburst and radio disappearing act of the magnetar 1E 1547.0–5408. Astrophys. J. 945, 153 (2023).
Article ADS Google Scholar
Bildsten, L. et al. Observations of accreting pulsars. Astrophys. J. Suppl. Ser. 113, 367–408 (1997).
Article ADS Google Scholar
Camilo, F., Reynolds, J., Johnston, S., Halpern, J. P. & Ransom, S. M. The magnetar 1E 1547.0–5408: radio spectrum, polarimetry, and timing. Astrophys. J. 679, 681–686 (2008).
Article ADS CAS Google Scholar
Zeng, S., Philippov, A., Juno, J., Beloborodov, A. M. & Popova, E. Origin of pulsed radio emission from magnetars. Astrophys. J. Lett. 996, L20 (2026).
Article ADS CAS Google Scholar
Stewart, R. et al. X-ray polarization of the magnetar 1E 1841-045. Astrophys. J. Lett. 985, L35 (2025).
Article ADS Google Scholar
Rigoselli, M. et al. IXPE detection of highly polarized X-rays from the magnetar 1E 1841-045. Astrophys. J. Lett. 985, L34 (2025).
Article ADS Google Scholar
Tiengo, A. et al. The dust-scattering X-ray Rings of the anomalous X-ray pulsar 1E 1547.0–5408. Astrophys. J. Lett. 710, 227–235 (2010).
Article Google Scholar
Heyl, J. et al. The detection of polarized X-ray emission from the magnetar 1E 2259+586. Mon. Not. R. Astron. Soc. 527, 12219–12231 (2024).
Article ADS Google Scholar
Medin, Z. & Lai, D. Condensed surfaces of magnetic neutron stars, thermal surface emission, and particle acceleration above pulsar polar caps. Mon. Not. R. Astron. Soc. 382, 1833–1852 (2007).
Article ADS CAS Google Scholar
Dinh Thi, H., Baring, M. G., Hu, K., Harding, A. K. & Barchas, J. A. Monte Carlo simulations of polarized radiative transfer in neutron star atmospheres. Astrophys. J. 992, 188 (2025).
Article ADS Google Scholar
Radhakrishnan, V. & Cooke, D. J. Magnetic poles and the polarization structure of pulsar radiation. Astrophys. Lett. 3, 225 (1969).
ADS Google Scholar
Li, B.-P., Gao, Z.-F., Ma, W.-Q. & Zhang, W.-F. Probing the origin of magnetar X-ray polarization diversity: a multi-wavelength geometrical study of 1E 1547.0-5408 and 1E 2259+586. Preprint at arxiv.org/abs/2604.10477 (2026).
Barchas, J. A., Hu, K. & Baring, M. G. Polarized Radiation Transfer in Neutron Star Surface Layers. Mon. Not. R. Astron. Soc. 500, 5369–5392 (2021).
Article ADS CAS Google Scholar
Hu, K., Baring, M. G., Barchas, J. A. & Younes, G. Intensity and polarization characteristics of extended neutron star surface regions. Astrophys. J. 928, 82 (2022).
Article ADS Google Scholar
Dinh Thi, H. et al. Pulsed, polarized X-ray emission from neutron star surfaces: the effects of vacuum birefringence in the magnetosphere. Astrophys. J. 1000, 73 (2026).
Article ADS CAS Google Scholar
Marshall, H. L. et al. The globe orbiting soft X-ray (GOSoX) polarimeter concept study. In Proc. SPIE 11822, Optics for EUV, X-Ray, and Gamma-Ray Astronomy X Vol. 11822 (SPIE, 2021).
Wilms, J., Allen, A. & McCray, R. On the absorption of X-rays in the interstellar medium. Astrophys. J. 542, 914–924 (2000).
Article ADS CAS Google Scholar
Verner, D. A., Ferland, G. J., Korista, K. T. & Yakovlev, D. G. Atomic data for astrophysics. II. New analytic FITS for photoionization cross sections of atoms and ions. Astrophys. J. 465, 487 (1996).
Article ADS CAS Google Scholar
Taverna, R. et al. The long quest for vacuum birefringence in magnetars: 1E 1547.0-5408 and the elusive smoking gun. Astrophys. J. 1002, 102 (2026).
Article ADS CAS Google Scholar
Everett, J. E. & Weisberg, J. M. Emission beam geometry of selected pulsars derived from average pulse polarization data. Astrophys. J. 553, 341–357 (2001).
Article ADS Google Scholar
Ashton, G. et al. BILBY: a user-friendly Bayesian inference library for gravitational-wave astronomy. Astrophys. J. Suppl. Ser. 241, 27 (2019).
Article ADS CAS Google Scholar
Speagle, J. S. DYNESTY: a dynamic nested sampling package for estimating Bayesian posteriors and evidences. Mon. Not. R. Astron. Soc. 493, 3132–3158 (2020).
Article ADS Google Scholar
Gil, J., Gronkowski, P. & Rudnicki, W. Geometry of the emission region of PSR 0950+08. Astron. Astrophys. 132, 312–316 (1984).
ADS Google Scholar
Rankin, J. M. Toward an empirical theory of pulsar emission. IV. geometry of the core emission region. Astrophys. J. 352, 247 (1990).
Article ADS Google Scholar
Johnston, S. et al. The thousand-pulsar-array programme on MeerKAT - XI. application of the rotating vector model. Mon. Not. R. Astron. Soc. 520, 4801–4814 (2023).
Article ADS CAS Google Scholar
Lower, M. E., Johnston, S., Shannon, R. M., Bailes, M. & Camilo, F. The dynamic magnetosphere of Swift J1818.0-1607. Mon. Not. R. Astron. Soc. 502, 127–139 (2021).
Article ADS CAS Google Scholar
Desvignes, G. et al. A freely precessing magnetar following an X-ray outburst. Nat. Astron. 8, 617–627 (2024).
Article ADS PubMed PubMed Central Google Scholar
Foreman-Mackey, D., Hogg, D. W., Lang, D. & Goodman, J. emcee: the MCMC mammer. Publ. Astron. Soc. Pac. 125, 306 (2013).
Article ADS Google Scholar
Doroshenko, V. et al. Determination of X-ray pulsar geometry with IXPE polarimetry. Nat. Astron. 6, 1433–1443 (2022).
Article ADS Google Scholar
Naghizadeh-Khouei, J. & Clarke, D. On the statistical behaviour of the position angle of linear polarization. Astron. Astrophys. 274, 968 (1993).
ADS Google Scholar
González-Caniulef, D., Caiazzo, I. & Heyl, J. Unbinned likelihood analysis for X-ray polarization. Mon. Not. R. Astron. Soc. 519, 5902–5912 (2023).
Article ADS Google Scholar
Hibschman, J. A. & Arons, J. Polarization sweeps in rotation-powered pulsars. Astrophys. J. 546, 382–393 (2001).
Article ADS Google Scholar
Tong, H., Wang, P. F., Wang, H. G. & Yan, Z. Rotating vector model for magnetars. Mon. Not. R. Astron. Soc. 502, 1549–1556 (2021).
Article ADS Google Scholar
Lai, D. & Ho, W. C. G. Resonant conversion of photon modes due to vacuum polarization in a magnetized plasma: implications for X-ray emission from magnetars. Astrophys. J. 566, 373–377 (2002).
Article ADS Google Scholar
van Adelsberg, M. & Lai, D. Atmosphere models of magnetized neutron stars: QED effects, radiation spectra and polarization signals. Mon. Not. R. Astron. Soc. 373, 1495–1522 (2006).
Article ADS Google Scholar
Deller, A. T., Camilo, F., Reynolds, J. E. & Halpern, J. P. The proper motion of PSR J1550-5418 measured with VLBI: a second magnetar velocity measurement. Astrophys. J. Lett. 748, L1 (2012).
Article ADS Google Scholar
Download references
Acknowledgements
This work reports observations obtained with the IXPE, a joint US (NASA) and Italian (ASI) mission, led by Marshall Space Flight Center (MSFC). The research uses data products provided by the IXPE Science Operations Center (MSFC), using algorithms developed by the IXPE Collaboration (MSFC, Istituto Nazionale di Astrofisica (INAF), Istituto Nazionale di Fisica Nucleare (INFN) and ASI Space Science Data Center (SSDC)) and distributed by the High-Energy Astrophysics Science Archive Research Center (HEASARC). Murriyang, the Parkes radio telescope of CSIRO, is part of the Australia Telescope National Facility (https://ror.org/05qajvd42), which is funded by the Australian Government for operation as a National Facility managed by CSIRO. We acknowledge the Wiradjuri people as the Traditional Owners of the Observatory site. This project was supported by resources and expertise provided by CSIRO IMT Scientific Computing and made use of the Ngarrgu Tindebeek supercomputer at the OzSTAR National Facility at Swinburne University of Technology. The OzSTAR programme receives funding in part from the Astronomy National Collaborative Research Infrastructure Strategy (NCRIS) allocation provided by the Australian Government, and from the Victorian Higher Education State Investment Fund (VHESIF) provided by the Victorian Government. G.Y. acknowledges constructive discussion with A. Philippov on radio emission from magnetars and P. Ray on radio/X-ray timing analysis.
Funding
The material is based on work supported by NASA under award no. 80GSFC24M0006. G.Y. acknowledges NASA support under grant nos. 80NSSC25K7257 and 80NSSC25K0283, through which R.E.S. and A.V.K. are partially supported. M.E.L. is supported by an Australian Research Council Discovery Early Career Research Award DE250100508. M.G.B. thanks NASA for generous support under grant nos. 80NSSC24K0589, 80NSSC25K7257 and 80NSSC25K0079. W.C.G.H. acknowledges support through grant no. 80NSSC23K0078 from NASA. J.B.C. acknowledges support under NASA award no. 80GSFC21M0006. F.C., A.K.H., T.E., C.P.H., P.K., M.N., P.S. and Z.W. do not declare relevant funding.
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature thanks the anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
Extended Data Fig. 1 Model-dependent and model-independent polarization characteristics comparison.
A comparative study between model-independent Stokes Q and U polarization characteristics obtained through using IXPEOBSSIM (Panel a) and through applying a polconst*bbodyrad spectral model in Xspec (Panel b) at four energy bins: 2–3 keV, 3–4 keV, 4–5 keV, and 2–8 keV. The resulting PD and PA for the two methods are comparable to each other. Notably, both show a non-linear trend in the energy-dependence of the PD as the 4–5 keV band approaches the value of the 2–3 keV band at ~ 1σ level. Future observations with higher count statistics are needed to complement this study to more directly probe the nature of the high-energy polarization.
Extended Data Fig. 2 Phase-averaged IXPE spectra.
The simultaneous spectral νFν models of the NICER+IXPE observations. Panel a displays the best fit model, a single absorbed blackbody with a linear polarization component: constant*tbabs(pollin*bbodyrad). NB: Only the IXPE spectra are displayed here for the sake of visual clarity. Panel b shows the data divided by the folded model for the absorbed BB model. The right-hand panels show the normalized Stokes Q/I (Panel c) and U/I (Panel d) spectra in linear space for the three IXPE DUs with the solid lines showing the best fit of the linear polarization component. The quasi-thermal BB is accompanied by a strong polarization signal that decreases as a function of energy.
Extended Data Fig. 3 Phase-averaged NICER spectra.
Simultaneous spectral modeling of the NICER+IXPE observations in the νFν representation (see also Extended Data Fig. 2). The top panel shows the best-fit model to the NICER spectra, described by a single absorbed blackbody. For visual clarity, we do not display the SCORPEON background model components. The bottom panel shows the ratio of the NICER data to the folded model.
Extended Data Fig. 4 X-ray spectro-polarimetric posterior distributions.
Corner plots of the one- and two-dimensional posterior distributions of the X-ray spectro-polarimetric parameters from the best-fit model shown in Table 1 produced by an MCMC chain (Panels a-u). The contours denote the 11.8%, 39.3%, 67.5%, 86.4% credible regions (corresponding to 0.5, 1, 1.5, and 2σ).
Extended Data Fig. 5 Radio RVM posterior distributions.
Corner plots of the one- and two-dimensional posterior distributions of the radio RVM fit parameters (Panels a-n). The contours denote the 39%, 87%, and 99% credible regions (corresponding to 1, 2, and 3σ).
Extended Data Fig. 6 Phase-resolved polarization characteristics in three X-ray bands.
Complementary plot to Fig. 2, with IXPE intensity, PD, and PA (shown in the lower three panels, respectively) binned according to 2–3 keV (a–e), 3–4 keV (f–j), and 4–8 keV (k–o).
Extended Data Fig. 7 RVM fit to the phase-resolved X-ray PA.
Panels a-n: Black contours represent the one- and two-dimensional posterior distributions for the RVM parameters α, ζ, β, ψ0, and ϕ0, obtained from a Bayesian fit to the 2-4 keV phase-resolved polarization angle profile using the55 wrapped-angle likelihood. Contours denote the 39%, 87%, and 99% credible regions (corresponding to 1, 2, and 3σ). The blue contours are those derived from the radio RVM fit overlaid for ease of comparison. Panel o: Phase-resolved X-ray PA measurements (black points) overlaid with the maximum-posterior RVM model (solid black line). Light gray curves show random posterior draws, illustrating the range of model realizations consistent with the data.
Extended Data Fig. 8 Stokes Q/I, U/I, and I modeling.
Simulated Stokes Q/I (panel a), U/I (panel b), and intensity (panel c) pulse profiles from MAGTHOMSCATT for a single hotspot wedge offset from the magnetic pole as a function of rotational phase (solid lines). The black dots represent the Stokes Q/I (panel a), U/I (panel b), and intensity (panel c) data extracted from IXPE in the 2–3 keV energy range. The case with the best statistical fit (lowest combined total χ2; see Methods) incorporates vacuum birefringence, and corresponds to a magnetic colatitude of θm = [0°, 17°] and longitude ϕm = [0°, 120°] (zero longitude contains the rotation and magnetic axes); it is displayed in blue. The orange solid curves on the center and right panel show the corresponding polarization profiles wherein VB is turned off - these fits are statistically worse than those with VB on (see Methods). Panel (d) displays the comparison in the Stokes Q-U space between the observed data (black dots) and the simulated results obtained with the best-fit wedge-shaped hotspot with (blue line) and without (orange line) including magnetospheric VB, plotted for one rotational cycle. The red solid line represents the best result among the VB-off cases (see Methods), corresponding to a pole-centered circular hotspot with a magnetic inclination of α = 2° and a viewing angle of ζ = 20. 5°.
Full size table
Full size table
Supplementary information
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
Reprints and permissions
About this article
Cite this article
Stewart, R.E., Dinh Thi, H., Younes, G. et al. Vacuum birefringence and the polarized X-ray emission from a radio magnetar. Nature (2026). https://doi.org/10.1038/s41586-026-10859-z
Download citation
Received:
Accepted:
Published:
Version of record:
DOI: https://doi.org/10.1038/s41586-026-10859-z