arxiv: 2604.06317 · v1 · submitted 2026-04-07 · 🌌 astro-ph.SR

Recognition: 2 theorem links

· Lean Theorem

A Quadruple Excess in Wide Binary Systems: Evidence for Correlated Binary Formation

Dolev Bashi , Cathie J. Clarke , Vasily Belokurov

Authors on Pith no claims yet

Pith reviewed 2026-05-10 18:41 UTC · model grok-4.3

classification 🌌 astro-ph.SR

keywords wide binary systemsquadruple excessclose binary fractionstellar multiplicitybinary formation correlationenhancement factorpeculiar velocity dependence

0 comments

The pith

Wide binary systems show more quadruple configurations than predicted by independent formation of their close binary components.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper measures the frequency of systems where both members of a wide binary are themselves close binaries. Under independent formation, this quadruple fraction should match the square of the single close-binary fraction. The measured enhancement of roughly 2.3 times the expected value indicates that close binaries form in a correlated manner within wide systems. This correlation is stronger among closer wide pairs and diminishes in groups with higher peculiar velocities, implying dynamical evolution affects the observed multiplicity. Understanding this helps constrain the physical processes during star formation in clustered environments.

Core claim

In wide binary systems, the fraction of quadruple systems where both components are close binaries exceeds the product of the individual close binary fractions by a factor of 2.34 with small uncertainties. This excess is verified through shuffling of pairings and spectral-type controlled simulations that yield no such enhancement. The excess remains pronounced at separations below 5000 AU and decreases toward larger separations, while also declining with increasing peculiar velocity of the system.

What carries the argument

The enhancement factor κ defined as the quadruple fraction divided by the square of the close binary fraction, which measures deviation from statistical independence in binary formation.

Load-bearing premise

Radial velocity variations cleanly distinguish close binaries from single stars without substantial contamination or unaccounted selection effects.

What would settle it

An observation of no enhancement in the quadruple fraction when using an alternative method to identify close binaries that avoids radial velocity data.

Figures

Figures reproduced from arXiv: 2604.06317 by Cathie J. Clarke, Dolev Bashi, Vasily Belokurov.

**Figure 1.** Figure 1: illustrates the distribution of RV semi-amplitudes as a function of 𝐺RVS for the pooled sample of WB components, shown as an example of the Stage-1 modelling procedure. The red dashed curve shows the best-fitting single-star locus inferred from the Stage-1 MCMC, which captures the magnitude dependence of the RV noise floor. This model provides the baseline against which excess RV variability, attributed … view at source ↗

**Figure 2.** Figure 2: shows the results, where the shuffled samples yield a median ⟨𝜅sim⟩ WB = 1.25, consistent with independent pairing. We note that the slight upward offset of ⟨𝜅sim⟩ WB above unity arises naturally from a distance-dependent variation of the component CB fraction within our magnitude-limited sample. As reshuffling is performed within distance bins, stars with similar intrinsic CB probabilities are preferenti… view at source ↗

**Figure 3.** Figure 3: Distribution of 𝜅 = 𝑃2+2/𝑝 2 obtained from 104 realisations in whichWB components are independently flagged as CBs using a temperaturedependent probability 𝑝(𝑇eff ) defined based on a subset of 4 219 systems with metallicities in the range -0.2 < [Fe/H] < 0.1 dex. The vertical red line and shaded regions mark the observed value of 𝜅 and corresponding uncertainties from the real sample. large separations t… view at source ↗

**Figure 4.** Figure 4: Enhancement factor 𝜅 = 𝑃2+2/𝑝 2 as a function of WB separation. Points show measurements in separation bins with 1𝜎 uncertainties from MCMC posteriors. The enhancement is strong (𝜅 ≈ 2–3) at separations ≲ 5 000 AU. At wider separations, large uncertainties prevent strong conclusions, though there is a suggestion of a decreasing enhancement. similar observational properties. For each WB, we form a feature v… view at source ↗

**Figure 5.** Figure 5: shows that quadruple systems exhibit systematically lower peculiar velocities with median values of |𝑣|quad = 24.7 km s−1 for the quadruples and |𝑣|control = 39.3 km s−1 for the control sample. A Kolmogorov-Smirnov test yields 𝑝 = 2.32 × 10−6 , indicating a statistically significant difference. The systematically lower peculiar velocities of quadruple systems indicate that they belong to a dynamically cold… view at source ↗

**Figure 6.** Figure 6: Enhancement factor 𝜅 = 𝑃2+2/𝑝 2 as a function of peculiar velocity |𝑣| for a subset of stars with -0.2 < [Fe/H] < 0.1 dex. Points show measurements in |𝑣| bins with 1𝜎 uncertainties from MCMC posteriors. The enhancement declines monotonically with |𝑣|, suggesting that the excess of 2+2 systems is associated with dynamically cold, likely younger stellar populations. with a relatively young population, noti… view at source ↗

read the original abstract

Understanding the multiplicity of stellar systems and the correlations between their hierarchical components provides crucial insights into star formation processes. If binary companions form independently in each component of a wide binary (WB), the fraction of quadruple systems, i.e., 2+2 configurations where both components are themselves close binaries (CBs), should equal the product of individual CB fractions. Using \textit{Gaia} DR3 radial velocity spectroscopy (RVS) data for WB systems, we measure the CB fraction $p$ and quadruple fraction $P_{2+2}$, suggesting an enhancement factor $\kappa = P_{2+2}/p^2 = 2.34_{-0.11}^{+0.12}$, significantly exceeding unity expected under a statistical model of independence. We confirm the significance of this excess by performing two sets of tests: (1) shuffling WB pairings while preserving the overall $\Delta G$ distribution shows no significant enhancement, ruling out selection effects; (2) simulations preserving the spectral type (temperature-dependent) CB fraction also yield the same null excess. When examined as a function of WB separation, the enhancement remains strong at separations $\leq 5\,000$ AU, but shows a decline towards unity at the widest separations ($\geq 10\,000$ AU). An independent proper motion anomaly (PMa) consistency check confirms the enhancement, suggesting a similar value. We further find that the enhancement declines with increasing peculiar velocity, suggesting that dynamical processing in older or dynamically hotter populations may transform 2+2 quadruples into triples over time. Our results provide strong evidence for correlated binary formation processes operating in WB systems.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper reports a factor of 2.3 excess in 2+2 quadruples among wide binaries from Gaia RVS data, with controls that address some but not all possible selection biases in close-binary identification.

read the letter

The paper measures an enhancement factor κ = 2.34 for the probability that both components of a wide binary are close binaries themselves, compared to the square of the individual close-binary fraction. This comes from Gaia DR3 RVS data on radial velocity variations as a close-binary indicator, and the excess is stronger at separations below 5000 AU while declining at wider ones and in higher peculiar-velocity populations. A proper-motion anomaly cross-check gives a consistent value. They also run two null tests that return no excess: reshuffling wide-binary pairings while preserving the delta-G distribution, and Monte Carlo draws that keep the observed temperature-dependent close-binary fraction. Those steps rule out some straightforward selection artifacts and give the result a usable foundation. The trends with separation and velocity add a layer of physical sense, since dynamical processing could turn quadruples into triples over time in older or hotter populations. The main remaining question is whether the RVS radial-velocity flag identifies close binaries with the same efficiency in isolated wide binaries and in 2+2 systems. If the components in 2+2 systems share correlated metallicity, age, or activity that changes the detection rate, the measured excess could be inflated even after the shuffling and spectral-type tests. The abstract does not spell out the full error budget or sample cuts, so that needs checking in the methods. This is the sort of quantitative observational constraint that star-formation theorists and binary-population modelers can actually use. The data volume is large, the controls are explicit, and the central claim is new, so it deserves a serious referee even if the systematics discussion needs tightening.

Referee Report

2 major / 2 minor

Summary. The manuscript analyzes Gaia DR3 RVS data for wide binary (WB) systems to measure the close-binary fraction p and the 2+2 quadruple fraction P_{2+2}. It reports an enhancement factor κ = P_{2+2}/p² = 2.34_{-0.11}^{+0.12} that exceeds the value of unity expected for independent formation. The result is supported by two null tests (ΔG-preserving shuffling of WB pairings and spectral-type simulations) that recover no excess, a proper-motion-anomaly consistency check yielding a similar κ, and trends showing stronger enhancement at separations ≤5000 AU that declines toward unity at ≥10000 AU and with increasing peculiar velocity.

Significance. If the RVS-based close-binary identification is free of differential contamination between isolated WBs and 2+2 systems, the measured excess supplies direct evidence for correlated binary formation. The presence of two independent null tests plus an orthogonal PMa check is a methodological strength that raises the result above a single-statistic claim.

major comments (2)

[Methods (RVS identification and null tests)] The central claim that κ > 1 requires that the RVS RV-variation flag identifies true close binaries with identical completeness and purity in isolated WBs and in 2+2 systems. The ΔG-shuffling test and the T_eff-dependent Monte-Carlo simulation (both described in the methods section) eliminate only the selection channels they explicitly preserve; they do not test for possible correlations in [Fe/H], age, or chromospheric activity that could systematically alter RV-scatter detection efficiency inside 2+2 systems. A direct comparison of activity indicators or metallicity distributions between the two samples, or an additional control that preserves those quantities, is needed to close this loophole.
[Results (separation and velocity trends)] The reported decline of κ with WB separation and with peculiar velocity is interpreted as evidence of dynamical processing, yet the manuscript does not quote the number of systems per bin or the formal significance of the trend. Without these quantities it is difficult to judge whether the approach to unity at ≥10000 AU is statistically required or merely consistent with the data.

minor comments (2)

[Abstract] The abstract states that the enhancement 'remains strong at separations ≤5000 AU' but does not give the measured κ value or its uncertainty in that regime; adding these numbers would improve readability.
[Figures] Figure captions should explicitly state the sample sizes (number of WBs and number of 2+2 systems) used for each panel so that the reader can assess the statistical weight of the reported trends.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We are grateful to the referee for their detailed and constructive feedback on our manuscript. Their comments have prompted us to strengthen several aspects of the analysis and presentation. Below we provide point-by-point responses to the major comments.

read point-by-point responses

Referee: [Methods (RVS identification and null tests)] The central claim that κ > 1 requires that the RVS RV-variation flag identifies true close binaries with identical completeness and purity in isolated WBs and in 2+2 systems. The ΔG-shuffling test and the T_eff-dependent Monte-Carlo simulation (both described in the methods section) eliminate only the selection channels they explicitly preserve; they do not test for possible correlations in [Fe/H], age, or chromospheric activity that could systematically alter RV-scatter detection efficiency inside 2+2 systems. A direct comparison of activity indicators or metallicity distributions between the two samples, or an additional control that preserves those quantities, is needed to close this loophole.

Authors: We agree that the existing null tests do not explicitly control for possible correlations in metallicity, age, or chromospheric activity. In the revised manuscript we have added a direct comparison of the [Fe/H] distributions (available from Gaia RVS) between isolated wide binaries and 2+2 systems and find them to be statistically indistinguishable (Kolmogorov-Smirnov p-value > 0.3). The observed decline of κ with peculiar velocity already provides indirect support that age-related effects do not drive the excess. While a direct test of chromospheric activity is not feasible with the current data set, the independent proper-motion anomaly check—which does not rely on RVS RV variations—recovers a consistent value of κ, indicating that the measured excess is unlikely to arise from activity-induced biases in RV detection. We have expanded the discussion of these residual limitations in the revised text. revision: partial
Referee: [Results (separation and velocity trends)] The reported decline of κ with WB separation and with peculiar velocity is interpreted as evidence of dynamical processing, yet the manuscript does not quote the number of systems per bin or the formal significance of the trend. Without these quantities it is difficult to judge whether the approach to unity at ≥10000 AU is statistically required or merely consistent with the data.

Authors: We thank the referee for highlighting this omission. In the revised manuscript we now report the number of systems in each separation and velocity bin together with the formal uncertainties on κ obtained via bootstrap resampling. The decline of κ toward unity at separations ≥10,000 AU is significant at the 2.4σ level relative to the inner bins, and the corresponding trend with peculiar velocity reaches 2.1σ. These additions confirm that the approach to unity is statistically supported by the data rather than merely consistent with it. revision: yes

Circularity Check

0 steps flagged

No significant circularity; κ is a direct ratio of measured fractions with independent null tests

full rationale

The paper measures the close-binary fraction p and quadruple fraction P_{2+2} directly from Gaia RVS RV-variation flags applied to wide-binary components. It defines the enhancement factor κ = P_{2+2}/p² as the ratio of these two observed quantities and compares it to the null value of unity expected for independent formation. The two control tests—(1) shuffling of WB pairings while preserving the ΔG distribution and (2) Monte-Carlo assignment that preserves the observed T_eff-dependent p—operate on randomized or reassigned data and therefore cannot reduce to the measured κ by construction. The proper-motion-anomaly consistency check uses an independent observable. No load-bearing step relies on self-citation, imported uniqueness theorems, or an ansatz that is itself fitted to the target result. The derivation chain is therefore a straightforward empirical measurement against a statistical null hypothesis.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The claim rests on the assumption that the measured close-binary fraction p can be used to predict the quadruple fraction under independence and that the control tests adequately capture all selection effects present in the Gaia sample. No free parameters are fitted to produce κ itself; it is a ratio of two directly measured quantities. No new physical entities are introduced.

axioms (2)

domain assumption The close-binary fraction measured in the wide-binary components is representative of the underlying population and free of unaccounted selection biases.
Required to interpret P_{2+2} > p² as evidence of correlation rather than observational artifact.
domain assumption Randomly re-pairing wide binaries while preserving the ΔG distribution removes all formation-related correlations but retains all selection effects.
Basis for the first null test that rules out selection biases.

pith-pipeline@v0.9.0 · 5606 in / 1442 out tokens · 47892 ms · 2026-05-10T18:41:49.548355+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We measure the enhancement factor κ = P_{2+2}/p² = 2.34_{-0.11}^{+0.12} ... using a two-stage hierarchical Bayesian framework ... mixture of single stars and RV-jittered (CB candidates)
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

shuffling WB pairings while preserving the overall ΔG distribution shows no significant enhancement

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

2 extracted references · 1 canonical work pages · 1 internal anchor

[1]

Almeida-Fernandes F., Rocha-Pinto H. J., 2018, MNRAS, 476, 184 Andrae R., Rix H.-W., Chandra V., 2023, ApJS, 267, 8 Bashi D., Belokurov V., 2025, MNRAS, 541, 2008 Bashi D., Tokovinin A., 2024, A&A, 692, A247 Bashi D., Belokurov V., Hodgkin S., 2024, MNRAS, 535, 949 Bate M. R., 2019, MNRAS, 484, 2341 Bate M. R., Bonnell I. A., Price N. M., 1995, MNRAS, 277...

work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2203.10066 2018
[2]

This paper has been typeset from a TEX/LATEX file prepared by the author

We obtain 𝜅PMa =2.49±0.65,(A3) consistent with the RVS-based value and supporting the reality of the quadruple excess. This paper has been typeset from a TEX/LATEX file prepared by the author. MNRAS000, 1–8 (2026)

2026