Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Maschberger (2013)

San Diego State University

Abstract (paraphrased)

Proposes the L3L_3 IMF, a three-parameter functional form that is a heavy-tailed approximation to the lognormal: a low-mass power law and a high-mass power law joined smoothly. The standard Kroupa (2001) and Chabrier (2003) IMFs are essentially indistinguishable from it. Its decisive practical advantage is that the cumulative distribution function and its inverse (the quantile / mass-generating function) are closed-form, so sampling needs no special functions and no Newton iteration.

The L3L_3 functional form (verified against the paper, Table 1)

With the auxiliary function G(m)G(m) (Table 1, Eq. 1),

G(m)=(1+(mμ)1α)1β,G(m) = \left(1 + \left(\tfrac{m}{\mu}\right)^{1-\alpha}\right)^{1-\beta},

the pdf, CDF, and quantile are (Table 1, Eqs. 2–4)

pL3(m)=A(mμ)α(1+(mμ)1α)β,A=(1α)(1β)μ1G(mu)G(ml),p_{L_3}(m) = A\left(\tfrac{m}{\mu}\right)^{-\alpha} \left(1 + \left(\tfrac{m}{\mu}\right)^{1-\alpha}\right)^{-\beta}, \qquad A = \frac{(1-\alpha)(1-\beta)}{\mu}\,\frac{1}{G(m_u)-G(m_l)},
PL3(m)=G(m)G(ml)G(mu)G(ml),m(u)=μ[(u[G(mu)G(ml)]+G(ml))11β1]11α.P_{L_3}(m) = \frac{G(m)-G(m_l)}{G(m_u)-G(m_l)}, \qquad m(u) = \mu\left[\Big(u\,[G(m_u)-G(m_l)] + G(m_l)\Big)^{\frac{1}{1-\beta}} - 1\right]^{\frac{1}{1-\alpha}}.

The canonical single-star parameters (Table 1; system/binary values in parentheses):

parametervaluemeaning
α\alpha2.3 (2.3)2.3\ (2.3)high-mass exponent
β\beta1.4 (2.0)1.4\ (2.0)low-mass turnover
μ\mu0.2 (0.2)M0.2\ (0.2)\,M_\odotscale parameter
mlm_l0.01M0.01\,M_\odotlower limit (normalization)
mum_u150M150\,M_\odotupper limit (normalization)

The effective low-mass exponent is γ=α+β(1α)=0.48\gamma = \alpha + \beta(1-\alpha) = 0.48 (single-star; Table 1, Eq. 5), i.e. p(m)m0.48p(m)\propto m^{-0.48} as m0m\to 0. The limits ml,mum_l, m_u “are only needed for the normalization” (Table 1 caption). Maschberger states the Salpeter exponent as α=+2.35\alpha = +2.35 in the linear convention (§2.1).

Use in progenax

Notes

progenax default IMF. The closed-form inverse CDF makes it the only canonical IMF that samples in O(1)\mathcal{O}(1) per particle without Newton iterations — important for fast, differentiable IC generation. The infamous “peak” of the L3L_3 pdf is at the scale parameter μ\mu, which is not where the two power laws cross (Table 1, Eq. 11).

References
  1. Maschberger, T. (2013). On the function describing the stellar initial mass function. Monthly Notices of the Royal Astronomical Society, 429, 1725–1733. 10.1093/mnras/sts479