---
title: 'GX 349+2: Sco-like Z Source in NS-LMXBs'
url: https://www.emergentmind.com/topics/gx-349-2
type: topic
---

# GX 349+2: Sco-like Z Source in NS-LMXBs

Searching arXiv for recent GX 349+2 papers to ground the article.
GX 349+2, also known as Sco X-2, is a bright Sco-like Z source in the class of neutron-star low-mass X-ray binaries (NS-LMXBs) accreting at near-Eddington luminosities. In hardness–intensity and color–color diagrams it traces predominantly the normal branch (NB) and flaring branch (FB), with a weak or absent horizontal branch (HB) in many observations, although HB intervals have been identified in at least one XMM-Newton study [1511.00599; 2602.05989]. Across recent work, GX 349+2 has emerged as a key laboratory for studying the coupling between the inner accretion disk, the neutron-star boundary or spreading layer, Comptonizing plasma, relativistic reflection, long and short X-ray lags, and weak but measurable X-ray polarization. These observables together indicate a near-surface accretion flow in which branch-dependent geometric and radiative changes occur on timescales from seconds to hundreds of seconds, while the source remains persistently luminous and spectrally complex [1708.01652; 2505.00813; 2602.05989].

## 1. Source class, nomenclature, and phenomenology

GX 349+2 is classified as a Sco-like Z source, that is, a weakly magnetized NS-LMXB accreting at near-Eddington rates [1511.00599; 2504.04601]. In this class, motion along the Z track is conventionally described in terms of the HB, NB, and FB, though GX 349+2 often shows a long NB and elongated FB with little or no HB in many campaigns [1511.00599; 1708.01652]. The source is therefore frequently compared with Sco X-1 rather than with Cyg-like Z sources.

Several studies adopt a distance of 9.2 kpc for converting fluxes to luminosities and emitting radii [2602.05989; 2505.00813; 2507.07163]. At this distance, the observed X-ray fluxes place GX 349+2 at luminosities of order \(10^{38}\,\mathrm{erg\,s^{-1}}\), consistent with its identification as a near-Eddington system [2505.00813]. NuSTAR observations covering the full Z track inferred Eddington fractions rising from \(\sim 0.70\) in the NB or vertex to \(\sim 1.41\) in the brightest FB state, using 3–30 keV fluxes and a canonical \(1.4\,M_\odot\) neutron star [1708.01652].

The source has long been known to exhibit strong flaring activity, broad iron-line emission, and characteristic Z-source timing behavior. Recent work extends that picture by connecting branch-resolved timing, reflection, and polarimetry to the structure of the innermost accretion flow. This suggests that GX 349+2 is especially informative for testing how the disk, boundary layer, corona, and reflection region reorganize as the source moves along the Z track [1708.01652; 2507.00302; 2602.05989].

## 2. Z-track morphology and observational campaigns

A systematic RXTE study analyzed 138 observations from 1996 to 2011 and found that the most extensive Z track was traced during 1998 January 9–29 [1511.00599]. In that epoch, the hardness–intensity diagram was partitioned into 23 regions, with 8 on the NB and 15 on the FB; no HB was observed [1511.00599]. The HID used hardness \(=\) count-rate ratio \((8.7\text{–}19.7~\mathrm{keV})/(6.2\text{–}8.7~\mathrm{keV})\) and intensity \(=\) count rate in \(2.0\text{–}19.7~\mathrm{keV}\) [1511.00599]. The source executed repeated FB-to-NB cycles during that interval.

A two-day NuSTAR observation beginning 2016-06-06 (ObsID 30201026002) covered approximately the full Z track and selected five spectral regions: NB, the NB–FB vertex (VX, the soft apex), FB1, FB2, and FB3 [1708.01652]. Hardness was defined as \((10\text{–}16~\mathrm{keV})/(6.4\text{–}10~\mathrm{keV})\), and no HB was detected, consistent with the Sco-like classification [1708.01652]. This campaign was central in establishing that the broad Fe K\(\alpha\) line persists across states.

An XMM-Newton EPIC-pn Timing-mode observation on 2008-03-19 (ObsID 0506110101) provided a different view by isolating eleven \(\approx 2000\) s timing segments and three \(\approx 2200\) s transition windows, using soft and hard bands of \(0.8\text{–}2.0\) keV and \(2.0\text{–}10.0\) keV, respectively [2602.05989]. In that dataset, two segments in the HID locus were identified as HB, while the remaining segments were assigned to NB or FB. This is significant because it enabled a direct comparison of cross-correlation morphology among all three canonical branches within a single observation [2602.05989].

The 2024 IXPE campaign, combined with NuSTAR and in some analyses Swift/XRT, added a polarization-based state diagnostic. One model-independent analysis found GX 349+2 mainly in the NB with brief FB excursions and no HB during the IXPE observation [2507.00302]. Another IXPE+NuSTAR analysis described the source as predominantly in NB during the first NuSTAR pointing and mainly in NB with short excursions into the soft apex (SA) during the second, with no deep FB coverage [2507.07163]. These differences reflect distinct state-selection strategies and are relevant when comparing branch-resolved polarization claims.

## 3. Spectral structure: disk, boundary layer, Comptonization, and iron-line reflection

Broadband and medium-band spectroscopy consistently require multiple continuum components. In IXPE+NuSTAR spectro-polarimetric modeling, the continuum was represented by an absorbed combination of disk blackbody, hotter blackbody, and Comptonization, plus a line or reflection component [2505.00813]. Specifically, one study used
\[
\mathrm{tbabs}\times(\mathrm{bbodyrad}+\mathrm{diskbb}+\mathrm{diskline}+\mathrm{nthcomp}),
\]
interpreting diskbb as the disk, bbodyrad as the neutron-star surface or boundary/spreading layer, nthcomp as thermal Comptonization, and diskline as the Fe-K reflection feature [2505.00813]. Branch-resolved fits gave \(kT_{\rm in}\simeq 0.62\) keV (SA), \(\simeq 0.70\) keV (FB), and \(\simeq 0.74\) keV (NB), while the hotter blackbody remained near \(kT_{\rm BB}\simeq 1.29\text{–}1.30\) keV, with an emitting radius \(\sim 15\text{–}19\) km at 9.2 kpc [2505.00813]. The Comptonized component had \(kT_e \simeq 2.42\text{–}2.62\) keV and \(\Gamma \simeq 1.0\text{–}1.33\), indicating a cool, optically thick scattering region [2505.00813].

A separate IXPE+NuSTAR+Swift/XRT study emphasized the “Western model,” in which soft disk photons are Comptonized in an optically thick, cool corona while a separate thermal component originates from the neutron-star surface or boundary layer [2504.04601]. Its preferred spectral model,
\[
\mathrm{const}\times \mathrm{tbabs}\times(\mathrm{compTT}+\mathrm{compbb}+\mathrm{Gaussian}),
\]
yielded \(N_H = (0.56 \pm 0.03)\times 10^{22}\,\mathrm{cm^{-2}}\), compTT seed temperature \(T_0 = 0.48 \pm 0.02\) keV, \(kT_e = 2.65 \pm 0.02\) keV, optical depth \(\tau = 6.96 \pm 0.16\), and a weak surface Comptonization layer with \(\tau = 0.05 \pm 0.01\) [2504.04601]. The corresponding Compton \(y\)-parameter was reported as \(y \approx 0.94\text{–}1.00\) for compTT and \(y \approx 0.02\) for compbb, supporting a picture of nearly saturated scattering in a slab-like configuration and negligible Comptonization in the surface layer [2504.04601].

NuSTAR spectroscopy across the Z track demonstrated that simple two-thermal-component fits are inadequate because of a broad asymmetric Fe K\(\alpha\) residual in all states [1708.01652]. Comptonized continua improve the fit but do not remove the need for reflection [1708.01652]. Full reflection modeling with relconv \(\otimes\) reflionx\(_{\rm BB}\) plus diskbb and bbodyrad showed that the inner disk radius remains broadly stable across NB through FB2, with a weighted mean \(R_{\rm in}\approx 17.5^{+1.5}_{-1.0}R_g\), corresponding to \(\approx 36.4\) km for a canonical \(1.4\,M_\odot\) neutron star [1708.01652]. In that analysis, the Fe line persisted in all states and the source did not show strong evidence for large monotonic disk motion along the track.

The line centroid depends on the model and dataset. NuSTAR diskline fits gave \(E_{\rm line}\approx 6.40\text{–}6.50\) keV across the Z track [1708.01652], whereas IXPE+NuSTAR spectro-polarimetry placed the line or diskline near \(6.66\text{–}6.70\) keV in some states [2505.00813], and the XMM-Newton transition study fixed Fe K\(\alpha\) at 6.7 or 6.76 keV depending on the continuum model [2602.05989]. A plausible implication is that line centroid differences partly reflect differing continuum prescriptions, ionization assumptions, and energy-band coverage rather than a single directly comparable physical measurement.

## 4. Timing and cross-correlation phenomenology

Cross-correlation analyses have revealed two distinct lag regimes in GX 349+2. The RXTE archival study used soft and hard bands of \(2\text{–}5\) keV and \(16\text{–}30\) keV, with 16 s bins, and computed discrete CCFs using XRONOS CROSSCOR [1511.00599]. Along the 1998 Z track, 18 of 23 regions showed positive correlations, 5 were ambiguous, and no anti-correlations were found [1511.00599]. Typical cross-correlation coefficients ranged from \(\sim 0.27\) to \(\sim 0.95\), and the corresponding positively correlated lags were short, usually \(|\tau| \lesssim 10\) s and generally within a few to a few tens of seconds, though one LNB region reached \(-36.4 \pm 9.2\) s [1511.00599]. Outside that canonical Z-track interval, anti-correlations appeared in 10 observations, positive correlations in 100, and ambiguous behavior in 5 [1511.00599]. The anti-correlated lags spanned tens to thousands of seconds, with examples from \(-6107 \pm 46\) s to \(+1177 \pm 57\) s [1511.00599].

The XMM-Newton study focused instead on \(0.8\text{–}2.0\) keV and \(2.0\text{–}10.0\) keV light curves binned at 50 s and computed a normalized CCF
\[
R_{xy}(\tau)=\frac{\sum_t (x(t)-\bar{x})(y(t+\tau)-\bar{y})}{(N-1)\sigma_x \sigma_y}.
\]
To quantify asymmetric peaks, each CCF was modeled with a Gaussian on a linear baseline, described as line + Gaussian + line, and the Gaussian centroid or peak provided the lag estimate [2602.05989]. This study found strongly branch-dependent behavior: HB segments showed markedly asymmetric CCFs with relatively low peak correlation coefficients and non-zero hard lags of order a few hundred seconds; NB and FB segments showed highly symmetric CCFs peaking at zero lag; and flux-transition windows exhibited asymmetric CCFs with lags of a few tens to a few hundreds of seconds [2602.05989].

The lag detections in the XMM-Newton work were tested with two independent Monte Carlo approaches. First, 10,000 data-driven pairs of synthetic light curves were generated by perturbing each data point within its measurement uncertainty while preserving variability; second, 10,000 Timmer–König red-noise simulations with \(\Gamma_{\rm PSD}=2.1\) and the observed mean were used to construct 95% confidence envelopes [2602.05989]. For all HB and transition intervals with non-zero lags, the observed CCFs lay outside the 95% confidence bands, and the lag distributions from the data-driven Monte Carlo were centered close to the measured values, leading to the conclusion that the lags are intrinsic at approximately 95% confidence [2602.05989].

These results are not identical to the RXTE phenomenology because the underlying selections differ. RXTE’s on-track analysis covered NB and FB only and emphasized short lags with positive correlations [1511.00599], while XMM-Newton isolated HB and transition intervals in which longer, asymmetric hard lags became visible [2602.05989]. This suggests that the lag behavior of GX 349+2 is highly state- and geometry-dependent, and that datasets lacking HB coverage may systematically miss the longer-lag regime.

## 5. Interpretations of lags and inner-flow dynamics

The long lags in GX 349+2 have been interpreted in two closely related but not identical ways. In the RXTE study, long-term lags in NS-LMXBs were argued to be inconsistent with short-lag mechanisms such as simple Comptonization or with black-hole-style truncated-disk models, and instead were attributed to an extended accretion-disk corona (ADC) [1511.00599]. In that framework, the corona extends to radii \(R \sim (2\text{–}70)\times 10^4\) km, Comptonizes disk seed photons to produce hard X-rays, and competes with viscous inflow to the neutron-star surface. Two timescales were emphasized: the diffusion or Comptonization escape time,
\[
t_{\rm diff}\approx \frac{(\tau+\tau^2)R}{c},
\]
and the viscous inflow time,
\[
t_{\rm visc}\sim \left(\frac{1}{\alpha}\right)\left(\frac{R}{H}\right)^2 \left(\frac{1}{\Omega_K}\right).
\]
Using \(\alpha \approx 0.1\) and \(H/R \approx 0.1\), the paper noted \(t_{\rm visc}\approx 4\times 10^2\) s at \(R=3\times 10^8\) cm and \(t_{\rm visc}\approx 2.3\times 10^3\) s at \(R=10^9\) cm, values compatible with the observed hundreds-to-thousands of seconds lags [1511.00599].

The XMM-Newton work instead located the relevant dynamics much closer to the star, in the boundary layer or corona near the inner edge of the disk [2602.05989]. It argued that hard lags of \(\sim 200\) to a few hundred seconds are too long for light-travel or Comptonization delays and are better explained by slow viscous readjustment of the boundary layer:
\[
t_{\rm visc} \sim \left(\frac{1}{\alpha}\right)\left(\frac{R}{H}\right)^2 \left(\frac{1}{\Omega_K}\right), \qquad
\Delta t \sim \left(\frac{3}{2}\right)\left(\frac{\Delta R}{R}\right)t_{\rm visc}.
\]
Adopting \(R \approx 15\) km, \(\Delta R \approx 2\) km, \(H/R \approx 0.1\), \(\Omega_K \approx 1.2\times 10^4\,\mathrm{s^{-1}}\), and \(\Delta t \approx 200\) s yields \(\alpha \approx 8\times 10^{-6}\), described as an extremely low effective viscosity [2602.05989]. The paper further proposed a depletion timescale for a low-viscosity boundary layer,
\[
t_{\rm dep} = \frac{\Omega-\Omega_{\rm NS}}{\alpha(\Omega_K^2-\Omega^2)},
\]
and showed that for \(\alpha \lesssim 10^{-6}\text{–}10^{-7}\) the resulting \(t_{\rm dep}\) values reproduce the tens-to-hundreds of seconds lags measured in the HB and during transitions [2602.05989].

The same study estimated the radial extent of the boundary layer from the Popham–Sunyaev scaling
\[
\log(R_{\rm BL}-R_{\rm NS}) \approx 5.02 + 0.245 \log(\dot{M}/10^{-9.85}M_\odot\,\mathrm{yr}^{-1}),
\]
obtaining \(R_{\rm BL}\approx 20\text{–}34\) km for one continuum parameterization and \(R_{\rm BL}\approx 13\text{–}15\) km for another [2602.05989]. It also noted that prior XMM-Newton work had suggested a large ADC spanning tens to hundreds of km, and showed that for a thick flow, \(\Delta R \approx \alpha(H/R)^2 v_K \Delta t\) gives \(\Delta R \sim 30\text{–}300\) km for \(H/R \approx 0.1\), \(\Delta t \approx 200\) s, and \(\alpha \approx 10^{-3}\text{–}10^{-4}\), consistent with an extended corona dynamically coupled to the boundary layer [2602.05989].

Taken together, these interpretations are not mutually exclusive. The RXTE results favor an extended ADC for long-term lags outside the canonical Z track [1511.00599], whereas the XMM-Newton study attributes branch-dependent HB and transition lags to slow mechanical readjustment or depletion in a low-viscosity boundary layer or inner coronal region [2602.05989]. A plausible implication is that different lag classes in GX 349+2 probe different radial zones of the accretion flow.

## 6. X-ray polarization and accretion geometry

IXPE established that GX 349+2 is weakly but significantly polarized in X-rays. A joint IXPE+NuSTAR study reported a full-exposure \(2\text{–}8\) keV polarization degree \(PD = 1.1 \pm 0.3\%\) and polarization angle \(PA = 32^\circ \pm 6^\circ\), with \(3.9\sigma\) significance [2505.00813]. The same work found energy-resolved values of \(PD = 1.1 \pm 0.2\%\) at \(PA = 17^\circ \pm 6^\circ\) in \(2\text{–}4\) keV, \(PD = 0.9 \pm 0.4\%\) at \(PA = 49^\circ \pm 12^\circ\) in \(4\text{–}6\) keV, and \(PD = 3.1 \pm 1.1\%\) at \(PA = 46^\circ \pm 9^\circ\) in \(6\text{–}8\) keV, suggesting a possible increase near the Fe-K band [2505.00813].

A model-independent survey of IXPE Z-source observations gave for GX 349+2 a whole-observation \(2\text{–}8\) keV polarization of \(PD = 1.1\% \pm 0.3\%\) and \(PA = 30^\circ \pm 7^\circ\), with \(3.8\sigma\) significance [2507.00302]. In that analysis, branch-resolved measurements yielded \(PD = 1.0\% \pm 0.3\%\), \(PA = 38^\circ \pm 8^\circ\) in the NB and \(PD = 2.1\% \pm 0.6\%\), \(PA = 7^\circ \pm 8^\circ\) in the FB [2507.00302]. The PD increase from NB to FB was stated to be consistent within uncertainties at 90% confidence, whereas the PA rotation of approximately \(30^\circ\) was reported as significant [2507.00302]. This is one of the clearest branch-dependent PA rotations reported among Z sources.

A different IXPE+NuSTAR+Swift/XRT analysis, also with branch-resolved selection, reported significant \(2\text{–}8\) keV polarization in both FB and NB: FB \(PD = 1.74 \pm 0.52\%\), \(PA = 19.4 \pm 8.9^\circ\); NB \(PD = 0.80 \pm 0.22\%\), \(PA = 35.4 \pm 7.9^\circ\) [2504.04601]. In the NB it further found an energy-dependent increase from \(PD = 0.78 \pm 0.22\%\), \(PA = 17.9 \pm 8.1^\circ\) in \(2\text{–}4\) keV to \(PD = 1.32 \pm 0.40\%\), \(PA = 53.2 \pm 8.6^\circ\) in \(4\text{–}8\) keV, implying a PA rotation of approximately \(35^\circ\text{–}40^\circ\) [2504.04601].

These low net polarization degrees are interpreted through Stokes-vector addition of multiple emission components:
\[
Q = \sum_i p_i I_i \cos(2\psi_i), \qquad
U = \sum_i p_i I_i \sin(2\psi_i), \qquad
PD_{\rm tot} = \frac{\sqrt{Q^2+U^2}}{I_{\rm tot}}.
\]
Because disk, boundary/spreading layer, Comptonization, and reflection may have different polarization angles, the net \(PD\) can remain near \(1\%\) even when individual components are polarized at a few percent [2505.00813]. This multi-component dilution is a central result of the IXPE literature on GX 349+2.

The most detailed normal-branch IXPE+NuSTAR spectropolarimetric study concluded that the spectrum is dominated below \(\sim 7\) keV by the disk and above \(\sim 7\) keV by a compact hotter layer, with energy-dependent polarization strengthening toward high energy [2507.07163]. In a component-resolved fit, the disk had \(P = 1.0\% \pm 0.7\%\) and \(\psi = 6^\circ(+24^\circ,-23^\circ)\), while the compact Comptonizing region had \(P = 3.7\% \pm 3.0\%\) and \(\psi = 68^\circ \pm 26^\circ\), with an angular difference of about \(60^\circ \pm 35^\circ\) [2507.07163]. That work stated that the geometry is slightly more favorable to a spreading layer than to a classic boundary layer, although it also emphasized model degeneracies and the difficulty of disentangling reflection polarization [2507.07163].

## 7. Synthesis, limitations, and outstanding issues

Across spectroscopy, timing, and polarimetry, GX 349+2 is consistently described as a near-Eddington neutron-star accretor in which the inner disk reaches close to the star while a substantial boundary or spreading layer and a cool Comptonizing region shape the X-ray output [1708.01652; 2505.00813; 2602.05989]. The disk does not appear to undergo large monotonic radial excursions through most of the Z track in the NuSTAR reflection analysis, remaining around \(R_{\rm in}\approx 17.5R_g\) on average [1708.01652]. At the same time, XMM-Newton transition spectra requiring a hot disk component with \(kT_{\rm in}\approx 1.0\text{–}1.5\) keV and relatively low \(N_{\rm dbb}\), together with prior detection of twin kHz QPOs at 712 and 978 Hz, were argued to support an inner flow extending very close to the neutron-star surface or ISCO [2602.05989]. These are not strictly identical radius inferences, and the literature treats them with appropriate caution.

A recurrent point of agreement is that branch changes are driven primarily by reorganization of the inner thermal and Comptonizing regions rather than by wholesale movement of the outer observable disk edge [1708.01652; 2602.05989]. In the timing domain, NB and FB often show tightly coupled soft–hard variability with symmetric or near-zero-lag CCFs, whereas HB segments and transitions reveal asymmetric CCFs with lags of tens to hundreds of seconds [2602.05989]. In polarimetry, GX 349+2 remains weakly polarized overall, but shows evidence for energy-dependent and branch-dependent PA rotation, implying that the dominant polarized component changes with spectral state and photon energy [2504.04601; 2507.00302].

Several limitations remain explicit in the literature. The XMM-Newton lag study is based on a single \(0.8\text{–}10\) keV observation with 50 s timing bins, and exact lag values and uncertainties are not tabulated despite Monte Carlo validation at approximately 95% confidence [2602.05989]. IXPE branch-resolved results depend strongly on state selection, exposure per branch, and whether brief SA or FB excursions are grouped with NB-dominated intervals [2505.00813; 2507.00302; 2507.07163]. Reflection constraints in some flaring states are weakened by high ionization, degeneracies with continuum shape, and spectral mixing across multiple flares [1708.01652]. More complex reflection models than simple diskline were attempted in polarimetric studies but often remained unconstrained at current signal-to-noise [2505.00813].

The present evidence nonetheless supports a coherent picture. GX 349+2 consists of a bright inner disk, a neutron-star boundary or spreading layer with characteristic radii from a few km to tens of km depending on model and definition, an optically thick low-temperature Comptonizing medium, and a persistent Fe-K reflection signature. Its lag behavior implies both fast coupled variability and slower structural readjustments, and its polarization implies weak intrinsic net asymmetry after vector summation of several components. This suggests that GX 349+2 is best understood not through any single diagnostic but through the combined evolution of CCF morphology, continuum decomposition, reflection, and Stokes geometry along the Z track [1511.00599; 1708.01652; 2505.00813; 2602.05989].

Source: https://www.emergentmind.com/topics/gx-349-2