4U 1700-377: High-Mass X-ray Binary Insights

• 4U 1700-377 is an eclipsing high-mass X-ray binary defined by a compact object accreting from an O6.5 Iaf+ supergiant, noted for its strong X-ray variability and pulse ambiguity.
• Eclipse timing and orbital studies reveal a short, decaying orbital period with asymmetrical ingress/egress features, indicative of accretion wakes and wind inhomogeneities.
• Spectral and polarimetric analyses identify candidate cyclotron features and high polarization degrees, supporting a neutron-star interpretation despite mass ambiguities.

Searching arXiv for recent and relevant papers on 4U 1700-377 to ground the article. arxiv_search(query="4U 1700-377 high mass X-ray binary timing spectral polarization cyclotron eclipse", max_results=10) arxiv_search(query="4U 1700-377", max_results=10) 4U 1700-37 is an eclipsing, wind-fed high-mass X-ray binary in which a compact object accretes from the O6.5 Iaf+ supergiant HD 153919. Across X-ray timing, spectroscopy, eclipse diagnostics, astrometry, and now polarimetry, it has emerged as a prototypical but unusually ambiguous system: its orbit is short, its X-ray variability is strong, coherent pulsations have not been firmly detected, and the nature of the compact object has remained debated because its inferred mass lies near the neutron-star/black-hole transition. Recent work nevertheless converges on a picture in which 4U 1700-37 is most likely a neutron-star HMXB embedded in a dense, structured stellar wind, with candidate cyclotron features implying a magnetic field of order $10^{12}\,\mathrm{G}$, strong eclipse reprocessing, and a runaway origin from NGC 6231 (West-Ocampo et al., 4 Feb 2026, Meij et al., 2021, Ninoyu et al., 18 Sep 2025).

1. System classification and evolutionary setting

4U 1700-37 is identified with HD 153919 and is classified as a high-mass X-ray binary with an O supergiant donor. Gaia-based kinematic analysis confirmed NGC 6231 as the parent cluster of the system, with a distance of $1.58 \pm 0.07$ kpc for 4U 1700-37 and $1.63 \pm 0.15$ kpc for NGC 6231. The system moves with a space velocity of $63 \pm 5$ km/s with respect to the cluster, and its trace-back indicates a kinematical age of $2.2 \pm 0.1$ Myr. Isochrone fitting gave an adopted age of $4.7 \pm 0.4$ Myr for NGC 6231, implying that the progenitor of the compact object exploded within $3.0 \pm 0.5$ Myr of cluster formation and therefore had an initial mass $>60\,M_\odot$ (Meij et al., 2021).

The same reconstruction gives current component masses of $60 \pm 11\,M_\odot$ for HD 153919 and $2.54 \pm 0.27\,M_\odot$ for the compact object, with orbital separation $1.58 \pm 0.07$0 and orbital period $1.58 \pm 0.07$1 days (Meij et al., 2021). This mass estimate places the compact object in a regime that has motivated both neutron-star and low-mass black-hole interpretations. The Gaia study argues that the current high space velocity and reconstructed evolutionary history make a large natal kick more likely, and therefore favor a neutron star over a black hole, while not absolutely ruling out a low-mass black hole (Meij et al., 2021).

A broader implication of this reconstruction is that 4U 1700-37 provides a rare constraint set for post-supernova binary evolution. The system has been discussed as a possible Galactic prototype for progenitors of gravitational-wave events such as GW190412, because its inferred history combines very high progenitor mass, mass transfer, supernova recoil, and survival as a compact HMXB (Meij et al., 2021).

2. Orbital ephemeris, eclipse geometry, and orbital evolution

Because detectable pulsations are absent, the long-term orbital evolution of 4U 1700-37 has been derived primarily from eclipse timing rather than pulse arrival-time analysis. A systematic eclipse-timing study used archival mid-eclipse measurements together with new determinations from RXTE-ASM, Swift-BAT, MAXI-GSC, and RXTE-PCA light curves, fitting the sequence of eclipse times with the quadratic ephemeris

$1.58 \pm 0.07$2

That analysis yielded $1.58 \pm 0.07$3 MJD, $1.58 \pm 0.07$4 days, and an orbital period decay rate $1.58 \pm 0.07$5, explicitly smaller than earlier estimates (Islam et al., 2016).

An earlier long-baseline ephemeris based on RXTE/ASM and INTEGRAL/ISGRI monitoring gave $1.58 \pm 0.07$6, $1.58 \pm 0.07$7 days, and $1.58 \pm 0.07$8 (Falanga et al., 2015). The later reanalysis attributed the larger earlier decay estimates partly to underestimated systematic uncertainties, especially in single-orbit observations affected by flares and dips near ingress and egress (Islam et al., 2016). The coexistence of these determinations is central to the orbital-evolution literature on the source.

The geometry of the eclipse has also been used to constrain the binary. The 2015 ephemeris study reported a semi-eclipse angle $1.58 \pm 0.07$9, ingress phase $1.63 \pm 0.15$0, ingress duration $1.63 \pm 0.15$1, egress phase $1.63 \pm 0.15$2, and egress duration $1.63 \pm 0.15$3, with sharp but asymmetric eclipse transitions. The asymmetry was interpreted primarily as the effect of accretion wakes trailing the compact object (Falanga et al., 2015).

Eccentricity and apsidal motion remain a distinct controversy. Using ten years of Swift-BAT eclipse timings and eclipse durations, one study found no evidence for the periodic modulation expected from significant apsidal motion. For apsidal motion rates greater than $1.63 \pm 0.15$4 degrees per year, it derived $1.63 \pm 0.15$5 from mid-eclipse times and $1.63 \pm 0.15$6 from eclipse-duration variations, in tension with earlier optical radial-velocity estimates of $1.63 \pm 0.15$7 (Islam et al., 2016). The timing data therefore favor a nearly circular orbit, although the discrepancy with optical/UV radial-velocity work has not been fully eliminated.

3. Timing phenomenology and broadband variability

The timing phenomenology of 4U 1700-37 is dominated by strong aperiodic variability and an unusually persistent non-detection of coherent pulsations. Insight-HXMT observations obtained in 2020 during out-of-eclipse found significant flux variations on kilo-second timescales while the hardness ratio between 10–30 keV and 2–10 keV remained relatively stable, and no evident pulsations were found over the frequency range $1.63 \pm 0.15$8–$1.63 \pm 0.15$9 Hz (Xiao et al., 2023).

A NuSTAR and ASTROSAT/LAXPC timing study similarly found no evidence for coherent pulsations or quasi-periodic oscillations in the frequency range $63 \pm 5$0 mHz to $63 \pm 5$1 Hz and did not confirm earlier claims of 97 min, 67 s, or other modulations (Bala et al., 2020). A later homogeneous NuSTAR analysis extended the search to periods from $63 \pm 5$2 s up to $63 \pm 5$3 s, again finding no coherent pulsations and constraining the pulsed fraction to $63 \pm 5$4 at 90% confidence in the 3–79 keV band (West-Ocampo et al., 4 Feb 2026).

Suzaku observations add one of the main exceptions to the otherwise null periodicity record. In that dataset, no coherent pulsations were found in the $63 \pm 5$5 Hz to 2 Hz interval, but a quasi-periodic oscillation at $63 \pm 5$6 mHz was detected at more than the $63 \pm 5$7 level (Jaisawal et al., 2015). This is consistent with an accretion flow that is highly structured but not stably modulated by a directly observed spin signal.

On much longer timescales, low-frequency power-density-spectrum analysis with RXTE/ASM identified a break frequency for 4U 1700-37 of $63 \pm 5$8. With $63 \pm 5$9 days, the normalized break frequency is $2.2 \pm 0.1$0, placing the source among wind-accreting HMXBs with short viscous times and only short-radius discs (İçdem et al., 2011). This timing diagnostic has been used to distinguish 4U 1700-37 from Roche-lobe-overflow systems, and suggests that any disc present is small and probably transient.

Eclipse timing behavior introduces an additional scale. Archival XMM-Newton and AstroSat eclipse observations revealed flares during eclipse with minimum count-rate doubling times of 100 s and 200 s, respectively (Rikame et al., 2024). The rapidity of these events indicates that the reprocessing region is larger than, but comparable to, the companion star rather than being an extremely extended halo (Rikame et al., 2024).

4. Spectral states, continuum modeling, and the wind-fed environment

Broadband spectroscopy has consistently shown that the X-ray continuum of 4U 1700-37 can be described by models commonly used for accreting pulsars, but with strong model dependence in the detailed interpretation of residuals. Insight-HXMT spectroscopy over 2–100 keV found that the spectrum is well described by phenomenological pulsar continua such as a power law with a high-energy cutoff. That analysis explicitly accounted, for the first time, for the differing detector orientations of the Insight-HXMT instruments, thereby enabling reliable results even in the presence of stable field-of-view contamination (Xiao et al., 2023).

Suzaku spectroscopy over 1–70 keV found that partially absorbed high-energy cutoff power-law and partially absorbed NPEX continua describe the source well. The same analysis detected Fe K$2.2 \pm 0.1$1 at $2.2 \pm 0.1$2 keV and Fe K$2.2 \pm 0.1$3 at $2.2 \pm 0.1$4 keV, and used time-resolved spectroscopy of 20 intervals to trace rapid changes in absorption and line behavior (Jaisawal et al., 2015). The column density increased sharply after orbital phase $2.2 \pm 0.1$5, coincident with an eclipse-like low-flux segment, and this was interpreted as the passage of an accretion wake across the line of sight (Jaisawal et al., 2015).

A different line of spectral analysis used BeppoSAX, Suzaku, and RXTE data to model the source with a low-temperature blackbody, two Comptonized components, and an iron line. In that framework, $2.2 \pm 0.1$6 for all spectral states, whereas $2.2 \pm 0.1$7 is quasi-constant near 2 when $2.2 \pm 0.1$8 keV and falls into the range $2.2 \pm 0.1$9 when $4.7 \pm 0.4$0 keV (Seifina et al., 2016). The authors interpreted this index behavior, together with the need for two Comptonized seed-photon components, as characteristic of neutron-star systems rather than black-hole systems (Seifina et al., 2016).

The observational picture that emerges is that of a source whose short-term flaring, dips, absorption changes, and iron-line variability are governed by an inhomogeneous stellar wind rather than by a steady, large disc. This is reinforced by the Suzaku result that rapid source-flux variations across the observation are consistent with accretion of inhomogeneously distributed matter in the wind of the supergiant companion (Jaisawal et al., 2015). The spectral behavior therefore links directly to wind structure, accretion wakes, and variable local obscuration.

5. Cyclotron features, magnetic-field estimates, and the compact-object question

Candidate cyclotron resonant scattering features have been reported at several energies, and reconciling these reports has become central to determining the compact object’s magnetic field and identity. A NuSTAR plus ASTROSAT/LAXPC study reported a possible cyclotron feature at $4.7 \pm 0.4$1 keV, consistently present in different continuum models with at least $4.7 \pm 0.4$2 confidence, while failing to detect the previously reported $4.7 \pm 0.4$3 keV feature (Bala et al., 2020). Interpreted as a fundamental electron cyclotron line, the $4.7 \pm 0.4$4 keV feature implies

$4.7 \pm 0.4$5

and thus a magnetic field of approximately $4.7 \pm 0.4$6 G in the emission region (Bala et al., 2020).

The same study also identified a rare Ni K$4.7 \pm 0.4$7 emission line at $4.7 \pm 0.4$8 keV with $4.7 \pm 0.4$9 significance, further emphasizing the dense and chemically informative wind environment (Bala et al., 2020). By contrast, Suzaku spectroscopy had earlier found an absorption-like residual at $3.0 \pm 0.5$0 keV, most notably with the partial covering NPEX continuum, and estimated a surface magnetic field of $3.0 \pm 0.5$1 G if that feature were a cyclotron line (Jaisawal et al., 2015). Insight-HXMT later reported hints of cyclotron absorption around $3.0 \pm 0.5$2 keV or/and $3.0 \pm 0.5$3 keV (Xiao et al., 2023).

A recent uniform NuSTAR reanalysis addressed these discrepancies directly. It found that absorbed blackbody plus cutoff power-law models fit the spectra, but left residuals around 20 keV and 40–50 keV; these features improve the fits but do not constitute firm cyclotron detections under the preferred continuum. The inferred magnetic field range is $3.0 \pm 0.5$4–$3.0 \pm 0.5$5 G, and intensity-resolved spectroscopy suggests possible shifts of the apparent line centroid (West-Ocampo et al., 4 Feb 2026). This places the line-identification problem in a model-dependent regime rather than establishing a single uncontested cyclotron energy.

The same NuSTAR study favors quasi-spherical subsonic accretion and an equilibrium spin period of $3.0 \pm 0.5$6 ks, a regime compatible with the lack of detected pulsations (West-Ocampo et al., 4 Feb 2026). This does not by itself prove the compact object is a neutron star, but it narrows the viable physical parameter space. Taken together, the CRSF candidates, the $3.0 \pm 0.5$7-scale field estimates, and the spectral-state diagnostics all weigh toward a neutron-star interpretation (Bala et al., 2020, Seifina et al., 2016, West-Ocampo et al., 4 Feb 2026).

That conclusion remains formally short of unanimity. The compact-object mass estimate of $3.0 \pm 0.5$8 straddles the canonical neutron-star/black-hole divide, and the lack of persistent pulsations has historically sustained the ambiguity (Meij et al., 2021). A balanced reading of the recent literature is therefore that a neutron star is more strongly supported than a black hole, but the case rests on convergent indirect evidence rather than on a single decisive observable.

6. Eclipse reprocessing, emission-line diagnostics, and X-ray polarization

Eclipse observations have become one of the most informative probes of the circumsource environment in 4U 1700-37 because the direct continuum is suppressed and reprocessed components dominate. A systematic search of archival data found flares during eclipse in 4U 1700-37, and spectral comparison of flare and non-flare intervals showed changes in the power-law photon index together with multiple emission lines. In 4U 1700-37, eleven lines were modeled, including Ne K$3.0 \pm 0.5$9, Ne IX, Ne X, Mg XI, Mg XII, Si XIII, S XV, Ca XIX, Fe K$>60\,M_\odot$0, Fe XXV, and Fe XXVI; the line fluxes rose with the overall eclipse-flare continuum, indicating formation in the binary environment rather than in the interstellar medium (Rikame et al., 2024).

The same eclipse-flare analysis found a soft excess below $>60\,M_\odot$1–3 keV that remained unchanged between flare and non-flare states. Because this component did not track the variable reprocessed continuum, it was interpreted as emission from an extremely thin shell of stellar wind immediately above the photosphere of the companion star, with an observed emission measure of $>60\,M_\odot$2 and shell thickness $>60\,M_\odot$3 in the bremsstrahlung interpretation (Rikame et al., 2024). This separates at least two reprocessing zones: a rapidly varying scattering region and a quasi-static near-photospheric soft emitter.

The polarimetric breakthrough came with IXPE, which provided the first statistically significant detection of X-ray polarization from 4U 1700-37. The average polarization degree was $>60\,M_\odot$4 in 2–8 keV, with polarization angle $>60\,M_\odot$5; above 5 keV, the polarization degree exceeds 10%, placing the source among the highest-polarization HMXBs observed to date (Ninoyu et al., 18 Sep 2025). The phase dependence is particularly striking: polarization peaks above 20% during short ingress and egress intervals just before and after eclipse, remains only $>60\,M_\odot$6 during eclipse, and does not peak at mid-eclipse (Ninoyu et al., 18 Sep 2025).

These results suggest that the dominant polarized component is produced by scattering in material spatially localized between the compact object and the O-type companion, likely associated with large-scale inhomogeneities in the stellar wind and its interaction with the compact object’s emission (Ninoyu et al., 18 Sep 2025). The IXPE study explored contributions from disk winds and orbital reflection, but concluded that neither mechanism alone explains the strong, phase-localized polarization increase around eclipse boundaries (Ninoyu et al., 18 Sep 2025). A plausible implication is that wind clumping, focused flows, or related asymmetric structures are not merely secondary complications but major determinants of the observable X-ray geometry.

In combination, eclipse spectroscopy and polarimetry recast 4U 1700-37 as a system in which line emission, obscuration, scattering, and polarization are all controlled by a structured wind on binary scales. This interpretation is consistent with the accretion-wake signatures seen in folded light curves and time-resolved spectroscopy, and it makes the source a particularly valuable laboratory for studying reprocessing in supergiant HMXBs (Falanga et al., 2015, Jaisawal et al., 2015, Ninoyu et al., 18 Sep 2025).