MAXI J1807+132: Neutron-Star X-Ray Binary

• MAXI J1807+132 is a neutron-star low-mass X-ray binary exhibiting recurrent outbursts, thermonuclear bursts, and truncated-disk accretion.
• Multiwavelength observations reveal detailed orbital dynamics, rapid reflare activity, and correlated spectral state transitions.
• Binary modeling and optical studies constrained its short 4.26 hr orbital period, high inclination, and significant donor contribution.

Searching arXiv for papers on MAXI J1807+132 to ground the article in the literature. MAXI J1807+132 is a neutron-star X-ray transient and neutron-star low-mass X-ray binary (NS LMXB) identified through thermonuclear Type-I X-ray bursts during its 2019 outburst (Albayati et al., 2020). It was first recognized by MAXI/GSC on 2017 March 13 and has since undergone outbursts in 2017, 2019, and 2023 (Shidatsu et al., 2017, Saavedra et al., 27 Aug 2025). The source is an atoll source (Rout et al., 4 Jun 2025) located at Galactic latitude $b = 15.501^\circ$ (Saavedra et al., 27 Aug 2025), and its quiescent optical behavior has provided the orbital solution: an orbital period of $P_{\rm orb} = 4.258 \pm 0.008$ hr, binary inclination $i = 72\pm5^\circ$, and mass ratio $q = 0.24^{+0.19}_{-0.14}$ (Saavedra et al., 27 Aug 2025). Current work therefore places MAXI J1807+132 among short-period neutron-star transients whose phenomenology spans thermonuclear bursting, truncated-disk outbursts, jet-related optical/infrared behavior, and rapid reflare activity (Albayati et al., 2020, Rout et al., 2024, Rout et al., 4 Jun 2025).

1. Discovery, classification, and observational setting

MAXI J1807+132 was first recognized with the MAXI/GSC nova-search system on 2017 March 13 (Shidatsu et al., 2017). The accurate Swift position was determined as

$$(\alpha^{2000}, \delta^{2000}) = (18^{\mathrm h}08^{\mathrm m}07^{\mathrm s}.549, +13^\circ15'05.''40)$$

with $(l,b)=(40.123127^\circ, 15.501653^\circ)$ and 90% uncertainty of $0.''16$ (Shidatsu et al., 2017). It lies at high Galactic latitude, about $15^\circ$ above the plane (Jiménez-Ibarra et al., 2018).

The source was initially discussed as a candidate neutron-star low-mass X-ray binary because its early Swift/XRT spectra in the 2017 decay could be described by a blackbody with a relatively low temperature $0.1\text{--}0.5$ keV plus a hard power-law component with photon index $\sim 2$, and because the optical/X-ray flux correlation was consistent with known neutron star LMXBs (Shidatsu et al., 2017). However, the 2017 data did not allow an unambiguous compact-object identification, and subsequent optical and X-ray work emphasized that both a black hole and a neutron star remained possible at that stage (Jiménez-Ibarra et al., 2018).

The classification issue was resolved by the 2019 outburst. NICER detected three thermonuclear Type-I X-ray bursts during a five-day interval in late October 2019, establishing that the accretor has a solid surface and is therefore a neutron star (Albayati et al., 2020). Later work explicitly describes the system as a neutron-star low-mass X-ray binary transient and an atoll source (Rout et al., 2024, Rout et al., 4 Jun 2025).

A notable aspect of the source’s observational history is recurrence. The system underwent outbursts in 2017, 2019, and 2023 (Saavedra et al., 27 Aug 2025). This repeated transient activity, together with the short orbital period measured in quiescence, places MAXI J1807+132 within the population of short-period NS transients whose accretion geometry and emission mechanisms can be studied across widely separated luminosity regimes (Saavedra et al., 27 Aug 2025).

2. Orbital period, ellipsoidal modulation, and binary geometry

The orbital solution for MAXI J1807+132 was obtained from quiescent optical photometry acquired with the 2.5 m Isaac Newton Telescope during three nights in 2022: 24–25 June and 28 July (Saavedra et al., 27 Aug 2025). The campaign used time-resolved $P_{\rm orb} = 4.258 \pm 0.008$0-band photometry, yielding 86 exposures of 600 s each with WFC/Chip 4 and the Harris $P_{\rm orb} = 4.258 \pm 0.008$1 filter (Saavedra et al., 27 Aug 2025). Because the source was faint in quiescence, the photometry was extracted with optimal photometry via the HiPERCAM pipeline and calibrated against local reference stars tied to the GSC 2.2 catalog (Saavedra et al., 27 Aug 2025).

The mean quiescent brightness was measured as $P_{\rm orb} = 4.258 \pm 0.008$2 mag, consistent with the faintest pre-discovery Pan-STARRS measurements transformed to $P_{\rm orb} = 4.258 \pm 0.008$3 mag (Saavedra et al., 27 Aug 2025). This consistency supports the conclusion that the source was in a deep quiescent state during the INT observations (Saavedra et al., 27 Aug 2025). The light curve showed clear variability on hour timescales, and a Lomb–Scargle periodogram applied to the first two nights revealed a strongest signal at

$P_{\rm orb} = 4.258 \pm 0.008$4

corresponding to

$P_{\rm orb} = 4.258 \pm 0.008$5

The peak lay above the $P_{\rm orb} = 4.258 \pm 0.008$6 level for white noise, and after simulations of $P_{\rm orb} = 4.258 \pm 0.008$7 synthetic light curves using a $P_{\rm orb} = 4.258 \pm 0.008$8 prescription with the Timmer & König method, none reproduced the observed power at $P_{\rm orb} = 4.258 \pm 0.008$9, giving $i = 72\pm5^\circ$0 (Saavedra et al., 27 Aug 2025).

The physically significant interpretation comes from folding the data at twice the photometric period. Folding at

$i = 72\pm5^\circ$1

produces a double-humped light curve with two unequal minima, the expected morphology of ellipsoidal modulation from a Roche-lobe-filling donor (Saavedra et al., 27 Aug 2025). In this picture, the companion is tidally distorted, so the projected area changes twice per orbit; the unequal minima arise because gravity darkening is strongest near the inner Lagrange point, with the deeper minimum expected around superior conjunction of the donor (Saavedra et al., 27 Aug 2025). The deeper minimum was defined as phase 0.5 at

$i = 72\pm5^\circ$2

The phase stability of this minimum across different nights secures the orbital-period identification (Saavedra et al., 27 Aug 2025).

Binary modeling used \texttt{XRBinary} with \texttt{emcee} MCMC fitting (Saavedra et al., 27 Aug 2025). The model assumed a Roche-lobe-filling donor co-rotating with the orbit plus a cool accretion disk contributing additional optical light. The donor spectrum was taken from Kurucz atmospheres, with limb darkening from Claret’s nonlinear law and gravity darkening following Claret’s prescription (Saavedra et al., 27 Aug 2025). The disk was assumed cylindrically symmetric, with height profile

$i = 72\pm5^\circ$3

outer radius

$i = 72\pm5^\circ$4

and steady-state viscous temperature law

$i = 72\pm5^\circ$5

To reduce degeneracy, the modeling fixed $i = 72\pm5^\circ$6 K, $i = 72\pm5^\circ$7, $i = 72\pm5^\circ$8, $i = 72\pm5^\circ$9, and donor albedo 0.5 (Saavedra et al., 27 Aug 2025).

The best-fit parameters were

$q = 0.24^{+0.19}_{-0.14}$0

$q = 0.24^{+0.19}_{-0.14}$1

$q = 0.24^{+0.19}_{-0.14}$2

$q = 0.24^{+0.19}_{-0.14}$3

and

$q = 0.24^{+0.19}_{-0.14}$4

with $q = 0.24^{+0.19}_{-0.14}$5 for 67 degrees of freedom, or $q = 0.24^{+0.19}_{-0.14}$6 (Saavedra et al., 27 Aug 2025). The posterior distributions imply that the companion contributes $q = 0.24^{+0.19}_{-0.14}$7 of the total $q = 0.24^{+0.19}_{-0.14}$8-band flux, summarized as roughly 30–50%, consistent with an early M-dwarf donor (Saavedra et al., 27 Aug 2025). The same analysis yielded a donor radial-velocity semi-amplitude $q = 0.24^{+0.19}_{-0.14}$9 km s$$(\alpha^{2000}, \delta^{2000}) = (18^{\mathrm h}08^{\mathrm m}07^{\mathrm s}.549, +13^\circ15'05.''40)$$0 and neutron-star semi-amplitude $$(\alpha^{2000}, \delta^{2000}) = (18^{\mathrm h}08^{\mathrm m}07^{\mathrm s}.549, +13^\circ15'05.''40)$$1 km s$$(\alpha^{2000}, \delta^{2000}) = (18^{\mathrm h}08^{\mathrm m}07^{\mathrm s}.549, +13^\circ15'05.''40)$$2 from

$$(\alpha^{2000}, \delta^{2000}) = (18^{\mathrm h}08^{\mathrm m}07^{\mathrm s}.549, +13^\circ15'05.''40)$$3

An H$$(\alpha^{2000}, \delta^{2000}) = (18^{\mathrm h}08^{\mathrm m}07^{\mathrm s}.549, +13^\circ15'05.''40)$$4 centroid shift reported in previous spectroscopy was noted to be broadly consistent with this $$(\alpha^{2000}, \delta^{2000}) = (18^{\mathrm h}08^{\mathrm m}07^{\mathrm s}.549, +13^\circ15'05.''40)$$5 scale (Saavedra et al., 27 Aug 2025).

This orbital solution supersedes earlier period inferences from the 2017 decay. In 2018, the quasi-periodic reflares of the discovery outburst had been used to estimate $$(\alpha^{2000}, \delta^{2000}) = (18^{\mathrm h}08^{\mathrm m}07^{\mathrm s}.549, +13^\circ15'05.''40)$$6 h for a neutron-star interpretation or $$(\alpha^{2000}, \delta^{2000}) = (18^{\mathrm h}08^{\mathrm m}07^{\mathrm s}.549, +13^\circ15'05.''40)$$7 h for a black-hole interpretation (Jiménez-Ibarra et al., 2018). The directly measured $$(\alpha^{2000}, \delta^{2000}) = (18^{\mathrm h}08^{\mathrm m}07^{\mathrm s}.549, +13^\circ15'05.''40)$$8 hr is close to the earlier black-hole-case estimate, but the later Type-I burst detections show that the compact object is in fact a neutron star (Albayati et al., 2020, Saavedra et al., 27 Aug 2025). This resolves one of the main early ambiguities in the source’s interpretation.

3. Optical and spectroscopic behavior from outburst to quiescence

The 2017 discovery outburst was already in decline when detailed optical monitoring began (Jiménez-Ibarra et al., 2018). The first optical point was about 3 magnitudes brighter than the Pan-STARRS quiescent level, and over the first $$(\alpha^{2000}, \delta^{2000}) = (18^{\mathrm h}08^{\mathrm m}07^{\mathrm s}.549, +13^\circ15'05.''40)$$9 days the source decayed from $(l,b)=(40.123127^\circ, 15.501653^\circ)$0 to roughly $(l,b)=(40.123127^\circ, 15.501653^\circ)$1 (Jiménez-Ibarra et al., 2018). After that, the light curve became highly non-monotonic, with up to 7 re-brightening events superposed on the overall decay and amplitudes of about $(l,b)=(40.123127^\circ, 15.501653^\circ)$2 mag above quiescence (Jiménez-Ibarra et al., 2018). The brightest four episodes appeared quasi-periodic with a recurrence time of $(l,b)=(40.123127^\circ, 15.501653^\circ)$3 days (Jiménez-Ibarra et al., 2018).

A WHT light curve obtained on 2017 July 22 in SDSS-$(l,b)=(40.123127^\circ, 15.501653^\circ)$4 had mean brightness $(l,b)=(40.123127^\circ, 15.501653^\circ)$5 mag, flickering up to $(l,b)=(40.123127^\circ, 15.501653^\circ)$6 mag, and no periodic modulation detected over $(l,b)=(40.123127^\circ, 15.501653^\circ)$7 hours (Jiménez-Ibarra et al., 2018). The authors argued that the source may not yet have reached true optical quiescence because even the faint states still showed flickering and emission-line activity (Jiménez-Ibarra et al., 2018). This is consistent with the later need for deeper quiescent observations in 2022 to reveal the stable orbital clock (Saavedra et al., 27 Aug 2025).

Optical spectroscopy during the 2017 event showed the canonical outburst phenomenology of a transient LMXB. Early spectra displayed Balmer emission lines up to H$(l,b)=(40.123127^\circ, 15.501653^\circ)$8, He II $(l,b)=(40.123127^\circ, 15.501653^\circ)$9, and weaker He I lines near $0.''16$0 Å and $0.''16$1 Å (Jiménez-Ibarra et al., 2018). A distinguishing feature was that these emission lines were embedded in deep, broad absorption troughs, especially in the bluer Balmer lines (Jiménez-Ibarra et al., 2018). Later spectra changed markedly: two were essentially featureless, while the faintest showed H$0.''16$2 in emission (Jiménez-Ibarra et al., 2018). No donor-star absorption features were detected (Jiménez-Ibarra et al., 2018), which is compatible with the later conclusion that even in quiescence the donor contributes only part of the optical flux and that disk light remains substantial (Saavedra et al., 27 Aug 2025).

The line centroids measured in 2017 were $0.''16$3, $0.''16$4, $0.''16$5, and $0.''16$6 km s$0.''16$7 in spectra 1, 2, 3, and 6, respectively, yielding an average systemic velocity

$0.''16$8

This large negative systemic velocity was argued to be difficult to explain by normal Galactic rotation unless the source were at implausibly large distance, leading to discussion of a natal kick (Jiménez-Ibarra et al., 2018). The high Galactic latitude and later distance estimate of $0.''16$9 kpc with $15^\circ$0 kpc reinforce the view that the system occupies a region well above the thin disk (Saavedra et al., 27 Aug 2025). A plausible implication is that the kinematic arguments raised from the early spectroscopy remain relevant in interpreting the system’s Galactic history.

During the 2023 outburst, the optical behavior was tracked in far greater detail with LCO/Faulkes monitoring, UVOT, infrared, radio, and polarimetry (Rout et al., 4 Jun 2025). The source exhibited a fast rise, a plateau or gradual decay phase, a brief flare near the end, and then a sequence of rapid, high-amplitude reflares (Rout et al., 4 Jun 2025). Color–magnitude behavior showed the source becoming bluer during the initial rise, a pronounced redward excursion near the outburst peak, and then a return toward bluer colors during decay (Rout et al., 4 Jun 2025). A uniform-temperature blackbody of constant area reproduced the diagonal track expected from an irradiated disk, while the redward outliers were interpreted as an additional component, likely the jet (Rout et al., 4 Jun 2025).

4. Thermonuclear bursting and confirmation of the neutron star

The decisive evidence that MAXI J1807+132 contains a neutron star came from the 2019 outburst, when NICER detected three thermonuclear Type-I X-ray bursts during observations spanning 27–31 October 2019 (Albayati et al., 2020). Type-I bursts require unstable nuclear burning on a solid surface, excluding a black hole accretor (Albayati et al., 2020).

The three bursts, labeled B1, B2, and B3, showed rapid rises, cooling-like decays, blackbody-like burst spectra, and energetics consistent with thermonuclear events (Albayati et al., 2020). Their common timing behavior included rise times of about 4 s and long decay tails $15^\circ$1 min, and the hardness ratio rose during the burst rise and dropped during decay (Albayati et al., 2020). The observed morphologies—slow rises and long decays—were interpreted as indicative of mixed H/He fuel rather than pure helium (Albayati et al., 2020).

Time-resolved spectroscopy used the variable persistent-flux method with

$15^\circ$2

while the persistent emission was modeled as

$15^\circ$3

with $15^\circ$4 (Albayati et al., 2020). Before the bursts, the persistent fluxes were $15^\circ$5, $15^\circ$6, and $15^\circ$7 erg cm$15^\circ$8 s$15^\circ$9 for B1, B2, and B3, respectively (Albayati et al., 2020).

Peak burst blackbody temperatures were comparatively low:

• B1: $0.1\text{--}0.5$0 keV
• B2: $0.1\text{--}0.5$1 keV
• B3: $0.1\text{--}0.5$2 keV (Albayati et al., 2020)

The reported peak bolometric fluxes, fluences, and timescales were:

Burst	Peak $0.1\text{--}0.5$3	Fluence	$0.1\text{--}0.5$4
B1	$0.1\text{--}0.5$5 erg s$0.1\text{--}0.5$6 cm$0.1\text{--}0.5$7	$0.1\text{--}0.5$8 erg cm$0.1\text{--}0.5$9	$\sim 2$0 s
B2	$\sim 2$1 erg s$\sim 2$2 cm$\sim 2$3	$\sim 2$4 erg cm$\sim 2$5	$\sim 2$6 s
B3	$\sim 2$7 erg s$\sim 2$8 cm$\sim 2$9	$P_{\rm orb} = 4.258 \pm 0.008$00 erg cm$P_{\rm orb} = 4.258 \pm 0.008$01	$P_{\rm orb} = 4.258 \pm 0.008$02 s

The bursts showed no strong evidence of reaching the Eddington luminosity, and the brightest burst did not show convincing photospheric radius expansion (Albayati et al., 2020). Under the assumption that B1 did reach $P_{\rm orb} = 4.258 \pm 0.008$03, with hydrogen-rich

$P_{\rm orb} = 4.258 \pm 0.008$04

the authors derived an upper limit

$P_{\rm orb} = 4.258 \pm 0.008$05

for the source distance (Albayati et al., 2020). Later photometric distance work produced a more specific estimate of $P_{\rm orb} = 4.258 \pm 0.008$06 kpc (Saavedra et al., 27 Aug 2025), well within this upper bound.

A distinctive result from the burst study was the $P_{\rm orb} = 4.258 \pm 0.008$07 s pause during the rise of B1 (Albayati et al., 2020). The burst rose from about 62 to 500 counts s$P_{\rm orb} = 4.258 \pm 0.008$08 in about 2 s, stalled for about 1.6 s, and then rapidly climbed to a peak of about 2250 counts s$P_{\rm orb} = 4.258 \pm 0.008$09 (Albayati et al., 2020). The paper compared this feature with a similar pause reported in SAX J1808.4–3658 and argued that the pause likely reflects a physical process that can occur in both pulsating and non-pulsating systems, and does not obviously depend on fuel composition, peak luminosity, or whether the system is an X-ray pulsar (Albayati et al., 2020). No burst oscillations were detected; the upper limit was approximately 10% fractional amplitude at 95% confidence (Albayati et al., 2020).

5. Accretion-flow evolution during the 2023 outburst

The 2023 outburst is the best-characterized event in the source’s history, with complementary X-ray and multiwavelength analyses (Rout et al., 2024, Rout et al., 4 Jun 2025). NICER monitoring showed that the outburst began on MJD 60132.1 with about 3 counts s$P_{\rm orb} = 4.258 \pm 0.008$10 in 1–10 keV after two earlier non-detections (Rout et al., 2024). The source then brightened extremely rapidly, from $P_{\rm orb} = 4.258 \pm 0.008$11 counts s$P_{\rm orb} = 4.258 \pm 0.008$12 to $P_{\rm orb} = 4.258 \pm 0.008$13 counts s$P_{\rm orb} = 4.258 \pm 0.008$14 in only a few days, corresponding to a rise by a factor of about 20 in one day near the start of the outburst evolution (Rout et al., 2024). A thermonuclear burst was detected at the first local maximum around MJD 60136.2 (Rout et al., 2024).

The subsequent evolution comprised a plateau phase from roughly MJD 60136–60146, a sharp drop to about 39 counts s$P_{\rm orb} = 4.258 \pm 0.008$15 in less than 3 days with increasing hardness, a short flare near MJD 60150.2 reaching about 50 counts s$P_{\rm orb} = 4.258 \pm 0.008$16, a final decay by MJD 60155.1, and a post-outburst reflare starting near MJD 60158.4 that rose again to $P_{\rm orb} = 4.258 \pm 0.008$17 counts s$P_{\rm orb} = 4.258 \pm 0.008$18 (Rout et al., 2024). The source underwent a full hard $P_{\rm orb} = 4.258 \pm 0.008$19 intermediate $P_{\rm orb} = 4.258 \pm 0.008$20 soft transition and then returned to the hard state during decay, tracing a clear q-shaped hysteresis loop in the hardness–intensity diagram (Rout et al., 2024). In the color–color diagram it occupied the island branch at the beginning and end and the banana branch for most of the rest of the outburst, consistent with standard atoll-source behavior at sub-Eddington accretion rates (Rout et al., 2024).

The hardness was defined as

$P_{\rm orb} = 4.258 \pm 0.008$21

with intensity taken as the 0.5–10 keV count rate (Rout et al., 2024). Early hard-state hardness was around $P_{\rm orb} = 4.258 \pm 0.008$22, dropping to $P_{\rm orb} = 4.258 \pm 0.008$23 in the rapid rise, then decreasing gradually from $P_{\rm orb} = 4.258 \pm 0.008$24 to $P_{\rm orb} = 4.258 \pm 0.008$25 in the plateau soft state, and later increasing again during decay; the reflare had hardness about $P_{\rm orb} = 4.258 \pm 0.008$26 (Rout et al., 2024).

Fast-timing analysis found power spectra dominated by band-limited noise with a cutoff around 1–10 Hz, fitted with a single Lorentzian whose centroid was consistent with zero (Rout et al., 2024). The fractional rms, computed over 0.1–10 Hz, evolved from $P_{\rm orb} = 4.258 \pm 0.008$27 early in the rise to $P_{\rm orb} = 4.258 \pm 0.008$28 by MJD 60146.2 and $P_{\rm orb} = 4.258 \pm 0.008$29 by MJD 60153.1, before increasing again during hard excursions (Rout et al., 2024). The rms and hardness were strongly positively correlated, with Spearman coefficient $P_{\rm orb} = 4.258 \pm 0.008$30 and $P_{\rm orb} = 4.258 \pm 0.008$31-value $P_{\rm orb} = 4.258 \pm 0.008$32; excluding the reflare, the coefficient was $P_{\rm orb} = 4.258 \pm 0.008$33 with $P_{\rm orb} = 4.258 \pm 0.008$34-value $P_{\rm orb} = 4.258 \pm 0.008$35 (Rout et al., 2024). The reflare is the clear exception, showing low hardness but still $P_{\rm orb} = 4.258 \pm 0.008$36 rms (Rout et al., 2024).

Spectral modeling of the NICER data employed the three-component model

$P_{\rm orb} = 4.258 \pm 0.008$37

where TBfeo accounts for interstellar absorption with variable O and Fe abundances, diskbb is a multicolor disk blackbody, bbodyrad is a single-temperature boundary-layer blackbody, and nthComp describes Comptonized emission from the corona (Rout et al., 2024). Best-fit absorption values were

$P_{\rm orb} = 4.258 \pm 0.008$38

O abundance $P_{\rm orb} = 4.258 \pm 0.008$39 solar, and Fe abundance $P_{\rm orb} = 4.258 \pm 0.008$40 solar at $P_{\rm orb} = 4.258 \pm 0.008$41 (Rout et al., 2024).

The inferred geometry was strongly state dependent. In hard states at the beginning and end of the outburst, and during the hard excursion near MJD 60147, the boundary-layer component was not required or became unconstrained, and the diskbb component dominated the thermal emission; this was interpreted as evidence that the disk was truncated far from the neutron star (Rout et al., 2024). In brighter intermediate and soft states, all three components were required and the disk radius stabilized at a much smaller value consistent with reaching the last stable orbit (Rout et al., 2024). Using correction factors $P_{\rm orb} = 4.258 \pm 0.008$42 and $P_{\rm orb} = 4.258 \pm 0.008$43, the inferred stable radius was approximately $P_{\rm orb} = 4.258 \pm 0.008$44 km for 5 kpc / $P_{\rm orb} = 4.258 \pm 0.008$45 and $P_{\rm orb} = 4.258 \pm 0.008$46 km for 1 kpc / $P_{\rm orb} = 4.258 \pm 0.008$47, corresponding roughly to $P_{\rm orb} = 4.258 \pm 0.008$48 and $P_{\rm orb} = 4.258 \pm 0.008$49, respectively (Rout et al., 2024). The disk moved inward by a factor of about 40 in $P_{\rm orb} = 4.258 \pm 0.008$50 days (Rout et al., 2024).

The characteristic timing frequency $P_{\rm orb} = 4.258 \pm 0.008$51 was strongly anti-correlated with the disk inner radius, with Spearman coefficient $P_{\rm orb} = 4.258 \pm 0.008$52 and $P_{\rm orb} = 4.258 \pm 0.008$53-value $P_{\rm orb} = 4.258 \pm 0.008$54 (Rout et al., 2024). This was interpreted as evidence that the corona shrinks as the source softens: in the hard state the corona fills the space between the neutron star and the truncated disk, whereas rising accretion rate drives the disk inward and contracts the corona (Rout et al., 2024). The authors further suggested that the corona may evolve from a more horizontal geometry in the hard state to a more vertical or compact structure in the intermediate state (Rout et al., 2024). This suggests that MAXI J1807+132 provides a relatively clean case of correlated disk truncation and fast-variability evolution in an atoll transient.

6. Multiwavelength emission, jets, polarization, and reflares

The 2023 outburst was covered by a broadband campaign including NICER, Swift/XRT + UVOT, LCO/Faulkes, VLT/FORS2 polarimetry, REM $P_{\rm orb} = 4.258 \pm 0.008$55-band, MeerKAT, and VLA follow-up (Rout et al., 4 Jun 2025). One of the central results is a measured delay between optical and X-ray rise times. A first optical brightening or precursor occurred on MJD 60106.2, about 26 days before the first significant NICER X-ray detection on MJD 60132.1; the main optical rise began between MJD 60120.2 and 60127.2, so the optical rise preceded the X-ray rise by about 4–12 days (Rout et al., 4 Jun 2025). This timing is interpreted as evidence for the disk instability model with a truncated inner disk: the heating wave reaches the truncated inner edge, after which the inner disk moves inward on a viscous timescale and X-rays rise later (Rout et al., 4 Jun 2025). This interpretation is explicitly linked to the separate NICER result that the disk was highly truncated at the beginning of the outburst (Rout et al., 2024, Rout et al., 4 Jun 2025).

The source also showed a slow long-term brightening during quiescence beginning about 400 days before the main outburst, with linear rise rates of $P_{\rm orb} = 4.258 \pm 0.008$56 mag/yr in $P_{\rm orb} = 4.258 \pm 0.008$57, $P_{\rm orb} = 4.258 \pm 0.008$58 mag/yr in $P_{\rm orb} = 4.258 \pm 0.008$59, $P_{\rm orb} = 4.258 \pm 0.008$60 mag/yr in $P_{\rm orb} = 4.258 \pm 0.008$61, and $P_{\rm orb} = 4.258 \pm 0.008$62 mag/yr in $P_{\rm orb} = 4.258 \pm 0.008$63 (Rout et al., 4 Jun 2025). This was interpreted as standard disk-instability-model disk mass buildup (Rout et al., 4 Jun 2025).

The optical/UV spectral energy distributions were fitted with

$P_{\rm orb} = 4.258 \pm 0.008$64

with $P_{\rm orb} = 4.258 \pm 0.008$65 varying from roughly $P_{\rm orb} = 4.258 \pm 0.008$66 to $P_{\rm orb} = 4.258 \pm 0.008$67 in quiescence, $P_{\rm orb} = 4.258 \pm 0.008$68 during the optical flare on MJD 60106.2, $P_{\rm orb} = 4.258 \pm 0.008$69 near the outburst peak on MJD 60134.4, and $P_{\rm orb} = 4.258 \pm 0.008$70 by MJD 60157.3 during decay (Rout et al., 4 Jun 2025). During reflare peaks, $P_{\rm orb} = 4.258 \pm 0.008$71 was typically $P_{\rm orb} = 4.258 \pm 0.008$72 (Rout et al., 4 Jun 2025). Steep or bluer SEDs during rise and flare peaks are consistent with disk emission, whereas flat or mildly inverted SEDs during the plateau and decay are consistent with either a viscous irradiated disk or a self-absorbed compact jet (Rout et al., 4 Jun 2025).

Optical/X-ray correlations were fit as

$P_{\rm orb} = 4.258 \pm 0.008$73

with $P_{\rm orb} = 4.258 \pm 0.008$74 in $P_{\rm orb} = 4.258 \pm 0.008$75 and $P_{\rm orb} = 4.258 \pm 0.008$76 in $P_{\rm orb} = 4.258 \pm 0.008$77, with intermediate values in other bands and a trend to steeper $P_{\rm orb} = 4.258 \pm 0.008$78 at longer wavelengths (Rout et al., 4 Jun 2025). The paper notes that $P_{\rm orb} = 4.258 \pm 0.008$79 is expected for a viscous disk, $P_{\rm orb} = 4.258 \pm 0.008$80 for an irradiated disk, and larger $P_{\rm orb} = 4.258 \pm 0.008$81 may indicate jet contribution; the measured correlations are positive but scattered, and the wavelength dependence is suggestive but not highly significant statistically (Rout et al., 4 Jun 2025).

The evidence for jet synchrotron emission is strongest during the intermediate and hard states of the plateau phase. MeerKAT detected the source significantly once, on MJD 60141, near the outburst peak, with flux about 0.17 mJy (Rout et al., 4 Jun 2025). REM $P_{\rm orb} = 4.258 \pm 0.008$82-band observations showed variable emission during the decay or plateau (Rout et al., 4 Jun 2025). Together with the reddening in optical colors and flat-to-mildly inverted optical SEDs, these observations support a jet contribution during the decay phase (Rout et al., 4 Jun 2025). The inferred optically thick synchrotron component likely contributes at OIR wavelengths; if the $P_{\rm orb} = 4.258 \pm 0.008$83-band is near the upper limit of the jet break, the optically thick spectral index is $P_{\rm orb} = 4.258 \pm 0.008$84, while the optically thin component has slope about $P_{\rm orb} = 4.258 \pm 0.008$85, suggesting that the jet break may lie near the IR/optical regime (Rout et al., 4 Jun 2025).

A particularly unusual finding was the optical polarimetry. The polarization was always low, $P_{\rm orb} = 4.258 \pm 0.008$86, but changed dramatically in the third VLT/FORS2 epoch at MJD 60149.0 (Rout et al., 4 Jun 2025). The formalism was given as

$P_{\rm orb} = 4.258 \pm 0.008$87

and

$P_{\rm orb} = 4.258 \pm 0.008$88

where $P_{\rm orb} = 4.258 \pm 0.008$89 is the polarization degree and $P_{\rm orb} = 4.258 \pm 0.008$90 the polarization angle (Rout et al., 4 Jun 2025). In the first epoch, $P_{\rm orb} = 4.258 \pm 0.008$91 in all bands with $P_{\rm orb} = 4.258 \pm 0.008$92; in the third epoch, $P_{\rm orb} = 4.258 \pm 0.008$93 became undetectable in $P_{\rm orb} = 4.258 \pm 0.008$94 and $P_{\rm orb} = 4.258 \pm 0.008$95, increased to $P_{\rm orb} = 4.258 \pm 0.008$96 in $P_{\rm orb} = 4.258 \pm 0.008$97 and $P_{\rm orb} = 4.258 \pm 0.008$98 in $P_{\rm orb} = 4.258 \pm 0.008$99, and the polarization angle rotated by about $i = 72\pm5^\circ$00 relative to earlier epochs (Rout et al., 4 Jun 2025). This occurred less than a day before a short X-ray flare on MJD 60150, after which the source entered the soft state (Rout et al., 4 Jun 2025). The interpretation offered is compression and reordering of magnetic fields in the jet, possibly a discrete ejection during the hard-to-soft transition, followed by quenching of the compact jet (Rout et al., 4 Jun 2025). The paper emphasizes that a $i = 72\pm5^\circ$01 polarization-angle swing is usually seen in radio during jet transitions, so a similar optical/IR event is unusual and possibly the first such event in an X-ray binary (Rout et al., 4 Jun 2025).

Rapid reflares are a recurrent theme in the source’s phenomenology. In 2017, the decay exhibited quasi-periodic re-brightenings with recurrence time $i = 72\pm5^\circ$02 days (Jiménez-Ibarra et al., 2018). In 2023, the main outburst was followed by at least six high-amplitude rapid reflares with approximate peaks at MJD 60163.4, 60175.4, 60189.2, 60200.8, 60218.9, and 60239.2, separated by $i = 72\pm5^\circ$03 days (Rout et al., 4 Jun 2025). Their amplitudes were about 2 orders of magnitude above quiescence, with rise times of $i = 72\pm5^\circ$04 days; the X-rays had the largest amplitude, followed by UV and optical (Rout et al., 4 Jun 2025). The reflares are described as high-amplitude and short-duration, which is unusual for NS X-ray binaries (Rout et al., 4 Jun 2025). Their SEDs are generally consistent with irradiated disk emission, colors remain on the hot branch, and the disk temperature stays above $i = 72\pm5^\circ$05 K even at faint optical levels, arguing against traveling heating/cooling fronts as the primary cause (Rout et al., 4 Jun 2025). Proposed explanations include weak propeller or trapped disk behavior, irradiation of the outer disk, enhanced mass transfer from the companion, or echoes of earlier accretion events, but the mechanism remains poorly constrained (Rout et al., 4 Jun 2025). The X-ray analysis likewise found that the post-outburst reflare was not a simple scaled-down repeat of the main outburst, because the spectrum appeared soft while retaining relatively high variability (Rout et al., 2024).

7. Distance, Galactic location, and broader significance

The quiescent optical study used the fitted orbital period to extend the empirical quiescent absolute-magnitude–orbital-period relation, originally developed for black hole X-ray transients, into the neutron-star regime (Saavedra et al., 27 Aug 2025). The adopted relation was

$i = 72\pm5^\circ$06

For $i = 72\pm5^\circ$07 d and faintest quiescent $i = 72\pm5^\circ$08, this gives

$i = 72\pm5^\circ$09

Using

$i = 72\pm5^\circ$10

with extinction $i = 72\pm5^\circ$11 mag from 3D dust maps, the inferred distance is

$i = 72\pm5^\circ$12

At Galactic latitude $i = 72\pm5^\circ$13, this corresponds to

$i = 72\pm5^\circ$14

placing MAXI J1807+132 well above the Galactic plane (Saavedra et al., 27 Aug 2025).

This distance estimate is significant in several respects. First, it is consistent with the upper limit $i = 72\pm5^\circ$15 kpc obtained from the Type-I burst analysis under the assumption of photospheric radius expansion in the brightest burst (Albayati et al., 2020). Second, the quiescent study states that the resulting distance agrees with independent estimates from X-ray spectral modeling and with expectations from donor-star brightness, strengthening both the binary interpretation and the early M-dwarf classification (Saavedra et al., 27 Aug 2025). Third, it bears directly on earlier controversies. The 2018 work noted that the optical/X-ray luminosity plane would place the source in the black-hole region at $i = 72\pm5^\circ$16 kpc, while a neutron star seemed more plausible at $i = 72\pm5^\circ$17 kpc or less (Jiménez-Ibarra et al., 2018). The later confirmation of the neutron star by thermonuclear bursts and the photometric distance near 6.3 kpc show that the optical/X-ray ratio was not, by itself, a reliable classifier for this object (Albayati et al., 2020, Saavedra et al., 27 Aug 2025). This is an instructive example of how source class, distance, and accretion state can complicate empirical luminosity-plane diagnostics.

The 2024 NICER spectral study also discussed a broad distance range of 1–5 kpc when interpreting the inner disk radius, boundary-layer radius, magnetic field, and luminosity (Rout et al., 2024). At 1 kpc, the blackbody or boundary-layer radius is only 2–8 km, while at 5 kpc it is 10–40 km; the larger distance was considered somewhat more natural if the blackbody is interpreted as the boundary layer (Rout et al., 2024). The magnetic field inferred under the assumption that the disk truncates at the magnetospheric or Alfvén radius was $i = 72\pm5^\circ$18 G for $i = 72\pm5^\circ$19 kpc and $i = 72\pm5^\circ$20 G for $i = 72\pm5^\circ$21 kpc (Rout et al., 2024). The high-distance case implied a rather strong field compared with the typical $i = 72\pm5^\circ$22 G found in many atoll sources (Rout et al., 2024). In light of the later $i = 72\pm5^\circ$23 kpc photometric distance, this tension becomes a salient open issue rather than a resolved inconsistency. This suggests that either the simple truncation-to-magnetosphere mapping is incomplete for this source or that some of the geometric assumptions entering the radius estimate require revision.

In the broader context of compact-binary studies, MAXI J1807+132 is important because it now has a relatively complete phenomenological portrait. It is a short-period NS LMXB with a quiescently visible ellipsoidal donor, a high inclination, a moderate mass ratio, repeated outbursts, Type-I bursts, state transitions with hysteresis, evidence for truncated-disk evolution, OIR jet signatures, a rare near-orthogonal optical polarization rotation, and unusually strong rapid reflares (Saavedra et al., 27 Aug 2025, Albayati et al., 2020, Rout et al., 2024, Rout et al., 4 Jun 2025). The source is also significant methodologically: the 2025 orbital paper explicitly extends the quiescent $i = 72\pm5^\circ$24–$i = 72\pm5^\circ$25 correlation beyond black hole transients into the neutron-star regime, proposing a purely photometric distance estimator for NS transients in quiescence (Saavedra et al., 27 Aug 2025). As a result, MAXI J1807+132 has become a reference case for connecting quiescent binary geometry to the multiwavelength physics of transient accretion.