Pith. Sign in

REVIEW

A census of massive eclipsing binaries in the Whirlpool Galaxy

Solar and Stellar Astrophysics · Submitted 16 Sep 2026 · 136 references

Pith's verdict: Impressive first extragalactic massive EB census, but the central demographic result has an unquantified systematic from a radius mismatch between the forward model and the data that could shift the close-binary fraction by ~30%.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on September 18, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 6 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2609.19268 (version 1)
Pith ID4SW3X45E, Pith's permanent record of this paper and its review
CitationAlessandro Savino, Charlie Conroy, Cheyanne Shariat, Daniel R. Weisz, Kareem El-Badry, Maude Gull. "A census of massive eclipsing binaries in the Whirlpool Galaxy." arXiv:2609.19268 (2026). https://pith.science/paper/4SW3X45E
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

A census of massive eclipsing binaries in the Whirlpool Galaxy (M51) finds that 27.4% ± 7.1% of massive stars are in close binaries with periods of 2–20 days, and the period distribution is consistent with being log-flat, suggesting that…

The reading

The paper aims to establish the intrinsic close-binary demographics of massive stars in M51, a Milky-Way-like spiral galaxy, by identifying 173 eclipsing binary candidates from HST light curves and correcting for observational biases through detailed injection-recovery simulations. It seeks to answer whether the close-binary fraction and period distribution vary with metallicity by comparing results to the Milky Way, LMC, and SMC. If correct, it provides the first such census beyond the Local Group and shows that close massive binary formation is largely independent of metallicity. This matters because binary interactions shape the evolution of massive stars and affect the production of transients, compact objects, and ionizing radiation.

What carries the argument

The core machinery is a two-step selection-function correction built on synthetic binary populations generated with the COSMIC code and propagated through the HST survey. First, photometric completeness is quantified via artificial star tests injected into the real images. Second, light-curve recovery is measured by injecting synthetic eclipsing binaries into real carrier-star light curves, running a Lomb–Scargle period search, and visually classifying the folded light curves. The observed EB count and period distribution are then compared to the forward-modeled predictions via a binomial likelihood to infer the intrinsic close-binary fraction and period-slope parameter π.

Extended reading notes

Core claim

By combining 34 epochs of HST/ACS red-optical photometry with single-epoch HST/WFC3 UV-optical photometry, the authors identify 173 massive eclipsing binary candidates in M51 with orbital periods 1–30 days. After correcting for selection effects using forward modeling and injection-recovery simulations, they infer that 27.4% ± 7.1% of massive stars reside in close binaries (2 ≤ P_orb/d < 20, mass ratio ≥ 0.1), and the intrinsic period distribution is approximately log-flat (π = 0.29 ± 0.32). These demographics are similar to those in the Milky Way, LMC, and SMC, despite a fivefold range in metallicity, indicating that close massive binary formation depends only weakly on metallicity.

Load-bearing premise

The entire inference of the intrinsic close-binary fraction and period distribution depends on the forward model of the binary population and selection function, which assumes a constant star formation rate, solar metallicity, the Moe & Di Stefano 2017 multiplicity distributions, and an exponential dust+spatial model; if any of these are wrong, the corrected demographics could be systematically biased.

Editorial extensions

If this is right

  • If the close-binary fraction in M51 is ~27%, then roughly 990 out of 3616 massive stars in the parent sample are in close binaries, of which only 50 are detected as EBs, implying a strong selection bias that must be accounted for in all extragalactic EB surveys.
  • The log-flat period distribution (π ≈ 0) means that close binaries are equally likely per logarithmic period interval, which matches the Milky Way and LMC but differs from the SMC early-B EB sample (π ≈ –0.9), hinting at a possible metallicity or age effect that warrants further study.
  • The weak metallicity dependence of close-binary formation across 0.2–1.1 Z_sun suggests that binary evolution models for high-redshift (low-metallicity) star formation can safely adopt local-Universe multiplicity distributions.
  • The catalog of ~100,000 luminous stars and variables in M51 provides a public reference for studies of stellar variability, Cepheids, and evolved massive stars beyond the Local Group.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Because the M51 EB sample is weighted toward luminous, partly evolved stars (median radii ~3 times larger than Local Group EBs), the inferred close-binary fraction might be a lower limit for the zero-age main-sequence population; a dedicated search for main-sequence dwarfs would test this.
  • The similarity of close-binary demographics across galaxies with different star-formation histories and environments implies that the physics of close binary formation – likely disk fragmentation and migration – is robust against global galaxy properties, a prediction that could be tested with radiation-hydrodynamic simulations of star formation.
  • Extending this method to very metal-poor Local Group galaxies (e.g., WLM, Z ~ 0.14 Z_sun) could reveal a threshold metallicity where binary properties change, which would directly inform models for Population III star binaries and their role in the first supernovae.
Original abstract

We present a census of massive eclipsing binary (EB) candidates in the Whirlpool Galaxy (M51), a face-on spiral galaxy at 7.5 Mpc with $Z\sim Z_\odot$. Using 34 epochs of red-optical Hubble Space Telescope (HST) photometry and a single epoch of broadband UV-optical photometry, we identify 173 massive EB candidates ($>5~{\rm M_\odot}$) with orbital periods $P_{\rm orb}$=1-30 days. The sample traces young stellar populations across the spiral arms of M51. Joint light-curve and spectral energy distribution fits for 136 systems yield binaries spanning $T_{\rm eff}=10{,}000$-$50{,}000$ K, $R=7$-$66~{\rm R_\odot}$, and masses of $5$-$120~{\rm M_\odot}$. We use detailed injection-recovery simulations to quantify the survey's completeness and infer intrinsic properties of the massive binary population. At the population level, we infer a close-binary fraction of $f_{\rm close}=27.4\pm7.1\%$ for $2\leq P_{\rm orb}/{\rm d}<20$ and $0.1\leq q\leq1$. Over the same period range, the intrinsic period distribution is consistent with being log-uniform, ${\rm d}N/{\rm d}\log P_{\rm orb}\propto P_{\rm orb}^{\pi}$ with $\pi=0.29\pm0.32$. These demographics broadly match those of massive-star populations in the Milky Way, LMC, and SMC, with tentative evidence for a lower close-binary fraction at lower metallicity. Across the $\sim0.2$-$1.1~Z_\odot$ range spanned by these galaxies, the overall similarities suggest that the formation of close massive binaries depends only weakly on metallicity. Our survey establishes EBs as a powerful tool for measuring massive-binary demographics beyond the Local Group and provides a rich panchromatic census of $\sim$100,000 luminous stars and variables in M51.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-23: 1 request, 8,163 input tokens.

What the abstract itself states, as probabilities: claims first 0.06 · states limitation 0.51 · human participants 0.01 · controlled comparison 0.16 · randomized assignment 0.01 · states quantitative result 0.99 · compares with prior results 0.93 · mentions code or data availability 0.03.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

We present a census of massive eclipsing binary (EB) candidates in the Whirlpool Galaxy (M51), a face-on spiral galaxy at 7.5 Mpc with $Z\sim Z_\odot$. Using 34 epochs of red-optical Hubble Space Telescope (HST) photometry and a single epoch of broadband UV-optical photometry, we identify 173 massive EB candidates ($>5~{\rm M_\odot}$) with orbital periods $P_{\rm orb}$=1-30 days. The sample traces young stellar populations across the spiral arms of M51. Joint light-curve and spectral energy distribution fits for 136 systems yield binaries spanning $T_{\rm eff}=10{,}000$-$50{,}000$ K, $R=7$-$66~{\rm R_\odot}$, and masses of $5$-$120~{\rm M_\odot}$. We use detailed injection-recovery simulations to quantify the survey's completeness and infer intrinsic properties of the massive binary population. At the population level, we infer a close-binary fraction of $f_{\rm close}=27.4\pm7.1\%$ for $2\leq P_{\rm orb}/{\rm d}<20$ and $0.1\leq q\leq1$. Over the same period range, the intrinsic period distribution is consistent with being log-uniform, ${\rm d}N/{\rm d}\log P_{\rm orb}\propto P_{\rm orb}^{\pi}$ with $\pi=0.29\pm0.32$. These demographics broadly match those of massive-star populations in the Milky Way, LMC, and SMC, with tentative evidence for a lower close-binary fraction at lower metallicity. Across the $\sim0.2$-$1.1~Z_\odot$ range spanned by these galaxies, the overall similarities suggest that the formation of close massive binaries depends only weakly on metallicity. Our survey establishes EBs as a powerful tool for measuring massive-binary demographics beyond the Local Group and provides a rich panchromatic census of $\sim$100,000 luminous stars and variables in M51.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 We present a census of massive eclipsing binary (EB) candidates in the Whirlpool Galaxy (M51), a face-on spiral galaxy at 7.5 Mpc with $Z\sim Z_\odot$. contribution 0.98 0.77 0.03 asserted 0.98
s1 Using 34 epochs of red-optical Hubble Space Telescope (HST) photometry and a single epoch of broadband UV-optical photometry, we identify 173 massive EB candidates ($>5~{\rm M_\odot}$) with orbital periods $P_{\rm orb}$=1-30 days. method 0.45 0.98 0.06 asserted 1.00
s2 The sample traces young stellar populations across the spiral arms of M51. result 0.65 0.76 0.04 asserted 1.00
s3 Joint light-curve and spectral energy distribution fits for 136 systems yield binaries spanning $T_{\rm eff}=10{,}000$-$50{,}000$ K, $R=7$-$66~{\rm R_\odot}$, and masses of $5$-$120~{\rm M_\odot}$. result 0.98 0.96 0.06 asserted 1.00
s4 We use detailed injection-recovery simulations to quantify the survey's completeness and infer intrinsic properties of the massive binary population. method 0.99 0.97 0.04 asserted 0.98
s5 At the population level, we infer a close-binary fraction of $f_{\rm close}=27.4\pm7.1\%$ for $2\leq P_{\rm orb}/{\rm d}<20$ and $0.1\leq q\leq1$. result 1.00 0.96 0.28 asserted 0.99
s6 Over the same period range, the intrinsic period distribution is consistent with being log-uniform, ${\rm d}N/{\rm d}\log P_{\rm orb}\propto P_{\rm orb}^{\pi}$ with $\pi=0.29\pm0.32$. result 1.00 0.95 0.09 hedged possibility 0.99
s7 These demographics broadly match those of massive-star populations in the Milky Way, LMC, and SMC, with tentative evidence for a lower close-binary fraction at lower metallicity. interpretation 0.57 0.81 0.21 hedged possibility 0.80
s8 Across the $\sim0.2$-$1.1~Z_\odot$ range spanned by these galaxies, the overall similarities suggest that the formation of close massive binaries depends only weakly on metallicity. interpretation 1.00 0.78 0.34 hedged possibility 0.90
s9 Our survey establishes EBs as a powerful tool for measuring massive-binary demographics beyond the Local Group and provides a rich panchromatic census of $\sim$100,000 luminous stars and variables in M51. interpretation 0.55 0.95 0.50 asserted 1.00

What kind of paper the classifier reports this is: empirical study 0.70 · evidence from observation or experiment 0.76.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. Optical color–absolute magnitude diagram com￾paring eclipsing binaries in M51 to those in the local Gaia sample. The upper region shows luminous HST/ACS sources in M51 (black) and the EB candidates detected among them (blue, this work). The lower region shows nearby Gaia sources (black) and all eclipsing binary candidates with ϖ/σϖ ≥ 10 (magenta; Gaia Collaboration et al. 2023; Mowlavi et al. 2023). M51 EBs are exclusively among the most massive and luminous binaries relative to local samples. high-mass tail of the IMF and ∼ 20× more of the most massive stars. Conroy et al. (2018) s
Figure 2. HST footprint of our M51 survey, shown atop an optical image of M51. Cyan outlines show the 34 ACS/WFC visits that provide the F606W and F814W time series, and the yellow outline shows the Cycle 30 WFC3/UVIS pointing. The background RGB image is reconstructed using F814W+F658N (R), F555W (G), and F435W (B) from the HST program of Mutchler et al. (2005). crowded field photometric reductions. From these de￾faults, we apply three main changes motivated by the new DOLPHOT 3.1 version and its specific ap￾plication to the M51 time domain. First, we use SubResRef=1, SecondPass=1, and…
Figure 3. Selection of point sources from HST/ACS photometry. Left: all DOLPHOT sources identified in M51, with median photometric errors shown in purple crosses. Middle: only sources that pass the point-source DOLPHOT quality cuts (outlined in Section 2.2). Right: the eclipsing binary candidates (blue) identified among these point sources. The gray shaded region shows the parent sample of blue sources from which we searched for EBs. We also overlay solar-metallicity MIST II evolutionary tracks for 8, 20, 40 M⊙ stars. All magnitudes are corrected for foreground extinction. atively faint in…
Figure 4. Completeness and photometric uncertainty for the M51 HST photometry computed from artificial star tests. The top panel shows the completeness function for each filter. The bottom panels show the difference between recovered and input magnitudes for ASTs that pass the same quality cuts applied to the science photometry. The 90th percentile magnitudes for EB candidates are marked in vertical lines. The ridge-like features are due to the in￾put stellar catalog (see text). The 50% completeness lim￾its are mF606W = 27.10 and mF814W = 25.96 for ACS, and mF275W = 24.76, mF336W = 25.56,…
Figure 5. Spatial completeness maps derived from artificial star tests. For each filter, we show the recovered fraction (i.e., completeness) in 40 × 40 spatial bins among artificial stars brighter than that filter’s 50% completeness magnitude (m50;
Figure 6. Examples of rejected and accepted EB candidate light curves. The top row shows HST/ACS light curves for three examples of rejected candidates. The first system shows sawtooth-like variability consistent with a pulsator, the second is likely an alias, and the third shows shallow sinusoidal variability. The bottom row shows three examples of selected EB candidate light curves at different periods. These show coherent, smooth variability with two apparent eclipses per orbital cycle. tice, we find that the orbital period is usually twice the highest-power Lomb-Scargle period. During…
Figure 7. Example light curve and SED joint fit. The left panel shows the ACS F606W/F814W light curves and phased WFC3 points folded on the fixed orbital period. Thin curves show posterior draws, and thick curves show the maximum-posterior model. The right panel shows the WFC3 SED evaluated with the same T1, T2, R1, R2, AV , inclination, and WFC3 phase offset used in the light curve model. tion, and phase of primary eclipse, while the WFC3 pho￾tometry constrains the temperatures, radii, and extinc￾tion through the UV–optical SED. The WFC3 points are evaluated at their observed phases after…
Figure 8. Sensitivity of our binary search, as determined by injection–recovery tests. The color shows the probabil￾ity that a close binary with a given orbital period, inclina￾tion, and properties similar to the observed EBs is recovered in our HST/ACS light curves (PLC; Equation (11)). Low￾inclination systems have small detection probabilities be￾cause they rarely eclipse, and long-period systems are harder to recover because the eclipse duty cycle and phase coverage decrease. 3.3.3. Summary
Figure 9. Spatial distribution of the massive eclipsing binary candidates. The grayscale image shows the den￾sity of all luminous sources mF814W,0 < 26, while the pink contours show the density of blue luminous sources mF606W,0 − mF814W,0 < 0.8, from which the EB candidates were selected. Blue points mark EB candidates; the larger white-filled symbols are EBs with only ACS photometry, and the filled symbols have WFC3 photometry as well. The dashed red outline shows the WFC3 footprint. EB candi￾dates trace the star-forming spiral arms of M51, which coin￾cide with the blue massive star…
Figure 10. Panchromatic color-magnitude diagrams of the M51 sample from HST/WFC3. Gray points show field sources that pass WFC3 star cuts, and blue points show the Gold EB candidates. From left to right, the panels show NUV–optical (left), blue optical–red optical (left middle), optical only (right middle), and red optical only (right). All observed magnitudes have been corrected for Milky Way foreground extinction and placed on the infinite aperture system. Orange curves show example MIST evolutionary tracks at the distance of M51 for solar-metallicity stars of mass 8, 12, 20, 30, 80 M⊙ and…
Figure 11. Observed orbital period distribution for the gold eclipsing binary sample. The blue histogram shows the EB candidates at their adopted orbital periods. The dashed ver￾tical line marks Porb = 2 d. The observed sample is strongly weighted toward short periods, consistent with geometric se￾lection biases favoring close binaries. identify a separate set of 10 short-period EB candidates with Porb = 1.3–2 d (red dashed line), which are not used for the close binary demographics in the following sections due to high alias contamination in this region from the sparse ACS light curves. The…
Figure 12. Stellar parameters for the EB candidates from the joint light curve+SED fits. Diagonal panels show the 1D distributions and off-diagonal panels show 2D parameter spaces. Each point represents the best-fit parameters of one eclipsing binary. We include a dust-corrected WFC3 blue-optical CMD in the top right, which corrects each EB by its fitted local AV . Gray points show the observed WFC3 field population, and orange curves show solar-metallicity MIST tracks for 8, 12, 30, 80 M⊙ stars at the distance of M51 without internal extinction. The fitted EB population hosts hot and…
Figure 13. HR diagram and apparent ages of massive binaries in M51. Panel A shows both components in the HR diagram, colored by the inferred MIST age, with solar-metallicity MIST tracks overplotted. The dashed black curve shows the classical Humphreys–Davidson (HD) luminosity limit (Humphreys & Davidson 1979; Sabhahit & Vink 2026). Panel B shows the difference between the two component ages for each binary, sorted by the mean Teff . Each point represents one binary with the posterior medians (blue) and 1σ range shown. Most systems are consistent with coeval components at the current level of…
Figure 14. 2D posterior for the intrinsic period slope (π) and close-binary fraction (fclose). Filled contours enclose the 68% and 95% density regions. The covariance arises be￾cause longer-period-weighted populations have lower selec￾tion probabilities and therefore require larger fclose. region: LN (π, fclose) =  Nparent NEB  [fcloseS(π)]NEB × [1 − fcloseS(π)]Nparent−NEB . (17) This likelihood takes the form of a binomial distribution where fcloseS(π) is the probability that a system in the parent sample contains a close binary that is detected as an EB. The observed periods primarily…
Figure 15. Observed and forward-modeled period distribu￾tions over 2 ≤ Porb/d < 20. The dashed orange curve shows the intrinsic COSMIC binary population, while the solid black curve shows the same population after applying the survey selection function. Blue points show the observed 50 Gold EBs in the primary main-sequence sample. The observed period distribution of massive binaries in M51 is consistent with a roughly log-flat (π = 0.29 ± 0.32) intrinsic period dis￾tribution. F336W CMD, corresponding to a raw EB fraction of 1.38%. The likelihood in Equation (15) constrains the fraction of…
Figure 16. Close-binary fractions (top) and period slopes (bottom) for massive stars in M51 and nearby galaxies, over Porb = 2– 20 d and q ≥ 0.1. The period slope (π) is defined by dN/d log Porb ∝ P π orb, with the dotted line marking a log-flat distribution (π = 0). ‘EB’ and ‘RV’ denote results from eclipsing-binary and radial-velocity surveys, whereas ‘O’ and ‘B’ refer to the approximate spectral types of each sample. Approximate metallicities are shown above each host galaxy (Croxall et al. 2015; Russell & Dopita 1992; Hunter et al. 2007). The early-B EB results and Galactic O-star…
Figure 17. Comparing M51 EB parameters to massive-binary samples in the MW, LMC, and SMC. Diagonal panels show 1D distributions; off-diagonal panels show 2D scatter plots. Blue points and histograms show the M51 parameters. Red markers show EB results from 30 Dor (Tarantula Nebula), LMC, and SMC samples (Mahy et al. 2020; Harries et al. 2003; Hilditch et al. 2005; North et al. 2010; Taormina et al. 2024). Green open markers show SB2 results from the MW, Cyg OB2, and 30 Dor (Sana et al. 2012; Kobulnicky et al. 2014; Almeida et al. 2017). For all samples, M1 is the more massive component and q…
Figure 18. Spatial distribution of the simulated close bina￾ries before (top) and after (bottom) photometric selection. The intrinsic systems follow an exponential disk + spiral-arm stellar distribution. using hz,young = 100 pc. These values are based on measurements of young stellar disks and 3D radiative transfer fits to edge-on galaxies (Bahcall & Soneira 1980; Xilouris et al. 1999; Bianchi 2007; Schechtman-Rook & Bershady 2013; De Geyter et al. 2014; De Looze et al. 2014). Lastly, we convolve this exponential profile with a Gaussian kernel density estimate (KDE) constructed using the blue…
Figure 19. Comparison between the forward-modeled binary population (blue) and the 50 Gold EBs used for the primary demographic inference (black). All systems satisfy 2 ≤ Porb/d < 20, F336W < 24, and lie blueward of the TAMS in the F336W − F606W vs. F336W CMD. Diagonal panels show normalized 1D distributions. Off-diagonal panels show 2D scatter plots. The model spans similar parameter ranges to the observed sample, although it favors somewhat smaller radii. are nearly identical (AV = 0.424 and 0.422 mag), al￾though the simulated distribution is slightly narrower. The largest difference is in…
Figure 20. Extinction as a function of the temperature of the hotter component for the M51 EB candidates. Points show the posterior medians, with error bars spanning the 16th–84th percentiles. The weak correlation (ρS = 0.49) and the presence of hot systems across a broad range of AV indi￾cate that the highest fitted temperatures are likely intrinsic, and not exclusively driven by large extinction. and AV is expected, the broad scatter and the pres￾ence of high-temperature systems at modest extinction indicate that the hottest inferred temperatures are not driven solely by the…
Figure 21. WFC3 CMD of point sources within 0.5 ′′ of a cataloged Chandra source. The grayscale background shows all WFC3 sources that pass the DOLPHOT point-source cuts in F475W and F814W. Colored points show the 345 ACS– Chandra matches that pass the same cuts, colored by the cataloged X-ray luminosity. Magnitudes are corrected for Milky Way foreground extinction. The broad CMD distribu￾tion illustrates that the X-ray matches have possible optical counterparts spanning a wide range of colors and evolution￾ary states. D.2. Magellanic Clouds For O stars, we use the published epoch RVs and…
Figure 22. ACS phase-folded light curves for Gold EB candidates, sorted by increasing orbital period. Black points show F606W and blue points show F814W, folded on the adopted orbital period and repeated over two cycles

Showing 22 of the paper's 25 figures. See all figures in the PDF ↗ Figures reproduced from arXiv:2609.19268 by the authors.

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.