The initial conditions and initial mass functions of Alpha Persei, Pleiades and Praesepe

arXiv:2607.17300 · astro-ph.GA, astro-ph.SR · Submitted 2026-07-19 · Read on arXiv

Listen

Transcript

Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.

Vera: Next we'll be talking about the paper "The initial conditions and initial mass functions of Alpha Persei, Pleiades and Praesepe".

Jocelyn: The paper was written by L. Hobart, H. Baumgardt and S. Sweet from School of Mathematics and Physics, University of Queensland, St. Lucia, QLD 4072, Australia and University of Queensland.

Vera: Stay tuned as we take you through the paper and discuss its implications.

The Scope and Implications: Vera: We've just touched on why these clusters are important, but let's go deeper into the scope of "The initial conditions and initial mass functions of Alpha Persei, Pleiades and Praesepe" as presented in its title.

Jocelyn: The title suggests a focus not only on what we see now—the present-day distribution—but also on what happened at birth: the initial conditions.

Subrahmanyan: That dual focus is crucial because it allows them to model the evolutionary path, connecting our current observations back to the physical processes of star formation itself.

Vera: It’s a way of reversing time, essentially seeing how the initial mass function shaped what we see today in these three distinct locations.

Jocelyn: The implications for finding a universal pattern are huge; if we find that these specific local clusters have different starting points, it challenges our assumption that star formation is uniform across the galaxy.

Subrahmanyan: Exactly, and this paper seems to be providing strong evidence that the IMF—the initial mass function—is not universally constant in relation to environmental factors.

Vera: We’re seeing a lot of detailed data from UKIDSS and HIPPARCOS helping define these initial conditions, which gives us a very robust picture of the local sky.

Jocelyn: It’s a powerful demonstration that the specific geometry and density of star formation are influencing how many stars we expect to find in certain mass bins.

Subrahmanyan: The goal is to provide constraints on the universality of the IMF, and this paper offers substantial data to achieve that goal.

Vera: So, setting us up for a big test: whether these local findings represent a real environmental dependence or just a random fluctuation.

Key Findings and Results: Vera: Moving on to the summary of "The initial conditions and initial mass functions of Alpha Persei, Pleiades and Praesepe," what are the most striking results they present?

Jocelyn: The researchers found a very specific pattern in how stars are distributed that tells us a lot about their star-forming environment.

Subrahmanyan: They determined that the best-fitting initial mass function is described by a three-stage broken power law, which explains why we see different slopes at different parts of the stellar mass range.

Vera: That's fascinating because it suggests there isn't one single mathematical description for how stars are born, but rather a complex interplay of stages and break masses.

Jocelyn: And when they factor in the resolved binaries—those stars we can see as pairs—the measured binary fractions fall right around twenty to twenty-four percent for each cluster.

Subrahmanyan: That percentage is important because it means that a significant portion of the total stellar population is actually hidden, which impacts how we calculate the overall mass function.

Vera: It’s a constant reminder that our observations are always incomplete; we're seeing only one part of the stellar story.

Jocelyn: The finding also shows that while the low-mass slope is sensitive to our mass-luminosity conversion, it’s not the high-mass slope that varies much, with those measurements showing a dispersion of around zero point three sigma.

Subrahmanyan: This lack of variation in the high mass IMF is a key finding, as it suggests that while star formation at low masses might be unpredictable, the massive stars are consistently behaving across different environments.

Vera: So, when comparing these findings to our expectation of a Salpeter-like distribution, they found a tendency to have fewer massive stars than we traditionally expected.

Jocelyn: It’s like finding out that the "rulebook" for star birth has some specific exceptions for every single cluster.

Methodological Improvements: Vera: The technical rigor in "The initial conditions and initial mass functions of Alpha Persei, Pleiades and Praesepe" is really what sets this work apart; let’s talk about the improvements they suggest for future methods.

Jocelyn: Their approach goes far beyond just counting stars; they use complex tools to ensure that their data collection is thorough and reliable.

Subrahmanyan: They are essentially bridging the gap between observational astronomy and advanced computational physics by training machine learning emulators on N-body simulations.

Vera: That's a huge leap forward, combining those detailed simulations with AI to be able to predict what the initial conditions must have been like.

Jocelyn: It allows us to account for the effects of dynamical evolution, which is something we can’t do simply by looking at current photos.

Subrahmanyan: Without modeling how gravity and stellar interactions pull stars out of the cluster over time, our current mass function would be fundamentally flawed.

Vera: They are using a Bayesian framework to determine the most probable initial number of stars, which is a very smart way to handle all those uncertainties in any measurement.

Jocelyn: It’s not just one single best guess; it’s a range of possibilities that makes the results much more reliable for us observers.

Subrahmanyan: This approach allows them to find initial mass function parameters that are consistent with the observed present-day clusters, which is a major theoretical win.

Conclusion and Wrap-Up: Vera: We've spent quite a lot of time going through the data and the methods in "The initial conditions and initial mass functions of Alpha Persei, Pleiades and Praesepe."

Jocelyn: It feels like we've reached a very clear understanding that these clusters are offering us a nuanced look at how stars start their lives.

Subrahmanyan: The core conclusion is that the IMF isn't just one simple power law; it has complex features and may vary based on the environment where it formed.

Vera: And I really appreciate that they have provided a clear path forward for future work, especially by noting how better mass-luminosity relations will help us reduce systematic errors.

Jocelyn: It gives us observers a strong baseline to compare against other deep surveys and helps define what we expect from the local Milky Way.

Subrahmanyan: The authors also leave us with the important realization that our current sample might not be complete, suggesting there could be more variability in star formation than we've observed so far.

Vera: That sense of ongoing discovery is always exciting, knowing that even with such detailed work, there' more to find out.

Jocelyn: It’s a great foundation for the next generation of cluster studies across the galaxy.

Final thoughts on "The initial conditions and initial mass functions of Alpha Persei, Pleiades and Praesepe": Subrahmanyan: Before we wrap up, I think it's worth reiterating that this paper is a tremendous effort in connecting our current observations back to the fundamental physics of stellar birth.

Vera: It really highlights the difference between what we *expect* to see based on old models and what these precise measurements of "The initial conditions and initial mass functions of Alpha Persei, Pleiades and Praesepe" show us today.

Jocelyn: We have a much more informed picture now of how stellar populations are distributed, which is a huge benefit for our survey work.

Subrahmanyan: The authors have provided the tools to quantify not just the stars we see, but the entire process of how they came to be.

Vera: Thank you all for walking us through this incredibly detailed work, and I look forward to discussing what other discoveries await us in our next segment.

School of Mathematics and Physics, University of Queensland, St. Lucia, QLD 4072, Australia · University of Queensland

astro-ph.GA, astro-ph.SR

Submitted: 2026-07-19

Updated: 2026-07-19

Comments: Accepted for publication in PASA. 27 pages, 17 figures, 7 tables

Project page: https://lachlanhobart.github.io/Data/References

License: http://creativecommons.org/licenses/by/4.0/

Importance score: 67/100

The gist: The following is a detailed summary of the scientific paper, quoting relevant sections where necessary: The paper investigates "The initial conditions and initial mass functions of Alpha Persei,

Key concepts

Initial Mass Function (IMF): The IMF is a measure of the distribution of stellar masses. The paper suggests it is not universally constant, but rather has complex features and may vary depending on the environment where the stars formed.
Initial Conditions: This refers to modeling the physical processes of star formation itself—the starting state of a cluster before it evolved. It allows researchers to connect current observations back to how the initial mass function shaped what we see today.
Bayesian Framework: The authors used this framework to determine the most probable initial number of stars. It handles uncertainties by providing a range of possibilities, making the results more reliable than relying on a single best guess.

Terminology

alpha — 58× in this paper · explained in 1 episode
Initial Mass Function — 24× in this paper
Metallicity — 16× in this paper · explained in 3 episodes
Mass Segregation — 12× in this paper · explained in 2 episodes
N-body simulations — 9× in this paper · explained in 2 episodes
Bayesian Framework — 6× in this paper · explained in 3 episodes
Initial Mass Function (IMF) — 5× in this paper · explained in 2 episodes
Velocity Dispersion — 5× in this paper · explained in 2 episodes
Core — 3× in this paper · explained in 2 episodes
Abundance — 2× in this paper · explained in 1 episode
Accuracy — 2× in this paper · explained in 1 episode
Calibration — 2× in this paper · explained in 7 episodes
Benchmark — 1× in this paper · explained in 2 episodes
Consistency — 1× in this paper · explained in 3 episodes
European Space Agency — 1× in this paper
Forward modeling — 1× in this paper · explained in 2 episodes
Population Synthesis — 1× in this paper · explained in 2 episodes
Systematic Bias — 1× in this paper · explained in 3 episodes

Summary

The following is a detailed summary of the scientific paper, quoting relevant sections where necessary:

The paper investigates The initial conditions and initial mass functions of Alpha Persei, Pleiades and Praesepe, aiming to determine the stellar Initial Mass Function (IMF) and structural properties of these open clusters. The authors acknowledge that studying the IMF is a fundamental concept in astrophysics, as it plays a central role in connecting star formation on small scales to the evolution of galaxies (Introduction).

Methodology and Data Acquisition

The study utilizes Gaia DR3 proper motions, parallaxes, and photometry to identify cluster members. Cluster membership was determined by filtering using the chi squared test based on proper motions and parallaxes (Equation 1). The sample was supplemented with near-infrared UKIDSS and optical HIPPARCOS survey data, allowing for stellar mass determination down to 0.10 – 0.17 M.

Addressing Observational Biases

The authors address significant observational challenges. A key complication is that "dynamical evolution, where processes such as two-body relaxation, mass segregation, and tidal stripping can preferentially remove low-mass stars, altering the present-day mass function (PDMF) from the initial distribution" (Introduction). Furthermore, unresolved binaries complicate mass estimates.

Binary Correction

To correct for these biases, a Bayesian framework was employed to measure and correct for unresolved binaries in each cluster using photometry and Monte Carlo simulations. The observed stick-out binary fraction (f SO) was found to be (18.3 plus or minus 1.4)%, (16.4 plus or minus 1.2)%, and (18.8 plus or minus 1.3)% for Alpha Persei, Pleiades, and Praesepe respectively, with the overall present-day unresolved binary fractions ranging between (20.0 plus or minus 0.8)% and (23.8 plus or minus 1.2)% (Abstract).

** Mass Estimation**

The stellar mass function (MF) was inferred using a forward modeling approach that accounts for measurement uncertainties on individual stellar masses, marginalizing over the unobserved true mass m true (Equation 6). The authors tested three methods for mass estimation:

Theoretical Isochrone Models: Using PARSEC v2.0 and MIST models, which incorporate updated input physics and pre-main sequence evolution tracks.
Empirical Mass-Luminosity Relation: Using a B-spline fit derived from eclipsing binary systems (Table 6).

** Modeling the IMF**

The stellar mass function is modeled as a multistage power law, which can be broken into two or three stages. The best-fitting initial mass function was found to be described by a three-stage broken power-law distribution with break masses in the ranges 0.24–0.50 M and 0.91–1.20 M and average slopes of alpha med = 1.72 plus or minus 0.09 and alpha high = 2.98 plus or minus 0.22 (Abstract).

** Accounting for Dynamical Evolution**

To determine the most probable initial conditions, the authors employed a novel approach combining N-body simulations with machine learning emulators and Hamiltonian Monte Carlo algorithms. This allowed them to account for the dynamical evolution by simulating cluster states within a King (1962) density profile.

** Key Findings and Conclusions**

The results indicate that the IMF of the open clusters is top-light compared to a Salpeter MF (Abstract). The observed high-mass IMF slope (alpha high = 2.98 plus or minus 0.22) is significantly steeper than field star estimates, suggesting that star formation in low-mass open clusters may be intrinsically deficient in the highest-mass stars relative to the broader galactic field and massive cluster populations (Discussion).

The study also found evidence of scatter: "We find evidence of some scatter in the high mass IMF between the individual clusters, as described by the measured dispersions in mass function slopes of sigma high = 0.29 plus or minus 0.16, whereas the intermediate mass slope shows no significant variation with sigma med = 0.00 plus or minus 0.14 " (Abstract).

In summary, the research provides constraints on the universality of the IMF and its dependence on environmental conditions, demonstrating that while an average open cluster IMF exists (l = 1.24 plus or minus 0.29, m = 1.72 plus or minus 0.09, and h = 2.98 plus or minus 0.22), any exact dependencies on environmental conditions... cannot be accurately determined or quantified from the current sample of open clusters (Conclusion).

Improvements for AI systems

[Disclaimer: The following suggestions assume access to high-performance computing resources and specialized astrophysical datasets. Given the high stakes, these improvements prioritize robust statistical rigor and interpretability.]

Based on the provided data (Stellar Luminosity Functions, parameter posterior distributions, multi-parameter fitting), the current methodology is heavily reliant on sophisticated statistical inference (likely Markov Chain Monte Carlo or similar Bayesian sampling). The primary opportunity for AI improvement lies in transitioning from purely iterative sampling methods to highly efficient, physically constrained deep learning architectures.

Here are three critical improvements and the resulting capabilities of the enhanced AI system:

The Improvement: Replace or augment traditional MCMC sampling for parameter estimation (alpha l, alpha m, alpha h, f bin, ini, etc.) with a deep Variational Inference (VI) architecture. This involves training a neural network to approximate the complex posterior distribution p(theta D) (where theta are the parameters and D is the observed data).

What the Improved AI System Can Do:

Accelerated Convergence: It can map out high-dimensional, correlated parameter spaces (like those involving alpha l, alpha m, alpha h) orders of magnitude faster than traditional MCMC methods. This allows researchers to test significantly larger parameter grids or incorporate more physical variables (e.g., metallicity gradients) without prohibitive computational cost.
Quantified Uncertainty Mapping: The VI framework provides a highly efficient estimate of the full posterior distribution and its associated credible intervals, allowing for rapid calculation of marginalized likelihoods for individual parameters while fully accounting for their interdependence (e.g., how uncertainty in alpha l impacts the derived value of f bin, ini).
Feature Importance: By analyzing the latent space structure of the VI encoder, the system can autonomously rank which input features (e.g., specific magnitude bins in the SLF) contribute most critically to constraining specific astrophysical parameters.

Abstract

We have determined the initial mass function (IMF) and structural properties of the open clusters Alpha Persei, Pleiades and Praesepe using Gaia DR3 astrometry and photometry. Cluster members were identified using primarily Gaia astrometry, supplemented with near-infrared UKIDSS and optical HIPPARCOS survey data with stellar masses down to 0.10-0.17 MSun. We measure and correct for unresolved binaries in each cluster using photometry and Monte Carlo simulations in a Bayesian framework, finding present-day fractions between (20.0 plus or minus0.8)% and (23.8 plus or minus1.2)%. Through a novel approach that combines N-body simulations with machine learning emulators and MCMC algorithms, we have also determined the most probable initial number of stars, binary fraction, half-mass radius and mass function under the assumption that early gas removal does not significantly influence the subsequent cluster evolution. We find a best-fitting initial mass function described by a three-stage broken power-law distribution with break masses in the ranges 0.24-0.50 MSun and 0.91-1.20 MSun and average slopes of alpha med=1.72 plus or minus0.09 and alpha high=2.98 plus or minus0.22. While the low-mass slope remains sensitive to the adopted mass-luminosity relation, this IMF reproduces the present-day clusters well once their dynamical evolution is taken into account. We find evidence of scatter in the high mass IMF between the individual clusters, as described by the measured dispersions in mass function slopes of sigma high=0.29 plus or minus0.16, whereas the intermediate mass slope shows no significant variation with sigma med=0.00 plus or minus0.14. Our findings provide additional evidence for a (compared to a Salpeter MF) top-light IMF and a cluster-to-cluster variation of the IMF. They therefore provide constraints on the universality of the IMF and its dependence on environmental conditions.