====== Population Synthesis ======

Population synthesis is a theoretical framework used to interpret the integrated spectra of galaxies as the sum of light emitted by their constituent stellar populations. Because the light of "normal" galaxies originates almost entirely from stars, and stellar evolution is relatively well-understood, astronomers can model a galaxy's spectral energy distribution (SED) by simulating its history of star formation and chemical enrichment.

  * **Fundamental Model Assumptions**: To synthesize a galaxy's population, several empirical and theoretical parameters must be defined:

//Initial Mass Function (IMF):// This specifies the mass distribution of stars at their time of birth. The fraction of stars in a mass interval $dm$ is denoted by $\phi(m)dm$, where the distribution is normalized as $$\int_{ml}^{mu} \ dm \ m \ \phi(m)=1M_{\odot}$$.

The integration limits are typically $m_L \approx 0.08 - 0.1 M_{\odot}$ and $m_U \approx 100 - 150 M_{\odot}$. 

Objects below $m_L$ threshold are unable to achieve the core temperatures and pressures necessary to ignite stable hydrogen fusion. These objects are classified as brown dwarfs rather than true stars.

Stars significantly more massive than $m_U$ are rarely observed for two primary reasons: (1) Observational Constraints: Their lifespans are extremely short (on the order of a few million years), making them statistically rare in any given survey. (2) Physical Instability: According to the theory of stellar structure, stars exceeding this mass generate excessive radiation pressure. This pressure eventually overcomes gravity (reaching the Eddington Limit), preventing the formation of a stable hydrostatic configuration.

A widely used standard is the //**Salpeter IMF**//, defined as $\phi(m) \propto m^{-2.35}$ for stars more massive than $1 M_{\odot}$.

//Star-Formation Rate (SFR):// This is the mass of gas converted into stars per unit time, denoted as $\psi(t)$. In many "standard models," the SFR is assumed to decrease exponentially over time: $\psi(t) = \tau^{-1} \exp[-(t - t_f)/\tau]$, where $\tau$ is the characteristic duration and $t_f$ is the onset of star formation.
 
//Chemical Evolution:// Stars are born with the metallicity ($Z$) of the interstellar medium (ISM) at that time. As stars evolve and die (via supernovae, planetary nebulae, or winds), they eject metal-enriched material back into the ISM, causing $Z(t)$ to increase over time.

  * **Mathematical Formulation**: 

The total spectral luminosity of a galaxy at time $t$, $F_{\lambda}(t)$, is calculated as a convolution of the star-formation history and the spectral energy distribution of the stellar population:

$$F_{\lambda}(t) = \int_{0}^{t} dt' \psi(t - t') S_{\lambda, Z(t-t')}(t')$$

In this equation, $S_{\lambda, Z}(t')$ represents the energy emitted per wavelength by a group of stars with initial metallicity $Z$ and age $t'$. This function accounts for the various evolutionary tracks stars follow in the Hertzsprung–Russell diagram (HRD) and their corresponding isochrones (positions of equal age in the HRD).

[{{ :courses:ast403:spectral_evolution.jpg?600 | Fig 1: Spectrum of a stellar population with solar metallicity that was instantaneously born a time $t$ ago; $t$ is given in units of $10^9$ years}}]

  * **Spectral and Color Evolution**:

The integrated spectrum of a stellar population changes dramatically as it ages:

//Young Populations:// Initially, the light is dominated by the most massive, hot blue stars, which emit intense ultraviolet (UV) radiation.

//Aging Populations:// After approximately $10^7$ years, UV flux diminishes significantly. After $10^8$ to $10^9$ years, massive stars evolve into red supergiants and giants (RGB), causing the flux in the near-infrared (NIR) to increase.

//Color Indices:// As a population ages without new star formation (a process called passive evolution), it becomes systematically redder. For example, the evolution is faster in the $B-V$ color index than in $V-K$.

//Mass-to-Light Ratio ($M/L$):// As stars age, the $M/L$ ratio increases because the total mass of the stars remains relatively constant while their total luminosity decreases. NIR filters (like the K-band) are often used as indicators of total stellar mass because the K-band $M/L$ ratio is less sensitive to age than blue bands.

[{{ :courses:ast403:color_evo.jpg?600 | Fig 2: (left) Evolution of colors between $0\le t \le 17 ×10^9$ yr
for a stellar population with star-formation rate given by the experssion of $\psi(t)$, (right) The dependence of colors and $M/L$ on the metallicity
of the population. The typical colors for four different morphological types of galaxies are plotted. For each $\tau$, the evolution begins at the lower left, i.e., as a blue population in both color indices. In the case of constant star-formation, the population never becomes redder than Irr’s; to achieve redder colors, $\tau$ has to be smaller. }}]



  * **Key Diagnostic Features:**

//The 4000 Â Break:// A prominent jump in the spectra of older stellar populations caused by metal line absorption in stellar atmospheres; it is a primary feature used for //photometric redshift// estimates.

//Age-Metallicity Degeneracy:// A significant challenge in population synthesis is that a galaxy's color can be made redder by either increasing its age or its metallicity. Quantitatively, an increase in age by a factor $X$ is nearly equivalent to an increase in metallicity by a factor $0.65X$.

//Nebular Emission:// While stellar light dominates the continuum, H II regions surrounding young hot stars contribute strong emission lines to the spectrum, which serve as direct tracers of current star-formation rates.

  * **Applications:**

Population synthesis allows astronomers to determine the star-formation histories of galaxies by comparing model predictions with observed colors and spectra. It is essential for interpreting high-redshift observations, tracing the cosmic star-formation history (Madau diagram), and estimating the ages of the oldest galaxies to provide lower limits on the age of the Universe.

  * **Summary:**

• A simple model of star-formation history reproduces
the colors of today’s galaxies fairly well.

• (Most of) the stars in elliptical and S0 galaxies are
old – the earlier the Hubble type, the older the stellar
population.

• Detailed models of population synthesis provide
information about the star-formation history, and
predictions by the models can be compared with
observations of galaxies at high redshift (and thus
smaller age).

  * **More** 
  * Astrobite article about SED fittings: [[https://astrobites.org/2026/01/30/sed-fitting/]]
  * Various SED fitting softwares: [[http://www.sedfitting.org/Welcome.html]]