---
title: '4U 1907+09: Wind-Fed X-ray Pulsar in HMXB'
url: https://www.emergentmind.com/topics/4u-1907-09
type: topic
---

# 4U 1907+09: Wind-Fed X-ray Pulsar in HMXB

4U 1907+09 is an accretion-powered X-ray pulsar in a high-mass X-ray binary, generally described in the recent literature as a classical wind-fed supergiant system, though some source summaries preserve earlier or alternative donor classifications. It is distinguished by a long spin period of about \(440\)–\(444\) s, an eccentric \(\sim 8.38\)-day orbit, recurrent dips and flares, a fundamental cyclotron resonance scattering feature (CRSF) near \(17\)–\(19\) keV with a harmonic reported near \(36\)–\(40\) keV, and a genuine mid-infrared bow shock indicating runaway motion through the interstellar medium. Modern studies with Suzaku, RXTE, INTEGRAL, AstroSat, NuSTAR, XMM-Newton, Swift, and IXPE have made 4U 1907+09 a reference source for phase-resolved cyclotron-line spectroscopy, wind-fed accretion variability, accretion-regime transitions, and, most recently, X-ray polarimetry [2504.02791] [2507.17873] [1102.2437].

## 1. Source classification and binary parameters

The published literature represented here is not uniform in its source summary, but it consistently identifies 4U 1907+09 as a high-mass X-ray binary containing an accreting neutron star. The donor has been listed as B2 III–IV in one Suzaku overview, as O8–O9 Ia in later timing and variability studies, and as O9.5 Iab in the Spitzer bow-shock analysis. The orbital period is reported as approximately \(8.38\) d or, more precisely, \(8.3753\) d, with eccentricity \(e \simeq 0.28\). The pulsar was discovered at approximately \(437.5\) s and is now observed at periods near \(442\)–\(444\) s, consistent with long-term spin evolution containing extended spin-down episodes and torque reversals [1111.1800] [1212.3433] [1102.2437] [2507.17873].

| Quantity | Reported value(s) | Context |
|---|---|---|
| Donor classification | B2 III–IV; O8–O9 Ia; O9.5 Iab | Different source summaries |
| Orbital period | \(8.3753\) d; \(\sim 8.38\) d | Timing and orbital studies |
| Eccentricity | \(e \simeq 0.28\) | Timing, wind-accretion, polarimetry studies |
| Spin period | \(\sim 437.5\) s at discovery; \(441.2088(2)\) s; \(442.92 \pm 0.03\) s; \(443.73^{+0.03}_{-0.01}\) s; \(443.99(4)\) s | Historical and modern timing measurements |
| Distance | \(4\) kpc adopted; minimum \(\approx 5\) kpc; Gaia EDR3 \(d \simeq 1.9\) kpc | Different analyses adopt different distances |

Distance estimates are likewise non-uniform. The bow-shock study adopted \(d=4\) kpc as a compromise among earlier estimates, the RXTE/ROTSE timing work described the source as lying at a minimum distance of \(\approx 5\) kpc, and recent IXPE and NuSTAR analyses used the Gaia EDR3-based distance \(d \simeq 1.9\) kpc for luminosity estimates [1102.2437] [1212.3433] [2507.17873] [2504.02791]. This suggests that luminosities quoted in different studies should be compared only after accounting for the assumed distance.

## 2. Orbit, stellar wind, and runaway environment

The system is a wind-fed binary with strong orbital modulation. One RXTE-based analysis reported two phase-locked flares per orbit, separated by \(\approx 0.5\) in orbital phase. The same work found that dips are most frequent between orbital phases \(\approx 0.1\)–\(0.6\) and are not observed between phases \(\approx 0.7\)–\(0.8\), while the hydrogen column density \(N_{\mathrm H}\) is maximal just after periastron and minimal during phases with frequent dips. That anti-correlation led to the conclusion that the dips are not caused by increased absorption but by decreased accretion when the neutron star traverses low-density wind regions [1212.3433].

Soft-X-ray monitoring with XMM-Newton and Swift extended this picture by resolving both clumps and larger-scale structures. In that campaign, the \(0.3\)–\(10\) keV flux peaked sharply at orbital phases \(0.5\)–\(0.8\), near periastron at phase \(\approx 0.7\), while spectral hardening at phases \(0.8\)–\(1.0\) was driven by both an \(N_{\mathrm H}\) increase by a factor \(\approx 2\) and a flatter photon index. This was interpreted as evidence for enhanced obscuration by a trailing stream or wake crossing the line of sight [2205.03023].

A distinctive property of 4U 1907+09 is its genuine mid-infrared bow shock. Spitzer/MIPS \(24\,\mu\mathrm m\) images revealed a clear arcuate bow shock with apex at \(\simeq 0.7\) arcmin from the star; at \(d=4\) kpc this corresponds to \(R_{\rm obs} \simeq 0.8\) pc. A local astrometric solution yielded \(\mu_\alpha \cos\delta=-7.7\pm 6.3\) mas yr\(^{-1}\) and \(\mu_\delta=+1.5\pm 5.3\) mas yr\(^{-1}\), implying a peculiar transverse velocity \(v_{\rm tr}\simeq 160\pm115\) km s\(^{-1}\) at the adopted distance. The direction of motion is consistent with the symmetry axis of the bow shock and shows that the system is moving away from the Galactic plane [1102.2437].

The bow-shock interpretation uses the standard stand-off relation
\[
R_0 = \sqrt{\frac{\dot{M}\,v_w}{4\pi\,\rho\,v_*^2}} \, .
\]
With the parameters adopted in the Spitzer analysis, \(\dot{M} \simeq 6\times10^{-6}\,M_\odot\,\mathrm{yr^{-1}}\), \(v_w \simeq 1690\) km s\(^{-1}\), \(v_* \simeq 160\) km s\(^{-1}\), and \(R_{\rm obs}\simeq0.8\) pc, the ambient density was estimated as \(n_a \simeq 2\) cm\(^{-3}\) [1102.2437]. In a broader HMXB mid-infrared survey, 4U 1907+09 and Vela X-1 were identified as the only two genuine bow shocks among the surveyed objects, emphasizing the rarity of such systems [1811.06781].

## 3. Spin evolution, torque states, and timing noise

The spin history of 4U 1907+09 is unusually structured for a wind-fed pulsar. Earlier measurements established long-term spin-down from the first Tenma period of \(\approx 437.5\) s, followed by a reduced spin-down rate around 2001, a torque reversal to spin-up reported after 2004 May, and then a return to spin-down before mid-2007. An RXTE phase-connected timing solution over MJD \(54280\)–\(55600\) gave \(P=441.2088(2)\) s at epoch \(t_0=54467.61(6)\) MJD, \(\nu=2.2665002(8)\times10^{-3}\) Hz, and \(\dot{\nu}=-3.672(1)\times10^{-14}\) Hz s\(^{-1}\), consistent with the resumed long-term spin-down state [1106.5957].

Several works summarize the secular torque history numerically. The pre-1998 spin-down rate was quoted as \(\dot{\nu}=-3.54\times10^{-14}\) Hz s\(^{-1}\), the post-2004 spin-up as \(\dot{\nu}=+2.58\times10^{-14}\) Hz s\(^{-1}\), and the recent spin-down as \(\dot{\nu}=-3.59\times10^{-14}\) Hz s\(^{-1}\). Near \(P\approx441\) s these correspond to \(\dot{P}\approx +6.9\times10^{-9}\) s s\(^{-1}\), \(-5.0\times10^{-9}\) s s\(^{-1}\), and \(+7.0\times10^{-9}\) s s\(^{-1}\), with characteristic spin-change timescales of order \(2000\)–\(2800\) years [1212.3433].

Recent observations show that the source remains in a spin-down state. NuSTAR measured \(P=442.92\pm0.03\) s in 2018 and derived a long-term spin-down rate of \(0.1971(4)\) s yr\(^{-1}\) across 2004–2018. A later NuSTAR observation in 2024 found \(P=443.99(4)\) s during the on-state, and IXPE measured \(P=443.73^{+0.03}_{-0.01}\) s with \(\dot{P}=(5^{+3}_{-6})\times10^{-8}\) s s\(^{-1}\) [2211.06038] [2504.02791] [2507.17873].

The source also exhibits stochastic torque fluctuations. Using Deeter’s polynomial estimator method, one RXTE/INTEGRAL analysis found a flat power spectrum of \(\dot{\nu}\) fluctuations between \(f \simeq 1/1300\) d\(^{-1}\) and \(1/75\) d\(^{-1}\), corresponding to white noise in \(\dot{\nu}\) and a random walk in \(\nu\), with noise strength
\[
S = 1.27\times10^{-21}\ \mathrm{Hz\ s^{-2}} \, .
\]
A later proceedings contribution compared 4U 1907+09 with magnetars and GX 1+4 and qualitatively concluded that noise strength scales up with spin-down rate, though it did not tabulate a source-specific value for 4U 1907+09 [1106.5957] [1701.04332].

Timing phenomenology beyond the coherent pulse includes transient QPO claims. A \( \sim 18 \) s QPO had been reported during a 1-hour flare on 1996 February 23; a Bayesian reanalysis recovered a QPO candidate around \(17.9\) s with false-alarm probability \(3.6\times10^{-3}\), corresponding to \(2.91\sigma\), while a classical power-density-spectrum estimate gave \(3.54\sigma\) [1701.04332]. Another review notes the source as showing occasional \(\approx 65\) mHz quasi-periodic oscillations [2205.03023].

## 4. Broad-band X-ray spectrum and cyclotron lines

The X-ray continuum of 4U 1907+09 has been modeled with the standard phenomenological continua used for accreting pulsars. Suzaku studies tested highecut\(\times\)power law, NPEX, and compTT; for phase-resolved CRSF work the NPEX and compTT fits were preferred because they gave consistent line parameters. INTEGRAL spectra were fit with POWERLAW\(\times\)HIGHECUT or CUTOFFPL modified by multiplicative Gaussian absorption lines. NuSTAR analyses used highecut\(\times\)power law, CUTOFFPL, compTT, and the Becker–Wolff bulk-plus-thermal Comptonization model. XMM-Newton spectra below \(10\) keV were described by \(\mathrm{TBabs}\times\mathrm{pcfabs}\times\)power law plus a Gaussian Fe K\(\alpha\) line [1304.6252] [1309.0875] [2211.06038] [2205.03023].

The defining spectral features are the CRSFs. Across missions, the fundamental is consistently found near \(17\)–\(19\) keV, while a harmonic is reported near \(36\)–\(40\) keV. Representative measurements include Suzaku phase-averaged values \(E_1=17.96^{+0.20}_{-0.19}\) keV with NPEX and \(18.07\pm0.18\) keV with compTT, both using a Lorentzian optical-depth profile; AstroSat/LAXPC measured \(E_{\rm cyc}=18.5\pm0.2\) keV; INTEGRAL found \(E_{c1}=18.5\pm0.7\) keV and \(E_{c2}=38^{+8}_{-5}\) keV with a Gaussian optical-depth line profile; NuSTAR in 2018 measured \(17.29\pm0.009\) keV and \(38.14\pm0.037\) keV; and NuSTAR in 2024 measured \(17.60^{+0.12}_{-0.24}\) keV and \(38.02^{+0.89}_{-0.82}\) keV [1304.6252] [1906.02917] [1309.0875] [2211.06038] [2504.02791].

Not every dataset detected the harmonic. The 2007 Suzaku phase-resolved analysis stated that the harmonic near \(\sim36\) keV was not detected with PIN, consistent with limited statistics in that band, and the LAXPC study likewise reported no evidence for a harmonic in the available statistics. These non-detections therefore coexist with secure detections in other broad-band observations [1304.6252] [1906.02917].

The CRSFs imply a surface magnetic field of order \(2\times10^{12}\) G through
\[
E_{\mathrm{cyc}} \approx 11.6\,\mathrm{keV}\,\frac{B_{12}}{1+z},
\qquad B_{12}=\frac{B}{10^{12}\,\mathrm G}.
\]
Using \(E_{\mathrm{cyc}}\approx18\) keV and a canonical neutron-star redshift \(1+z\approx1.30\), Suzaku and AstroSat analyses obtained \(B\approx 2.0\times10^{12}\) G, while the 2024 NuSTAR study quoted \(B \simeq 2\times10^{12}\) G and the 2018 NuSTAR work found \(B \approx (1.94\)–\(2.12)\times10^{12}\) G depending on model and assumed redshift [1304.6252] [1906.02917] [2211.06038] [2504.02791].

Iron fluorescence is also established. Suzaku phase-averaged fits included narrow Fe K\(\alpha\) and K\(\beta\) lines at \(6.42\) and \(7.10\) keV with equivalent widths \(\approx 40\)–\(52\) eV and \(\approx 8\)–\(11\) eV, respectively. XMM-Newton measured \(E_{\rm Fe}=6.412^{+0.009}_{-0.005}\) keV, while the 2024 NuSTAR observation found \(E_{\rm Fe}=6.31^{+0.05}_{-0.06}\) keV with \(EW\simeq46\) eV [1304.6252] [2205.03023] [2504.02791].

A more unusual result is a broad absorption feature near \(8\) keV reported from NuSTAR. The 2018 NuSTAR analysis found \(E=7.97\pm0.07\) keV, \(\sigma=1.26\pm0.11\) keV, and depth \(0.33\pm0.07\), and argued that its energy and breadth favor an origin related to Fe XXV K-\(\beta\) and/or blended Ni K\(\alpha\)/K\(\beta\) absorption rather than a cyclotron interpretation [2211.06038]. IXPE discussed this feature as having been reported previously but did not constrain it within its \(2\)–\(8\) keV band [2507.17873].

## 5. Pulse-phase dependence and luminosity dependence of the CRSF

4U 1907+09 is among the best-studied examples of pulse-phase–dependent cyclotron spectroscopy. Suzaku phase-resolved analysis used \(25\) overlapping phase bins, of which \(8\) were independent, and found that the fundamental CRSF energy varies by \(\sim19\%\) over the pulse. In the 2007 observation, \(E_1\) ranged from a minimum \(\sim15\) keV near the second pulse peak at phase \(\approx0.7\)–\(0.8\) to maxima \(\sim20\) keV near phase \(\approx0.3\) and again near phase \(\approx0.6\). The line depth showed a double-peaked pattern, ranging from \(\approx0.2\) to \(1.4\), and the width tracked the energy variation [1304.6252].

The same phase-dependence is robust against continuum choice and, within a modest luminosity interval, against luminosity changes. A Suzaku proceedings study compared two observations differing by a factor of \(\approx 2\) in luminosity and found the same phase-resolved CRSF pattern in both epochs. It also showed that NPEX and CompTT reproduce essentially identical phase trends, leading to the conclusion that the \(\approx19\%\) modulation is intrinsic and not an artifact of continuum parameterization or luminosity drift [1309.5790].

AstroSat/LAXPC confirmed the phase dependence with \(10\) independent phase bins. In that work the CRSF centroid varied by \(\approx12\%\) around \(18.5\) keV, reaching its minimum at the main peak and increasing on its trailing side and around the secondary peak. The line strength was nearly constant over most phases but increased sharply during the rise of the secondary peak; the authors emphasized two features, namely different energy dependence of the two pulse peaks and a strong CRSF only around the secondary peak, as indications of a deviation from a dipole geometry of the magnetic field or a complex beaming pattern from the two poles [1906.02917].

NuSTAR phase-resolved spectroscopy in 2024 likewise found that the spectral parameters vary systematically with pulse phase. The cutoff energy varies in phase with the pulse profile, while photon index and e-folding energy vary out of phase. The CRSF centroid varies with phase and is broadly correlated with the pulse shape. The pulse profiles themselves are energy dependent: the 2024 NuSTAR light curve showed an asymmetric double-peaked structure with phase separation \(0.47\), the peak near phase \(\approx0.25\) fading above \(\approx15\) keV, and pulse fraction increasing with energy [2504.02791].

The luminosity dependence of the CRSF is more nuanced. The 2024 NuSTAR observation found that the fundamental CRSF energy remained consistent with being constant across a \(25.8\pm0.9\)-fold flux swing, with \(E_{\rm cyc,1}\simeq18.1\)–\(18.6\) keV in off-state, low-flux on-state, and high-flux on-state spectra within uncertainties. The 2017 AstroSat flare analysis similarly found no significant CRSF change across a factor \(\approx2.6\) in flux. By contrast, an INTEGRAL literature synthesis found evidence for a positive correlation between CRSF energy and luminosity, with slope \(4\pm3\) keV per \(10^{37}\) erg s\(^{-1}\) and Pearson coefficient \(r=0.92\), though it emphasized the tentative nature of that conclusion and its sensitivity to individual data points [2504.02791] [1906.02917] [1309.0875].

Interpretation of the pulse-phase behavior is not unique. One Suzaku study argued that the pattern of deepest and widest lines near pulse peaks and shallowest and narrowest lines near off-pulse favors a fan-beam emission pattern, whereas the 2024 NuSTAR study described the measured luminosity \(L_{3-50\,\mathrm{keV}}\simeq2.85\times10^{35}\) erg s\(^{-1}\) as consistent with a “pencil” beam radiation pattern expected from a collisionless gas-mediated shock [1304.6252] [2504.02791]. This suggests that the observational phenomenology is being interpreted in terms of geometry-dependent beam mixtures rather than a single immutable beam configuration.

## 6. Dips, off-states, and accretion-regime transitions

Short-timescale variability is a central property of 4U 1907+09. Suzaku observations in 2006 and 2007 showed dips and flares, including a “deep dip” spanning approximately \(60\)–\(100\) ks into the 2006 observation. Color-color analysis demonstrated that some intervals are consistent with variable absorption, but the deep dip is not consistent with absorption-only scenarios and instead points to a temporary reduction or cessation of accretion onto the neutron star. The broader interpretation was a clumpy wind from the supergiant donor, producing stochastic changes in both absorbing column and mass-accretion rate [1111.1800].

Subsequent analyses emphasized that pulsations persist in the low state. A dedicated Suzaku study of the dipping activity concluded that the source continues to pulsate in the “off” state, that transitions between “on” and “off” may be either dip-like or flare-like, and that off-states can account for up to \(\sim60\%\) of the observing time. Using \(d=4\) kpc, that work quoted representative luminosities \(L_X\approx8.6\times10^{35}\) erg s\(^{-1}\) in the on-state and \(L_X\approx8.6\times10^{34}\) erg s\(^{-1}\) in the off-state, with an on/off ratio \(\approx10\), and proposed that 4U 1907+09 may be a missing link between supergiant fast X-ray transients and ordinary accreting pulsars [1210.4428].

The longer RXTE/INTEGRAL campaign revised the fraction of exposure time spent in dips from \(28\%\) to \(24\%\), again finding dips concentrated at orbital phases \(0.1\)–\(0.6\) and absent at \(0.7\)–\(0.8\). Because the dip occurrence is anti-correlated with \(N_{\mathrm H}\), that analysis favored decreases in accretion rate over enhanced absorption. It wrote the standard scalings
\[
r_{\mathrm m}\approx \xi \left(\frac{\mu^4}{2GM\dot{M}^2}\right)^{1/7}, \qquad
r_{\mathrm{co}}=\left(\frac{GM}{\Omega^2}\right)^{1/3},
\]
and argued that when \(r_{\mathrm m}\) approaches or exceeds \(r_{\mathrm{co}}\), centrifugal inhibition can temporarily reduce accretion and pulsed flux [1212.3433].

This line of argument was developed further in the gated-accretion interpretation. For 4U 1907+09, one study proposed that normal states correspond to Rayleigh–Taylor-instability-dominated entry through the magnetosphere, whereas off-states correspond to Kelvin–Helmholtz leakage at a reduced luminosity. Using the observed low-state luminosity, it argued that routine switching between these regimes is easier to accommodate if the surface dipole is stronger than the field directly inferred from the \(\sim19\) keV line, specifically \(B\sim10^{13}\) G, with the observed CRSF then forming at a greater height where the local field is weaker [1210.4428].

An alternative but related framework is quasi-spherical settling accretion. In this picture, appropriate for \(L_X\lesssim\) a few \(\times10^{36}\) erg s\(^{-1}\), a hot quasi-static shell forms above the magnetosphere and plasma entry is regulated by the cooling time near the Alfvén radius. Two subsonic regimes are possible: Compton-cooled and radiatively cooled. The characteristic transition luminosity was written as
\[
L_{\dag}\sim 3\times10^{35}\ \mathrm{erg\ s^{-1}}\ \mu_{30}^{-3/10},
\]
with an off-state luminosity scale
\[
L_{x,\mathrm{rad}}\approx10^{35}\ \mathrm{erg\ s^{-1}}\ \mu_{30}^{7/33}.
\]
The model proposed that 4U 1907+09 enters off-states when the accretion-column optical depth drops, the beam switches from fan to pencil, equatorial Compton cooling is suppressed, and the flow transitions to the radiative-cooling regime [1209.4962] [1407.3163].

A different semi-analytical treatment emphasized photoionization feedback in eccentric wind accretion. Applied qualitatively to 4U 1907+09 with \(e=0.28\), it suggested that strong photoionization can inhibit wind acceleration and produce off-states with durations of several hundreds of seconds and flux drops by factors \(\sim10\)–\(20\), while the bright flare at periastron is more likely due to a temporary switch to disk accretion, outside the range of validity of a pure wind-accretion treatment [2012.00392].

Recent NuSTAR data demonstrate that the phenomenology remains active. The 2024 observation captured deep dips and flares with a \(25.8\pm0.9\)-fold flux swing in the \(3\)–\(25\) keV band. Despite that range, the fundamental and harmonic CRSFs persisted across all flux states, and the pulse shape remained similar in low- and high-flux on-states, supporting the interpretation that the short-term variability reflects changes in mass inflow at the magnetospheric boundary rather than large structural changes in the line-forming region [2504.02791].

## 7. X-ray polarimetry and the contemporary picture

IXPE added a new diagnostic by providing the first high-quality polarization measurements of 4U 1907+09. Two observations were obtained near periastron in 2024 November. The first yielded a phase-averaged polarization degree \( {\rm PD}=6.0\pm1.6\% \) and polarization angle \( {\rm PA}=69^\circ\pm8^\circ \); the second gave \( {\rm PD}=2.2\pm1.6\% \) and \( {\rm PA}=46^\circ\pm23^\circ \); the combined dataset gave \( {\rm PD}=3.7\pm1.1\% \) and \( {\rm PA}=63^\circ\pm9^\circ \) in the \(2\)–\(8\) keV band [2507.17873].

The polarimetric results are notable for their energy dependence. In phase-averaged analysis, the null hypothesis of energy-independent PA was rejected at \(p=0.049\), corresponding to \(\approx1.7\sigma\). More strikingly, the phase-resolved analysis found a probable \(\approx90^\circ\) PA rotation between the adjacent \(4\)–\(5\) and \(5\)–\(6\) keV bands within phase \(0.25\)–\(0.375\). The largest phase-resolved polarization degree occurred at phase \(0.15\)–\(0.25\), with \( {\rm PD}=13.5\pm3.3\% \) and \( {\rm PA}=71^\circ\pm7^\circ \) from XSPEC fitting; model-independent pcube analysis gave \(17.1\pm4.7\%\) and \(60^\circ\pm8^\circ\) for the same phase interval [2507.17873].

The same study examined three short flares, which contributed \(\approx30\%\) of the total photons, and found that flare-only and non-flare polarimetric properties are consistent within uncertainties. It therefore concluded that flares do not significantly affect the energy-phase-dependent PA and that the global accretion geometry remains stable during these events [2507.17873].

In the IXPE interpretation, a simple rotating-vector model with energy-independent PA is disfavored. The observations instead motivate either superposition of at least two polarized spectral components with nearly orthogonal PAs and different energy dependences, or magnetized-plasma transfer effects including partial mode conversion. This modern polarimetric view is consistent with the broader observational record: 4U 1907+09 combines a stable cyclotron field scale of order \(2\times10^{12}\) G, strongly phase-dependent line formation, structured wind-fed variability, and evidence that pulse morphology and polarization are set by a complex, phase-dependent emission geometry rather than by a single-axis, single-component beam [2507.17873].

Source: https://www.emergentmind.com/topics/4u-1907-09