Galaxy-Wide Stellar Initial Mass Function

• Galaxy-wide Stellar Initial Mass Function (gwIMF) is a composite stellar mass distribution built over finite star-formation epochs that reflects local physical conditions.
• IGIMF theory integrates cluster-scale IMFs over space and time, linking molecular-cloud star formation with galaxy-wide diagnostics and evolution.
• Empirical evidence indicates gwIMF variations with SFR and metallicity, impacting stellar mass estimates, supernova rates, and chemical enrichment models.

The galaxy-wide stellar initial mass function (gwIMF) is the initial mass function of all newly formed stars in a whole galaxy, considered over a finite star-formation epoch rather than within a single star-forming cloud or embedded cluster. In this usage, the gwIMF is a composite, time-dependent, and potentially environment-dependent stellar mass distribution that links molecular-cloud-scale star formation to galaxy evolution, stellar population synthesis, chemical enrichment, remnant production, and the cosmological matter cycle. A central theme of modern gwIMF work is that the gwIMF need not equal the stellar IMF of any one star-forming region, and need not be universal across galaxies or across cosmic time (Kroupa et al., 2024, Kroupa et al., 2021).

1. Definitions and conceptual status

The basic IMF notation is

$dN=\xi(m)\,dm,$

with logarithmic form

$$dN=\xi(m)\,dm=\xi_{\rm L}(m)\, d\log_{10}m,\qquad \xi_{\rm L}(m)=\left(m\ln 10\right)\xi(m).$$

Within this terminology, the stellar IMF refers to the IMF of stars born in one molecular cloud clump or embedded cluster, the composite IMF is the sum of more than one stellar IMF in a larger region, and the gwIMF is the composite IMF of a whole galaxy; it can vary with time (Kroupa et al., 2024).

A core conceptual distinction is that the gwIMF is not merely the canonical stellar IMF applied to a larger system. Kroupa and Jeřábková formulate the limiting case in which the IMF is interpreted as a universal invariant probability density distribution function: then

${\rm IMF}={\rm cIMF}={\rm gwIMF}.$

If, instead, the IMF varies and/or is not a PDF, then in general

${\rm IMF}\ne{\rm cIMF}\ne{\rm gwIMF}.$

This distinction is central to IGIMF theory and to most modern discussions of galaxy-integrated IMF variation (Kroupa et al., 2021).

The gwIMF matters because it controls the galaxy-wide production of ionizing photons, core-collapse supernovae, white dwarfs, neutron stars, black holes, metals, and long-lived low-mass stars. It therefore enters directly into star-formation-rate calibrations, stellar mass estimates, chemical evolution models, dynamical mass inferences, and the interpretation of unresolved galactic light (Kroupa et al., 2024).

2. Formal constructions of the gwIMF

A general composite-IMF expression is

$$\xi_{\rm cIMF}(m) = \int_V \, \int_P \, \int_t \, \xi_{\rm dens}(m: {\vec r}, \{P_{\rm i}\}^{N_{\rm P}},t) \, {\rm d}V \, {\rm d}^{N_{\rm P}}P_{\rm i} \, {\rm d}t,$$

where the local stellar-mass-dependent star-formation rate density depends on position, time, and physical parameters such as metallicity, temperature, pressure, angular momentum, and cosmic-ray field. This is the most general statement of the gwIMF problem: the galaxy-wide IMF emerges by integrating local star formation over space, physical conditions, and time (Kroupa et al., 2024).

In IGIMF theory, the gwIMF is built from embedded clusters. For a galaxy treated over one star-formation epoch $\delta t$,

$$\xi_{\rm IGIMF}(m;t) = \int_{M_{\rm ecl,min}}^{M_{\rm ecl,max}(\psi(t))} \xi\!\left(m\le m_{\rm max}\left(M_{\rm ecl}\right):Z\right)\,\xi_{\rm ecl}(M_{\rm ecl})\,dM_{\rm ecl}.$$

Here $M_{\rm ecl,min}\approx 5\,M_\odot$, $M_{\rm ecl,max}$ depends on the galaxy-wide SFR $\psi$, $$dN=\xi(m)\,dm=\xi_{\rm L}(m)\, d\log_{10}m,\qquad \xi_{\rm L}(m)=\left(m\ln 10\right)\xi(m).$$0 is the embedded-cluster mass function (ECMF), and the stellar IMF in each cluster is truncated at the cluster-dependent maximum stellar mass $$dN=\xi(m)\,dm=\xi_{\rm L}(m)\, d\log_{10}m,\qquad \xi_{\rm L}(m)=\left(m\ln 10\right)\xi(m).$$1 (Kroupa et al., 2024, Kroupa et al., 2021).

The ECMF is written as

$$dN=\xi(m)\,dm=\xi_{\rm L}(m)\, d\log_{10}m,\qquad \xi_{\rm L}(m)=\left(m\ln 10\right)\xi(m).$$2

with typically $$dN=\xi(m)\,dm=\xi_{\rm L}(m)\, d\log_{10}m,\qquad \xi_{\rm L}(m)=\left(m\ln 10\right)\xi(m).$$3, and in one empirical parameterization

$$dN=\xi(m)\,dm=\xi_{\rm L}(m)\, d\log_{10}m,\qquad \xi_{\rm L}(m)=\left(m\ln 10\right)\xi(m).$$4

A second bridge from local to global scales is the relation between the maximum cluster mass and the galaxy-wide SFR,

$$dN=\xi(m)\,dm=\xi_{\rm L}(m)\, d\log_{10}m,\qquad \xi_{\rm L}(m)=\left(m\ln 10\right)\xi(m).$$5

Low-SFR galaxies therefore do not form very massive clusters and become top-light; high-SFR galaxies form more massive clusters and can become top-heavy (Kroupa et al., 2024, Kroupa et al., 2021).

The canonical stellar IMF used as benchmark is a piecewise power law,

$$dN=\xi(m)\,dm=\xi_{\rm L}(m)\, d\log_{10}m,\qquad \xi_{\rm L}(m)=\left(m\ln 10\right)\xi(m).$$6

with canonical slopes

$$dN=\xi(m)\,dm=\xi_{\rm L}(m)\, d\log_{10}m,\qquad \xi_{\rm L}(m)=\left(m\ln 10\right)\xi(m).$$7

$$dN=\xi(m)\,dm=\xi_{\rm L}(m)\, d\log_{10}m,\qquad \xi_{\rm L}(m)=\left(m\ln 10\right)\xi(m).$$8

$$dN=\xi(m)\,dm=\xi_{\rm L}(m)\, d\log_{10}m,\qquad \xi_{\rm L}(m)=\left(m\ln 10\right)\xi(m).$$9

IGIMF-based work then allows these slopes, especially ${\rm IMF}={\rm cIMF}={\rm gwIMF}.$0, to vary with cluster mass, density, and metallicity. One widely used high-mass prescription is

${\rm IMF}={\rm cIMF}={\rm gwIMF}.$1

with

${\rm IMF}={\rm cIMF}={\rm gwIMF}.$2

while low-mass slopes can vary with metallicity as

${\rm IMF}={\rm cIMF}={\rm gwIMF}.$3

with ${\rm IMF}={\rm cIMF}={\rm gwIMF}.$4 and ${\rm IMF}={\rm cIMF}={\rm gwIMF}.$5 in the IGIMF-2021 implementation (Haslbauer et al., 2024, Zonoozi et al., 5 Feb 2025).

3. Environmental and temporal variation

A systematic IGIMF grid shows that the gwIMF varies over wide ranges of SFR and metallicity. For

${\rm IMF}={\rm cIMF}={\rm gwIMF}.$6

and

${\rm IMF}={\rm cIMF}={\rm gwIMF}.$7

the gwIMF is top-heavy relative to the canonical IMF; for

${\rm IMF}={\rm cIMF}={\rm gwIMF}.$8

the gwIMF becomes top-light regardless of metallicity; and for

${\rm IMF}={\rm cIMF}={\rm gwIMF}.$9

the gwIMF can become bottom-heavy regardless of SFR. The same framework predicts that massive ellipticals should have formed with a gwIMF that was top-heavy within the first few hundred Myr and then evolved into a bottom-heavy gwIMF in the metal-enriched galactic center (Jerabkova et al., 2018).

The GalIMF chemical-evolution framework makes this explicitly time dependent by recalculating the gwIMF every

${\rm IMF}\ne{\rm cIMF}\ne{\rm gwIMF}.$0

with the gwIMF at each step set by the current galaxy-wide SFR and metallicity. In that formalism, the resulting gwIMF, stellar yields, supernova rates, and metal abundance co-evolve, so that dwarf galaxies naturally become top-light at low SFR while massive early-type galaxies experience top-heavy phases at high SFR and bottom-heavy low-mass behavior after metal enrichment (Yan, 2021).

Time dependence also appears in chemically motivated toy models for giant ellipticals. Ferreras et al. argue that a time-independent bottom-heavy gwIMF cannot reproduce, simultaneously, the high fraction of low-mass stars, the high stellar metallicity, the ${\rm IMF}\ne{\rm cIMF}\ne{\rm gwIMF}.$1-enhancement, and the remnant-related constraints. Their illustrative two-stage models use an early phase with

${\rm IMF}\ne{\rm cIMF}\ne{\rm gwIMF}.$2

for

${\rm IMF}\ne{\rm cIMF}\ne{\rm gwIMF}.$3

followed by a later bottom-heavy phase with

${\rm IMF}\ne{\rm cIMF}\ne{\rm gwIMF}.$4

or a less extreme early phase with

${\rm IMF}\ne{\rm cIMF}\ne{\rm gwIMF}.$5

showing that the gwIMF of a massive elliptical is more plausibly a time-dependent history than a single immutable function (Weidner et al., 2013).

A related Galactic fossil argument comes from binary population synthesis of carbon-enhanced metal-poor stars. That work concludes that the effective IMF of the low-metallicity Milky Way halo was more high-mass dominated than the present-day IMF and likely transitioned near

${\rm IMF}\ne{\rm cIMF}\ne{\rm gwIMF}.$6

supporting a metallicity- and epoch-dependent Galactic gwIMF rather than a universal one (Suda et al., 2013).

4. Empirical evidence across galaxy populations

One of the foundational galaxy-scale results is the ATLAS${\rm IMF}\ne{\rm cIMF}\ne{\rm gwIMF}.$7 dynamical study of 260 nearby early-type galaxies by Cappellari et al. Using two-dimensional SAURON stellar kinematics and axisymmetric Jeans Anisotropic Models, the paper compares a dynamical stellar mass-to-light ratio to a population-synthesis value based on a Salpeter reference IMF,

${\rm IMF}\ne{\rm cIMF}\ne{\rm gwIMF}.$8

through the ratio

${\rm IMF}\ne{\rm cIMF}\ne{\rm gwIMF}.$9

The inferred IMF normalization varies systematically from Kroupa/Chabrier-like at low $$\xi_{\rm cIMF}(m) = \int_V \, \int_P \, \int_t \, \xi_{\rm dens}(m: {\vec r}, \{P_{\rm i}\}^{N_{\rm P}},t) \, {\rm d}V \, {\rm d}^{N_{\rm P}}P_{\rm i} \, {\rm d}t,$$0 to Salpeter-like and heavier-than-Salpeter at high $$\xi_{\rm cIMF}(m) = \int_V \, \int_P \, \int_t \, \xi_{\rm dens}(m: {\vec r}, \{P_{\rm i}\}^{N_{\rm P}},t) \, {\rm d}V \, {\rm d}^{N_{\rm P}}P_{\rm i} \, {\rm d}t,$$1, producing differences “up to a factor of three in mass.” This is a galaxy-integrated, but centrally weighted, result: the kinematics typically extend to about one projected half-light radius $$\xi_{\rm cIMF}(m) = \int_V \, \int_P \, \int_t \, \xi_{\rm dens}(m: {\vec r}, \{P_{\rm i}\}^{N_{\rm P}},t) \, {\rm d}V \, {\rm d}^{N_{\rm P}}P_{\rm i} \, {\rm d}t,$$2, so it constrains an effective gwIMF normalization for the inner $$\xi_{\rm cIMF}(m) = \int_V \, \int_P \, \int_t \, \xi_{\rm dens}(m: {\vec r}, \{P_{\rm i}\}^{N_{\rm P}},t) \, {\rm d}V \, {\rm d}^{N_{\rm P}}P_{\rm i} \, {\rm d}t,$$3, not the full virialized galaxy (Cappellari et al., 2012).

A second, independent central-aperture test uses stellar chromospheric $$\xi_{\rm cIMF}(m) = \int_V \, \int_P \, \int_t \, \xi_{\rm dens}(m: {\vec r}, \{P_{\rm i}\}^{N_{\rm P}},t) \, {\rm d}V \, {\rm d}^{N_{\rm P}}P_{\rm i} \, {\rm d}t,$$4 emission in the centers of NGC 1407 and NGC 2695. Because stars below $$\xi_{\rm cIMF}(m) = \int_V \, \int_P \, \int_t \, \xi_{\rm dens}(m: {\vec r}, \{P_{\rm i}\}^{N_{\rm P}},t) \, {\rm d}V \, {\rm d}^{N_{\rm P}}P_{\rm i} \, {\rm d}t,$$5 contribute only a few percent of the integrated $$\xi_{\rm cIMF}(m) = \int_V \, \int_P \, \int_t \, \xi_{\rm dens}(m: {\vec r}, \{P_{\rm i}\}^{N_{\rm P}},t) \, {\rm d}V \, {\rm d}^{N_{\rm P}}P_{\rm i} \, {\rm d}t,$$6-band light but about $$\xi_{\rm cIMF}(m) = \int_V \, \int_P \, \int_t \, \xi_{\rm dens}(m: {\vec r}, \{P_{\rm i}\}^{N_{\rm P}},t) \, {\rm d}V \, {\rm d}^{N_{\rm P}}P_{\rm i} \, {\rm d}t,$$7 of the integrated $$\xi_{\rm cIMF}(m) = \int_V \, \int_P \, \int_t \, \xi_{\rm dens}(m: {\vec r}, \{P_{\rm i}\}^{N_{\rm P}},t) \, {\rm d}V \, {\rm d}^{N_{\rm P}}P_{\rm i} \, {\rm d}t,$$8 luminosity for a Kroupa IMF, the $$\xi_{\rm cIMF}(m) = \int_V \, \int_P \, \int_t \, \xi_{\rm dens}(m: {\vec r}, \{P_{\rm i}\}^{N_{\rm P}},t) \, {\rm d}V \, {\rm d}^{N_{\rm P}}P_{\rm i} \, {\rm d}t,$$9 ratio is a low-mass-star diagnostic. The measured $\delta t$0-to-$\delta t$1-band flux ratio is a factor of

$\delta t$2

higher in NGC 1407 than in NGC 2695, in agreement with previously inferred IMF differences from gravity-sensitive optical absorption lines. However, the measurement is explicitly confined to the central $\delta t$3 COS aperture and is therefore a probe of aperture-scale IMF variation rather than a direct gwIMF measurement (Dokkum et al., 2021).

For star-forming galaxies, Gunawardhana et al. use the GAMA survey to infer the effective high-mass IMF slope from dust-corrected H$\delta t$4 equivalent width and broadband color. In their highest-redshift volume-limited sample, the inferred high-mass slope obeys

$\delta t$5

and, as a function of SFR surface density,

$\delta t$6

In that analysis, low-SFR galaxies are consistent with Salpeter-like or steeper high-mass slopes, whereas highly star-forming galaxies require flatter slopes, implying a more top-heavy effective galaxy-integrated IMF at the high-mass end (Gunawardhana et al., 2011).

Resolved Milky Way populations provide a complementary route. Gaia DR3 open-cluster mass functions fitted with a Kroupa-like broken power law show that the break mass increases systematically with cluster age, with Spearman

$\delta t$7

strengthening after distance-bias correction to

$\delta t$8

while the intermediate- and high-mass slopes do not show significant age dependence. This is interpreted as direct evidence that the IMF is not universal across Galactic stellar populations and implies that the Milky Way’s effective gwIMF should evolve as the mixture of cluster-scale IMFs changes with epoch (Steinhardt et al., 24 Mar 2026).

Field-star star counts sharpen the low-mass side of the argument. In the Solar neighborhood, M-dwarf counts in the range $\delta t$9 show an effective low-mass slope rising from about

$$\xi_{\rm IGIMF}(m;t) = \int_{M_{\rm ecl,min}}^{M_{\rm ecl,max}(\psi(t))} \xi\!\left(m\le m_{\rm max}\left(M_{\rm ecl}\right):Z\right)\,\xi_{\rm ecl}(M_{\rm ecl})\,dM_{\rm ecl}.$$0

to about

$$\xi_{\rm IGIMF}(m;t) = \int_{M_{\rm ecl,min}}^{M_{\rm ecl,max}(\psi(t))} \xi\!\left(m\le m_{\rm max}\left(M_{\rm ecl}\right):Z\right)\,\xi_{\rm ecl}(M_{\rm ecl})\,dM_{\rm ecl}.$$1

For dynamically cold stars over $$\xi_{\rm IGIMF}(m;t) = \int_{M_{\rm ecl,min}}^{M_{\rm ecl,max}(\psi(t))} \xi\!\left(m\le m_{\rm max}\left(M_{\rm ecl}\right):Z\right)\,\xi_{\rm ecl}(M_{\rm ecl})\,dM_{\rm ecl}.$$2,

$$\xi_{\rm IGIMF}(m;t) = \int_{M_{\rm ecl,min}}^{M_{\rm ecl,max}(\psi(t))} \xi\!\left(m\le m_{\rm max}\left(M_{\rm ecl}\right):Z\right)\,\xi_{\rm ecl}(M_{\rm ecl})\,dM_{\rm ecl}.$$3

whereas dynamically hot stars over the same interval show

$$\xi_{\rm IGIMF}(m;t) = \int_{M_{\rm ecl,min}}^{M_{\rm ecl,max}(\psi(t))} \xi\!\left(m\le m_{\rm max}\left(M_{\rm ecl}\right):Z\right)\,\xi_{\rm ecl}(M_{\rm ecl})\,dM_{\rm ecl}.$$4

This indicates that the low-mass IMF depends on both metallicity and formation epoch or environment, providing an empirical ingredient for time-dependent gwIMF models (Li et al., 2023).

A broader synthesis of low-mass gwIMF constraints across ultra-faint dwarfs, the SMC, Milky Way field stars, and metal-rich early-type galaxies finds that a metallicity-dependent low-mass IMF model is preferred over an invariant one. In an AIC comparison on a heterogeneous $$\xi_{\rm IGIMF}(m;t) = \int_{M_{\rm ecl,min}}^{M_{\rm ecl,max}(\psi(t))} \xi\!\left(m\le m_{\rm max}\left(M_{\rm ecl}\right):Z\right)\,\xi_{\rm ecl}(M_{\rm ecl})\,dM_{\rm ecl}.$$5 compilation, the invariant model gives

$$\xi_{\rm IGIMF}(m;t) = \int_{M_{\rm ecl,min}}^{M_{\rm ecl,max}(\psi(t))} \xi\!\left(m\le m_{\rm max}\left(M_{\rm ecl}\right):Z\right)\,\xi_{\rm ecl}(M_{\rm ecl})\,dM_{\rm ecl}.$$6

a logarithmic metallicity law gives

$$\xi_{\rm IGIMF}(m;t) = \int_{M_{\rm ecl,min}}^{M_{\rm ecl,max}(\psi(t))} \xi\!\left(m\le m_{\rm max}\left(M_{\rm ecl}\right):Z\right)\,\xi_{\rm ecl}(M_{\rm ecl})\,dM_{\rm ecl}.$$7

and a $$\xi_{\rm IGIMF}(m;t) = \int_{M_{\rm ecl,min}}^{M_{\rm ecl,max}(\psi(t))} \xi\!\left(m\le m_{\rm max}\left(M_{\rm ecl}\right):Z\right)\,\xi_{\rm ecl}(M_{\rm ecl})\,dM_{\rm ecl}.$$8-linear metallicity law gives

$$\xi_{\rm IGIMF}(m;t) = \int_{M_{\rm ecl,min}}^{M_{\rm ecl,max}(\psi(t))} \xi\!\left(m\le m_{\rm max}\left(M_{\rm ecl}\right):Z\right)\,\xi_{\rm ecl}(M_{\rm ecl})\,dM_{\rm ecl}.$$9

supporting the view that metal-poor systems are bottom-light or near canonical while metal-rich systems are bottom-heavy in the low-mass regime (Yan et al., 2024).

5. Modeling frameworks and galaxy-scale consequences

Integrated-light IMF inference is itself model dependent. A hierarchical Bayesian framework for composite stellar populations shows that fitting a realistic composite galaxy spectrum with a single SSP biases the inferred IMF slope to a higher value than the true one, whereas using two SSPs removes most of the bias and allows recovery of the imposed IMF–velocity-dispersion trend. The same framework reconstructs the variable IMF of the mock spectra for signal-to-noise ratios exceeding $M_{\rm ecl,min}\approx 5\,M_\odot$0, which is directly relevant to spectroscopic gwIMF work in unresolved galaxies (Dries et al., 2017).

For galaxy photometry and star-formation diagnostics, the photGalIMF extension of GalIMF couples environment-dependent gwIMFs to chemical evolution and stellar isochrones. In local star-forming galaxies, it finds that stellar masses inferred from near-IR light differ from canonical estimates by about a factor of $M_{\rm ecl,min}\approx 5\,M_\odot$1, while H$M_{\rm ecl,min}\approx 5\,M_\odot$2-based SFRs differ by factors of about $M_{\rm ecl,min}\approx 5\,M_\odot$3 on average and in some cases up to $M_{\rm ecl,min}\approx 5\,M_\odot$4. The resulting IGIMF-corrected local main sequence is fit as

$M_{\rm ecl,min}\approx 5\,M_\odot$5

for

$M_{\rm ecl,min}\approx 5\,M_\odot$6

and the gas-depletion time is fit by

$M_{\rm ecl,min}\approx 5\,M_\odot$7

These results operationalize gwIMF variation in observationally accessible quantities (Haslbauer et al., 2024).

A complementary SPS implementation, SPS-VarIMF, computes spectra, luminosities, colors, and remnant populations for time-dependent gwIMFs under the IGIMF theory. It finds that late-type galaxies can, in principle, have their underlying gwIMF constrained from UV and optical colors and luminosities, whereas early-type galaxies remain photometrically degenerate under different gwIMF assumptions even though their gwIMF-dependent $M_{\rm ecl,min}\approx 5\,M_\odot$8 ratios differ by up to an order of magnitude. In that framework, massive present-day elliptical galaxies would have been

$M_{\rm ecl,min}\approx 5\,M_\odot$9

times as bright as at present when they were forming (Zonoozi et al., 5 Feb 2025).

gwIMF variation also propagates into galaxy scaling relations. In an IGIMF-based reinterpretation of the baryonic Tully–Fisher relation, the true baryonic mass is

$M_{\rm ecl,max}$0

with

$M_{\rm ecl,max}$1

Assuming instead a constant

$M_{\rm ecl,max}$2

causes massive galaxies to have their masses in stars and remnants underestimated, while low-mass gas-dominated galaxies are less affected. In this picture, some apparent high-mass BTFR offsets are reinterpreted as stellar-population modeling effects induced by a varying gwIMF rather than purely dynamical anomalies (Zonoozi et al., 16 Jul 2025).

Historically, semi-analytic galaxy formation work also explored IMF variation at the whole-galaxy level. One precursor implemented a redshift-dependent modified Chabrier IMF with characteristic mass evolving from

$M_{\rm ecl,max}$3

to

$M_{\rm ecl,max}$4

finding that the evolving IMF changes the gas recycle fraction $M_{\rm ecl,max}$5 in

$M_{\rm ecl,max}$6

and thereby affects the SFR–$M_{\rm ecl,max}$7 relation and stellar-mass-density accounting. This model was explicitly phenomenological and universal at fixed redshift rather than a self-consistent IGIMF, but it illustrates how strongly galaxy-scale inference depends on gwIMF assumptions (Kang et al., 2010).

6. Scope, caveats, and active controversies

A recurring limitation is that many methods constrain IMF normalization more robustly than IMF shape. The ATLAS$M_{\rm ecl,max}$8 dynamical analysis demonstrates systematic variation of the effective IMF normalization, but it cannot distinguish uniquely between a dwarf-rich and a remnant-rich heavy IMF; the same $M_{\rm ecl,max}$9 can arise from bottom-heavy and remnant-heavy power-law cases. In gwIMF language, this is strong evidence for a non-universal galaxy-scale mass normalization, not a unique recovery of the full functional form (Cappellari et al., 2012).

A second limitation is spatial scope. Some of the strongest early-type evidence is centrally weighted: Cappellari et al. constrain the inner $\psi$0, and the chromospheric $\psi$1 test is confined to the $\psi$2 central aperture. This suggests that central IMF gradients can bias attempts to identify a single galaxy-wide IMF, especially in massive ellipticals where radial variation is itself part of the observational picture [(Cappellari et al., 2012); (Dokkum et al., 2021)].

Method dependence is also substantial. The low-mass gwIMF synthesis literature emphasizes that star counts, stellar population synthesis, and galaxy chemical evolution probe different mass ranges, use different IMF priors, and carry different systematics. In that setting, metallicity appears to capture the dominant trend across most current data, but the same review stresses systematic discrepancies between methods, IMF functional-form dependence, binary-star systematics, SFH degeneracies, and the likelihood that metallicity is not the only relevant variable (Yan et al., 2024).

The broader IGIMF framework remains physically motivated rather than universally accepted. Current debates include whether IMF variation with density and metallicity is universally established, whether temperature or pressure must enter explicitly beyond $\psi$3 and $\psi$4, whether stochastic sampling can still explain some cluster data, how the IMF of Population III stars should be treated, and whether the brown-dwarf IMF has an environmental dependence. The theory therefore functions as an active program of multi-scale star formation rather than a closed empirical law (Kroupa et al., 2024, Kroupa et al., 2021).

Taken together, the literature supports a precise but limited conclusion: the gwIMF is best understood as an emergent, time-dependent, and environment-dependent composite IMF, distinct from the stellar IMF of a single embedded cluster. Low-SFR galaxies tend toward top-light gwIMFs, high-SFR systems toward top-heavy gwIMFs, and high-metallicity populations toward bottom-heavy low-mass behavior, while both central-aperture and population-level evidence argue against strict IMF universality. What remains unsettled is the exact mapping from local cloud physics to the full radial, temporal, and functional form of the gwIMF in real galaxies (Jerabkova et al., 2018, Kroupa et al., 2024, Yan, 2021).