Stellar initial mass function in the 100-pc solar neighbourhood

Chao Liu; Jiadong Li; Yu-Ting Wang

arxiv: 2506.12987 · v2 · submitted 2025-06-15 · 🌌 astro-ph.GA · astro-ph.SR

Stellar initial mass function in the 100-pc solar neighbourhood

Yu-Ting Wang , Chao Liu , Jiadong Li This is my paper

Pith reviewed 2026-05-19 09:12 UTC · model grok-4.3

classification 🌌 astro-ph.GA astro-ph.SR

keywords stellar initial mass functionGaia DR3solar neighbourhoodunresolved binariescolor-magnitude diagrammass function slopesbinary fraction

0 comments

The pith

Gaia DR3 data determine the stellar initial mass function near the Sun with a break at 0.4 solar masses and slopes 0.75 and 2.07.

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

The paper fits a parametrized initial mass function to the color-magnitude diagram of stars inside 100 parsecs by forward-modeling Gaia DR3 photometry and astrometry. It corrects for Malmquist and Lutz-Kelker biases, metallicity-dependent mass-luminosity relations, and unresolved binaries whose components are drawn from a single IMF after a dynamical-evolution mimic that reproduces the observed mass-ratio distribution. The resulting IMF over 0.25 to 1 solar mass has a low-mass slope of 0.75, a high-mass slope of 2.07, and a break mass of 0.40 solar masses, all with uncertainties of a few percent. A reader cares because this distribution sets the number of low-mass stars, brown dwarfs, and the total stellar mass budget that shapes the local galactic disk and the rate of planetary and compact-object formation.

Core claim

The authors determine an averaged stellar IMF over 0.25 to 1.0 solar masses with power-law indices α1 = 0.75 and α2 = 2.07 separated by a break at 0.40 solar masses. The same fit yields a mean binary fraction of 26 percent and an angular-resolution constraint on Gaia DR3 of 1.11 arcseconds, together with a ξ-parameter value of 0.5070 for direct comparison with other IMF determinations.

What carries the argument

Forward modeling of the observed color-magnitude diagram that synthesizes unresolved binaries by drawing all components from the same IMF after a dynamical-evolution process that recovers the measured present-day mass-ratio distribution.

If this is right

The IMF parameters align with canonical forms but carry substantially smaller uncertainties than previous local determinations.
The average binary fraction over the fitted mass range is 26 percent.
The Gaia DR3 angular resolution is constrained to 1.11 arcseconds.
The ξ-parameter of the IMF is fixed at 0.5070 for comparison with other studies.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

If the derived break mass and slopes apply to other low-density environments, the local stellar mass budget would be lower than in denser star-forming regions.
Future astrometric missions could repeat the same forward-modeling exercise on a larger volume to test whether the 0.4 solar-mass break is universal.
Refinements to the binary synthesis step would allow the method to be applied to open clusters where dynamical processing is stronger.

Load-bearing premise

The binary population is synthesized by imitating dynamical evolution in star clusters so that all components are drawn from the same IMF while recovering the observed present-day mass-ratio distribution.

What would settle it

A volume-limited census of resolved binary mass ratios in the same 100-pc volume that deviates systematically from the distribution recovered by the dynamical-evolution mimic would falsify the derived IMF slopes and break mass.

read the original abstract

The stellar initial mass function (IMF) is among the most fundamental distributions in astrophysics, defined as the mass spectrum of stars produced in a single star-formation event. Even in the solar neighbourhood, where measurements can be conducted via star counting, disentangling the IMF from observational effects remains challenging. In this work we introduce a new parametrisation of the stellar IMF in the 100-pc solar neighbourhood, leveraging the high-precision astrometric and photometric data from \textsl{Gaia} DR3: we model the colour-magnitude diagram of the field star population while accounting for observational uncertainties, Malmquist bias, Lutz-Kelker bias, variations in the mass-luminosity relation arising from metallicity differences, and the effects of unresolved binaries. In particular, we synthesise the binary population with a process imitating the dynamical evolution observed in star clusters to enforce that all components are drawn from the same IMF, while simultaneously recovering the observed present-day mass-ratio distribution. We determine an averaged stellar IMF over $0.25<m<1.0~M_{\odot}$ that aligns with canonical IMFs but achieves significantly tighter constraints: $\alpha_1=0.75^{+0.06}_{-0.04}$, $\alpha_2=2.07^{+0.04}_{-0.03}$, and a break point at $m_{\mathrm{break}}=0.40^{+0.01}_{-0.01}$ $\mathrm{M_{\odot}}$. Our inference also yields an averaged binary fraction over $0.25<m<1.0~M_{\odot}$ of approximately 26\%, and constrains the \textsl{Gaia} DR3 angular resolution to $1.11^{+0.11}_{-0.08}$ arcsec. We also provide the $\xi$-parameter for our IMF, which is $0.5070_{-0.0096}^{+0.0068}$, to facilitate direct comparison with other IMF determinations.

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note private letter to a colleague

Gaia DR3 forward modeling yields tighter local IMF slopes and break mass than prior work, but the cluster-dynamical binary synthesis risks coupling into those constraints for a mixed field sample.

read the letter

The paper fits a broken power-law IMF to the Gaia DR3 color-magnitude diagram for stars within 100 pc, jointly modeling unresolved binaries, metallicity effects on the mass-luminosity relation, and the main selection biases. They report α1 = 0.75 with asymmetric errors of +0.06/-0.04, α2 = 2.07 +0.04/-0.03, and a break at 0.40 solar masses, plus a binary fraction near 26% and a Gaia resolution constraint around 1.11 arcsec. The numbers line up with canonical forms but come with smaller uncertainties than most earlier local determinations. The forward-modeling setup that includes Malmquist, Lutz-Kelker, and photometric scatter is a clear step forward from simple star-counting approaches. The binary component is constructed to draw all stars from the same IMF while reproducing the observed mass-ratio distribution through an imitation of cluster dynamical evolution. That choice lets them separate binary effects from the underlying IMF shape in a single fit. The main soft spot is exactly there. The solar-neighborhood population is a mix of stars born in different clusters and processed over Gyr timescales, so an imitation tuned to cluster conditions may not match the actual present-day binary statistics in the field. Any mismatch would shift the inferred slopes and break mass, and the reported precision could be optimistic. The abstract gives no validation plots or sensitivity tests on that assumption, which leaves the central result harder to assess. This work is aimed at researchers who need updated local IMF parameters for star-formation theory or chemical-evolution calculations. Readers who already use Gaia CMDs or binary population synthesis will find the joint constraints useful even if they adjust the binary prescription later. The approach is coherent enough and the data quality high enough that it deserves a serious referee rather than a desk reject. I would send it out for review with a request for explicit checks on how the binary model affects the IMF posteriors.

Referee Report

2 major / 2 minor

Summary. The paper presents a forward-modeling analysis of the stellar initial mass function (IMF) in the 100-pc solar neighborhood using Gaia DR3 astrometric and photometric data. The authors model the observed color-magnitude diagram while incorporating Malmquist bias, Lutz-Kelker bias, metallicity-dependent mass-luminosity variations, and unresolved binaries. Binaries are synthesized via a process that imitates dynamical evolution in star clusters, ensuring all components are drawn from the same IMF while reproducing the observed present-day mass-ratio distribution. They report an averaged IMF over 0.25 < m < 1.0 M⊙ with parameters α1 = 0.75^{+0.06}_{-0.04}, α2 = 2.07^{+0.04}_{-0.03}, m_break = 0.40^{+0.01}_{-0.01} M⊙, an average binary fraction of ~26%, a Gaia DR3 resolution of 1.11^{+0.11}_{-0.08} arcsec, and ξ = 0.5070^{+0.0068}_{-0.0096}.

Significance. If the binary synthesis accurately captures field-star statistics, the work delivers substantially tighter constraints on the local IMF than most prior determinations, along with a directly comparable ξ-parameter. The forward-modeling framework that simultaneously fits multiple observational effects is a methodological strength.

major comments (2)

[Binary modeling (methods)] Binary synthesis section: the assumption that a cluster-dynamical imitation process (which draws all components from the same IMF and recovers the observed mass-ratio distribution) fully represents the mixed field population in the 100-pc sample is load-bearing for the reported precision on α1, α2, and m_break. Stars in the solar neighborhood originate from diverse birth environments with potentially different initial binary fractions and Gyr-scale processing; any mismatch would couple directly into the fitted slopes and break mass. An explicit test using alternative field-appropriate binary prescriptions or a sensitivity analysis would be required to support the claimed uncertainties.
[Results and fitting procedure] Fitting and validation: the abstract states that the model accounts for Malmquist and Lutz-Kelker biases and yields tight posteriors, yet the manuscript provides insufficient detail on MCMC convergence diagnostics, prior sensitivity tests, or checks against post-hoc sample cuts. Without these, it is not possible to confirm that the reported 1σ uncertainties (e.g., ±0.01 M⊙ on m_break) are not underestimated.

minor comments (2)

[Abstract] Notation for asymmetric uncertainties is inconsistent between the abstract and the ξ-parameter reporting; standardize the superscript/subscript order.
[Discussion] Add a brief comparison table of the derived IMF parameters against at least three recent solar-neighborhood determinations to quantify the claimed improvement in precision.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive comments and positive evaluation of the work's significance. We address each major comment below and will revise the manuscript to incorporate additional tests and details as outlined.

read point-by-point responses

Referee: [Binary modeling (methods)] Binary synthesis section: the assumption that a cluster-dynamical imitation process (which draws all components from the same IMF and recovers the observed mass-ratio distribution) fully represents the mixed field population in the 100-pc sample is load-bearing for the reported precision on α1, α2, and m_break. Stars in the solar neighborhood originate from diverse birth environments with potentially different initial binary fractions and Gyr-scale processing; any mismatch would couple directly into the fitted slopes and break mass. An explicit test using alternative field-appropriate binary prescriptions or a sensitivity analysis would be required to support the claimed uncertainties.

Authors: We agree that the binary modeling choice is important for the precision of the reported IMF parameters. Our synthesis method was selected specifically to ensure all binary components are drawn from a single IMF while matching the observed present-day mass-ratio distribution in the field. Although the solar neighborhood includes stars from varied birth sites, the local volume is dominated by populations with comparable dynamical histories to those in clusters. To directly address the referee's concern, we will add an explicit sensitivity analysis in the revised manuscript, testing variations in initial binary fractions and alternative mass-ratio distributions to quantify their impact on the fitted α1, α2, and m_break. revision: yes
Referee: [Results and fitting procedure] Fitting and validation: the abstract states that the model accounts for Malmquist and Lutz-Kelker biases and yields tight posteriors, yet the manuscript provides insufficient detail on MCMC convergence diagnostics, prior sensitivity tests, or checks against post-hoc sample cuts. Without these, it is not possible to confirm that the reported 1σ uncertainties (e.g., ±0.01 M⊙ on m_break) are not underestimated.

Authors: We acknowledge that the current manuscript lacks sufficient detail on the statistical validation of the fitting procedure. In the revised version, we will expand the relevant section to include MCMC convergence diagnostics (e.g., trace plots, autocorrelation times, and Gelman-Rubin statistics), results from prior sensitivity tests, and robustness checks using alternative post-hoc sample selections. These additions will demonstrate that the reported uncertainties, including the tight constraint on m_break, are not underestimated. revision: yes

Circularity Check

0 steps flagged

No significant circularity; IMF parameters are direct outputs of fit to Gaia data

full rationale

The paper's derivation consists of a statistical inference of IMF parameters (α1, α2, m_break) by forward-modeling the observed color-magnitude diagram of the 100-pc solar-neighborhood field population in Gaia DR3. The model incorporates external inputs for observational biases, metallicity-dependent mass-luminosity relations, and a binary synthesis prescription that imitates cluster dynamical evolution to reproduce the present-day mass-ratio distribution. These modeling choices are stated as inputs to the fit rather than quantities derived from the target IMF parameters. No equations or steps reduce the reported slopes or break mass to the inputs by construction, and no load-bearing self-citations or uniqueness theorems are invoked in the abstract or described method. The result is therefore a standard parameter estimation against external data and remains self-contained.

Axiom & Free-Parameter Ledger

5 free parameters · 2 axioms · 0 invented entities

The central claim rests on the assumed broken-power-law form, the cluster-dynamical binary model, and the adopted mass-luminosity-metallicity relation; these are not derived inside the paper.

free parameters (5)

α1 = 0.75
Low-mass slope of the IMF, fitted to the data
α2 = 2.07
High-mass slope of the IMF, fitted to the data
m_break = 0.40
Mass at which the power-law index changes, fitted to the data
binary_fraction = 0.26
Average fraction of stars in binaries, fitted to the data
gaia_resolution = 1.11
Effective angular resolution of Gaia DR3, fitted to the data

axioms (2)

domain assumption The stellar IMF follows a broken power-law form between 0.25 and 1.0 solar masses
Stated as the model being fitted; not derived from first principles in the abstract
domain assumption Binary components are drawn from the same IMF and their present-day mass-ratio distribution can be recovered by imitating cluster dynamical evolution
Explicitly invoked in the binary-synthesis step

pith-pipeline@v0.9.0 · 5891 in / 1627 out tokens · 40922 ms · 2026-05-19T09:12:49.644930+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

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

IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

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

Our best estimation for the averaged field star IMF is α1 = 0.81+0.06−0.05, α2 = 2.12+0.04−0.04 and mbreak = 0.41+0.01−0.01

What do these tags mean?

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