{"id":"193bd2d5-a75b-48c4-a734-5bcd5b373ad6","arxiv_id":"2506.22673","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Open clusters produce white dwarf-white dwarf merger rates of 46 to 1,400 per cubic gigaparsec per year, including super-Chandrasekhar mergers that may be fast radio burst progenitors.","lead":"This paper simulates open star clusters with masses of 100, 1,000, and 10,000 suns and counts how often their dead stars, such as white dwarfs, neutron stars, and black holes, collide and merge. It finds that white dwarf mergers from these common clusters may happen more often than previously thought and could explain some fast radio bursts.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quoted WD-WD and super-Chandrasekhar rates scale linearly with the assumed 100% initial binary fraction; realistic binary fractions and mass-function-weighted f_cluster may put the lower bound near the globular-cluster rate.","rationale":"I read the paper as a solid, well-scoped N-body study of open clusters, with public code and an explicit limitation section. The WD-WD and WD-NS channels are genuinely underexplored, and the comparison to isolated BSE runs is a useful control. The strongest claim is quantitative: open clusters produce 70–780 Gpc^-3 yr^-1 super-Chandrasekhar WD-WD mergers, larger than globular clusters. The most load-bearing condition is that the simulated efficiency is representative of real open clusters, which in turn requires the 100% binary fraction and the f_cluster normalization to be jointly realistic. Neither alone would flip the claim (halving the binary fraction still leaves the lower bound above 10), but together they can place the lower bound at the globular-cluster level. The reader's conditional verdict already flags the binary fraction; I agree and sharpen it by pointing to the joint normalization. A simple rescaling with observational binary fractions and mass-function-integrated f_cluster would settle whether the claim survives. I do not see an internal logical error; the concern is about external calibration. Therefore the verdict remains CONDITIONAL, and my read does not change the reader's recommendation.","tokens_in":20881,"tokens_out":19255,"duration_ms":219506,"concrete_test":"Rescale the Table 3 super-Chandrasekhar efficiencies (3.5e-4 and 3.8e-4 M_sun^-1 for 10^3 and 10^4 M_sun) by binary-fraction factors 0.5 and 0.3 and by f_cluster values 0.01, 0.03, and 0.1 for the combined 10^3–10^4 M_sun mass range, using psi(0) ≈ 1.85e-3 M_sun yr^-1 Mpc^-3. If the resulting rates remain above ~10 Gpc^-3 yr^-1 for f_cluster ≥ 0.03 and binary fraction ≥ 0.5, the qualitative dominance over globular clusters stands; if they fall to ~10 or below, the abstract's rate range should be presented as an upper envelope rather than a fiducial estimate.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline rate range for super-Chandrasekhar WD-WD mergers (70–780 Gpc^-3 yr^-1, exceeding the ~10 Gpc^-3 yr^-1 globular-cluster value) is the product of three factors: merger efficiency per solar mass from the N-body runs, the local star-formation rate density, and f_cluster. The efficiency is computed for a 100% initial binary fraction (Section 2). For the intermediate-mass primaries that dominate WD-WD mergers (roughly 2–8 M_sun), observed binary fractions are typically 0.3–0.7, so the primordial contribution to eta is overestimated by a factor of 1.4–3. The quoted ranges then assign f_cluster = 0.1–1 separately to each cluster mass (Section 3.3), although f_cluster ≈ 0.1 is the total cluster formation efficiency for a Milky Way-like galaxy; once the cluster mass function (dN/dM ∝ M^-2) is integrated, the fraction of star formation in 10^3–10^4 M_sun clusters alone is plausibly only ~0.02–0.05. Taking binary fraction 0.5 and f_cluster = 0.03 for the relevant mass range reduces the super-Chandrasekhar lower bound from 70 to roughly 10 Gpc^-3 yr^-1, at the globular-cluster level. This is not an internal inconsistency, but it means the 'larger than globular clusters' claim is not robust to jointly realistic normalization choices. The paper's own BH-BH range at f_cluster=1 (20–200 Gpc^-3 yr^-1) already violates the LIGO-Virgo upper limit (17.9–44), confirming that the f_cluster=1 end of the quoted ranges is not physically viable for at least one channel.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"Using the PETAR N-body code with SSE/BSE stellar and binary evolution, the authors simulate open clusters of initial masses 10^2, 10^3, and 10^4 M_sun (596, 729, and 17 realizations, respectively), assuming a 100% primordial binary fraction, solar metallicity, and a Milky Way potential. They compare these cluster runs with isolated binary evolution (\"BSE\") and report merger efficiencies per solar mass for WD-WD, WD-NS, BH-BH, BH-NS, and NS-NS systems. These efficiencies are converted into local volumetric merger rates using the Madau-Fragos star formation rate density and f_cluster = 0.1-1, yielding WD-WD rates of 46-460, 130-1300, and 140-1400 Gpc^-3 yr^-1 for the three cluster masses, super-Chandrasekhar WD-WD rates of 70-780 Gpc^-3 yr^-1 (claimed to exceed globular-cluster estimates), CO WD-WD rates of 0.14%-2.6% of the observed Type Ia supernova rate, and BH-NS mergers only in the 10^3 M_sun clusters.","tokens_in":21202,"tokens_out":9735,"duration_ms":105565,"significance":"The paper fills a genuine gap by simulating WD-WD and WD-NS mergers in low-mass open clusters and by providing many realizations of stochastic low-mass clusters. The use of well-tested public codes, the explicit comparison with isolated binary evolution, and the authors' candid caveats about small samples are strengths. If the normalization assumptions were secure, the claimed open-cluster contribution to super-Chandrasekhar WD-WD mergers would be an important result for FRB progenitor studies. However, the absolute rate claims, and in particular the \"larger than globular clusters\" conclusion, are not robust to jointly realistic choices of the primordial binary fraction and the cluster formation efficiency, as detailed below.","major_comments":[{"comment":"The load-bearing assumption of a 100% primordial binary fraction is not tested or propagated into the quoted rates. Since only 2%-6% of WD-WD mergers are dynamically formed (Section 3.2.1), the WD-WD, WD-NS, and BH-NS efficiencies in Table 2 scale almost linearly with the assumed binary fraction. Observed binary fractions for the intermediate-mass stars that dominate WD-WD progenitors are typically 0.3-0.7, implying an overestimate by factors of roughly 1.4-3. The manuscript itself uses this assumption to explain why its BH-BH efficiency is about ten times larger than that of Kumamoto et al. (2019) (Section 3.2.3). A sensitivity study or a rate quoted for a fiducial binary fraction is needed before the headline \"larger than globular clusters\" statement can be considered robust.","section":"Section 2 / Table 2 / Section 3.2.1"},{"comment":"The rate ranges use f_cluster = 0.1 and f_cluster = 1 as lower and upper bounds, applied separately to each cluster mass. f_cluster = 1 is not a physically meaningful upper bound: it would require that all star formation occur in clusters of a single mass, and the paper's own BH-BH rate at this end (20-200 Gpc^-3 yr^-1) violates the LIGO-Virgo constraint of 17.9-44 Gpc^-3 yr^-1 (Section 3.3). The mass-dependent normalization behind R_Local is likewise not an observationally grounded cluster formation efficiency: setting f_cluster,100 + f_cluster,1000 + f_cluster,10000 = 1 with f_cluster,M proportional to M^-1 assigns about 90% of all star formation to 10^2 M_sun clusters. The authors should derive rates from a continuous, observationally motivated cluster mass function and state explicitly what fraction of star formation passes through the simulated mass range.","section":"Section 3.3, Eq. (2), and Table 2"},{"comment":"The headline super-Chandrasekhar WD-WD rate range (70-780 Gpc^-3 yr^-1) combines the two uncalibrated normalizations above with an efficiency that only counts mergers occurring within 500 Myr or 1 Gyr, while the merger-time distributions in Figure 3 are still rising at the end of the simulations. The reported rates are therefore truncated estimates for the adopted model, not full-lifetime yields. With a more realistic binary fraction of about 0.5 and a mass-function-weighted f_cluster of order 0.03 for the relevant cluster mass range, the lower end of the quoted range drops to roughly the globular-cluster value, so the qualitative conclusion that open clusters dominate the super-Chandrasekhar WD-WD channel is not yet established.","section":"Section 3.3 and Section 3.4.1"}],"minor_comments":[{"comment":"The displayed formula for the rate density is garbled: the text reads \"R = ... z f_cluster\" and omits both the efficiency eta and the star formation rate density psi. Please restate the equation with all symbols defined.","section":"Eq. (2)"},{"comment":"Several entries in Table 2 are based on fewer than ten events (bracketed values), including the BH-NS and NS-NS channels. The text acknowledges limited sample sizes, but the abstract and conclusions quote these rates without uncertainties; Poisson confidence intervals should be reported and propagated into the quoted ranges.","section":"Table 2 and Section 3.3"},{"comment":"The conclusion that \"BH-NS mergers only occur in 10^3 M_sun clusters\" is based on nine events in that model and zero events in the 10^4 M_sun model; given the smaller total sample mass of the latter, the Poisson upper limit is not strongly inconsistent with the 10^3 M_sun efficiency. This claim should be softened or accompanied by a statistical test.","section":"Section 3.2.4"},{"comment":"The text should clarify that the mass-dependent normalization f_cluster,M proportional to M^-1 is normalized to unity only over the three discrete cluster masses, which is not equivalent to integrating a continuous cluster mass function with dN/dM proportional to M^-2.","section":"Section 3.3, R_Local"},{"comment":"The figures use vertical lines to indicate single events, but the event counts per channel and per cluster mass are never listed explicitly; adding the raw counts to Table 2 would make the statistical weight of each efficiency transparent.","section":"Figure 3 / Figure 4"}],"recommendation":"major_revision","confidential_remarks":"The simulations themselves are a useful resource, and the internal comparison between cluster dynamics and isolated binary evolution is credible. My main concern is that the abstract and conclusions present absolute rate ranges that are dominated by two uncalibrated assumptions (100% binary fraction and f_cluster = 0.1-1), and the paper's own BH-BH rate at f_cluster = 1 already violates LIGO-Virgo limits. I request a revision that either rescales the central rates to fiducial observed binary fractions and a physically motivated cluster formation efficiency, or explicitly labels the headline numbers as maximal estimates under optimistic assumptions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read the Cary et al. paper. The useful new thing is the first systematic N-body treatment of WD-WD and WD-NS mergers in 10^2-10^4 M_sun open clusters, with a clean comparison to isolated BSE runs. Previous open cluster work focused on BH-BH and BH-NS, so this fills a real gap. The qualitative trends look solid: BH-BH are dynamically enhanced only in 10^4 M_sun clusters, BH-NS show up only in 10^3 M_sun clusters, and WD-WD mergers are slightly suppressed by dynamics relative to isolated evolution.\n\nThe soft spot is the normalization. The headline rates scale linearly with the assumed 100% primordial binary fraction, and the authors admit this produces a BH-BH efficiency about ten times Kumamoto et al. 2019. Observed binary fractions for the relevant intermediate-mass stars are more like 0.3-0.7, so the efficiencies are likely overestimated by a factor of 1.4-3. The rate ranges also use f_cluster=0.1-1, but f_cluster=1 is already inconsistent with LIGO-Virgo BH-BH upper limits (the paper notes this), and f_cluster=0.1 is the total cluster formation efficiency, not the fraction in one mass bin. A mass-function-weighted f_cluster for 10^3-10^4 M_sun clusters is more like 0.02-0.05, which would put the super-Chandrasekhar WD-WD lower bound near the globular cluster value. The BH-NS rate for 10^3 M_sun clusters rests on nine events; several other channels have even less.\n\nSo the qualitative story and the methodology are worth keeping, but the abstract's 'higher than globular clusters' claim is not robust to jointly realistic choices. The paper needs a binary fraction sensitivity test, or at least a clearly stated caveat that the rates are upper-ish estimates.\n\nWho is this for: people working on compact merger rate models, LISA predictions, and FRB progenitor channels. It deserves a serious referee. After a revision that addresses the normalization, it can be a solid contribution. I would send it to review.","headline":"First systematic N-body look at WD-WD and WD-NS mergers in low-mass open clusters, with interesting FRB implications, but headline rates are not robust to binary fraction and cluster formation efficiency choices.","tokens_in":21843,"tokens_out":7829,"would_cite":true,"duration_ms":73873,"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":"Open clusters, though short-lived, contribute compact-object merger rates that rival or exceed globular clusters, especially for white-dwarf-white-dwarf mergers that may power fast radio bursts.","keywords":["open star clusters","compact binary mergers","white dwarf mergers","merger rate density","fast radio bursts","Type Ia supernovae","N-body simulations","black hole-neutron star mergers"],"falsifier":"Measure the primordial binary fraction in young open clusters with masses near $10^2$ to $10^4$ solar masses, for example through spectroscopic binary surveys of embedded clusters; if the fraction is substantially below 100%, the reported WD-WD, WD-NS, and BH-NS rates should be scaled down proportionally, and a fraction near 50% would pull the super-Chandrasekhar WD-WD rate toward the globular-cluster comparison level.","tokens_in":20604,"feed_emoji":"🌌","tokens_out":7030,"duration_ms":66940,"temperature":0.7,"pith_summary":"The paper argues that previous merger-rate censuses missed a major channel: open clusters, which are far more numerous than globular clusters and form throughout cosmic time. From N-body simulations of $10^2$, $10^3$, and $10^4$ solar-mass open clusters with full stellar evolution, the authors derive merger efficiencies per solar mass and convert them into local-universe rate densities. They find white-dwarf-white-dwarf mergers at 46-1400 Gpc$^{-3}$ yr$^{-1}$, higher than globular-cluster estimates, with super-Chandrasekhar WD-WD mergers at 70-780 Gpc$^{-3}$ yr$^{-1}$ as candidate fast-radio-burst progenitors. They also find BH-NS mergers only in intermediate-mass $10^3$ solar-mass clusters and dynamically formed, sometimes eccentric BH-BH mergers only in $10^4$ solar-mass clusters. A sympathetic reader would care because the results reposition open clusters as nonnegligible contributors to the transient sky.","feed_headline":"Open clusters rival globulars as white-dwarf merger factories","feed_subtitle":"White-dwarf mergers from simulated open clusters hit 46-1400 per cubic gigaparsec per year, above globular-cluster estimates.","key_machinery":"The carrying object is the merger efficiency $\\eta$, the number of mergers per solar mass of cluster, extracted from N-body simulations and converted into a cosmic rate density $R = \\eta \\, \\dot{\\rho}_*(z) \\, f_{\\rm cluster}$ using a cosmic star formation rate density and a cluster formation efficiency $f_{\\rm cluster}$ bracketed between 0.1 and 1. The simulations evolve hundreds of open clusters of $10^2$, $10^3$, and $10^4$ solar masses with a tree-based N-body integrator, single and binary stellar evolution, natal kicks, and a Milky Way tidal field; matched isolated-binary runs provide the control that isolates the dynamical effect of the cluster environment.","core_discovery":"The central claim is that the compact-binary merger rate density of the local universe receives a substantial contribution from open clusters, with the WD-WD channel dominating. Using a large suite of N-body cluster simulations plus matched isolated-binary control runs, the authors report WD-WD merger efficiencies of $2.5\\times10^{-4}$ to $7.4\\times10^{-4}$ per solar mass, which translate to rate densities of 46-460, 130-1300, and 140-1400 Gpc$^{-3}$ yr$^{-1}$ for $10^2$, $10^3$, and $10^4$ solar-mass clusters, respectively. These exceed previous globular-cluster estimates. Super-Chandrasekhar WD-WD mergers, whose total mass exceeds the Chandrasekhar limit, occur at 70-780 Gpc$^{-3}$ yr$^{-1}$, above the roughly $10$ Gpc$^{-3}$ yr$^{-1}$ estimated for globular clusters, and are proposed as magnetar and fast-radio-burst progenitors; carbon-oxygen WD-WD mergers that may produce Type Ia supernovae account for only 0.14%-2.6% of the observed local Type Ia rate. The paper also establishes mass-dependent behavior: BH-BH mergers are dynamically formed, some with high eccentricity, only in $10^4$ solar-mass clusters, while BH-NS mergers appear only in $10^3$ solar-mass clusters, with a local rate of 2.3-23 Gpc$^{-3}$ yr$^{-1}$ compatible with gravitational-wave detector limits.","pith_inferences":["Inference: because roughly 70% of WD-WD mergers in the $10^3$ solar-mass models occur outside the cluster, searches for repeating fast radio burst counterparts should also look in old tidal streams and field populations descended from dissolved open clusters, not only inside intact clusters.","Inference: if real open clusters have a primordial binary fraction near 50% rather than 100%, the reported rate ranges would shrink by about a factor of two but would likely remain within an order of magnitude of globular-cluster estimates, so the qualitative conclusion would survive.","Inference: the mass-dependence map found here suggests each merger channel has a preferred host-cluster mass, which could be tested by correlating gravitational-wave event host environments or fast radio burst environments with estimated cluster masses.","Inference: a direct check of the fast radio burst hypothesis would compare the predicted super-Chandrasekhar WD-WD rate of 70-780 Gpc$^{-3}$ yr$^{-1}$ with the volumetric rate of repeating fast radio bursts in old stellar environments; order-of-magnitude agreement would support the magnetar-from-merger channel."],"forward_implications":["Open clusters with masses near $10^3$ solar masses produce WD-WD merger rate densities of 130-1300 Gpc$^{-3}$ yr$^{-1}$, larger than globular-cluster estimates, so any census of white-dwarf mergers must include open clusters.","Super-Chandrasekhar WD-WD mergers in open clusters occur at 70-780 Gpc$^{-3}$ yr$^{-1}$, making old open clusters and their dissolved remnants candidate birth sites for magnetars and fast radio bursts.","Carbon-oxygen WD-WD mergers from open clusters supply only 0.14%-2.6% of the observed Type Ia supernova rate, so open clusters are a minor but not negligible Type Ia channel.","Dynamically formed and eccentric BH-BH mergers appear only in the most massive $10^4$ solar-mass open clusters, placing the active-dynamics boundary between $10^3$ and $10^4$ solar masses.","BH-NS mergers occur only in $10^3$ solar-mass clusters at 2.3-23 Gpc$^{-3}$ yr$^{-1}$, a rate within the range measured by gravitational-wave observatories.","Most WD-WD mergers from low- and intermediate-mass clusters happen after the host cluster is tidally disrupted, so the cluster's dynamical influence outlives the cluster itself."],"supporting_citations":[{"why":"Supplies the N-body simulation code used to evolve the open clusters with regularized treatment of close binaries and encounters.","marker":"L. Wang et al. 2020a"},{"why":"Supplies the single and binary stellar evolution prescriptions that turn stars into white dwarfs, neutron stars, and black holes.","marker":"J. R. Hurley et al. 2002"},{"why":"Provides the binary evolution model, including common envelope evolution and natal kick treatment, used in the simulations.","marker":"S. Banerjee et al. 2020"},{"why":"Provides the initial orbital-period, mass-ratio, and eccentricity distributions for binaries with primary masses above five solar masses.","marker":"H. Sana et al. 2012"},{"why":"Provides the initial binary distributions for lower-mass primaries and the initial mass function adopted for the cluster populations.","marker":"P. Kroupa 1995a, 1995b"},{"why":"Supplies the cosmic star formation rate density as a function of redshift used to convert merger efficiencies into volumetric rates.","marker":"P. Madau & T. Fragos 2017"},{"why":"Supplies the abundance, mass function, and cluster formation efficiency of open clusters that set the assumed $f_{\\rm cluster}$ range and the mass-dependent local rate.","marker":"S. F. Portegies Zwart et al. 2010"},{"why":"Provides the globular-cluster super-Chandrasekhar WD-WD merger rate, roughly 10 Gpc$^{-3}$ yr$^{-1}$, that the paper's open-cluster rates are compared against and exceed.","marker":"K. Kremer et al. 2021a"},{"why":"Provides prior N-body open-cluster BH-BH merger efficiencies without primordial binaries, which the authors use to attribute their roughly tenfold higher efficiency to the 100% binary fraction.","marker":"J. Kumamoto et al. 2019"},{"why":"Provides the observed local Type Ia supernova rate used to compute the 0.14%-2.6% contribution from open-cluster carbon-oxygen WD-WD mergers.","marker":"W. Li et al. 2011"}],"fun_headline_variants":["Open clusters outpace globulars in white-dwarf mergers","White-dwarf merger rate from open clusters tops globulars","Open clusters may power fast radio bursts via white-dwarf mergers","Open clusters exceed globulars in white-dwarf merger density","Open clusters: white-dwarf mergers up to 1400 per Gpc^3/yr"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The rate numbers scale linearly with the assumed primordial binary fraction, and the paper sets that fraction to 100% for every cluster; if real open clusters have lower binary fractions, all quoted merger rates are correspondingly too high.","fun_headline_variants_meta":{"raw":{"variants":["Open clusters outpace globulars in white-dwarf mergers","White-dwarf merger rate from open clusters tops globulars","Open clusters may power fast radio bursts via white-dwarf mergers","Open clusters exceed globulars in white-dwarf merger density","Open clusters: white-dwarf mergers up to 1400 per Gpc^3/yr"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.002169,"raw_usage":{"total_tokens":8578,"prompt_tokens":1290,"completion_tokens":7288,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":906,"completion_tokens_details":{"reasoning_tokens":7206}},"tokens_in":906,"tokens_out":7288,"duration_ms":63326,"temperature":1.0,"reasoning_tokens":7206,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T22:00:52.011329+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the primordial binary fraction in young open clusters with masses near $10^2$ to $10^4$ solar masses, for example through spectroscopic binary surveys of embedded clusters; if the fraction is substantially below 100%, the reported WD-WD, WD-NS, and BH-NS rates should be scaled down proportionally, and a fraction near 50% would pull the super-Chandrasekhar WD-WD rate toward the globular-cluster comparison level.","supporting_citations":[{"cited_title":"2017, ApJ, 840, 39","cited_arxiv_id":null,"evidence_quote":"Supplies the cosmic star formation rate density as a function of redshift used to convert merger efficiencies into volumetric rates."},{"cited_title":"S., & Tanikawa, A","cited_arxiv_id":null,"evidence_quote":"Provides prior N-body open-cluster BH-BH merger efficiencies without primordial binaries, which the authors use to attribute their roughly tenfold higher efficiency to the 100% binary fraction."}],"review_version":1}