The IXPE View of Neutron Star Low-Mass X-ray Binaries

Low-mass X-ray binaries hosting weakly magnetized neutron stars (NS-LMXBs) are among the brightest sources in the X-ray sky. Since 2021, the Imaging X-ray Polarimetry Explorer (IXPE) has provided new measurements of the X-ray polarization of these sources. IXPE observations have revealed that most NS-LMXBs are significantly polarized in the X-rays, providing unprecedented insight into the geometry of their accretion flow. In this review paper, we summarize the first results obtained by IXPE on NS-LMXBs, the emerging trends within each class of sources (atoll/Z), and possible physical interpretations.


Introduction
Accreting, weakly magnetized NS-LMXBs are very bright X-ray sources and represent excellent laboratories for studying the physics of accretion on compact objects.They are divided into two main classes, Z and atoll, according to the pattern they trace on X-ray colorcolor diagrams [1].Atoll sources are less luminous and have lower accretion rates compared with Z sources.The X-ray spectrum of weakly magnetized NS-LMXBs is generally well described by a soft thermal component and a hard Comptonization component, but their physical origin is still uncertain.In principle, the soft component could be due either to the accretion disk or to the neutron star, while the hard component could originate in a hot corona, in a boundary layer (BL) between the disk and the NS [2], or a more extended spreading layer (SL) around the NS [3,4].Variability studies showed that the spectrum of the most rapidly variable component is hard, and is likely associated with the BL rather than the disk [5,6].Furthermore, the primary X-ray emission can be reprocessed by the accretion disk, resulting in a reflection component, which typically consists of a Compton hump at 20-30 keV and fluorescent emission lines (see [7] and the references therein).
X-ray spectroscopy and timing techniques allow us to study the physical properties of NS-LMXBs in depth (see [8] for a recent review), but their geometry remains elusive.X-ray polarimetry crucially adds two more observables, namely the polarization degree (PD) and angle (PA), which depend on the geometry of the accretion flow.By measuring X-ray polarization, we can distinguish between radially extended regions such as the accretion disk, and vertically or spherically extended structures such as the BL or SL [9,10].Radiation reflected off the disk is also expected to be significantly polarized [11].
Before the launch of the Imaging X-ray Polarimetry Explorer (IXPE) [12,13], the Xray polarization properties of NS-LMXBs were almost unknown.The only source with a significant measurement was Sco X-1, from a 15-day observation with OSO-8 [14] and, more recently, a 322-day observation with PolarLight [15].A marginal detection of Cyg X-2 was also obtained by OSO-8 [16].As of December 2021, IXPE has re-opened the X-ray polarimetric window.Thanks to its greatly improved sensitivity, IXPE has measured the X-ray polarization in the 2-8 keV band of several NS-LMXBs, with only ∼1-day exposures.Moreover, IXPE is capable of performing energy-and time-resolved spectropolarimetry.
This paper summarizes the first X-ray spectro-polarimetric results obtained with IXPE on NS-LMXBs during the first two-year observational campaign.Sections 2 and 3 focus on atoll and Z sources, respectively.Section 4 presents the most peculiar objects.Section 5 concludes our review with a discussion of lessons learned and open questions.
1.4% ± 0.3% strongly variable 3  [27] 1 However, a PD of 6% is attributed to the hard component. 2The PA, however, varies with time. 3The PD varies both with time and energy (see text for details).
GS 1826−238 and GX 9+9 were observed by IXPE for 100 ks each in March and October 2022, respectively.Both were jointly observed by NICER, NuSTAR, and/or INTEGRAL, ensuring broad-band spectral coverage [17,18].The two sources show similar spectral properties, the most notable difference being the presence of reflection features: these are seen in GX 9+9, but not in GS 1826−238 [17,18] (see Figure 1a).In GS 1826−238, IXPE finds only an upper limit of PD < 1.3%, slightly lower than the PD measured for GX 9+9 [17].The polarization of GX 9+9 also shows a hint of an increasing trend with energy, reaching 2-3% in the 4-8 keV band [18,19] (Figure 1b).These results can be explained by the combination of Comptonization in a BL and reflection [18].In fact, the BL emission is not expected to be highly polarized [9,10], while the reflected radiation can provide a significant contribution to polarization at high energies [11,28].The higher polarization of GX 9+9 compared with GS 1826−238 can, thus, be related to the stronger reflection component [29].
In all sources, a low polarization is generally found below 3-4 keV, where the thermal emission dominates.So far, 4U 1820−303, a bright ultracompact NS-LMXB, offers the most dramatic example of this behavior [20].This source was observed by IXPE twice, in October 2022 (16 ks) and April 2023 (86 ks), in the same spectral state.The PD strongly increases in energy, from < 1% below 4 keV up to the surprisingly high value of 10% in the 7-8 keV band [20].The measurements of PA indicate that the soft and hard components are polarized orthogonally to each other [20] (Figure 1c).4U 1624−49 is another intriguing case.It was observed by IXPE in August 2023 for 200 ks, and it shows a PD of 3% in the 2-8 keV band, with an indication for an increase with energy, up to 6% in the 6-8 keV band [21] (Figure 1d).4U 1624−49 is a dipping source seen at a high inclination, but even in this case, a BL alone cannot explain the high PD [9,10].Possible scenarios to explain the high polarization include the presence of significant reflection and/or a radially extended X-ray emitting region, such as a corona [21] or scattering in a wind [30].Another possible scenario is that Comptonization takes place in an outfowing plasma with a mildly relativistic velocity, which can produce a PD up to ∼ 10% [31].This could, in principle, explain the high PD of 4U 1820−303, however, the question of why the PD is high only in the 7-8 keV band is still difficult to solve [20].As for 4U 1624−49, the spectropolarimetric fit favors a model including soft disk emission, a hard Comptonized component from a BL/SL, and disk reflection [32].[18]).(b) PD−PA contour plots of GX 9+9 (68, 90, and 99% confidence levels; adapted from [18]), (c) 4U 1820−303 (50, 95, and 99% confidence levels; adapted from [20]), and (d) 4U 1624−49 (68 and 95% confidence levels; adapted from [21]).
IXPE observations enabled polarization measurements along the whole Z track, typically consisting of a Horizontal Branch (HB), a Normal Branch (NB) and a Flaring Branch (FB); HB is the hardest spectral state, while FB is the softest (see Figure 2a).IXPE revealed that the polarization is correlated with the position along the Z track: PD is high (up to 4-5%) for HB [23,24], and low (1-2%) for NB and FB [22][23][24][25].This correlation was first observed in XTE J1701−462, a transient NS-LMXB that was targeted by IXPE twice for 50 ks in September/October 2022, three weeks after the onset of an outburst.During the first pointing, the source was in HB and exhibited a PD of 4.6% in the 2-8 keV band [23,33].Ten days later, the source was observed in NB and the PD in the 2-8 keV band had dropped to < 1.5% [23] (Figure 2b).A similar behavior was observed also in GX 5−1, which was targeted by IXPE twice for 50 ks, in March and April 2023.As in the case of XTE J1701−462, the energy-integrated PD is significantly higher when the source is in the HB (Figure 2c).GX 5−1 also displays a puzzling behavior for PA.In fact, the PA in the 2-3 keV band shows a ∼ 20 • difference compared with the PA at higher energies [24].Because in systems hosting a NS radiative processes occur at distances of at least 3 Schwarzschild radii, energy-dependent PA behavior would be difficult to attribute to general and special relativistic effects (which actually become effective at smaller distances), but must originate in geometrical configurations where components are not strictly orthogonal each other.Disk precession or warping with respect to the NS spin axis is a very intriguing explanation, albeit reproducing such effects via detailed Monte Carlo simulation is a challenging task.
IXPE measured a lower polarization degree (1-2%) in Cyg X-2 and Sco X-1, both observed in the soft branches only (NB and FB, respectively).These two sources are among the few NS-LMXBs known to host resolved radio jets [34,35].This enables the comparison between the X-ray PA and the direction of the radio jet, which is assumed to be the same as the symmetry axis of the system.Cyg X-2 was targeted by IXPE in April/May 2022 for 76 ks.Remarkably, in this source, PA is consistent with the orientation of the radio jet, which is expected to be parallel to the disk axis [22].The PA is also consistent with the previous measurement by OSO-8 [16], but with greater significance.Additionally, the spectropolarimetric fit provided a hint of a 90 • swap between the PA of the soft and hard components, the latter being aligned with the radio jet.These results were interpreted in the framework of a spectropolarimetric model, consisting of a standard disk for the soft component, and Comptonization in a BL/SL for the hard component [22].The PA of the soft component is indeed consistent with direct disk emission, while the PA of the hard component indicates that Comptonization takes place in a vertically extended BL/SL rather than in a region coplanar with the disk [22].Disk reflection can also be significantly polarized; unfortunately, this component was unconstrained owing to the limited statistic [22].However, as Comptonization in a BL/SL is expected to produce a low PD (not exceeding 1%), it is reasonable to infer that disk reflection contributes to the net polarization signal in the IXPE band [10].
Sco X-1, the brightest persistent source in the X-ray sky, was observed by IXPE in August 2023 for 24 ks, while it was mostly in the Soft Apex (SA) that connects NB to FB.The previous measurement with PolarLight produced a PD of 4% and PA parallel to the radio jet, at least in the high flux state, and a non-detection during a low flux state [15].IXPE instead measured a PD of 1%, with no evidence of a significant trend with energy, and a PA rotated by 40 • -50 • with respect to the jet [25].The most likely explanation for such PA misalignment is the presence of a preceding disk [25].Following the trend observed in XTE J1701−462 and GX 5−1, we can reasonably infer that the higher PD previously measured by PolarLight was due to the source being predominantly in HB, albeit the relatively narrow band of PolarLight did not allow for branch identification [15].[24]).In the latter plot, solid and dashed lines indicate contours obtained with XSPEC and IXPEOBSSIM [36], respectively.

Peculiar Sources
IXPE observed two peculiar objects that exhibit properties of both atoll and Z sources, namely Cir X-1 [26] and GX 13+1 [27].Cir X-1 was targeted during August 2023 in two pointings separated by six days, to cover two different parts of the orbit, for a total exposure time of 263 ks.GX 13+1 was observed continuously for 100 ks in October 2023.One of the most intriguing results is the observation of a significant variation of PA with time [26,27].Along the orbit of Cir X-1, the polarization is well constrained in the high-flux, soft phases (P2 and P3 in Figure 3a).The PA is found to change by ∼50 • between different phase intervals, while PD is 1-2% and is consistent with the remaining constant [26] (Figure 3b).As for GX 13+1, PA is found to rotate by ∼70 • during the two-day IXPE exposure; in this case, PD also varies in time, ranging from non-detectable during a dip in flux seen in the IXPE light curve, and up to 5% in the phase after the dip [27].Differently from Cir X-1, these variations in polarization properties are not related to clearly separate spectral states [27].These results challenge theoretical models, however, they broadly indicate the coexistence of two variable components with misaligned symmetry axes [26,27].[26]).(b) PD−PA contour plots of Cir X-1 (68, 95, and 99% confidence levels; adapted from [26]).

Summary
The IXPE satellite enabled a big step forward in our knowledge of accretion processes in weakly magnetized NS-LMXBs.In particular, it revealed that the hard emission component dominates the polarization signal, while the soft emission has a lower polarization [18,20,22,23,25].This trend is apparent in both atoll and Z sources.The spectropolarimetric data are generally well described by the combination of soft, thermal emission from the accretion disk and Comptonization in a BL or SL, with a possible contribution from reflected radiation [17,18,22,24].This geometry naturally produces a PA parallel to the jet [22].On the other hand, the unexpectedly high PD measured in some sources is not easy to explain; scattering in an extended corona, or perhaps in a wind, seen at high inclination, cannot be ruled out [20,21].Another puzzling result is the energy and time dependence of PA observed in some sources, challenging simple models based on axially symmetric geometries [24][25][26][27].A misalignment between the NS spin and the orbital axis has been suggested as a possible scenario [26,27].In any case, the IXPE polarimetric measurements provide constraints to the accretion geometry and represent a powerful test for theoretical predictions [9,10,30], paving the way for further developments.

igure 1 .
(a) IXPE and NuSTAR spectral data of GX 9+9 fitted with a model consisting of disk emission, Comptonization, and reflection (adapted from

Figure 2 .
Figure 2. (a) Typical Z track in the hardness−intensity diagram of GX 5−1, based on NuSTAR data simultaneous to IXPE observations (adapted from [24]).(b) PD−PA (denoted here as P and Ψ) contour plots in the 2-8 keV band of XTE J1701−462, showing a correlation with the position along the Z track (adapted from[23]).(c) PD-PA contour plots of GX 5−1 in the 2-8 keV band, showing a similar behavior (adapted from[24]).In the latter plot, solid and dashed lines indicate contours obtained with XSPEC and IXPEOBSSIM[36], respectively.

Table 1 .
The polarization degree of the atoll (top part of the table), Z sources (middle part of the table), and peculiar sources (bottom part of the table).