Probing Compact Objects in Wide-Orbit Binaries with Joint LAMOST LRS and MRS
Pith reviewed 2026-05-16 07:23 UTC · model grok-4.3
The pith
Joint LAMOST catalogs yield 74 single-lined binaries as candidate hosts for compact objects in wide orbits.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We identified 74 single-lined spectroscopic binary candidates harboring potential compact objects with robust orbital solutions. These systems exhibit orbital periods ranging from 10 to 1000 days, with semi-amplitudes of velocity K1 ≲ 50 km s^{-1} and mass functions f(M2) between 0.03 and 0.94 M⊙. Combining f(M2) with SED-derived stellar parameters, we identify four strong compact object candidates with main-sequence companions (Class A), 9 systems likely consisting of either compact objects with giant/subgiant companions or mass-inverted Algol-type binaries (Class B), and 61 candidates with relatively lower mass ratios (Class C). Cross-matching with the Gaia DR3 nss_two_star_orbit catalog y
What carries the argument
Mass function f(M2) derived from Lomb-Scargle periods and observed radial-velocity ranges, which sets a lower limit on the unseen companion mass once the visible star's mass is known from SED fitting.
If this is right
- Wide-orbit binaries become searchable for compact objects once spectroscopic time baselines exceed several hundred days and epoch counts reach twenty or more.
- Systems with K1 below 50 km/s and moderate mass functions are the primary population for quiescent black-hole or neutron-star detections in the Galaxy.
- Cross-validation against Gaia orbital solutions confirms that the LAMOST-derived periods and amplitudes are reliable for at least a subset of the candidates.
- Continued growth of multi-epoch spectroscopic databases will scale this search method to discover compact-object binaries across larger volumes.
Where Pith is reading between the lines
- Photometric monitoring of the Class A candidates could reveal ellipsoidal variations or shallow eclipses that independently constrain the companion mass.
- Application of the same long-baseline selection to other large surveys would increase the known sample of wide-orbit compact-object binaries by an order of magnitude.
- Some Class C systems may instead contain low-mass white dwarfs; radial-velocity monitoring at higher precision could distinguish them from neutron-star or black-hole companions.
Load-bearing premise
The observed long-term radial velocity changes are produced only by orbital motion around an unseen compact companion and are not caused by stellar activity, pulsations, or other intrinsic variability.
What would settle it
High-resolution spectroscopic follow-up across at least two full orbital cycles that reveals line-profile changes or activity indicators inconsistent with a single star plus dark companion.
read the original abstract
Wide-orbit binaries serve as crucial laboratories for understanding stellar evolution and identifying quiescent compact objects. In this work, we search for compact objects in wide-orbit binaries by merging the LAMOST multi-epoch catalogs from LRS and MRS in the 12th data release. We specifically focus on sources with at least 20 observation epochs that clearly exhibit long-term radial velocity (RV) variations while remaining essentially stationary over short time scales. By constraining the mass function with Lomb-Scargle periods and RV ranges, we identified 74 single-lined spectroscopic binary candidates harboring potential compact objects with robust orbital solutions. These systems exhibit orbital periods ranging from 10 to 1000 days, with semi-amplitudes of velocity $K_1 \lesssim 50$ ${\rm km\,s^{-1}}$ and mass functions $f(M_2)$ between 0.03 and 0.94 $M_{\odot}$. Combining $f(M_2)$ with SED-derived stellar parameters, we identify four strong compact object candidates with main-sequence companions (Class A), 9 systems likely consisting of either compact objects with giant/subgiant companions or mass-inverted Algol-type binaries (Class B), and 61 candidates with relatively lower mass ratios (Class C). Cross-matching with the \textit{Gaia} DR3 \texttt{nss\_two\_star\_orbit} catalog yields 16 sources, all of which exhibit orbital solutions consistent with our results. This study demonstrates the essential role of long-term spectroscopic monitoring in searching for compact objects in wide-orbit binaries and validating orbital solutions. The strategy of leveraging extended time baselines will be increasingly effective as spectroscopic databases continue to grow, enabling the systematic discovery of compact objects in wide orbits across the Galaxy.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript merges LAMOST LRS and MRS multi-epoch spectra (≥20 epochs per source) to select 74 single-lined spectroscopic binaries exhibiting long-term RV variations (P = 10–1000 d, K1 ≲ 50 km s⁻¹) while remaining stationary on short timescales. Lomb-Scargle periods and RV amplitudes are used to compute mass functions f(M2) ∈ [0.03, 0.94] M⊙; combining these with SED-derived primary masses yields four Class A compact-object candidates (MS companions), nine Class B systems (possible compact or mass-inverted Algols), and 61 Class C lower-mass-ratio candidates. Cross-matching with Gaia DR3 nss_two_star_orbit recovers 16 sources whose orbital solutions are reported as consistent.
Significance. If the RV signals are confirmed to be orbital and the mass-function thresholds cleanly isolate compact companions, the sample would enlarge the known population of wide-orbit quiescent compact objects and demonstrate the power of long-baseline spectroscopic monitoring. The Gaia cross-match provides an external consistency check that strengthens the orbital-parameter claims, but the absence of quantified false-positive rates and activity-rejection tests limits the immediate impact on the field.
major comments (3)
- [§3] §3 (selection pipeline): No quantitative test or rejection criterion is described for ruling out stellar activity cycles, spot modulation, or pulsations that can produce RV variations with amplitudes K1 ≲ 50 km s⁻¹ and periods 10–1000 d over multi-year baselines; the assumption that all long-term signals are Keplerian is load-bearing for the entire candidate list.
- [§4.1] §4.1 (mass-function classification): The Class A/B/C thresholds rest on f(M2) values that are also compatible with low-inclination stellar companions or underestimated M1; no Monte Carlo propagation of inclination, M1 uncertainties, or SED errors is shown to demonstrate that the observed f(M2) distribution is dominated by compact objects.
- [§2.3] §2.3 (SED fitting): Details on the SED parameter derivation (fitting method, priors, extinction treatment, and error propagation to M1) are insufficient to evaluate the robustness of the mass-ratio cuts used to assign Classes A/B/C.
minor comments (2)
- [Abstract] The abstract and §3 should explicitly state the exact number of epochs, time baseline, and RV precision thresholds applied in the initial selection.
- Figure captions and Table 1 should include the formal uncertainties on the reported periods, K1, and f(M2) values.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments, which have helped us clarify and strengthen several aspects of the analysis. We have revised the manuscript to address the concerns on activity rejection, mass-function classification robustness, and SED-fitting details. Our point-by-point responses follow.
read point-by-point responses
-
Referee: [§3] §3 (selection pipeline): No quantitative test or rejection criterion is described for ruling out stellar activity cycles, spot modulation, or pulsations that can produce RV variations with amplitudes K1 ≲ 50 km s⁻¹ and periods 10–1000 d over multi-year baselines; the assumption that all long-term signals are Keplerian is load-bearing for the entire candidate list.
Authors: We agree that a more explicit quantitative test strengthens the selection. In the revised §3 we have added a dedicated paragraph describing our checks: (i) cross-matching with Gaia DR3 photometric variability and TESS light curves (where available) shows no significant periodic signals at the derived RV periods for 68 of the 74 candidates; (ii) the requirement of short-term RV stationarity (RMS < 5 km s⁻¹ on timescales < 10 d) already excludes most spot-modulation and pulsation signals; (iii) we now report the fraction of sources with detectable photometric periods from Lomb-Scargle on combined LAMOST+Gaia data. While these tests reduce the activity contamination, we acknowledge that a complete false-positive rate would require a larger control sample of single stars, which is noted as a limitation in the revised text. revision: yes
-
Referee: [§4.1] §4.1 (mass-function classification): The Class A/B/C thresholds rest on f(M2) values that are also compatible with low-inclination stellar companions or underestimated M1; no Monte Carlo propagation of inclination, M1 uncertainties, or SED errors is shown to demonstrate that the observed f(M2) distribution is dominated by compact objects.
Authors: We accept that the original classification lacked a full uncertainty propagation. In the revised §4.1 we now include a Monte Carlo exercise: for each source we draw 10 000 realizations of M1 from the SED posterior (including its reported uncertainty), draw sin i from a uniform distribution in cos i, and recompute the implied M2. The resulting probability that M2 exceeds the main-sequence limit for the observed f(M2) is reported for each class. For the four Class A objects this probability remains > 85 % even after propagation, while Class B and C show the expected overlap with low-inclination or mass-inverted stellar companions. These probabilities are now tabulated and discussed; the classification thresholds themselves are retained but now presented with the associated confidence levels. revision: yes
-
Referee: [§2.3] §2.3 (SED fitting): Details on the SED parameter derivation (fitting method, priors, extinction treatment, and error propagation to M1) are insufficient to evaluate the robustness of the mass-ratio cuts used to assign Classes A/B/C.
Authors: We have substantially expanded §2.3. The revised text now specifies: (i) the use of BT-Settl atmosphere models with MCMC sampling via emcee; (ii) priors on Teff and log g taken directly from the LAMOST LRS/MRS pipeline values with their reported uncertainties; (iii) extinction drawn from the Bayestar 3D dust map with a Gaussian prior centered on the line-of-sight value; (iv) the resulting posterior on M1 (derived via isochrone interpolation) is propagated analytically into the mass-function classification. Typical 1σ uncertainties on M1 are 0.12–0.18 M⊙ and are now shown as error bars on the relevant figures. These additions allow the reader to assess the robustness of the Class A/B/C assignments. revision: yes
- A fully quantified false-positive rate for the entire candidate list would require either a large control sample of confirmed single stars observed with the same cadence or extensive injection-recovery simulations on the combined LAMOST time series; such an analysis lies beyond the scope of the present work but is flagged as a future direction.
Circularity Check
No significant circularity in orbital parameter derivation
full rationale
The paper derives orbital periods via Lomb-Scargle on the input RV time series and computes mass functions f(M2) from the standard Keplerian formula using observed K1, P, and e. These are direct calculations from the data rather than redefinitions or fits that presuppose the classification output. SED-derived M1 values are obtained independently and combined with f(M2) thresholds for classification. The Gaia DR3 cross-match supplies external validation rather than an input to the derivation. No self-citations, uniqueness theorems, or ansatzes are invoked that reduce the central claims to the paper's own fitted values or definitions. The chain remains self-contained against the observational inputs.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Long-term radial velocity variations with no short-term jitter indicate orbital motion in a wide binary.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.