A Bayesian approach to the analysis of time symmetry in light curves: Reconsidering Scorpius X-1 occultations

We present a new approach to the analysis of time symmetry in light curves, such as those in the x-ray at the center of the Scorpius X-1 occultation debate. Our method uses a new parameterization for such events (the bilogistic event profile) and provides a clear, physically relevant characterization of each event’s key features. We also demonstrate a Markov Chain Monte Carlo algorithm to carry out this analysis, including a novel independence chain configuration for the estimation of each event’s location in the light curve. These tools are applied to the Scorpius X-1 light curves presented in Chang et al. (2007), providing additional evidence based on the time series that the events detected thus far are most likely not occultations by TNOs.

💡 Research Summary

The paper introduces a Bayesian framework for testing time‑symmetry in astronomical light curves, with a focus on the controversial X‑ray dips observed in the Scorpius X‑1 data set. The authors argue that previous analyses relied on overly simplistic symmetric models and did not capture the detailed shape of each dip. To address this, they propose a “bilogistic event profile,” a parametric form that models the rise and fall of an event with two separate logistic‑type functions. This yields five interpretable parameters: event centre (τ), pre‑ and post‑event time scales (α₁, α₂), an asymmetry index (β), background level (μ), and depth (δ). The asymmetry index directly quantifies how much the rise time differs from the fall time, a key diagnostic for distinguishing true occultations by distant trans‑Neptunian objects (TNOs) from instrumental or statistical artifacts.

Parameter inference is performed using a Markov Chain Monte Carlo (MCMC) algorithm. A novel independence‑chain proposal is used for τ, drawing candidate positions from a prior distribution rather than a random‑walk step. This dramatically reduces autocorrelation and speeds convergence, even when the posterior is multimodal. The remaining parameters are updated with standard Metropolis–Hastings steps, each equipped with appropriate priors (log‑normal for time scales, uniform for asymmetry, etc.). Posterior samples are summarized with histograms and 95 % credible intervals, allowing the authors to assess both the magnitude and uncertainty of asymmetry for each dip.

The method is applied to the 58 dip events reported by Chang et al. (2007). For the majority of events, the posterior distribution of β is significantly different from zero, indicating pronounced asymmetry. In many cases the rise time is markedly shorter than the fall time, a pattern that contradicts the expectation for a TNO occultation, which should produce a symmetric, box‑like dip. Bayes factor calculations further favor the asymmetric bilogistic model over a symmetric alternative. Consequently, the authors conclude that the observed dips are far more likely to be caused by instrumental glitches, background fluctuations, or data‑processing artifacts rather than genuine occultations.

Beyond the specific Scorpius X‑1 case, the paper highlights the broader utility of the bilogistic profile and the independence‑chain MCMC scheme. The approach can be transferred to other high‑energy astrophysical time series—such as gamma‑ray bursts, pulsar timing irregularities, or fast radio burst profiles—where precise shape characterization is essential for physical interpretation. Future work is suggested to extend the model to simultaneous multi‑event fitting and to incorporate hierarchical priors that could capture population‑level asymmetry trends.

In summary, the study provides a rigorous Bayesian toolset for quantifying time‑symmetry in light curves, demonstrates its effectiveness on a contentious data set, and delivers compelling statistical evidence that the Scorpius X‑1 dips are not TNO occultations but rather manifestations of non‑astrophysical variability. This work therefore reshapes the debate on small‑body detection via X‑ray occultations and sets a new standard for statistical analysis in time‑domain astronomy.