## Discovery of polarized X-ray emission from the accreting millisecond pulsar SRGA J144459.2 -604207

Alessandro Papitto 1 /star , Alessandro Di Marco 2 , Juri Poutanen 3 4 , , Tuomo Salmi 5 6 , , Giulia Illiano 1 7 8 , , , Fabio La Monaca 2 7 8 , , , Filippo Ambrosino 1 , Anna Bobrikova 3 , Maria Cristina Baglio 9 , Caterina Ballocco 1 8 , , Luciano Burderi 10 11 , , Sergio Campana 9 , Francesco Coti Zelati 12 13 , , Tiziana Di Salvo 14 , Riccardo La Placa 1 , Vladislav Loktev 3 , Sinan Long 15 , Christian Malacaria 1 , Arianna Miraval Zanon 16 , Mason Ng 17 , Maura Pilia 18 , Andrea Sanna 11 , Luigi Stella 1 , Tod Strohmayer 19 20 , , and Silvia Zane 15

- 1 INAF Osservatorio Astronomico di Roma, Via Frascati 33, 00078 Monte Porzio Catone (RM), Italy
- 3 Department of Physics and Astronomy, 20014 University of Turku, Finland
- 2 INAF Istituto di Astrofisica e Planetologia Spaziali, Via del Fosso del Cavaliere 100, 00133 Roma, Italy
- 4 Space Research Institute, Russian Academy of Sciences, Profsoyuznaya 84 32, 117997 Moscow, Russia /
- 6 Department of Physics, P.O. Box 64, 00014 University of Helsinki, Finland
- 5 Anton Pannekoek Institute for Astronomy, University of Amsterdam, Science Park 904, 1098XH Amsterdam, the Netherlands
- 7 Dipartimento di Fisica, Università degli Studi di Roma 'Tor Vergata', Via della Ricerca Scientifica 1, 00133 Roma, Italy
- 9 INAF Osservatorio Astronomico di Brera, Via E. Bianchi 46, 23807 Merate (LC), Italia
- 8 Dipartimento di Fisica, Università degli Studi di Roma 'La Sapienza', Piazzale Aldo Moro 5, 00185 Roma, Italy
- 10 INAF IASF Palermo, Via Ugo La Malfa 153, 90146 Palermo, Italy /
- 12 Institute of Space Sciences (ICE, CSIC), Campus UAB, Carrer de Can Magrans s n, 08193 Barcelona, Spain /
- 11 Dipartimento di Fisica, Università degli Studi di Cagliari, SP Monserrato-Sestu, KM 0.7, Monserrato 09042, Italy
- 13 Institut d'Estudis Espacials de Catalunya (IEEC), 08860 Castelldefels (Barcelona), Spain
- 15 Mullard Space Science Laboratory, University College London, Holmbury St Mary, Dorking, Surrey RH5 6NT, UK
- 14 Dipartimento di Fisica e Chimica - Emilio Segrè, Università di Palermo, Via Archirafi 36, 90123 Palermo, Italy
- 16 ASI - Agenzia Spaziale Italiana, Via del Politecnico snc, 00133 Rome (RM), Italy
- 18 INAF Osservatorio Astronomico di Cagliari, Via della Scienza 5, 09047 Selargius (CA), Italy
- 17 MIT Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
- 19 Astrophysics Science Division, NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA
- 20 Joint Space-Science Institute, NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA

Received August 2, 2024; accepted December 7, 2024

## ABSTRACT

We report on the discovery of polarized X-ray emission from an accreting millisecond pulsar. During a 10 day-long coverage of the February 2024 outburst of SRGA J144459.2 -604207, the Imaging X-ray Polarimetry Explorer (IXPE) detected an average polarization degree of the 2-8 keV emission of 2 3% . ± 0 4% at an angle of 59 . ◦ ± 6 ◦ (East of North; uncertainties quoted at the 1 σ confidence level). The polarized signal shows a significant energy dependence with a degree of 4 0% . ± 0 5% between 3 and 6 keV and . &lt; 1 5% . (90% c.l.) in the 2-3 keV range. We used NICER, XMM-Newton, and NuSTAR observations to obtain an accurate pulse timing solution and perform a phase-resolved polarimetric analysis of IXPE data. We did not detect any significant variability of the Stokes parameters Q and U with the spin and the orbital phases. We used the relativistic rotating vector model to show that a moderately fan-beam emission from two point-like spots at a small magnetic obliquity ( /similarequal 10 ) is compatible with the observed pulse profile and ◦ polarization properties. IXPE also detected 52 type-I X-ray bursts, with a recurrence time ∆ t rec increasing from 2 to 8 h as a function of the observed count rate C as ∆ t rec ∝ C -0 8 . . We stacked the emission observed during all the bursts and obtained an upper limit on the polarization degree of 8 5% (90% c.l.). .

Key words. magnetic fields - methods: observational - polarization - pulsars: individual: SRGA J144459.2 -604207 - stars: neutron - X-rays: binaries

## 1. Introduction

The simultaneous measurement of the mass M and the equatorial radius R eq of a few neutron stars (NSs) is of paramount importance for constraining the equation of state of ultra-dense matter because of its one-to-one mapping to the NS mass-radius dependence (Lattimer &amp; Prakash 2016; Özel &amp; Freire 2016; Baym et al. 2018). Relativistic e ff ects bend the light trajectory

/star alessandro.papitto@inaf.it

and modify energy of the photons emitted by hot spots on the surface of millisecond pulsars (MSPs; see, e.g., Watts et al. 2016, and references therein). As a result, modeling the Xray pulse profiles of MSPs is one of the most powerful techniques to obtain the desired information on M and R eq (Poutanen &amp; Gierli´ nski 2003; Miller &amp; Lamb 2016). The pulse shape also depends on several geometrical parameters (e.g., the binary inclination i , the spot co-latitude θ , and its angular size ρ ) as well as spectral parameters (e.g., the angular pattern of the emitted radiation; see, e.g., Poutanen &amp; Beloborodov 2006 and

Section ). Breaking this degeneracy with X-ray spectral-timing data requires a very large number of counts (Lo et al. 2013; Miller &amp; Lamb 2015). A few megasecond-long observations of rotation-powered MSPs with the Neutron Star Interior Composition Explorer Mission (NICER) met this requirement and led to constraints on M and R eq with a relative uncertainty of ∼ 10%for three of them (Riley et al. 2019, 2021; Miller et al. 2019, 2021; Choudhury et al. 2024; Dittmann et al. 2024; Salmi et al. 2024; Vinciguerra et al. 2024).

Accreting millisecond pulsars (AMSPs; Patruno &amp; Watts 2021; Di Salvo &amp; Sanna 2022) represent an intriguing alternative. During accretion outbursts, they attain an X-ray luminosity of L X /similarequal 10 36 -10 37 erg s -1 , i.e. 5-6 orders of magnitude brighter than rotation-powered MSPs. The degeneracy between M and R eq and the spot geometrical and spectral parameters (see, e.g., Viironen &amp; Poutanen 2004) prevented to get meaningful constraints from the data taken with, e.g., the Rossi Xray Timing Explorer ( M = 1.2-1.6 M /circledot , R eq = 6-13 km by Poutanen &amp; Gierli´ nski 2003; M = 0.8-1.7 M /circledot , R eq = 5-13 km by Morsink &amp; Leahy 2011). However, Salmi et al. (2018) showed that if i and θ were known a priori with an accuracy of a few degrees, those data would have allowed measurements with a relative uncertainty comparable to or slightly smaller than that for rotation-powered MSPs.

X-ray polarimetry might provide key information on the geometrical parameters of the spots. Soft photons from the NS surface are Compton up-scattered by hot electrons in the accretion shock (Poutanen &amp; Gierli´ nski 2003) achieving polarization degree (PD) from a few per cent up to 1020% (Nagirner &amp; Poutanen 1994; Viironen &amp; Poutanen 2004; Salmi et al. 2021; Bobrikova et al. 2023). For the magnetic field strength of AMSPs (10 8 -10 9 G), the typical energy of X-ray photons is orders of magnitude larger than the fundamental cyclotron energy ( /similarequal 1-10 eV), and scattering is not a ff ected by the magnetic field. As a result of azimuthal symmetry, the polarization vector of the X-rays emitted from the accretion shock is expected to lie either in the meridional plane, formed by the photon momentum and the normal to the surface, or perpendicular to it (Viironen &amp; Poutanen 2004; Poutanen 2020). As the NS rotates, the polarization angle (PA) should thus swing as a function of the pulsar spin phase in accordance with the rotating vector model (RVM; Radhakrishnan &amp; Cooke 1969; Meszaros et al. 1988). In AMSPs, the relativistic motion of the emission regions causes additional rotation of the polarization plane (Poutanen 2020); also the oblate shape of a quickly rotating NS has to be accounted for (Loktev et al. 2020).

Here, we report on the discovery of polarized X-ray emission from the AMSP SRGA J144459.2 -604207 (SRGA J1444 in the following) by the Imaging X-ray Polarimetry Explorer (IXPE; Weisskopf et al. 2022) in the context of a campaign which included also NICER, X-ray Multi-Mirror Mission (XMMNewton), and Nuclear Spectroscopic Telescope Array (NuSTAR) observations. Discovered in outburst on 2024 February 21 by the Mikhail Pavlinsky ART-XC telescope on-board the Spectrum-Roentgen-Gamma observatory (SRG, Molkov et al. 2024) at a relatively bright 4-12 keV X-ray flux of 2 × 10 -9 erg s -1 cm -2 (see Fig. 1), SRGA J1444 showed coherent pulsations at 447.9 Hz first detected by NICER, coming from a NS in a 5.2 h orbit around a donor with a mass in the range 0.2-0.7 M /circledot (Ng et al. 2024). SRGA J1444 also showed type-I Xray bursts (Mariani et al. 2024) quite regularly, with a recurrence time increasing from 1.6 to 2.2 h as the X-ray flux was observed to decrease (Molkov et al. 2024). No optical counterpart could be identified (Baglio et al. 2024) at the most accurate X-ray po-

Fig. 1. Light curve of the 2024 outburst of SRGA J1444 observed by MAXI (Matsuoka et al. 2009). We converted 2-20 keV observed count rates into 2-10 keV flux values assuming that the spectrum is described by a power law with a photon index Γ = 1 9 absorbed by an equivalent . hydrogen column of N H = 2 9 . × 10 22 cm 2 -(Ng et al. 2024). Horizontal bars indicate the time intervals covered by observations of the instruments discussed in this paper.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00608v2-figure-452f8d87-3233-46a2-b2f6-aef821bcc74f.jpg)

sition obtained with High Resolution Camera observations on board the Chandra X-ray observatory (CXO, Illiano et al. 2024), most likely due to the extinction of A V /similarequal 7 mag (estimated 1 using the maps of Doroshenko 2024, at a distance of 8 kpc) in that direction. Russell et al. (2024) reported a radio counterpart with a flux density of ∼ 200 µ Jy at 5.5 GHz and a flat spectral shape, consistent with a compact jet or discrete ejecta.

This paper is organized as follows. Section 2 details the analysis of the di ff erent X-ray datasets considered. Section 3 presents the properties of the X-ray persistent 2 emission. The measurement of polarized emission is presented in Sect. 3.1.1 and its interpretation is discussed in Sect. 3.1.2. Section 3.2.1 reports the pulsar timing solution obtained from the analysis of the coherent X-ray pulsations observed by the various instruments involved in this analysis, Sect. 3.2.2 the observed phase-resolved variability of the polarized emission, and Sect. 3.2.3 its modeling. Section 4 outlines the properties of recurrence of the type-I X-ray bursts detected by IXPE, Sect. 4.1 its the X-ray polarization and Sect. 4.2 the physical implications. The main results of this paper are summarized in Sect. 5.

## 2. Observations

Table 1 lists the observations analyzed in this paper, using HEAS oft package version 6.33.2. IXPE started observing SRGA J1444 less than a week after the first detection of the source, in response to the trigger of the General Observer (GO) program 03250101 (PI: A. Papitto). Observations lasted ∼ 10 d with visibility gaps giving a duty cycle exceeding 60%. We used the pcube algorithm in the xpbin tool of the ixpeobssim package version 30.6.4 (Baldini et al. 2022), along with the IXPE calibration database (CALDB) released on 2024 February 28, to extract polarimetric information using the formalism from Kislat et al. (2015). We selected source photons from a circular region with a radius of 100 ′′ centered at the source posi-

1 see the tool available at http://astro.uni-tuebingen.de/nh3d/nhtool

2 Here, persistent refers to the accretion-powered X-ray emission seen during the outburst, as opposed to the burst emission.

Table 1. Observations analyzed in this paper.

| Observatory   | Obs. ID       | Start date (MJD)    | End Date (MJD)      | Instrument   | Exposure (s)   | Energy band (keV)   |
|---------------|---------------|---------------------|---------------------|--------------|----------------|---------------------|
| IXPE          | 03250101      | 60367.556           | 60377.676           | DU1 DU2      | 552930 553541  | 2-8                 |
| NICER         | 620419 663908 | 60361.831 60364.219 | 60363.908 60392.971 | XTI          | 9536 24805     | 0.3-10              |
| NuSTAR        | 80901307002   | 60366.459           | 60369.883           | FPMA FPMB    | 157675 157440  | 3-79                |
| XMM-Newton    | 0923171501    | 60368.102           | 60369.376           | EPIC-pn      | 110077         | 0.3-10              |

tion. We use the unweighted analysis implemented in ixpeob -ssim (Di Marco et al. 2022) to measure the Stokes parameters ( I , Q U , ). Figure 2 shows the 2-8 keV IXPE count rate, the polarization degree PD = √ q 2 + u 2 and angle PA = 1 2 arctan( u q / ) measured in the 3-6 keV band, where a significant polarization is detected (see Sect. 3.1.1 and Fig. 4). Here, q = Q I / and u = U I / . As the source count rate is relatively high ( ∼ 2-8 count s -1 ), the background contributes at most 5% of the total rate and we did not subtract it following the recommendation of Di Marco et al. (2023). For the IXPE spectral analysis, we used the response matrices 20240101\_ alpha 075present in the CALDB, computed ancillary response files using ixpecalcarf weighting the Stokes parameters following the analysis presented by Di Marco et al. (2022) and binned the I spectra to have at least 30 counts per bin. To ensure statistics high enough for a spectro-polarimetric analysis, we binned data to a constant resolution of 200 eV.

Ng et al. (2024) reported the analysis of NICER observations during the first three days of the outburst using publicly available data. Here, we also include data taken afterwards in the GO program ID 663908 (PI: Papitto). Visibility constraints and high background prevented NICER to collect data simultaneously with the IXPE (see Fig. 1). The monitoring briefly resumed on 2024 March 14 and more regularly after March 21, when the source count rate had already decreased close to the background value. We reduced the observations using calibration version xti 20240206 and the same filtering criteria adopted by Ng et al. (2024), retaining photons with energy ranging between 0.3 and 10 keV for a pulse timing analysis.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00608v2-figure-a23da604-3882-43f9-bf5f-83ace68ce628.jpg)

The European Photon Imaging Camera (EPIC) pn-CCD camera on-board XMM-Newton observed SRGA J1444 uninterruptedly for ∼ 110 ks. We retained 0.3-10 keV source photons extracted from a 21 pixel-wide stripe (equivalent to 86 ′′ . 1) around the position of the pulsar to perform a pulse timing analysis.

NuSTAR provided the longest overlap with IXPE data starting on 2024 February 26. We reduced the data using the calibration version indx 20240325, with the fine clock correction file ver. 192. A circular region with a wide radius of 140 ′′ was adopted to extract source photons in the 3-79 keV band.

The times of arrival of the photons recorded by all the instruments were transformed into the inertial reference frame centred at the Solar System barycenter using the Chandra coordinates reported by Illiano et al. (2024), R.A. = 14 h 44 m 58 s . 94, Dec = -60 41 55 ◦ ′ ′′ . 3 (J2000) and the DE 405 ephemeris. We detected 52, 5, 13 and 23 type-I X-ray bursts during IXPE, NICER, XMM-Newton and NuSTAR exposures, respectively. When analysing the persistent emission, we discarded intervals between 5-s before and 60-s after the onset of each burst.

Fig. 2. IXPE 2-8 keV light curve in 200 s time bins (top panel). Type-I X-ray bursts are easily recognized. PD (middle panel) and PA (bottom panel) of the total (i.e., persistent and bursting) emission measured in the 3-6 keV band and evaluated over ∼ 2 d-long intervals. The blue shaded region indicates the interval after the pulse amplitude increase observed at MJD 60376.5 in which a 90% c.l. upper limit on the PD is reported.

## 3. Persistent X-ray emission

## 3.1. Phase-averaged emission

## 3.1.1. Average polarization

Figure 3 shows the average normalized Stokes parameters observed during the persistent emission by the three IXPE detector units. The average Stokes parameters are q = -1 1% . ± 0 4% . and u = 2 0% . ± 0 4%, which translate into an average PD of . 2 3% . ± 0 4%, and an angle of . PA 59 ◦ ± 6 ◦ (measured East of North, errors are quoted at the 1 σ c.l. unless di ff erently specified; see Fig. 3). The probability of detecting such values of the normalized Stokes parameters from noise fluctuations is 1 1 . × 10 -7 (corresponding to a 5.3 σ c.l.). Figure 4

Fig. 3. Average normalized Stokes parameters observed during the 28 keV persistent emission for the single IXPE detector units (DUs) and for the sum.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00608v2-figure-cf6c472f-0b0b-4648-b910-95e558ce38da.jpg)

Fig. 4. PD (top) and PA (bottom) as a function of energy. Above 6 keV, polarization is compatible with zero and we report upper limits at 90% c.l., while the PA is unconstrained.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00608v2-figure-d8a54ab2-8f88-4794-9956-f455543d2401.jpg)

shows the results of an energy-resolved analysis. Polarized emission is detected at the highest statistical significance in the 36 keV energy range with an average PD3 -6 = 4 0% . ± 0 5% and . PA3 6 -= 57 ◦ ± 3 . ◦ Between 2 and 3 keV, the PD is not significant (0 7% . ± 0 6%) with an upper limit at 90% confidence . level of PD2 -3 &lt; 1 5%. Above 6 keV, the 5% upper limit on . the PD becomes larger than both PD3 -6 and the minimum detectable polarization at 99% confidence level in the 3-6 keV band (MDP99 = 1.4%). We did not detect any significant variation or rotation of the 3-6 keV polarization vector over 2 d-long intervals, even though the PD slightly increases halfway through the observation (see bottom panels of Fig. 2). We checked that no significant variation of the spectral parameters occurs in the time interval characterized by a slightly larger PD compared to the emission observed elsewhere.

Similar to Ng et al. (2024), we used the sum of a black-body ( bbodyrad in xspec ) and thermal Comptonization ( nthcomp ) to satisfactorily model ( χ 2 = 483 9 . for 441 d.o.f.) the 2-8 keV IXPE energy spectrum. We modelled absorption with the Tuebingen-Boulder model (Wilms et al. 2000) with the absorption column fixed to N H = 2 9 . × 10 22 cm -2 and quote uncertainties on the spectral parameters at 90% c.l. The thermal component has a temperature of kT bb = 0 52 . ± 0 01 keV and an appar-. ent size of R bb = (14 0 . ± 0 5) . d 8 km, where d 8 is the distance to the source in units of 8 kpc. Comptonization is described by an asymptotic power law with index Γ = 1 70 . ± 0 03 and a soft pho-. ton input with temperature &lt; 0 15 keV. The electron temperature . was fixed at a value beyond the IXPE band-pass kT e = 15 keV. The Comptonization component accounts for ∼ 95% of the 28 keV observed flux of (5 6 . ± 0 1) . × 10 -10 erg s -1 cm -2 . Spectral analysis of the NICER observation (Obs. ID 6639080103) performed ∼ 12 h before the start of the IXPE coverage indicates a smaller black-body size ( kT bb = 0 68 . ± 0 05 keV; . R bb = (5 7 . ± 1 0) . d 8 km) and a steeper slope of the Comptonized emission ( Γ = 1 90 . ± 0 01). .

The decrease of the PD below 3 keV, i.e. where the soft thermal component gives its highest contribution to the 2-8 keV spectrum, motivated us to perform a spectro-polarimetric analysis. By fitting the spectra of the three Stokes parameters with a constant polarization component to multiply the model presented above, tbabs*(bbodyrad+nthcomp)*polconst , gave PD = 2 6% . ± 0 5% and PA . = 59 ◦ ± 6 ◦ with a χ 2 / d.o.f. of 340 257, / compatible with the pcube results. By using separate constant polarization for the two spectral components gave a slight improvement of the fit χ 2 / d.o.f. = 328 255with a / bbodyrad component polarized with the PD = 8% ± 5% at PA = -43 ◦ ± 19 ◦ and the nthcomp with the PD = 4 6% . ± 1 1% and PA . = 55 ◦ ± 7 . ◦

## 3.1.2. Interpretation

For the first time, IXPE observations of SRGA J1444 allowed us to significantly detect polarization from an AMSP with an average 2-8 keV degree of 2 3% . ± 0 4%. Such a value increases . to 4 0% . ± 0 5% if the analysis is restricted to the 3-6 keV en-. ergy band. Comptonization by hot ( kT e ∼ 20-60 keV) electrons in slabs of moderate optical depth ( τ ≈ 1-2) located above hot spots ( kT bb ∼ 1 keV) on the NS surface accounts for more than 90% of the emission in the IXPE bandpass of AMSPs such as SRGA J1444 (Poutanen 2006) and represents the most natural explanation of the polarization we detected. Values of the PD as high as 10-20% were predicted by modelling the scattering in Thomson approximation (Viironen &amp; Poutanen 2004; Salmi et al. 2021). However, the actual degree of polarization is a function of the physical properties of the scattering medium (Nagirner &amp; Poutanen 1994). Recently, Bobrikova et al. (2023) applied the formalism for Compton scattering in a hot slab (Nagirner &amp; Poutanen 1993; Poutanen &amp; Svensson 1996) to describe the Stokes parameters as a function of energy and angle in terms of the electron temperature and optical depth of the slab and the seed photon temperature. They found values of the PD of a few per cent at the most, decreasing to values undetectable to IXPE at high ( &gt; 50 keV) electron temperatures. They concluded that chances for IXPE to detect polarization from an AMSP significantly increased for systems with a larger optical depth ( τ &gt; 2) and lower electron temperature ( kT e &lt; 30 keV). Molkov et al. (2024) fitted the 4-30 keV ART-XC spectrum with Comptonization from a medium with kT e ∼ 25 keV and τ ∼ 3 2 (evaluated from the photon index of . Γ = 1 8 as in, . e.g., Lightman &amp; Zdziarski 1987). Preliminary fitting of the 3-

79 keV NuSTAR spectrum (Malacaria et al., in prep.) suggests an even lower temperature ( kT e = 10 9 . + 0 2 . -0 1 . keV) and higher optical depth ( Γ = 2 2 . ± 0 1, . τ = 4 0 . ± 0 3). In the framework devel-. oped by Bobrikova et al. (2023), SRGA J1444 may thus present favorable spectral conditions to produce a relatively high degree of X-ray polarization compared to other AMSPs which might show even lower values. In addition, Bobrikova et al. (2023) showed that the observed polarization degree increased when an observer mostly sees the spot under high angles (i.e., | i -θ | close to 90 ). Such a configuration is supported by the results ◦ presented in Sect. 3.2.2.

We detected a significant decrease in the PD below 3 keV. In the IXPE bandpass, that energy band is where the nonupscattered blackbody emission contributes the most to the total observed flux. These photons are not expected to be polarized and might explain the decrease of the total PD. However, the soft thermal component accounts for just 5-10% of the total emission (7 7% . ± 1 4% between 2 and 3 keV, according to the . NICER spectral analysis reported in Sect. 3.1.1). As a result, the influence of such an unpolarized component at low energies does not readily explain the observed decrease of the PD. The analysis presented in Sect. 3.1.1 shows that, to explain the decrease of the total polarization, the soft thermal component should be polarized to a degree of 8% ± 5% and its polarization angle should di ff er by ∼ 90 ◦ from the polarization angle compared to the Comptonization component. However, the improvement in the fit description compared to a model with the same polarization properties for the two spectral components is modest. As thermal unscattered emission is not expected to be polarized, an intrinsic decrease of the polarization degree of the Compton scattered photons is a more likely explanation of the observed energy dependence. Bobrikova et al. (2023) predicted a swing by ∼ 90 ◦ of the polarization vector of the Comptonized emission at an energy of /similarequal 5 keV, only slightly higher than the turn-over energy of /similarequal 3 keV we observed. A comprehensive analysis of the spectral properties of the source, including also XMM-Newton data, is in preparation.

## 3.2. X-ray pulsed variability

## 3.2.1. Pulse phase timing solution

We used the pulsar timing solutions based on the first few days of NICER observations (Ng et al. 2024) to correct the arrival time of X-ray photons for the light travel time delays due to the pulsar orbit. We used two harmonic components to model the pulse profiles obtained by folding data taken in ∼ 500 slong intervals and sampled by 16 phase bins. Table 2 shows the timing solutions based on the observed evolution of the phase computed over the fundamental harmonic for the various instruments considered here, plotted in Fig. 5 with the best-fitting linear trend). We attribute the residual random phase variability by ∼ 0.05-0.1 cycles to pulse timing noise, as observed from other AMSPs with an even larger magnitude (Patruno &amp; Watts 2021; Di Salvo &amp; Sanna 2022). The two values of the pulse phase measured by NICER at MJD 60384, after an interruption of the monitoring due to visibility constraints, line up with the trend observed earlier and allow us to put a tight upper limit on the spin frequency derivative (˙ ν &lt; 1 2 . × 10 -13 Hz s -1 ). However, we caution that such an alignment could be fortuitous, as the pulse phase of AMSPs often shows erratic variability, especially in the late stages of an outburst (see, e.g. Illiano et al. 2023).

We used these ephemerides to fold 5 ks-long intervals of IXPE data and recover the pulsed signal despite the lower count-

Fig. 5. Pulse phase computed over the first harmonic of the pulses observed by NICER (magenta), NuSTAR (light blue), XMM-Newton (green) and IXPE (yellow) by folding 500 s-long intervals (5 ks for IXPE) around ν fold (see Table 2) with a reference epoch 60361 8 MJD. . Solid lines mark the best-fit linear trends.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00608v2-figure-22577a0d-ff91-46ab-afc2-757dda1f69d6.jpg)

Fig. 6. Pulse fractional amplitude and phase of the first harmonic component used to model the pulse profile obtained by folding IXPE data into 5 ks-long interval around the timing solution reported in the rightmost column of Table 2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00608v2-figure-288bdb3d-0eeb-400f-8de0-77ec249afaa2.jpg)

ing statistics. Figure 6 shows the fractional amplitude and phase of the fundamental harmonic observed by IXPE. An increase of the amplitude from ∼ 4 to 15% is observed after MJD 60376.5.

The frequencies measured by NuSTAR, IXPE and XMMNewton di ff er by δν/ν /similarequal 0 3 . × 10 -9 , 0 6 . × 10 -9 and 2 × 10 -9 with respect to the NICER measurement, respectively. The latter is compatible with the value quoted in the latest calibration technical note for the XMM-Newton clock relative time accuracy 3 ( δν/ν /similarequal 10 -8 ). The pulse phase recorded by NuSTAR and IXPE at epoch T 0 precedes those observed by NICER by 0.07-0.1 cycles. The NuSTAR shift can be understood in terms of the dependence of the pulse phase on energy already noticed by SRG ART-XC (Molkov et al. 2024); the reason for the IXPE / mismatch is less clear and will be addressed in a future paper. For consistency, the phase-resolved analysis presented in Sect. 3.2.2 is based on the IXPE timing solution, even if slightly less accurate.

3 https://xmmweb.esac.esa.int/docs/documents/CAL-TN-0018.pdf

Table 2. Pulse phase timing solutions of SRGA J1444.

|                                                                                                     | NICER                                                                                                                       | XMM-Newton                                                                                                                   | NuSTAR                                                                                                                         | IXPE                                                                                                                           |
|-----------------------------------------------------------------------------------------------------|-----------------------------------------------------------------------------------------------------------------------------|------------------------------------------------------------------------------------------------------------------------------|--------------------------------------------------------------------------------------------------------------------------------|--------------------------------------------------------------------------------------------------------------------------------|
| Interval (MJD) φ ( T 0 ) ν ( T 0 ) (Hz) ˙ ν (Hzs - 1 ) a 1 sin i / c (lt-s) P orb (s) T asc (MJD) e | 60361.837-60383.544 0 . 47(1) 447 . 87156112(2) < 1 . 2 × 10 - 13 0 . 65052(2) 18803 . 65(1) 60361 . 641306(3) < 9 × 10 - 5 | 60368.106-60369.374 0 . 500(7) 447 . 8715620(3) < 7 . 0 × 10 - 12 0 . 650491(5) 18803 . 65(1) 60361 . 641305(6) < 9 × 10 - 5 | 60366.465-60369.881 0 . 400(5) 447 . 87156127(2) < 1 . 5 × 10 - 12 0 . 650498(3) 18803 . 673(4) 60361 . 641297(2) < 3 × 10 - 5 | 60367.615-60377.619 0 . 37(1) 447 . 87156138(8) < 3 . 0 × 10 - 13 0.650486(13) 18803 . 665(4) 60361.641310(5) < 1 . 0 × 10 - 4 |

Notes. All the solutions presented in the table were obtained by folding data around ν fold = 447 87156116029 Hz with a reference epoch MJD . 60361.8. To facilitate the comparison, the frequency estimates were evaluated at T 0 = MJD 60368.0. All the epochs reported here are in the Barycentric Dynamical Time (TDB) system. Uncertainties at the 1 σ confidence level are reported in parentheses. Upper limits are given at a 90% confidence level.

## 3.2.2. Phase-resolved polarization variability

To perform a phase-resolved analysis of the polarimetric properties of SRGA J1444, we considered IXPE data corrected with the orbital parameters listed in Table 2. We used the pcube algorithm to evaluate in n = 16 phase bins the Stokes parameters Q and U , which are expected to be normally distributed (González-Caniulef et al. 2023; Suleimanov et al. 2023). We restricted the analysis to the 3-6 keV energy band in order to maximize the strength of the polarized signal. We used the pcube algorithm without acceptance corrections to measure the Stokes parameters in counts per second.

Figure 7 shows the phase-resolved Stokes parameters, the PD, and the PA observed before the pulse amplitude increase which occurred at MJD 60376.5. The Stokes parameters Q and U did not show significant variability. A fit with a constant function gave a χ 2 of 16.6 and 15.6 for 15 d.o.f., respectively, with average values of q = Q I / = -1 6% . ± 0 5% and . u = U I / = 3 5% . ± 0 4%. .

We also analyzed separately data taken after the pulse amplitude increase which occurred at MJD 60376.5 (see Fig. 8). Fitting with a constant gives q ′ = 0 5% . ± 1 4%and . u ′ = 2 4% . ± 2 7% . for χ 2 = 10 2 and 18.7 for 15 d.o.f. Both Stokes parameters are . compatible with zero within the uncertainties, even though the value of u ′ measured at phase 0.4 significantly deviated from the trend, leading to an increase of the PD up to ∼ 30%.

We also measured the Stokes parameters over 15 orbital phase bins to search for orbital variability of the polarization properties but found none significant. Constant values of Q and U modeled well the dependence of these values on the orbital phase ( χ 2 = 9 1 and 12.5 for 14 d.o.f., respectively). .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00608v2-figure-b274280f-f9a3-4b17-b7a1-db2f10af1588.jpg)

## 3.2.3. Modelling of the phase-resolved polarization properties

Assuming the accretion columns as the source of the polarized X-ray photons we observed, it is necessary to ask what spot geometrical configuration can explain at the same time the observed pulse profile and the non-significant phase-resolved variability of the Stokes parameters Q and U .

Unlike the case of highly magnetized, slowly rotating X-ray pulsars (Doroshenko et al. 2022; Tsygankov et al. 2022; Mushtukov et al. 2023; Marshall et al. 2022; Forsblom et al. 2023; Tsygankov et al. 2023; Malacaria et al. 2023; Doroshenko et al. 2023; Suleimanov et al. 2023; Poutanen et al. 2024; Forsblom et al. 2024), the polarization properties of

1

Fig. 7. Pulse phase dependence of polarimetric properties. From top to bottom: total flux normalized to the average, normalized Stokes parameters q = Q I / and u = U I / , the PD and PA observed by IXPE during the interval MJD 60367.6-60376.5, as a function of the spin phase. The solid red line indicates the model obtained by fitting the absolute Stokes parameters assuming two point-like spots and the parameters listed in text.

AMSPs are defined by the geometry of the spots, relativistic effects (Viironen &amp; Poutanen 2004; Poutanen 2020; Loktev et al. 2020), and possibly eclipses by the inner parts of the accretion disks (Poutanen et al. 2009; Ibragimov &amp; Poutanen 2009). This does not allow us to use a simple rotating vector model to fit the observed phase-resolved PD and PA, but requires a fit of the total Stokes vector ( I , Q U , ) which is the sum of the Stokes vectors of the two spots.

The spots' Stokes parameters depend on how the total and polarized radiation vary with the spin phase, photon energy

Fig. 8. Same as Fig. 7 but for the interval MJD 60376.5-60377.7 only.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00608v2-figure-98df6986-325a-4464-b3e6-14d791295ca1.jpg)

and zenith angle (Viironen &amp; Poutanen 2004; Poutanen 2020; Salmi et al. 2021; Bobrikova et al. 2023). A modeling that includes an accurate description of all these e ff ects is beyond the scope of this letter and will be presented elsewhere. Here, we draw considerations on the system's geometry by using the approximate formulae given by Poutanen &amp; Beloborodov (2006) and Poutanen (2020) for the phase-resolved Stokes parameters. We assume two point-like spots on the surface of a rapidly rotating NS (see Fig. 9 for a cartoon of the assumed geometry) which emit X-rays with a power-law spectrum with I ( α , ′ E ) = I 0 (1 + h cos α ′ ) E -( Γ -1) . Here, Γ = 1 7 is the spectral photon in-. dex, α ′ is the zenith angle (i.e., the angle between the normal to the spot and emission direction) in the co-rotating frame and h is an anisotropy parameter ( h = 0 for black-body emission, -1 &lt; h &lt; 0 for a fan-beam pattern such that expected from Comptonization in a thin slab; Poutanen &amp; Gierli´ski 2003). We n used Eq. (20-22) of Poutanen &amp; Beloborodov (2006) to determine the phase intervals in which the photons are visible to the observer and Eqs. (28-34) of the same paper to express the flux observed from each spot, Fi ( φ i + φ , 0 i , θ , i h ) ( i = 1 2) at a given , binary inclination i and colatitude θ i of the spot, and the resulting total flux F tot ( φ ) = N F ( 1 + F 2 ). Here, N is a normalization constant, φ 0 is the spin phase of the primary spot at the folding reference epoch. We assumed the phase of the secondary spot to be φ 2 = φ 1 + π and used the relativistic rotating vector model of Poutanen (2020) to express the phase-resolved Stokes parameters of each of the two spots (Eqs. 58-59), assuming that NFi gives the observed flux and that their polarization degree p 1 = p 2 does not depend on energy and zenith angle. We then used (Eqs. 21-23) of the same paper to compute the total expected Stokes vector, Q = Q 1 + Q 2 and U = U 1 + U 2 .

We fitted the phase-resolved Stokes parameters observed before the amplitude increase (see Fig. 7) with our simple model,

Fig. 9. Assumed geometry of the problem, adapted from Poutanen &amp; Beloborodov (2006). The binary inclination i , the co-latitutde of the two spots θ 1 and θ 2 , their phase φ 1 and φ 2 and the angle made by the direction of the emitted photon with respect to the normal to the spot in the co-moving frame α ′ are explictly indicated.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00608v2-figure-76923e20-b4cc-4c0e-9a54-0b5a09e7f8c5.jpg)

fixing the NS mass and radius to M = 1 4 M . /circledot and R eq = 10 km as the fit was insensitive to the choice of these parameters within a range of reasonable values for a NS. The fit (see Fig. 7) suggests a configuration with two visible spots as the fit χ 2 (38.8 for 41 d.o.f) improves compared to a fit with only one spot ( χ 2 = 63 0 . for 42 d.o.f.; the probability of chance improvement is /similarequal 10 -5 according to an F-test). A reasonably good fit can be obtained for any value of h ranging between 0 and -1; here we show the results obtained for a randomly chosen intermediate value of h = -0 3. We find a best-fit for . i = (74 1 . + 5 8 . -6 3 . ) ◦ , θ 1 = (11 8 . + 2 5 . -3 5 . ) ◦ , θ 2 = (172 6 . + 2 0 . -1 0 . ) ◦ , φ 0 = 0 57(4), . N = 6 90(6) . × 10 4 cnt, p 1 = 4 0% . ± 2 0% and . χ p = 57 2 ◦ . ± 0 5, where the latter is ◦ . the position angle of the pulsar angular momentum. We evaluated uncertainties at 1 σ confidence level from the variation of the fit χ 2 for the number of interesting parameters measured ( ∆ χ 2 = 8 15 for 7 parameters Lampton et al. 1976). Such un-. certainties have to be taken with caution, as they were obtained for fixed values of the NS mass, radius and emission anisotropy parameter. In addition, the large correlation between the parameters warrants a study of the posterior distribution of the parameters fitted using the likelihood function for I , Q and U , which will be presented in a forthcoming paper. The red line in Fig. 7 marks the predicted model of the normalized Stokes parameters obtained for these parameters. Bottom panels of the same figure show the total PD = √ q 2 + u 2 and PA = χ p + 1 2 tan -1 ( u q / ), respectively. The absence of significant variability with the pulse phase of the phase-resolved observed normalized Stokes parameters q and u can be explained by two almost antipodal spots if their magnetic inclination is small. Yet, the observed roughly constant phase-resolved trend gives little information on the other parameters which are essentially set by the total pulse profile alone. Molkov et al. (2024) used a modeling similar to ours (even though they assumed extended spots) and found qualita-

tively similar results ( i = 58 ◦ and θ 1 = 14 , ◦ θ 2 = π -θ 1 ). They modeled a pulse with a higher amplitude (9-12%) than ours ( A 1 = 5 2% . ± 0 1%) and a plateau close to the pulse mini-. mum, which they interpreted in terms of occultation by the inner accretion disk extending down to ∼ 25 km. The flattening of the pulse profile observed by IXPE is less pronounced and partial occultation is not required to obtain qualitative modeling of the profile. Because the pulse profile Molkov et al. (2024) modeled was observed when the X-ray flux was roughly double that observed by IXPE, this could indicate a recession of the inner accretion disk at a decreasing mass accretion rate.

After MJD 60376.5, the pulsed fraction observed by IXPE increased to A 1 = 18 7% . ± 0 8% (see Fig. 6). The same model . described above with the binary inclination and the anisotropy parameter fixed to i = 74 1 and ◦ . h = -0 3, respectively, provides . a good fit of the profile observed after that epoch, with a spot magnetic inclination and longitude increased by /similarequal 10 ◦ and ∼ 20 , ◦ respectively, and a polarization degree decreased to 2 0% . ± 1 5%, . compared to the observations performed previously. The observed change in the pulse profile shape could then reflect motions of the hot spots, as often assumed to explain the timing noise observed from AMSPs (see, e.g., Lamb et al. 2009). The profile predicted using our model is plotted as a red line in Fig. 8.

## 4. Type-I X-ray bursts

During the observations reported here, IXPE detected 52 typeI X-ray bursts which shared similar properties (see the inset of Fig. 10 for an example). The sample of peak count rates had an average of 38.7 cnt s -1 and a standard deviation of 2.9 cnt s -1 . After subtracting the persistent emission, the distribution of the fluence of the observed bursts had an average of 649 cnt with a standard deviation of 40; the fluence of the di ff erent bursts was very similar and only slightly increased for bursts occurring when the persistent count rate had decreased below 3 cnt s -1 . The relatively low number of counts observed by IXPE prevented us from performing a time-resolved spectral analysis of the single bursts. By stacking all the observed bursts, and fitting the observed spectra with the sum of a thermal component and the two-component model used to fit the persistent spectrum (see Sect. 3.1.1), we get 〈 kT burst 〉 = 1 83 . ± 0 05 keV . and R burst = (6 2 . ± 0 2) . d 8 km. Defining the burst duration as τ = F / F peak , and neglecting variations of the spectrum throughout the burst, gives τ = 16 8 . ± 1 6 s. Figure 10 shows the recur-. rence times ∆ t rec between the observed bursts 4 as a function of the average persistent count rate C measured during the longest possible contiguous interval before the burst. For 14 bursts, we scaled the time elapsed since the previous burst by a factor of 2 as the immediately preceding burst was likely missed due to interruption in the IXPE coverage (green points in Fig. 10). For one burst (red point), we used a scaling factor of 3.

## 4.1. X-ray polarimetry of burst emission

To search for a polarized signal during the X-ray bursts, we stacked the emission observed during all the events, retaining 32 s after the onset of each of the bursts for a total exposure of 1664 s. We measured q b = Q b / I b = -0 8% . ± 3 3% and . u b = U b / I b = 2 8% . ± 3 3%, resulting in . an upper limit on the PD of the sum of the burst and the persistent emission of

4 Consecutive bursts without data gaps could not be observed in IXPE data, but the uninterrupted XMM coverage allowed us to define the recurrence time unambiguously.

Fig. 10. Recurrence time of the type-I X-ray bursts observed by IXPE evaluated as the interval between consecutive bursts scaled by a factor 1 (blue points), 2 (green) and 3 (red) as a function of the persistent count rate observed before the burst C . The dashed line shows the relation ∆ t rec = KC -0 8 . with K = 12 75 h. The inset shows the profile of the . burst observed starting at MJD 60377.05969.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00608v2-figure-f2f3955e-90dd-4ed3-9d10-8d37368853dd.jpg)

PD UL = 7 2% (90% c.l.). We then subtracted the absolute Stokes . vectors observed during the persistent emission after normalizing them to the exposure time of 1664 s over which the burst emission was integrated, summing in quadrature the uncertainties, and obtained normalized and background subtracted Stokes parameters ˜b q = -0 7% . ± 3 9% and ˜b . u = 3 0% . ± 3 9%, corre-. sponding to an upper limit on the PD of 8.5% (90% c.l.).

## 4.2. Implications

To our knowledge, the observations presented here are the first in which IXPE caught thermonuclear X-ray bursts from an accreting NS. The ∼ 10 d coverage of the outburst of SRGA J1444 allowed us to detect 52 bursts and stack them when searching for a polarized signal. The upper limit we put on the PD of the burst emission (8.5%) is higher than the maximum value estimated by Lapidus &amp; Sunyaev (1985) assuming reflection of the burst emission by the accretion disk (3.7% for an inclination of 72 ). ◦ With a recurrence time variable between /similarequal 1 6 hr at the begin-. ning of the outburst (Molkov et al. 2024) and /similarequal 8 hr at the end of the IXPE observation presented here, SRGA J1444 was immediately recognized as a very prolific burster. The 11 Hz pulsar IGR J17480-2446 (which showed a recurrence time decreasing from /similarequal 30 to /similarequal 3 minutes; Motta et al. 2011; Chakraborty et al. 2011; Linares et al. 2012) and 4U 1636-536 (with a recurrence time of /similarequal 30 minutes and even less energetic short recurrence time bursts; Beri et al. 2019) are similar cases. Among accreting MSPs, SAX J17498-2021 showed bursts between 1 and 2 hr (Li et al. 2018), IGR J17511-3057 between 7 and 21 hr (Falanga et al. 2011) and IGR J17498-2921 every 16-18 hr (Falanga et al. 2012), in all cases with a recurrence time increasing as the X-ray flux declined. The ≈ 15 s duration of the bursts seen from SRGA J1444 suggests that they ignite in a mixed H He environment. Hot CNO cycle is expected to deplete the ac-/ creted hydrogen in 9 8( . X / 0 7)( . Z / 0 02) . -1 hr, where H and Z are the hydrogen and metals abundances (see Galloway et al. 2022, and references therein). The range of observed recurrence times ( &lt; 8 hr) suggests that complete depletion of hydrogen could only be achieved if the H abundance is sub-solar. For bursts igniting

in a mixed H He environment, the burst ignition depth is not ex-/ pected to depend much on the mass accretion rate ˙ M and the recurrence time scale should increase as ∆ t ∝ ˙ M -1 . Assuming that the IXPE count rate is a good tracer of the mass accretion rate, the observed dependence (see Fig. 10) is broadly compatible with such an expectation. Malacaria et al. (in prep.) will present a detailed analysis of the burst ignition process on this fast rotating pulsar allowed by the long and high duty cycle allowed by IXPE.

## 5. Conclusions

We presented the first X-ray polarimetric IXPE observation of an accreting millisecond pulsar in outburst, SRGA J1444, in the context of an observational campaign involving also NICER, NuSTAR and XMM-Newton observations. The main results are the following:

- -The 2-8 keV emission is significantly polarized, with an average degree of 2 3% . ± 0 4% at an angle of 59 . ◦ ± 6 ◦ (East of North).
- -We observed a significant change of the polarization properties with energies. The PD is maximum between 3 and 6 keV (4 0% . ± 0 5%) and decreases to . &lt; 2% (90% c.l.) in the 2-3 keV range.
- -The pulse phase of SRGA J1444 is rather stable throughout the observations presented here, with an upper limit on the spin frequency derivative of ˙ ν &lt; 1 2 . × 10 -13 Hz s -1 . During IXPE observations the pulse has an amplitude of ∼ 5%, before abruptly increasing to almost 20%, roughly one day before the end of the coverage.
- -The observed phase-resolved Stokes parameters are compatible with a constant function. A simple approximate model of the Stokes vector expected from the sum of two almost antipodal spots on the NS surface reproduces the observed values for a binary inclination of i = (74 1 . + 5 8 . -6 3 . ) ◦ and relatively small magnetic inclination of the spots θ 1 = (11 8 . + 2 5 . -3 5 . ) ◦ , θ 2 = (172 6 . + 2 0 . -1 0 . ) ◦ . The little phase-resolved variability shown by the Stokes parameters, if any, prevented us from obtaining constraints on the NS mass and radius based on our simple modeling.
- -IXPE observed 52 type-I X-ray bursts sharing similar properties (e.g. burst peak fluence, duration, fluence) with a recurrence time increasing as a function of the observed count rate, ∆ t rec ∝ C -0 8 . . After subtracting the persistent emission, the upper limit on the PD observed in the stacked burst emission is 8 5% (90% c.l.). .

The first significant detection of polarization from an accreting millisecond pulsar confirms theoretical expectations and strengthens the prospects of using this technique to measure the geometrical parameters of accreting MSPs and obtain constraints on the mass and radius of the NS from pulse profile modeling. Future IXPE observations of accreting millisecond pulsars outbursts will hopefully detect a larger phase-resolved variability and enable more accurate determination than those allowed by the present dataset.

Acknowledgements. We warmly thank the Directors and the Science Operations Team of IXPE, NICER, NuSTAR, and XMM for promptly scheduling the observations reported here and Bas Dorsman, Matteo Bachetti and Anna Watts for useful discussions.. The Imaging X-ray Polarimetry Explorer (IXPE) is a joint US and Italian mission. The US contribution is supported by the National Aeronautics and Space Administration (NASA) and led and managed by its Marshall Space Flight Center (MSFC), with industry partner Ball Aerospace (contract NNM15AA18C). The Italian contribution is supported by the Italian

Space Agency (Agenzia Spaziale Italiana, ASI) through contract ASI-OHBI2022-13-I.0, agreements ASI-INAF-2022-19-HH.0 and ASI-INFN-2017.13-H0, and its Space Science Data Center (SSDC) with agreements ASI-INAF-2022-14HH.0 and ASI-INFN 2021-43-HH.0, and by the Istituto Nazionale di Astrofisica (INAF) and the Istituto Nazionale di Fisica Nucleare (INFN) in Italy. This research used data products provided by the IXPE Team (MSFC, SSDC, INAF, and INFN) and distributed with additional software tools by the High-Energy Astrophysics Science Archive Research Center (HEASARC), at NASA Goddard Space Flight Center (GSFC). NICER is a 0.2-12 keV X-ray telescope operating on the International Space Station, funded by NASA. The NuSTAR mission is a project led by the California Institute of Technology, managed by the Jet Propulsion Laboratory, and funded by the National Aeronautics and Space Administration. Data analysis was performed using the NuSTAR Data Analysis Software (NuSTARDAS), jointly developed by the ASI Science Data Center (SSDC, Italy) and the California Institute of Technology (USA). XMM-Newton is an ESA science mission with instruments and contributions directly funded by ESA Member States and NASA. MN and this work were supported by NASA under grant 80NSSC24K1175. AP, GI, FA, RLP, CM and LS are supported by INAF (Research Grant 'Uncovering the optical beat of the fastest magnetised neutron stars (FANS)') and the Italian Ministry of University and Research (MUR) (PRIN 2020, Grant 2020BRP57Z, 'Gravitational and Electromagnetic-wave Sources in the Universe with current and next-generation detectors (GEMS)'). AP acknowledges support from the Fondazione Cariplo Cassa Depositi e Prestiti, grant / no. 2023-2560. JP thanks the Ministry of Science and Higher Education grant 075-15-2024-647 for support. AB acknowledges support from the Finnish Cultural Foundation grant 00240328. MCB acknowledges support from the INAFAstrofit fellowship. FCZ is supported by a Ramón y Cajal fellowship (grant agreement RYC2021-030888-I). TS acknowledges support from ERC Consolidator Grant No. 865768 AEONS (PI: Watts).

## References

- Baglio, M. C., Russell, D. M., Saikia, P., et al. 2024, The Astronomer's Telegram, 16487, 1
- Baldini, L., Bucciantini, N., Di Lalla, N., et al. 2022, SoftwareX, 19, 101194 Baym, G., Hatsuda, T., Kojo, T., et al. 2018, Reports on Progress in Physics, 81, 056902
- Beri, A., Paul, B., Yadav, J. S., et al. 2019, MNRAS, 482, 4397
- Bobrikova, A., Loktev, V., Salmi, T., &amp; Poutanen, J. 2023, A&amp;A, 678, A99
- Chakraborty, M., Bhattacharyya, S., &amp; Mukherjee, A. 2011, MNRAS, 418, 490
- Choudhury, D., Salmi, T., Vinciguerra, S., et al. 2024, ApJ, 971, L20
- Di Marco, A., Costa, E., Muleri, F., et al. 2022, AJ, 163, 170
- Di Marco, A., So ffi tta, P., Costa, E., et al. 2023, AJ, 165, 143
- Di Salvo, T. &amp; Sanna, A. 2022, in ASSL, Vol. 465, Millisecond Pulsars, ed. S. Bhattacharyya, A. Papitto, &amp; D. Bhattacharya, 87-124
- Dittmann, A. J., Miller, M. C., Lamb, F. K., et al. 2024, ApJ, 974, 295
- Doroshenko, V. 2024, A&amp;A, submitted, arXiv:2403.03127
- Doroshenko, V., Poutanen, J., Heyl, J., et al. 2023, A&amp;A, 677, A57
- Doroshenko, V., Poutanen, J., Tsygankov, S. S., et al. 2022, Nature Astronomy, 6, 1433
- Falanga, M., Kuiper, L., Poutanen, J., et al. 2011, A&amp;A, 529, A68
- Falanga, M., Kuiper, L., Poutanen, J., et al. 2012, A&amp;A, 545, A26
- Forsblom, S. V., Poutanen, J., Tsygankov, S. S., et al. 2023, ApJ, 947, L20
- Forsblom, S. V., Tsygankov, S. S., Poutanen, J., et al. 2024, A&amp;A, 691, A216
- Galloway, D. K., Johnston, Z., Goodwin, A., &amp; He, C.-C. 2022, ApJS, 263, 30
- González-Caniulef, D., Caiazzo, I., &amp; Heyl, J. 2023, MNRAS, 519, 5902
- Ibragimov, A. &amp; Poutanen, J. 2009, MNRAS, 400, 492
- Illiano, G., Coti Zelati, F., Marino, A., et al. 2024, The Astronomer's Telegram, 16510, 1
- Illiano, G., Papitto, A., Sanna, A., et al. 2023, ApJ, 942, L40
- Kislat, F., Clark, B., Beilicke, M., &amp; Krawczynski, H. 2015, Astroparticle Physics, 68, 45
- Lamb, F. K., Boutloukos, S., Van Wassenhove, S., et al. 2009, ApJ, 706, 417

Lampton, M., Margon, B., &amp; Bowyer, S. 1976, ApJ, 208, 177 Lapidus, I. I. &amp; Sunyaev, R. A. 1985, MNRAS, 217, 291 Lattimer, J. M. &amp; Prakash, M. 2016, Phys. Rep., 621, 127 Li, Z., De Falco, V ., Falanga, M., et al. 2018, A&amp;A, 620, A114 Lightman, A. P. &amp; Zdziarski, A. A. 1987, ApJ, 319, 643

- Linares, M., Altamirano, D., Chakrabarty, D., Cumming, A., &amp; Keek, L. 2012, ApJ, 748, 82
- Lo, K. H., Miller, M. C., Bhattacharyya, S., &amp; Lamb, F. K. 2013, ApJ, 776, 19 Loktev, V., Salmi, T., Nättilä, J., &amp; Poutanen, J. 2020, A&amp;A, 643, A84
- Malacaria, C., Heyl, J., Doroshenko, V., et al. 2023, A&amp;A, 675, A29 Mariani, I., Motta, S., Baglio, M. C., et al. 2024, The Astronomer's Telegram,
- 16475, 1
- Marshall, H. L., Ng, M., Rogantini, D., et al. 2022, ApJ, 940, 70
- Matsuoka, M., Kawasaki, K., Ueno, S., et al. 2009, PASJ, 61, 999

Article number, page 9 of 9

- Meszaros, P., Novick, R., Szentgyorgyi, A., Chanan, G. A., &amp; Weisskopf, M. C. 1988, ApJ, 324, 1056
- Miller, M. C. &amp; Lamb, F. K. 2015, ApJ, 808, 31
- Miller, M. C. &amp; Lamb, F. K. 2016, European Physical Journal A, 52, 63
- Miller, M. C., Lamb, F. K., Dittmann, A. J., et al. 2019, ApJ, 887, L24
- Miller, M. C., Lamb, F. K., Dittmann, A. J., et al. 2021, ApJ, 918, L28
- Molkov, S. V., Lutovinov, A. A., Tsygankov, S. S., et al. 2024, A&amp;A, 690, A353
- Morsink, S. M. &amp; Leahy, D. A. 2011, ApJ, 726, 56
- Motta, S., D'Aì, A., Papitto, A., et al. 2011, MNRAS, 414, 1508
- Mushtukov, A. A., Tsygankov, S. S., Poutanen, J., et al. 2023, MNRAS, 524, 2004
- Nagirner, D. I. &amp; Poutanen, J. 1993, A&amp;A, 275, 325
- Nagirner, D. I. &amp; Poutanen, J. 1994, Astrophysics and Space Physics Reviews, 9, 1
- Ng, M., Ray, P. S., Sanna, A., et al. 2024, ApJ, 968, L7

Özel, F. &amp; Freire, P. 2016, ARA&amp;A, 54, 401

- Patruno, A. &amp; Watts, A. L. 2021, in ASSL, Vol. 461, Timing Neutron Stars: Pulsations, Oscillations and Explosions, ed. T. M. Belloni, M. Méndez, &amp; C. Zhang, 143-208
- Poutanen, J. 2006, Advances in Space Research, 38, 2697
- Poutanen, J. 2020, A&amp;A, 641, A166
- Poutanen, J. &amp; Beloborodov, A. M. 2006, MNRAS, 373, 836
- Poutanen, J. &amp; Gierli´ski, M. 2003, MNRAS, 343, 1301 n
- Poutanen, J., Ibragimov, A., &amp; Annala, M. 2009, ApJ, 706, L129
- Poutanen, J. &amp; Svensson, R. 1996, ApJ, 470, 249
- Poutanen, J., Tsygankov, S. S., Doroshenko, V., et al. 2024, A&amp;A, 691, A123
- Radhakrishnan, V. &amp; Cooke, D. J. 1969, Astrophys. Lett., 3, 225
- Riley, T. E., Watts, A. L., Bogdanov, S., et al. 2019, ApJ, 887, L21
- Riley, T. E., Watts, A. L., Ray, P. S., et al. 2021, ApJ, 918, L27
- Russell, T. D., Carotenuto, F., van den Eijnden, J., et al. 2024, The Astronomer's Telegram, 16511, 1
- Salmi, T., Choudhury, D., Kini, Y., et al. 2024, ApJ, 974, 294
- Salmi, T., Loktev, V ., Korsman, K., et al. 2021, A&amp;A, 646, A23
- Salmi, T., Nättilä, J., &amp; Poutanen, J. 2018, A&amp;A, 618, A161
- Suleimanov, V. F., Forsblom, S. V., Tsygankov, S. S., et al. 2023, A&amp;A, 678, A119
- Tsygankov, S. S., Doroshenko, V., Mushtukov, A. A., et al. 2023, A&amp;A, 675, A48
- Tsygankov, S. S., Doroshenko, V., Poutanen, J., et al. 2022, ApJ, 941, L14 Viironen, K. &amp; Poutanen, J. 2004, A&amp;A, 426, 985
- Vinciguerra, S., Salmi, T., Watts, A. L., et al. 2024, ApJ, 961, 62
- Watts, A. L., Andersson, N., Chakrabarty, D., et al. 2016, Reviews of Modern Physics, 88, 021001
- Weisskopf, M. C., So ffi tta, P., Baldini, L., et al. 2022, JATIS, 8, 026002 Wilms, J., Allen, A., &amp; McCray, R. 2000, ApJ, 542, 914
Alessandro Papitto Alessandro Di Marco Juri Poutanen Tuomo Salmi Giulia Illiano Fabio La Monaca Filippo Ambrosino Anna Bobrikova Maria Cristina Baglio Caterina Ballocco Luciano Burderi Sergio Campana Francesco Coti Zelati Tiziana Di Salvo Riccardo La Placa Vladislav Loktev Sinan Long Christian Malacaria Arianna Miraval Zanon Mason Ng Maura Pilia Andrea Sanna Luigi Stella Tod Strohmayer Silvia Zane
Discovery of Polarized X-Ray Emission from the Accreting Millisecond Pulsar SRGA J144459.2-604207
We report on the discovery of polarized X-ray emission from an accreting
millisecond pulsar. During a 10-day-long coverage of the February 2024 outburst
of SRGA J144459.2-604207, the Imaging X-ray Polarimetry Explorer (IXPE)
detected an average polarization degree of the 2-8 keV emission of 2.3% +/-
0.4% at an angle of 59{\deg} +/- 6{\deg} (East of North; uncertainties quoted
at the 1${\sigma}$ confidence level). The polarized signal shows a significant
energy dependence with a degree of 4.0% +/- 0.5% between 3 and 6 keV and < 1.5%
(90% c.l.) in the 2-3 keV range. We used NICER, XMM-Newton, and NuSTAR
observations to obtain an accurate pulse timing solution and perform a
phase-resolved polarimetric analysis of IXPE data. We did not detect any
significant variability of the Stokes parameters Q and U with the spin and the
orbital phases. We used the relativistic rotating vector model to show that a
moderately fan-beam emission from two point-like spots at a small magnetic
obliquity ($\simeq$ 10{\deg}) is compatible with the observed pulse profile and
polarization properties. IXPE also detected 52 type-I X-ray bursts, with a
recurrence time $\Delta t_{rec}$ increasing from 2 to 8 h as a function of the
observed count rate $C$ as as $\Delta t_{rec} \simeq C^{-0.8}$ We stacked the
emission observed during all the bursts and obtained an upper limit on the
polarization degree of 8.5% (90% c.l.).