{"id":"4b864518-3f37-46b8-8654-e98f7a77864a","arxiv_id":"2501.01617","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"The Pleiades binary fraction is 0.34 +/- 0.02, with a mass-ratio distribution that rises, falls, and rises again across three power-law segments.","lead":"This paper measures the binary fraction and mass-ratio distribution of unresolved main-sequence binaries in the Pleiades using Gaia and 2MASS photometry. Its Bayesian multiband method reaches mass ratios down to 0.09 divided by the primary mass, yielding a binary fraction of 0.34 +/- 0.02.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed three-segment mass-ratio distribution is not tested against a single power law convolved with the mass-dependent qlim completeness; a forward model is needed before the low-q break can be considered real.","rationale":"The paper's central quantitative claims are fb=0.34 +/- 0.02 and a three-segment mass-ratio distribution. The RV comparison in Figure 8 (17 stars, all with Pb~1) validates only high-q systems, so the low-q break at q=0.32 rests entirely on the stacked PDFs and has no external anchor. The qlim line in Figure 4 is a hard cutoff at secondary mass 0.09 Msun, but a companion at 0.09 Msun is not photometrically detectable; the true completeness is a smooth function of flux ratio, photometric errors (0.01 mag Gaia, 0.02-0.03 mag 2MASS), and error floors. The paper uses qlim as a sharp boundary in Equations (1), (5), and (7), so posterior mass just above qlim is counted as binary whenever the model offers a modest chi-square improvement, but this is also where the likelihood is nearly flat in q. The stacked P(M1,q) therefore inherits the prior and the selection geometry, and the observed low-q deficiency relative to an extrapolated gamma2=-1 may simply reflect the scarcity of massive primaries that can probe low q and the incompleteness even for those primaries. The paper does not rule out the null hypothesis that the intrinsic q distribution is a single power law; this is precisely the test that should be run. This concern is load-bearing because if the three-segment shape is an artifact, the claimed detection of a low-q deficiency and high-q excess, which is the paper's main astrophysical interpretation, is not established. The reader's weakest_assumption centered on the empirical ridgeline, but the selection-function and forward-model issue is more directly relevant to the headline q-distribution claim, so agreement is only partial.","tokens_in":23666,"tokens_out":11661,"duration_ms":125009,"concrete_test":"Run an injection-recovery simulation using the authors' own empirical model and photometric errors. Generate mock samples with the observed M1 distribution and binary fraction, draw q from a single power law with index gamma=-1 (and also gamma=-0.5 and -1.5), add photometric errors as in Section 2.2, and run the full Bayesian stacking pipeline of Sections 3.2 and 4.4, including the qlim cutoff and uniform q prior. Fit the recovered nb(q) with both a single power law and the three-segment model. If the three-segment model is preferred by |Delta BIC|<10 or if recovered breakpoints near 0.32 and 0.79 appear under a single-power-law input, the claimed three-segment structure is a selection artifact and the central claim is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline claim (Section 4.4) that the mass-ratio distribution follows a three-segment power law with breakpoints at q=0.32 and q=0.79 rests on the stacked posterior P(M1,q) shown in Figure 12. However, the detectability of a binary is not a hard cutoff at qlim(M1)=0.09/M1: a secondary near the hydrogen-burning limit contributes only about 1-2% of the flux in the 2MASS bands for a 0.5 Msun primary, comparable to the photometric errors. The effective completeness is therefore a smooth, mass-dependent function of q that rises from zero near qlim to unity only at substantially larger q. The low-q segment (gamma1=0.8) and the break at qb1~0.32 may be produced entirely by this completeness function acting on a single power-law intrinsic distribution, especially because the number of primaries massive enough to probe q<0.3 is small and weighted by the steep mass function. The paper neither forward-models this selection nor compares the three-segment fit against a single power law convolved with an estimated completeness function. Consequently, the claimed deficiency at lower q and excess at higher q could be artifacts of the method rather than properties of the Pleiades binary population.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a multiband photometric fitting method for unresolved main-sequence binaries in the Pleiades, combining Gaia DR3 and 2MASS photometry with an empirical correction to the PARSEC model. The authors derive a Bayesian PDF for primary mass and mass ratio for each of 1154 members, stack the PDFs, and report a binary fraction fb = 0.34 ± 0.02, a rising binary fraction with primary mass, and a three-segment power-law mass-ratio distribution with exponents γ1 = 0.8 ± 0.2, γ2 = −1.0 ± 0.3, γ3 = 3.1 ± 1.0 and breakpoints at q = 0.32 ± 0.04 and q = 0.79 ± 0.05. The method is validated against 17 radial-velocity binaries from Torres et al. (2021) and is compared with several earlier determinations of the Pleiades binary fraction.","tokens_in":23998,"tokens_out":4535,"duration_ms":47377,"significance":"If the central claims hold, the paper would provide the broadest photometric census of unresolved Pleiades binaries to date, reaching mass ratios close to the theoretical hydrogen-burning limit, and would furnish a clear empirical target for models of binary formation and dynamical evolution. The Bayesian stacking approach is a genuine methodological step beyond point-estimate analyses, and the external comparison with the Torres et al. (2021) dynamical mass ratios is a useful validation for moderate- and high-q binaries. However, the headline fb uncertainty is only Poisson, the empirical model is calibrated on the same cluster it is used to analyze, and the claimed three-segment mass-ratio distribution is not tested against a forward model of the smooth completeness or against simpler distributions. These issues are load-bearing for the main conclusions, so the paper needs major revision before the claims can be accepted.","major_comments":[{"comment":"The empirical photometric model is calibrated on the same cluster it is then used to analyze, and the first calibration step assumes all members are single stars. At low masses the NIR corrections are large, up to about 0.5 mag (Table 3, low-Ms rows; Section 3.1.2), so if unresolved binaries or model errors bias the ITGP ridgeline, the derived M1 and q posteriors and hence the low-q part of the stacked PDFs are systematically biased. The external check against Torres et al. (2021) covers only 17 binaries and does not validate the q<0.3 regime in which the claimed low-q rise and the qb1 break occur. I request a simulation test: inject synthetic binaries with a known intrinsic q distribution into the sample and verify that the ITGP procedure recovers the true single-star ridgeline, or provide an independent low-mass NIR calibration from another cluster or from theoretical models.","section":"Section 3.1, Fig. 3, Table 3"},{"comment":"The three-segment power-law claim is not tested against a single power law convolved with the actual, smooth completeness function. Equation (1) and the hard cutoff qlim(M1) = 0.09/M1 in Section 3.2.1 define a detectability boundary, but a secondary near qlim contributes only about 1–2% of the flux in the 2MASS bands for a 0.5 Msun primary, which is comparable to the adopted photometric error floors in Section 2.2; detectability is therefore a smooth, mass-dependent function of q. Because the number of primaries massive enough to probe q<0.3 is small and weighted by the steep mass function, a single intrinsic power law combined with this completeness can plausibly produce the apparent low-q rise and the qb1 break. The paper should forward-model the selection by simulating the full fitting and stacking procedure with a single power-law intrinsic q distribution, and should compare the three-segment fit against one- and two-segment models using an information criterion. Without this test, the existence of the three-segment distribution is not established.","section":"Section 4.4, Fig. 12"},{"comment":"The reported uncertainty fb = 0.34 ± 0.02 is only the Poisson counting term nb^1/2 / n from Equation (7). It does not propagate the uncertainties from the empirical model corrections (Section 3.1), the adopted priors (log-uniform M1 and uniform q; Section 3.2.2), the membership probability cut Pk ≥ 0.5, the exclusion criterion χj^2 > 25 (Section 3.2.3), the Gaia error floors, or the assumed cluster parameters (age, AV, [Fe/H]) and individual distance moduli. Each of these can shift fb by several percent, and the large spread of previous measurements in Table 5 suggests that such systematics are non-negligible. I recommend a bootstrap over the sample and a marginalization over the calibration and nuisance parameters to produce a realistic systematic uncertainty alongside the Poisson term.","section":"Section 4.1–4.2, Eq. (7)"}],"minor_comments":[{"comment":"The title contains an obvious typo: 'Binary F raction' should be 'Binary Fraction'.","section":"Title/Abstract"},{"comment":"The phrase 'complete mass-ratio distribution' overstates the scope: the sample is limited to G<19, excludes MS+WD and higher-order multiples, and is restricted to q > qlim(M1). Suggest wording such as 'photometric mass-ratio distribution of unresolved MS+MS binaries with M1 > 0.11 Msun.'","section":"Abstract, Section 4.4"},{"comment":"The MiMO fit for age and extinction assumes a flat mass-ratio distribution (γq = 0), which is later contradicted by the measured three-segment q distribution; the text should state explicitly why this prior choice does not feed back into the empirical model and the final results.","section":"Section 2.3"},{"comment":"The gray histogram of nb(q) does not specify the bin width or the error prescription. Please state the binning and whether the error bars are Poisson errors on the integrated counts.","section":"Figure 12"},{"comment":"The interpretive statement that 'the mass ratio of real binary systems, whether or not having an illuminated companion, might be monotonically decreasing' is speculative and is not directly tested by the data; it should be clearly labeled as an interpretation rather than a measurement, or supported by a quantitative model.","section":"Section 4.5"}],"recommendation":"major_revision","confidential_remarks":"I see no reason to question the authors' integrity. The main methodological risk is the self-calibration of the empirical model on the same cluster and the absence of a forward-model test for the completeness-corrected q distribution; both are fixable with simulations. I would not reject the paper, but the central claims about the low-q slope and the three-segment shape should not be accepted until these tests are performed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should read this if you care about binary fractions in open clusters. The paper does something real: it pushes photometric detection of unresolved binaries down to q ~ 0.09/M1 over a mass range of 0.11–4.9 M_sun in the Pleiades, using a Bayesian multiband fit with an empirically corrected model. That is a step beyond previous optical-only or narrower mass-range studies.\n\nThe empirical ridgeline correction is a sensible way to handle model–observation mismatch, and the RV comparison for 17 systems with Torres et al. (2021) provides a genuine external check. The comparison with earlier binary fractions is honest, and the claim that optical-only studies miss low-q binaries is well supported.\n\nThe soft spot is the mass-ratio distribution. The three-segment power law is fit to the stacked PDF, but detectability is not a hard cutoff at qlim = 0.09/M1. A secondary near the hydrogen-burning limit adds only ~1–2% of the flux, comparable to photometric errors, so completeness is a smooth function of q and M1. The paper never forward-models this selection or tests the three-segment model against a single power law convolved with an estimated completeness. Given the small number of massive primaries probing q < 0.3, the low-q segment and the break near q = 0.32 could be artifacts of completeness rather than intrinsic features. The quoted uncertainties also ignore systematics from the empirical model, priors, and sample cuts; the f_b error bar is only Poisson.\n\nNone of this kills the central result. The binary fraction and its mass dependence are less vulnerable to the completeness critique, and the method itself is a clear improvement for cluster binary studies. I'd send this to a referee who will ask for a forward-model completeness test and a formal comparison against simpler alternatives—that is a revision, not a rejection. I'd bring it to reading group with a 'maybe'; the low-q question is worth arguing over.","headline":"A genuinely useful photometric binary census of the Pleiades, but the three-segment mass-ratio distribution needs a forward-model completeness test before being trusted.","tokens_in":24514,"tokens_out":2724,"would_cite":true,"duration_ms":25644,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A complete photometric census of unresolved binaries in the Pleiades gives a binary fraction of 0.34 and a mass-ratio distribution that rises, falls, and rises again across three power-law segments.","keywords":["open star clusters","binary stars","mass ratio distribution","binary fraction","Pleiades","photometric binaries","Bayesian statistics","main-sequence stars"],"falsifier":"Measure radial velocities or take high-resolution images of the Pleiades stars classified as low-mass single stars, especially those with inferred masses near the bottom of the main sequence where the near-infrared corrections are about 0.5 mag; if a significant fraction show companions with $q\\simeq0.2$--$0.4$, then the empirical ridgeline is biased and the claimed low-$q$ shape of the mass-ratio distribution is not physical.","tokens_in":23496,"feed_emoji":"✨","tokens_out":8395,"duration_ms":77078,"temperature":0.7,"pith_summary":"The paper attempts to establish a complete census of unresolved main-sequence binaries in the Pleiades: not just binaries with a massive companion, but all binaries down to the lowest stellar mass that can add light to the system. It combines optical and near-infrared magnitudes in a Bayesian six-band fit, using an empirical photometric model built by fitting the cluster's own single-star main sequence. The headline results are a binary fraction of $f_{\\rm b}=0.34\\pm0.02$ and a mass-ratio distribution that is not a single power law: it rises to $q\\simeq0.32$, falls to $q\\simeq0.79$, and rises again, with exponents $0.8\\pm0.2$, $-1.0\\pm0.3$, and $3.1\\pm1.0$. A complete low-mass-ratio census matters because low-$q$ binaries are the ones most sensitive to dynamical disruption, so the shape of the $q$ distribution is a direct test of how star formation and cluster dynamics shape binary populations.","feed_headline":"Pleiades stars: 34% harbor unresolved binary companions","feed_subtitle":"A combined optical and infrared fit reaches the lowest detectable companion masses and finds a three-part mass-ratio curve.","key_machinery":"The carrying mechanism is an empirical photometric model, $M_{\\rm emp}(M_s)=M_{\\rm model}(M_s)+\\Delta m(M_s)$, built for each band by fitting the ridgeline of the $\\delta m$--$M_s$ plane with a robust regression that discards outliers; this calibrates the theoretical isochrone to the actual single-star main sequence. A binary's combined magnitude is then $-2.5\\log_{10}(10^{-0.4M_{\\rm emp}(M_1)}+10^{-0.4M_{\\rm emp}(qM_1)})$, leaving only the primary mass $M_1$ and mass ratio $q$ as free parameters for each member. Individual Bayesian posterior densities $p_i(M_1,q)$ are stacked into a cluster-wide density $P(M_1,q)$, and all binary statistics, including $f_{\\rm b}$, the mass-ratio distribution, and their primary-mass dependence, are computed by integrating this stacked density rather than by thresholding best-fit values.","core_discovery":"The central claim is that unresolved main-sequence binaries in the Pleiades can be detected photometrically down to the theoretical lower mass limit, defined by a companion mass of $0.09\\,M_\\odot$, if optical and near-infrared magnitudes are fitted together against an empirical model calibrated on the cluster's own main sequence. For the 1154 main-sequence members, the paper reports a binary fraction $f_{\\rm b}=0.34\\pm0.02$ for all mass ratios above $q_{\\rm lim}(M_1)=0.09\\,M_\\odot/M_1$. Stacking the per-star posterior probability densities gives a mass-ratio distribution that rises as $q^{0.8\\pm0.2}$ below $q\\approx0.32$, falls as $q^{-1.0\\pm0.3}$ between $q\\approx0.32$ and $q\\approx0.79$, and rises again as $q^{3.1\\pm1.0}$ above $q\\approx0.79$: a deficiency of low-$q$ binaries and an excess of high-$q$ binaries relative to a single $q^{-1}$ power law. The paper also claims that binary fraction increases from about 0.18 to 0.94 with primary mass, that high-$q$ binaries are mostly low-mass primaries, and that using the Gaia RUWE value as a binary indicator catches only about one-third of photometric binaries.","pith_inferences":["Editorial inference: rebuilding the empirical model using only spectroscopically confirmed single stars would provide an independent check of whether the low-$q$ upturn and the $0.34$ fraction are robust; this is a test the paper does not perform.","Editorial inference: because the empirical corrections are calibrated on the Pleiades itself, applying the same six-band procedure to other clusters requires re-deriving each cluster's ridgeline, so the method's conclusions should be read as cluster-specific until recalibrated.","Editorial inference: if the excess of high-$q$ binaries is dominated by low-mass primaries, as the paper's stacked densities indicate, deeper surveys of very-low-mass companions in the Pleiades should find a higher companion frequency around low-mass stars than around solar-mass stars."],"forward_implications":["Optical-plus-infrared photometry pushes the detectable mass ratio to the theoretical stellar-mass limit, roughly doubling the low-$q$ reach of optical-only color-magnitude analyses.","Within the same $q>0.6$ selection used by optical-only studies, the paper reproduces previously published binary fractions, so earlier discrepancies are largely a mass-ratio coverage effect rather than a contradiction.","The mass-ratio distribution requires a three-segment power law rather than a single exponent, described as a fiducial $q^{-1}$ law with a low-$q$ deficit and a high-$q$ excess.","Binary fraction rises monotonically with primary mass, and the rise sits mainly in low- and intermediate-mass-ratio binaries, consistent with dynamical disruption of weakly bound systems plus the $q_{\\rm lim}$ selection effect.","The Gaia RUWE value flags only about one-third of the photometric binaries in this sample, so RUWE-based binary censuses are strongly incomplete."],"supporting_citations":[{"why":"Supplies the PARSEC theoretical model magnitudes and the 0.09 Msun lower mass limit that defines the detectable mass-ratio boundary.","marker":"Bressan et al. 2012"},{"why":"Provides the earlier empirical-isochrone correction strategy and the known rising binary-fraction trend with primary mass that this paper extends.","marker":"Li et al. 2020"},{"why":"Supplies the robust iterative-trimming Gaussian-process regression used to fit the single-star ridgeline in the delta-m versus M_s plane.","marker":"Li et al. 2021"},{"why":"Provides the independent spectroscopic binary sample with dynamical mass ratios used to validate the photometric q measurements and compare binary fractions.","marker":"Torres et al. 2021"},{"why":"Gives the optical-to-infrared pseudo-color binary fraction of the Pleiades that this paper reproduces within shared mass ranges and explains as including multiple systems and single-star contamination.","marker":"Malofeeva et al. 2023"},{"why":"Supplies the general mass-ratio distribution model whose three-part shape the Pleiades results are compared to.","marker":"Moe & Di Stefano 2017"},{"why":"Provides an optical-only CMD binary fraction for q>0.6 that the paper matches within the same selection.","marker":"Jadhav et al. 2021"},{"why":"Provides a mixture-model CMD measurement giving a roughly flat q distribution, used as a comparison target for the new three-segment result.","marker":"Niu et al. 2020"}],"fun_headline_variants":["Pleiades binary fraction 34%, down to 0.09 solar-mass companions","Three-part mass-ratio curve in Pleiades unresolved binaries","Photometric fit reveals Pleiades binary fraction and mass-ratio shape","Pleiades: 34% harbor companions, mass ratios show 3-segment pattern","Bayesian multiband analysis of Pleiades binaries: 34% fraction"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the ridgeline used to define single-star magnitudes is a true picture of single stars at every mass, including the low-mass range where the near-infrared corrections reach about 0.5 mag; if that ridgeline is pulled by unresolved binaries or model errors, the binary fraction and the low-mass-ratio rise would be biased without an independent check in that range.","fun_headline_variants_meta":{"raw":{"variants":["Pleiades binary fraction 34%, down to 0.09 solar-mass companions","Three-part mass-ratio curve in Pleiades unresolved binaries","Photometric fit reveals Pleiades binary fraction and mass-ratio shape","Pleiades: 34% harbor companions, mass ratios show 3-segment pattern","Bayesian multiband analysis of Pleiades binaries: 34% fraction"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000223,"raw_usage":{"total_tokens":1594,"prompt_tokens":1217,"completion_tokens":377,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":833,"completion_tokens_details":{"reasoning_tokens":274}},"tokens_in":833,"tokens_out":377,"duration_ms":4687,"temperature":1.0,"reasoning_tokens":274,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T22:23:40.136393+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure radial velocities or take high-resolution images of the Pleiades stars classified as low-mass single stars, especially those with inferred masses near the bottom of the main sequence where the near-infrared corrections are about 0.5 mag; if a significant fraction show companions with $q\\simeq0.2$--$0.4$, then the empirical ridgeline is biased and the claimed low-$q$ shape of the mass-ratio distribution is not physical.","supporting_citations":[],"review_version":1}