Stellar Population Models Based on the Milky Way Interstellar Medium-Corrected MaStar Stellar Library

This file contains simple stellar population (SSP) model spectra constructed from a version of the SDSS-IV MaNGA Stellar Library (MaStar; Yan et al. 2019, Abdurro'uf et al. 2022) that has been corrected for the effects of absorption in the CaII 3934, 3969 and NaI D 5891, 5897 transitions arising in the Milky Way's interstellar medium (ISM).  These corrections are described in full in Rubin et al. (2025), and our approach to constructing these SSP models is described in Maraston et al. (2020) and Rubin et al. (2025).  In brief, our models are calculated with the evolutionary population synthesis code of Maraston (2005), which is based on the fuel consumption theorem for the evaluation of the energetics of post-Main Sequence phases.  We use the calibrated median values of the stellar parameters calculated for the MaStar sample to generate representative stellar spectra as functions of effective temperature, surface gravity, and chemical composition.  These representative spectra are then used as input for the stellar population models.  The stellar parameter estimates are described in R. Yan et al. (2025, in preparation) and at https://www.sdss4.org/dr17/mastar/mastar-stellar-parameters/.  

We calculate SSPs using stars in metallicity bins centered at [Z/H] = -1.35, -0.33, 0.0, and +0.35 with an approximate bin width of 0.1 dex assuming a Salpeter IMF.  The SSP ages span 3 Myr to 15 Gyr and are calculated at 51 gridpoints.  For comparison, we also calculate the equivalent SSPs using the uncorrected MaStar spectra.  The datamodel is described below.

HDU1: 51 x 4 x 1 x 3 matrix describing the parameters of each SSP spectrum.  Each gridpoint (i,j,k) contains a 3-element array listing the age (in Gyr), metallicity, and IMF slope (in linear mass units)

HDU2: 2 x 4563 array containing the vacuum wavelength and spectral resolution (R) grids for models constructed from the uncorrected (original) stellar library.  The wavelength sampling is logarithmic and the wavelengths have units of Angstroms.  R = wave / (FWHM dwave)

HDU3: 51 x 4 x 1 x 4563 matrix containing the SSPs constructed from the uncorrected (original) stellar library in units of erg/s/Ang/Msun  

HDU4: 2 x 4542 array containing the vacuum wavelength and spectral resolution (R) grids for models constructed from the corrected (cleaned) stellar library.  The wavelength sampling is logarithmic and the wavelengths have units of Angstroms.  R = wave / (FWHM dwave)

HDU5: 51 x 4 x 1 x 4542 matrix containing the SSPs constructed from the corrected (cleaned) stellar library in units of erg/s/Ang/Msun