Superorbital Phase Evolution and a Soft–Hard X-Ray Phase Shift in LMC X-4

Abstract The superorbital period of LMC X-4 is among the most stable known in Roche-lobe overflow, high-mass X-ray binaries. We analyzed 33 yr of monitoring data from the Compton Gamma Ray Observatory Burst and Transient Source Experiment, the Rossi X-ray Timing Explorer All-Sky Monitor, the Neil Gehrels Swift Burst Alert Telescope, the Monitor of All-sky X-ray Image Gas Slit Camera, and the Fermi Gamma-ray Burst Monitor. The measured phases show a smooth long-term trend with superposed systematic fluctuations. Fits with cubic, quartic, and sinusoidal models indicate that the quartic and sinusoidal forms provide significantly better descriptions, with the sinusoidal model yielding an <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mn>890</mml:mn> <mml:msubsup> <mml:mrow> <mml:mn>0</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>−</mml:mo> <mml:mn>230</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>+</mml:mo> <mml:mn>210</mml:mn> </mml:mrow> </mml:msubsup> </mml:math> days modulation. Such a long timescale is unlikely to arise from orbital motion around a tertiary companion. The fluctuations resemble stochastic, glitch-like events on several-hundred-day timescales. Their rms period variation exceeds that of the smooth trend, yet the total rms period variation over 33 yr remains only 0.55%, demonstrating the exceptional stability of the superorbital period. During MJD 57000-60461, we detect a phase offset of 0.044 ± 0.010 cycles between the soft and hard X-ray bands. This offset can be reproduced by including a higher-harmonic term in the azimuthal disk model, allowing a transition from antisymmetric to asymmetric structure. A contemporaneous decline in the hard X-ray flux suggests a partial obscuration of the emission region, similar to the anomalous low state in Her X-1. This evolving-disk scenario may also explain the superorbital phase shift previously reported in Her X-1.