🔭 astrophysics

Optical Super-orbital Modulation of SMC X-1: Disk Precession and a Revised Pulsar Mass

By analyzing synchronized optical and X-ray variations in SMC X-1 using a modified ellipsoidal modulation model that accounts for disk precession and irradiation-induced emission shifts, researchers corrected a 20% underestimation in the donor's radial velocity to revise the neutron star's mass to approximately 1.35 $M_\odot$ .

Original authors: Masafumi Niwano, Nobuyuki Kawai, Michael Fausnaugh

Published 2026-05-21

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Masafumi Niwano, Nobuyuki Kawai, Michael Fausnaugh

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Cosmic Detective Story: Weighing a Dead Star

Imagine you are trying to weigh a ghost. In the universe, that "ghost" is a neutron star—the incredibly dense, dead core of a star that exploded. Scientists have a hard time weighing them because they are so small and far away.

The star in this story is called SMC X-1. It's a cosmic dance partner: a neutron star (the ghost) and a massive, living star (the donor) orbiting each other very closely.

For years, scientists thought the neutron star in this pair was very light—so light that it might be too small to exist according to the laws of physics. This paper says, "Wait a minute, we think we've been weighing it wrong."

The Problem: A Shifty Light Curve

When we look at this binary star system from Earth, the living star doesn't just look like a perfect ball. Because the neutron star is so heavy, it pulls on the living star like a giant magnet, stretching it into a teardrop shape.

As this teardrop spins, we see it from different angles:

When we see the "pointy" end, it looks smaller.
When we see the "flat" side, it looks bigger.

This creates a predictable "wobble" in the brightness, like a lighthouse beam. Scientists use this wobble to calculate how heavy the neutron star is.

But here's the glitch: The light curve wasn't behaving like a perfect, static wobble. Sometimes the "dip" in brightness was deep, and sometimes it was shallow. Sometimes the two bright peaks were equal, and sometimes one was taller than the other. It was as if the teardrop star was changing its shape or lighting up differently every few weeks.

The Culprit: A Spinning, Tilted Umbrella

The authors realized the culprit was an accretion disk—a swirling pizza-dough-shaped ring of gas and dust orbiting the neutron star.

Think of this disk like a giant, tilted umbrella that is slowly wobbling (precessing) as it spins.

The Umbrella Effect: As this umbrella tilts, it blocks some of the neutron star's X-ray light from hitting the living star.
The Flashlight Effect: When the umbrella tilts one way, it shines a bright "flashlight" (X-rays) onto the side of the living star facing us, making it glow brighter. When it tilts the other way, it blocks that light, making the star look dimmer.
The Shadow: The umbrella also casts a shadow on itself and the star, changing the shape of the light curve we see.

The authors built a new computer model (a "modified ellipsoidal modulation model") that accounts for this wobbling umbrella. They fed in data from two space telescopes: TESS (which watches the visible light) and MAXI (which watches the X-rays).

The Big Discovery: The "Center of Gravity" Trick

This is the most important part of the paper.

When the neutron star blasts the living star with intense X-rays, it heats up the side of the star facing the neutron star. This makes that side glow much brighter than the dark side.

Imagine you are trying to find the center of a spinning top by watching a sticker on it.

Normal Star: The sticker is right in the middle. You track the sticker, and you find the true center.
Irradiated Star (SMC X-1): The "sticker" (the bright spot) is glued to the side of the star facing the neutron star. It's not in the middle anymore!

Because the bright spot is shifted to the side, when scientists measured how fast the star was moving back and forth (its radial velocity), they were tracking the bright spot, not the true center of the star.

The Result: They were measuring a "fake" speed that was about 20% slower than the star's actual speed.

The New Weight

In physics, if you think a star is moving slower than it actually is, you calculate that the star pulling on it (the neutron star) must be lighter.

Old Calculation: Because they thought the speed was slow, they calculated the neutron star weighed about 1.06 times the mass of our Sun. This was dangerously close to the theoretical minimum, making it a "lightweight" anomaly.
New Calculation: Once the authors corrected for the "bright spot" trick and realized the star was actually moving 20% faster, the math changed. The neutron star is actually heavier.

The New Weight: The neutron star weighs about 1.35 times the mass of our Sun.

Why This Matters

This new weight is a "normal" weight for a neutron star. It fits perfectly with what we expect from the physics of how stars die and explode.

In summary:
The authors used high-tech space cameras to watch a wobbling, teardrop-shaped star. They realized a spinning disk of gas was acting like a shifting spotlight, tricking scientists into thinking the star was moving slower than it really was. Once they fixed this "spotlight illusion," they found the neutron star is actually a healthy, normal weight, solving a mystery that had scientists worried the laws of physics might be broken.

Technical Summary: Optical Super-orbital Modulation of SMC X-1: Disk Precession and a Revised Pulsar Mass

Problem Statement
The observational determination of the lower mass limit for neutron stars is critical for constraining the physics of core-collapse supernovae (CCSNe). The high-mass X-ray binary SMC X-1 is a pivotal object in this context because previous mass estimates placed its pulsar mass near or potentially below the theoretical lower limit of $\approx 1.1 M_\odot$ . These prior estimates relied on determining the binary mass ratio and orbital inclination, primarily through spectroscopic radial velocity (RV) measurements of the donor star and modeling of the optical ellipsoidal modulation. However, SMC X-1 exhibits complex optical super-orbital modulations (variations synchronized with the $\sim 40\text{--}60$ day X-ray super-orbital cycle) that were not fully accounted for in previous static models. Ground-based observations lacked the photometric precision and temporal sampling to resolve these variations within individual orbital cycles, leading to potential systematic errors in the derived binary parameters.

Methodology
The authors analyzed high-cadence optical light curves of SMC X-1 obtained by the Transiting Exoplanet Survey Satellite (TESS) across six sectors, combined with long-term X-ray monitoring data from the Monitor of All-sky X-ray Image (MAXI).

Observational Analysis: The study identified systematic variations in the optical orbital light curves synchronized with the X-ray super-orbital modulation. Specifically, the depth of the minimum at inferior conjunction and the asymmetry of the double-peaked light curve varied depending on the X-ray state (rising, decaying, on, or off).
Model Development: A modified ellipsoidal modulation code was developed in Python (utilizing NumPy, SciPy, Astropy, and CuPy for parallelization). Unlike static models, this code incorporates a precessing accretion disk that alters:
- The geometry of X-ray irradiation on the donor star's surface.
- The geometry of optical/UV (OUV) reprocessing on the disk surface.
- The mutual occultation between the donor, disk, and pulsar.
  The model calculates the temperature and blackbody radiation of surface elements (SEs) on both the donor and the disk, accounting for gravity darkening, X-ray heating, and OUV irradiation.
Parameter Estimation: Markov Chain Monte Carlo (MCMC) sampling was performed to fit the model to the TESS optical and MAXI X-ray light curves. The analysis tested two scenarios for the mass ratio ( $q = M_{NS}/M_s$ ): one based on observed RVs ( $\log q = -1.17$ ) and a second where $q$ was increased by 20% ( $\log q = -1.09$ ) to test the hypothesis that the observed donor RV is underestimated.

Key Results

Correlations: The study confirmed strong correlations between X-ray super-orbital flux and optical features: the inferior conjunction minimum is shallow during X-ray "off" states and deep during "on" states; the double-peak asymmetry correlates with the X-ray flux derivative (e.g., the left peak is higher during the X-ray rising phase).
Model Validation: The modified ellipsoidal model successfully reproduced the observed optical and X-ray light curves, providing observational evidence that the super-orbital modulation is driven by a precessing accretion disk. The analysis suggests the disk precession is retrograde.
Radial Velocity Bias: The intense X-ray irradiation of the donor star shifts the center of optical emission away from the gravitational center. The model predicts that this effect causes the observed radial velocity amplitude to be underestimated by approximately 20% (the apparent RV curve appears "clipped" compared to a true sine wave).
Revised Mass Estimate:
- Using the standard mass ratio ( $\log q = -1.17$ ), the derived pulsar mass is $M_{NS} = 1.11 \pm 0.07 M_\odot$ , consistent with previous studies.
- Correcting for the 20% RV underestimation (increasing $q$ to $\log q = -1.09$ ), the revised pulsar mass is $M_{NS} = 1.35 \pm 0.10 M_\odot$ .
Discrepancies: The model required an X-ray luminosity ( $L_X \approx 4 \times 10^{39}$ erg s $^{-1}$ ) significantly higher than typical observed values to reproduce the optical modulation intensity. The authors suggest this may imply anisotropic X-ray emission (beaming toward the donor) or spin-orbit misalignment, though the latter is difficult to reconcile with the system's circular orbit unless mass accretion during a common-envelope phase was inefficient.

Significance
The paper claims that the systematic variations in the optical light curve, driven by the precessing disk, were previously unmodeled in a way that affected mass determination. By developing a model that accounts for the changing irradiation geometry, the authors demonstrate that the donor's radial velocity has likely been underestimated due to irradiation-induced shifts in the emission center. Correcting for this effect yields a revised pulsar mass of $\approx 1.35 M_\odot$ . This value places the SMC X-1 pulsar well within the typical range for neutron stars, removing the tension with theoretical lower limits for core-collapse supernovae and suggesting that previous mass estimates were biased by uncorrected irradiation effects.