A method for constructing the joint mass function of binary stars

Amery Gration; Payel Das; Robert G. Izzard

REVIEW 2 major objections 4 minor 3 references

A method for constructing the joint mass function of binary stars

T0 review · 2 major / 4 minor · reviewed 2026-08-03 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read This paper derives the joint mass function of binary stars for any pairing rule by solving a Fredholm integral equation for the primary-mass function.

desk verdict A clean, useful method for building the joint mass function of binaries, with a reproducible-computation problem in the one new numerical result. read the letter →

arxiv 2602.10186 v2 pith:ZZVNGS7D submitted 2026-02-10 astro-ph.SR astro-ph.GA

Amery Gration , Robert G. Izzard , Payel Das This is my paper

classification astro-ph.SRastro-ph.GA

keywords binarystarsinitialmassfunctionjointFredholmintegralequationprimaryratiodistributionrandompairinguniform

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper presents a general method to construct the joint distribution of primary and secondary stellar masses in a binary population, given only the stellar mass function and the conditional distribution of secondary mass given primary mass. The key is an integral equation whose solution is the primary-mass function; the joint distribution then follows as a product. The authors recover the known random-pairing result and, for the first time, derive the joint distribution for uniform pairing, where the mass-ratio distribution is flat. They also show that not every combination of mass function and pairing rule is self-consistent; a valid joint distribution exists only if the integral equation has a non-negative solution. This gives a practical way to build and test models of binary-star populations.

What carries the argument

The central object is an inhomogeneous Fredholm integral equation of the second kind (Equation 6), whose kernel is the conditional secondary-mass distribution and whose known term is twice the stellar mass function. It follows directly from the mixture identity f_M = (f_M1 + f_M2)/2, which states that in a pure binary population each star is equally likely to be the primary or secondary, combined with the marginal relation f_M2(m2) = ∫ f_{M2|M1}(m2|m1) f_M1(m1) dm1. Solving for the primary-mass function f_M1 yields the joint mass function as a product, and the equation's invertibility provides a self-consistency test for any proposed mass function and pairing rule.

What would settle it

Simulate a binary population with a known pairing rule and known IMF, measure the mass function and conditional secondary-mass distribution, solve Equation 6 numerically, and compare the reconstructed joint mass function with the true one; a mismatch beyond numerical error would falsify the claim that the equation always yields the correct JMF.

Watch

Extended reading notes

Core claim

The primary-mass function of a binary population must satisfy an inhomogeneous Fredholm integral equation of the second kind: f_M1(m) = 2 f_M(m) − ∫ f_{M2|M1}(m|m1) f_M1(m1) dm1. Once solved, the joint mass function is f(M1,M2) = f_{M2|M1}(m2|m1) f_M1(m1). The authors derive this equation from the mixture identity f_M = (f_M1 + f_M2)/2 and the marginal relation for the secondary-mass function. They illustrate the method by recovering the known random-pairing joint mass function and by computing the previously unknown uniform-pairing joint mass function, using the canonical Kroupa IMF and a minimum mass ratio qmin=0.1. They also prove limiting behaviours: in main-sequence binaries the seconda

Load-bearing premise

The input mass function must describe the masses of all stars in a pure binary population, with every star equally likely to be the primary or secondary; if single stars are mixed in or the binary fraction is below one, the mixture identity f_M = (f_M1 + f_M2)/2 fails and the construction no longer holds.

Editorial extensions

If this is right

The joint mass function can now be constructed for any pairing scheme, not just random pairing, provided the mass function and conditional secondary-mass distribution are specified.
The uniform-pairing joint mass function is derived for the first time, showing a primary-mass peak at 0.11 solar masses (versus 0.16 for random pairing) and a corresponding shift toward more low-mass primaries and fewer high-mass primaries.
The method doubles as an existence test: if the integral equation yields a negative solution, the chosen mass function and conditional secondary-mass distribution cannot both describe a real binary population.
The limiting behaviour at the mass extremes is fixed for main-sequence binaries: no secondary can be less massive than the least-massive primary, and no primary can be more massive than the most-massive secondary, with the relevant component function equal to twice the mass function at each limit.
Population-synthesis codes can replace approximations (such as primary-constrained pairing) with the exact joint distribution for any chosen pairing kernel.

Reading between the lines

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

The inverse problem becomes plausible: given observed primary and secondary mass functions, one could solve Equation 6 for the pairing kernel, thereby empirically constraining the binary star formation mechanism.
For populations with a binary fraction below 100%, the mixture identity f_M = (f_M1 + f_M2)/2 must be generalised to a weighted mixture including single stars; the method could be extended to estimate the binary fraction alongside the joint distribution.
The Fredholm formulation may extend to higher-order multiples, such as triples, by writing a similar integral equation for the distribution of the most-massive component in an n-tuple system.
A data-driven test would compare the method's predictions against observed mass-ratio distributions: if a measured CMRF from a survey fails to yield a non-negative PMF, that pair of empirical inputs is internally inconsistent.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Dear [name],

The short version: this paper gives a general method for constructing the joint mass function of binaries — solve a Fredholm equation for the primary mass function — and it actually works. The random-pairing case is recovered cleanly, and the uniform-pairing result is new. Worth engaging with. But the appendix has an algebraic error that needs fixing before the new result is reproducible from the paper.

The derivation of Eq 6 is straightforward and the authors write it up well: Eq 4 (the MF is a fifty-fifty mixture of the primary and secondary mass functions) plus marginalization of the joint distribution gives an inhomogeneous Fredholm equation of the second kind. They are candid that non-negativity is the only real consistency condition. The random-pairing check is a nice touch — it reproduces Malkov & Zinnecker's PMF. The uniform-pairing solution, for a Kroupa-like IMF with q_min=0.1, is the first exact construction of that JMF I know of. The discussion of limiting behaviour at the low- and high-mass ends is careful and correct.

The soft spot is Appendix A, and the stress-test note is right about it. Eq 35 is correct: \hat f = 2 f_M - A \hat f. But the matrix form drops the factor of 2. The paper then solves (I+A)\hat f = f_M, which gives exactly half the right PMF. Because the equation is linear, the shape is right but the normalization is off by a factor of two. This is a trivial fix — replace f_M with 2 f_M on the RHS of Eqs 36 and 38 — but as printed the algorithm does not produce a properly normalized PDF. There is also no code or error estimate, so the figures can't be reproduced without guesswork. That's the main thing I'd ask the authors to address.

The other limitation, which the reader flagged, is that Eq 4 assumes a pure binary population. The paper states this explicitly in Section 2, but the implications for applying the method to observed IMFs are not explored. If the observed IMF is for a mixture of single and binary stars, or if the binary fraction is not 100%, Eq 4 doesn't hold. That's a genuine constraint, though it's also a scope choice rather than a flaw in the math.

Bottom line: the central idea is sound and useful for binary population synthesis and IMF inference. The appendix error is real but easily fixed. I'd send it to a serious referee. If the authors correct the normalization and post the code, I'll cite it.

Referee Report

2 major / 4 minor

Summary. The paper presents a method for constructing the joint mass function (JMF) of binary stars from a given mass function (MF) and a given conditional secondary-mass distribution (CSMF). The central result is a Fredholm integral equation of the second kind for the primary-mass function (PMF), Eq. (6), derived from the mixture identity f_M = (f_{M1}+f_{M2})/2. The JMF then follows as f_{M2|M1} f_{M1}. The method is illustrated in two cases: random pairing, where the known result is recovered, and uniform pairing, where a novel PMF/JMF is computed numerically. The authors also derive limiting behaviours for main-sequence binaries and discuss the sensitivity of the results to the mass-function support.

Significance. If the numerical issue is corrected, the paper offers a clean, general framework for constructing binary-star JMFs beyond random pairing. The integral-equation formulation is elegant and the random-pairing recovery is a valuable self-consistency check. The method's ability to handle arbitrary pairing schemes, including the previously untreated uniform-pairing case, is a genuine contribution. However, the paper's central new numerical result (uniform pairing, Figs 3–4) is currently unreproducible because of an inconsistency in the discretized equations in Appendix A. The discussion of the method's domain of validity (pure binary populations) also needs to be made explicit.

major comments (2)

[Appendix A, Eqs (35)–(38)] There is an internal inconsistency in the discretized Fredholm equation. Eq. (35) correctly retains the factor 2 on the right-hand side, \(\hat f_{M1}(m_j)=2f_M(m_j)-\sum_i w_i f_{M2|M1}(m_j|m_i)\hat f_{M1}(m_i)\), but Eq. (36) drops this factor, writing \(\hat f_{M1}=f_M-A\hat f_{M1}\), and Eqs (37)–(38) solve \((I+A)\hat f_{M1}=f_M\). Since the kernel integrates to 1, the printed scheme yields a PMF that integrates to 1/2, which is not a probability density. The figures therefore cannot have been produced by the algorithm as stated. The fix is to replace \(f_M\) with \(2f_M\) in Eqs (36) and (38). In addition, the quadrature rule, node count, and convergence tests are not reported, so the numerical result is not reproducible even after correcting the factor. Please provide these details, or release code.
[Section 2.1, Eq. (4)] The mixture identity \(f_M=(f_{M1}+f_{M2})/2\) is load-bearing and holds only for a pure binary population with a 100% binary fraction, where the MF is the distribution of the union of primary and secondary masses. If the supplied \(f_M\) is the standard field IMF, which includes single stars, or if the binary fraction is less than unity, Eq. (4) fails and the integral equation no longer constructs the JMF of the full population. The paper uses the canonical Kroupa IMF in both examples without explicitly stating this condition at the point of derivation. Section 4.1 qualitatively acknowledges that the IMF of binaries may differ from the canonical IMF, but the method's scope should be stated precisely in Section 2.1 so that readers do not apply it to an IMF that includes single stars.

minor comments (4)

[Section 2.2, first paragraph] The sentence 'In general the mass of the secondary star can be greater than the mass of the primary star' contradicts the definition of primary as the more massive component used later in Section 3.1 ('the more-massive star designated the primary'). Please revise for consistency; for main-sequence binaries the secondary is lighter by construction.
[Appendix A, Eq. (36) and matrix definitions] The notation around the matrix A is inconsistent: the indices in Eq. (35) run over i=0,...,n, whereas the vectors in Eq. (36) are indexed from 1 to n. Also, A should be defined such that its (j,i) entry is \(w_i f_{M2|M1}(m_j|m_i)\), not \(f_{M2|M1}(m_i,m_j)\) as written. Please make the indexing explicit.
[Section 3.2 and Appendix A] Since the uniform-pairing solution is numerical and new, the authors should state explicitly that the computed PMF is non-negative (as required by Eq. 6 to be a valid PDF) and provide some measure of numerical uncertainty, e.g., convergence with the number of quadrature nodes. A brief note in the figure captions or appendix would suffice.
[General presentation] The paper would benefit from a short summary of the conditions under which Eq. (6) admits a non-negative solution. The Fredholm theory guarantees a unique solution for many kernels, but the non-negativity constraint is problem-specific and is only checked after the fact. A short discussion or example of an inconsistent (MF, CSMF) pair would clarify the method's limitations.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the integral-equation construction is self-contained; the few misprints are correctness/reproducibility issues, not circular reductions.

full rationale

This is a self-contained mathematical construction. Equation 6 is obtained algebraically from the mixture identity f_M=(f_M1+f_M2)/2 (Eq 4) and the factorization f(M1,M2)=f_{M2|M1} f_{M1} (Eq 3); these are stated assumptions about the binary population, not results derived from the target PMF. Solving the Fredholm equation for f_M1 then determines the JMF without adjusting parameters. The random-pairing section deliberately starts from the known JMF (Eq 18) to form the CSMF, so recovering Equation 22 is a consistency check against an external benchmark, not a prediction forced by a fitted value. The uniform-pairing section chooses the CSMF explicitly (Eqs 29-31) with q_min=0.1 as an illustrative input and computes the PMF numerically; no quantity is fitted to the plotted output. Self-citations (Izzard & Halabi 2019; Izzard 2023) appear only in application/population-synthesis remarks and are not load-bearing. Two manuscript defects are correctness issues, not circularity: the displayed Eq 7 states the integral of the conditional density times f_M1 over the secondary mass as 1 rather than f_M1(m1) (the normalization conclusion still follows after correcting this), and in Appendix A the factor 2 is dropped between Eq 35 and Eq 36, so the printed matrix system (I+A) f_hat = f_M is not a faithful discretization of Eq 6. These affect reproducibility of the uniform-pairing figures but do not make the derivation equivalent to its inputs.

Assumptions & free parameters 1 free parameters · 5 assumptions · 0 invented entities

The method has essentially no fitted parameters beyond the illustrative q_min choice; the main axioms are domain assumptions about how the mass function relates to primary/secondary distributions. No new physical entities are introduced.

free parameters (1)

q_min = 0.1
Lower bound on the mass ratio in the uniform-pairing example; chosen by hand to illustrate the method, not fitted to data. The resulting PMF/JMF depend on this choice.

assumptions (5)

domain assumption The mass function f_M of stars in the binary population is the average of the primary and secondary mass functions: f_M = (f_M1 + f_M2)/2 (Eq 4).
This treats every star in a binary as equally likely to be drawn from the primary or secondary population. It is the central structural assumption that converts the problem into a Fredholm equation; if the input f_M is measured from single stars only, Eq 4 need not hold.
domain assumption A joint mass function exists for the chosen f_M and CSMF; the paper notes this may fail (non-negative solution required).
Eq 6 may have no PDF solution for inconsistent inputs; the method assumes one exists for the examples considered.
domain assumption For main-sequence binaries, secondary mass never exceeds primary mass (m2 <= m1) in Section 2.2.
Used to derive boundary behaviors of PMF and SMF. Not required for the general method but used for the illustrative limits.
domain assumption The canonical IMF of Kroupa (2001) with slopes -1.3 and -2.3 and break at 0.5 Msun is used as f_M.
An empirical input from the literature; the results are specific to this choice.
standard math The Nystrom method with product integration converges for this problem (Appendix A).
Background on integral equation numerics; no convergence proof supplied for the specific kernel.

how reviews work

0 comments

Cite this review

Pith. "Pith review of A method for constructing the joint mass function of binary stars." pith.science (2026). https://pith.science/paper/ZZVNGS7D

@misc{pith2026260210186,
  author       = {Pith},
  title        = {Pith review of: A method for constructing the joint mass function of binary stars},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ZZVNGS7D}},
  note         = {Machine review of arXiv:2602.10186}
}

read the original abstract

The initial mass function (IMF) describes the distribution of stellar masses in a population of newly born stars and is amongst the most fundamental concepts in astrophysics. It is not only the direct result of the star formation process but it also explains the evolution of galaxies' luminosities, metal yields, star-formation efficiencies, and supernova production rates. Because most stars exist in binary systems, however, a full statistical account of stellar mass requires not the IMF but rather the joint distribution of a binary population's primary- and secondary-star masses. This joint distribution must respect the IMF of the stars from which the population has been assembled as well as the distribution of mass ratios that results from the assembly mechanism. Despite its importance, this joint distribution is known only in the case of random pairing. Here we present a method for constructing it in the general case. We also illustrate the use of our method by recovering the known result for random pairing and by finding the previously unknown result for uniform pairing.

Figures

Figures reproduced from arXiv: 2602.10186 by the authors.

**Figure 1.** The primary-mass function, fM1 (Equation 22), and secondary-mass function, fM2 (Equation 27), alongside the mass function, fM (Equation 16). Here we have used the IMF of Kroupa (2001) and random pairing (Section 3.1) [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗

**Figure 4.** The joint mass function, f(M1,M2) , constructed using the IMF of Kroupa (2001) and uniform pairing (Section 3.2, Figure 3). It is non-zero for m2 ≤ m1 and m2 ≥ qminm1 where qmin = 0.1. As is the case for random pairing ( [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

3 extracted references · 1 linked inside Pith

[1]

E., 1997, The Numerical Solution of Integral Equa- tions of the Second Kind

Atkinson K. E., 1997, The Numerical Solution of Integral Equa- tions of the Second Kind. Cambridge Monographs on Applied and Computational Mathematics, Cambridge University Press Bate M. R., 2015, in Rucinski S. M., Torres G., Zejda M., eds, Astronomical Society of the Paciﬁc Conference Series Vol. 496, Living Together: Planets, Host Stars and Binaries. p...

1997
[31]

This equation does not have closed- form solution but we can instead use a numerical method, speciﬁcally the Nyström method ( Hackbusch 1995)

In this case, fM1 (m) = 2 fM (m) − ∫ mmax mmin 1B(m1)fM1 (m1) m1(1 − max(qmin, mmin/m1)) dm1 (33) where 1B is the indicator function of the interval B := [m, min(m/qmin, mmax)]. This equation does not have closed- form solution but we can instead use a numerical method, speciﬁcally the Nyström method ( Hackbusch 1995). When using the Nyström method we app...

1995
[120]

Birkhäuser Hennebelle P., Grudić M. Y., 2024, Annual Review of Astronomy and Astrophysics, 62, 63 Hypki A., Giersz M., 2013, Monthly Notices of the Royal Astro- nomical Society, 429, 1221 Izzard R., 2023, Astrophysics Source Code Library, p. ascl:2307.035 Izzard R. G., Halabi G. M., 2019, in Beccari G., Boﬃn H. M. J., eds, Cambridge Astrophysics, The Impa...

arXiv 2024

Pith tools

Reviewed August 3, 2026 · model on record in the stance chip above.

[1] [1]

E., 1997, The Numerical Solution of Integral Equa- tions of the Second Kind

Atkinson K. E., 1997, The Numerical Solution of Integral Equa- tions of the Second Kind. Cambridge Monographs on Applied and Computational Mathematics, Cambridge University Press Bate M. R., 2015, in Rucinski S. M., Torres G., Zejda M., eds, Astronomical Society of the Paciﬁc Conference Series Vol. 496, Living Together: Planets, Host Stars and Binaries. p...

1997

[2] [31]

This equation does not have closed- form solution but we can instead use a numerical method, speciﬁcally the Nyström method ( Hackbusch 1995)

In this case, fM1 (m) = 2 fM (m) − ∫ mmax mmin 1B(m1)fM1 (m1) m1(1 − max(qmin, mmin/m1)) dm1 (33) where 1B is the indicator function of the interval B := [m, min(m/qmin, mmax)]. This equation does not have closed- form solution but we can instead use a numerical method, speciﬁcally the Nyström method ( Hackbusch 1995). When using the Nyström method we app...

1995

[3] [120]

Birkhäuser Hennebelle P., Grudić M. Y., 2024, Annual Review of Astronomy and Astrophysics, 62, 63 Hypki A., Giersz M., 2013, Monthly Notices of the Royal Astro- nomical Society, 429, 1221 Izzard R., 2023, Astrophysics Source Code Library, p. ascl:2307.035 Izzard R. G., Halabi G. M., 2019, in Beccari G., Boﬃn H. M. J., eds, Cambridge Astrophysics, The Impa...

arXiv 2024