Stellar Population Synthesis Models Research Articles

Quiescent Galaxies at z ≳ 2.5: Observations versus Models

Abstract The presence of massive quiescent galaxies at high redshifts is still a challenge for most models of galaxy formation. The aim of this work is to compare the observed number density and properties of these galaxies with the predictions of state-of-the-art models. The sample of massive quiescent galaxies has been selected from the COSMOS2015 photometric catalog with z phot ≥ 2.5, , and a specific star formation rate (sSFR) of . The photometric spectral energy distributions (SEDs) of the selected galaxies have been thoroughly analyzed based on different stellar population synthesis models. The final sample includes only those galaxies qualified as quiescent in all SED fitting runs. The observed properties have been compared to theoretical models: the number density of quiescent galaxies with is reproduced by some models, although there is a large scatter in their predictions. Instead, very massive are underpredicted by most of the current models of galaxy formation: some of them, built on the CARNage simulation, are consistent with data up to z ∼ 4, while at higher redshifts the volume of the considered simulation is too small to find such rare objects. Simulated galaxies that match the observed properties in the sSFR–M * plane at z ∼ 3 have been analyzed by reconstructing their evolutionary paths: their merger trees suggest that AGN feedback could be the key process allowing for a rapid quenching of the star formation at z ≳ 4 and that its treatment should be improved in models.

New Constraints on Early-type Galaxy Assembly from Spectroscopic Metallicities of Globular Clusters in M87

Abstract The observed characteristics of globular cluster (GC) systems, such as metallicity distributions, are commonly used to place constraints on galaxy formation models. However, obtaining reliable metallicity values is particularly difficult because of our limited means to obtain high quality spectroscopy of extragalactic GCs. Often, “color–metallicity relations” are invoked to convert easier-to-obtain photometric measurements into metallicities, but there is no consensus on what form these relations should take. In this paper we make use of multiple photometric data sets and iron metallicity values derived from applying full-spectrum stellar population synthesis models to deep Keck/LRIS spectra of 177 GCs centrally located around M87 to obtain a new color–metallicity relation. Our new relation differs substantially from previous relations in the blue, and we present evidence that the M87 relation differs from that of the Milky Way GCs, suggesting environmental dependence of GC properties. We use our color–metallicity relation to derive a new GC metallicity-host galaxy luminosity relation for red and blue GCs and find a shallower relation for the blue GCs than what previous work has found and that the metal-poor GCs are more enriched than what was previously found. This could indicate that the progenitor satellite galaxies that now make up the stellar halos of early-type galaxies are more massive and formed later than previously thought, or that the properties of metal-poor GCs are less dependent on their present-day host, indicating a common origin.

New insight on the far-UV SED and He ii emission from low metallicity galaxies

AbstractUnderstanding the ionizing spectrum of low-metallicity galaxies is of great importance for modeling and interpreting emission line observations of early/distant galaxies.Although a wide suite of stellar evolution, atmosphere, population synthesis, and photoionization models, taking many physical processes into account now exist, all models face a common problem: the inability to explain the presence of nebular Heii emission, which is observed in many low metallicity galaxies, both in UV and optical spectra. Several possible explanations have been proposed in the literature, including Wolf-Rayet (WR) stars, binaries, very massive stars, X-ray sources, or shocks. However, none has so far been able to explain the major observations.We briefly discuss the He ii problem, available empirical data, and observed trends combining X-ray, optical and other studies. We present a simple and consistent physical model showing that X-ray binaries could explain the long-standing nebular Heii problem. Our model, described in Schaerer et al. (2019), successfully explains the observed trends and strength of nebular He ii emission in large samples of low metallicity galaxies and in individual galaxies, which have been studied in detail and with multi-wavelength observations. Our results have in particular important implications for the interpretation of galaxy spectra in the early Universe, which will be obtained with upcoming and future facilities.

Initial Mass Function Variation in Two Elliptical Galaxies Using Near-infrared Tracers

Abstract Using integral field spectroscopy, we demonstrate that gravity-sensitive absorption features in the zJ band (0.9–1.35 μm) can constrain the low-mass stellar initial mass function (IMF) in the cores of two elliptical galaxies, M85 and M87. Compared to the visible bands, the near-infrared (NIR) is more sensitive to light from low-mass dwarf stars, whose relative importance is the primary subject of the debate over IMF variations in nearby galaxies. Our analysis compares the observed spectra to the latest stellar population synthesis models by employing two different methods: equivalent widths and spectral fitting. We find that the IMF slopes in M85 are similar to the canonical Milky Way IMF with a median IMF-mismatch parameter α K = 1.26. In contrast, we find that the IMF in M87 is steeper than a Salpeter IMF with α K = 2.77. The derived stellar population parameters, including the IMF slopes, are consistent with those from recent results in the visible bands based on spectroscopic and kinematic techniques. Certain elemental abundances, e.g., Na and Fe, have dramatic effects on the IMF-sensitive features and therefore the derived IMF slopes. We show evidence for a high [Na/H] ∼ 0.65 dex in the core of M85 from two independent Na i absorption features. The high Na abundance may be the result of a recent galactic merger involving M85. This suggests that including [Na/H] in the stellar population model parameters is critical for constraining the IMF slopes in M85. These results confirm the viability of using NIR absorption features to investigate IMF variation in nearby galaxies.

Globular clusters in the outer halo of M 31

In this paper, we present photometry of 53 globular clusters (GCs) in the M 31 outer halo, including the GALEX far-ultraviolet (FUV) and near-ultraviolet (NUV), SDSS ugriz, 15 intermediate-band filters of BATC, and 2MASS JHKs bands. By comparing the multicolour photometry with stellar population synthesis models, we determine the metallicities, ages, and masses for these GCs, aiming to probe the merging/accretion history of M 31. We find no clear trend of metallicity and mass with the de-projected radius. The halo GCs younger than ∼8 Gyr are mostly located at the de-projected radii around 100 kpc, but this may be due to a selection effect. We also find that the halo GCs have consistent metallicities with their spatially associated substructures, which provides further evidence of the physical association between them. Both the disc and halo GCs in M 31 show a bimodal luminosity distribution. However, we should emphasise that there are more faint halo GCs which are not seen in the disc. The bimodal luminosity function of the halo GCs may reflect a different origin or evolution environment in their original hosts. The M 31 halo GCs include one intermediate metallicity group (−1.5 < [Fe/H] < −0.4) and one metal-poor group ([Fe/H] < −1.5), while the disc GCs have one metal-rich group more. There are considerable differences between the halo GCs in M 31 and the Milky Way (MW). The total number of GCs in M 31 is approximately three times greater than in the MW, however M 31 has about six times more halo GCs than the MW. Compared to the halo GCs of M 31, those of the MW are mostly metal-poor. Both the numerous halo GCs and the higher-metallicity component are suggestive of an active merger history of M 31.

Numerical models for the diffuse ionized gas in galaxies

Context. The diffuse ionized gas (DIG) constitutes the largest fraction of the total ionized interstellar matter in star-forming galaxies, but it is still unclear whether the ionization is driven predominantly by the ionizing radiation of hot massive stars, as in H II regions, or whether additional sources of ionization have to be considered. Key to understanding the ionization mechanisms in the DIG is the line emission by the ionized gas. Aims. We systematically explore a plausible subset of the parameter space involving effective temperatures and metallicities of the ionizing sources, the effects of the hardening of their radiation by surrounding “leaky” H II regions with different escape fractions, as well as different scenarios for the clumpiness of the DIG, and compute the resulting line strength ratios for a number of diagnostic optical emission lines. Methods. For the ionizing fluxes we computed a grid of stellar spectral energy distributions (SEDs) from detailed, fully non-LTE model atmospheres that include the effects of stellar winds and line blocking and blanketing. To calculate the ionization and temperature structure in the interstellar gas we used spherically symmetric photoionization models and state-of-the-art three-dimensional (3D) non-LTE radiative transfer simulations, considering hydrogen, helium, and the most abundant metals. We first applied these methods to classical H II regions around hot stars, using the model SEDs at different metallicities and effective temperatures as ionizing fluxes, and compute the SEDs of the escaping radiation for different escape fractions of hydrogen-ionizing photons. In a second step, we studied the effects of the escaping radiation on the more dilute and extended DIG. Using 3D models simulating a section of a galactic spiral arm, we computed the ionization structure in the DIG for different scenarios for the inhomogeneity of the gas, assuming ionization by a stellar population SED based on plausible parameters. Results. We provide quantitative predictions of how the line ratios from H II regions and the DIG vary as a function of metallicity Z, stellar effective temperature Teff, and escape fraction fesc from the H II region. The range of predicted line ratios reinforces the hypothesis that the DIG is ionized by (filtered) radiation from hot stars. At one-tenth solar metallicity, radiation hardening is mostly due to hydrogen and helium, whereas at solar metallicity absorption by metals plays a significant role. The effects of hardening are seen primarily in the increase in the emission line ratios of the most important cooling lines of the gas, [N II]∕Hβ and [O II]∕Hβ at lower Teff, and [O III]∕Hβ at higher Teff. For low Teff nearly the entire He I-ionizing radiation is absorbed in the H II regions, thereby preventing the formation of high ionization stages such as O III in the DIG. The ionization structure of the DIG depends strongly on both the clumping factor fcl = 〈nH2〉/〈nH2〉 and the large-scale distribution of the gas. In our simulations about 10% of the ionizing radiation produced by hot massive stars in a spiral arm is sufficient to ionize the DIG up to a height of approximately 1 kpc above the galactic plane for a clumping factor close to the observed value of fcl ~ 5. Even small changes in simulation parameters such as the clumping factor can lead to considerable variation in the ionized volume. Both for a more homogeneous gas and a very inhomogeneous gas containing both dense clumps and channels with low gas density, the ionized region in the dilute gas above the galactic plane can cease to be radiation-bounded, allowing the ionizing radiation to leak into the intergalactic medium. Comparison of observed and predicted line ratios indicates that the DIG is typically ionized with a softer SED than predicted by the chosen stellar population synthesis model.

SNITCH: seeking a simple, informative star formation history inference tool

Deriving a simple, analytic galaxy star formation history (SFH) using observational data is a complex task without the proper tool to hand. We therefore present SNITCH, an open source code written in Python, developed to quickly (~2 minutes) infer the parameters describing an analytic SFH model from the emission and absorption features of a galaxy spectrum dominated by star formation gas ionisation. SNITCH uses the Flexible Stellar Population Synthesis models of Conroy et al. (2009), the MaNGA Data Analysis Pipeline and a Markov Chain Monte Carlo method in order to infer three parameters (time of quenching, rate of quenching and model metallicity) which best describe an exponentially declining quenching history. This code was written for use on the MaNGA spectral data cubes but is customisable by a user so that it can be used for any scenario where a galaxy spectrum has been obtained, and adapted to infer a user-defined analytic SFH model for specific science cases. Herein we outline the rigorous testing applied to SNITCH and show that it is both accurate and precise at deriving the SFH of a galaxy spectra. The tests suggest that SNITCH is sensitive to the most recent epoch of star formation but can also trace the quenching of star formation even if the true decline does not occur at an exponential rate. With the use of both an analytical SFH and only five spectral features, we advocate that this code be used as a comparative tool across a large population of spectra, either for integral field unit data cubes or across a population of galaxy spectra.

The cosmic spectral energy distribution in the EAGLE simulation

The cosmic spectral energy distribution (CSED) is the total emissivity as a function of wavelength of galaxies in a given cosmic volume. We compare the observed CSED from the UV to the submm to that computed from the EAGLE cosmological hydrodynamical simulation, post-processed with stellar population synthesis models and including dust radiative transfer using the SKIRT code. The agreement with the data is better than 0.15 dex over the entire wavelength range at redshift $z=0$, except at UV wavelengths where the EAGLE model overestimates the observed CSED by up to a factor 2. Global properties of the CSED as inferred from CIGALE fits, such as the stellar mass density, mean star formation density, and mean dust-to-stellar-mass ratio, agree to within better than 20 per cent. At higher redshift, EAGLE increasingly underestimates the CSED at optical-NIR wavelengths with the FIR/submm emissivity underestimated by more than a factor of 5 by redshift $z=1$. We believe that these differences are due to a combination of incompleteness of the EAGLE-SKIRT database, the small simulation volume and the consequent lack of luminous galaxies, and our lack of knowledge on the evolution of the characteristics of the interstellar dust in galaxies. The impressive agreement between the simulated and observed CSED at lower $z$ confirms that the combination of EAGLE and SKIRT dust processing yields a fairly realistic representation of the local Universe.

Initial mass function variations cannot explain the ionizing spectrum of low metallicity starbursts

Aims. Observations of both galaxies in the distant Universe and local starbursts are showing increasing evidence for very hard ionizing spectra that stellar population synthesis models struggle to reproduce. Here we explore the effects of the assumed stellar initial mass function (IMF) on the ionizing photon output of young populations at wavelengths below key ionization energy thresholds. Methods. We use a custom set of binary population and spectral synthesis (BPASS) models to explore the effects of IMF assumptions as a function of metallicity, IMF slope, upper mass limit, IMF power law break mass and sampling. Results. We find that while the flux capable of ionizing hydrogen is only weakly dependent on IMF parameters, the photon flux responsible for the He II and O VI lines is far more sensitive to assumptions. In our current models this flux arises primarily from helium and Wolf-Rayet stars which have partially or fully lost their hydrogen envelopes. The timescales for formation and evolution of both Wolf Rayet stars and helium dwarfs, and hence inferred population age, are affected by choice of model IMF. Even the most extreme IMFs cannot reproduce the He II ionizing flux observed in some high redshift galaxies, suggesting a source other than stellar photospheres. Conclusions. We caution that detailed interpretation of features in an individual galaxy spectrum is inevitably going to be subject to uncertainties in the IMF of its contributing starbursts. We remind the community that the IMF is fundamentally a statistical construct, and that stellar population synthesis models are most effective when considering entire galaxy populations rather than individual objects.

Dark matter in massive early-type galaxies: 100 years of the jeans equations

In 1919 James Jeans published the book Problems of Cosmogony and Stellar Dynamics in which he summarized his work on dynamics of stellar systems based on his papers published from 1915 onwards. To mark the centenary of this publication we present here one application of his work relevant for contemporary research of galaxies: we analyze the problem of dark matter in massive early-type galaxies (ellipticals and lenticulars) using various available observational data. After discussing the basics of the Jeans equations we study their application on the integrated stellar spectra of galaxies, planetary nebula data and, especially, globular cluster data. We rely on both Newtonian and MOND frameworks and show their advantages and drawbacks. To infer the contribution of the dark component in early-type galaxies we rely on several stellar population synthesis models. It is shown that dark matter does not dominate in the inner regions of early-type galaxies, but becomes more important beyond three effective radii.

Spectroscopic Fundamental Parameters of 16 Young Large Magellanic Cloud Star Clusters

Flux-calibrated integrated spectra have been obtained for a sample of 16 blue compact and high surface brightness Large Magellanic Cloud (LMC) star clusters. The spectra cover the wavelength range 3700–6800 Å, with a typical resolution and dispersion of 14 Å and 3.5 Å/pixel, respectively. We simultaneously derive foreground E(B − V) reddening values and ages for the clusters by comparing the line strengths and continuum distribution of their integrated spectra with those of LMC cluster spectra templates. As a first approach to obtain cluster ages, we also use the equivalent width (EW) of Balmer absorption lines and some diagnostic diagrams involving the sum of the EWs of selected spectral lines. Through evolutionary synthesis models for simple stellar populations and using a full-spectrum fitting code, we also derive reddenings, ages, and metallicities (Z/Z⊙) for the cluster sample. No metallicity determinations have been previously performed for nine out of the 16 clusters here studied. We find good agreement between the cluster ages derived from the template-matching method and those ages estimated from theoretical evolutionary synthesis models. All of the studied clusters are found to be young, with ages ranging from 15 Myr (NGC 1693) to 80 Myr (NGC 2160), affected by reddening values typical for the LMC: . The derived metallicities Z/Z⊙ range from 0.2 up to 1.0 (solar metal content). We compare our results with those reported in the literature, and discuss the distribution of metallicity values. In an effort to create a spectral library at LMC metallicity level with several clusters per age range, the present cluster sample stands out as a useful complement to previous ones.

Mapping Lyman Continuum Escape in Tololo 1247–232

Abstract Low-redshift, spatially resolved Lyman continuum (LyC) emitters allow us to clarify the processes for LyC escape from these starburst galaxies. We use Hubble Space Telescope (HST) WFC3 and ACS imaging of the confirmed low-redshift LyC emitter Tol 1247−232 to study the ionization structure of the gas and its relation to the ionizing star clusters. We perform ionization parameter mapping (IPM) using [O iii] λλ4959, 5007 and [O ii] λ3727 imaging as the high- and low-ionization tracers, revealing broad, large-scale, optically thin regions originating from the center and reaching the outskirts of the galaxy, consistent with LyC escape. We carry out stellar population synthesis modeling of the 26 brightest clusters using our HST photometry. Combining these data with the nebular photometry, we find a global LyC escape fraction of f esc = 0.12, with uncertainties also consistent with zero escape and all measured f esc values for this galaxy. Our analysis suggests that, similar to other candidate LyC emitters, a two-stage starburst has taken place in this galaxy, with a 12 Myr old, massive central cluster likely having precleared regions in and around the center and the second generation of 2–4 Myr old clusters dominating the current ionization, including some escape from the galaxy.

Mid-IR colors and surface brightness fluctuations as tracers of stellar mass-loss in the TP-AGB

AbstractI present integrated colors and surface brightness fluctuation magnitudes in the mid-IR, derived from stellar population synthesis models that include the effects of the dusty envelopes around TP-AGB stars. The models are based on the Bruzual & Charlot CB* isochrones; they are single-burst, range in age from a few Myr to 14 Gyr, and comprise metallicities between Z = 0.0001 and Z = 0.04. I compare these models to mid-IR data of AGB stars and star clusters in the Magellanic Clouds, and study the effects of varying self-consistently the mass-loss rate, the stellar parameters, and the output spectra of the stars plus their dusty envelopes.

A new method to identify subclasses among AGB stars using Gaia and 2MASS photometry

Aims. We explore the wealth of high-quality photometric data provided by data release 2 (DR2) of the Gaia mission for long-period variables (LPVs) in the Large Magellanic Cloud (LMC). Our goal is to identify stars of various types and masses along the asymptotic giant branch. Methods. For this endeavour, we developed a new multi-band approach combining Wesenheit functions WRP,BP−RP and WKs,J−Ks in the Gaia BP, RP, and 2MASS J, Ks spectral ranges, respectively, and use a new diagram, (WRP,BP−RP − WKs,J−Ks) versus Ks, to distinguish between different kinds of stars in our sample of LPVs. We used stellar population synthesis models to validate our approach. Results. We demonstrate the ability of the new diagram to discriminate between O- and C-rich objects, and to identify low-mass, intermediate-mass, and massive O-rich red giants, as well as extreme C-rich stars. Stellar evolution and population synthesis models guide the interpretation of the results, highlighting the diagnostic power of the new tool to discriminate between stellar initial masses, chemical properties, and evolutionary stages.

Clocking the Evolution of Post-starburst Galaxies: Methods and First Results

Abstract Detailed modeling of the recent star formation histories (SFHs) of post-starburst (or “E+A”) galaxies is impeded by the degeneracy between the time elapsed since the starburst ended (post-burst age), the fraction of stellar mass produced in the burst (burst strength), and the burst duration. To resolve this issue, we combine GALEX ultraviolet photometry, SDSS photometry and spectra, and new stellar population synthesis models to fit the SFHs of 532 post-starburst galaxies. In addition to an old stellar population and a recent starburst, 48% of the galaxies are best fit with a second recent burst. Lower stellar mass galaxies (log M ⋆/M ☉ < 10.5) are more likely to experience two recent bursts, and the fraction of their young stellar mass is more strongly anticorrelated with their total stellar mass. Applying our methodology to other, younger post-starburst samples, we identify likely progenitors to our sample and examine the evolutionary trends of molecular gas and dust content with post-burst age. We discover a significant (4σ) decline, with a 117–230 Myr characteristic depletion time, in the molecular gas to stellar mass fraction with the post-burst age. The implied rapid gas depletion rate of 2–150 M ☉ yr−1 cannot be due to current star formation, given the upper limits on the current star formation rates in these post-starbursts. Nor are stellar winds or supernova feedback likely to explain this decline. Instead, the decline points to the expulsion or destruction of molecular gas in outflows, a possible smoking gun for active galactic nucleus feedback.

The stellar population and initial mass function of NGC 1399 with MUSE

We present spatially resolved measurements of the stellar initial mass function (IMF) in NGC 1399, the largest elliptical galaxy in the Fornax Cluster. Using data from the Multi Unit Spectroscopic Explorer (MUSE) and updated state-of-the-art stellar population synthesis models from Conroy et. al 2018, we use full spectral fitting to measure the low-mass IMF, as well as a number of individual elemental abundances, as a function of radius in this object. We find that the IMF in NGC 1399 is heavier than the Milky Way in its centre and remains radially constant at a super-salpeter slope out to 0.7 R$_{e}$. At radii larger than this, the IMF slope decreases to become marginally consistent with a Milky Way IMF just beyond R$_{e}$. The inferred central V-band M/L ratio is in excellent agreement with the previously reported dynamical M/L measurement from Houghton et. al 2006. The measured radial form of the M/L ratio may be evidence for a two-phase formation in this object, with the central regions forming differently to the outskirts. We also report measurements of a spatially resolved filament of ionised gas extending 4 arcseconds (404 pc at $D_{\mathrm{L}}=21.1$ Mpc) from the centre of NGC 1399, with very narrow equivalent width and low velocity dispersion ($65\pm14$ \kms). The location of the emission, combined with an analysis of the emission line ratios, leads us to conclude that NGC 1399's AGN is the source of ionising radiation.

DirtyGrid I: 3D Dust Radiative Transfer Modeling of Spectral Energy Distributions of Dusty Stellar Populations

Abstract Understanding the properties of stellar populations and interstellar dust has important implications for galaxy evolution. In normal star-forming galaxies, stars and the interstellar medium dominate the radiation from ultraviolet (UV) to infrared (IR). In particular, interstellar dust absorbs and scatters UV and optical light, re-emitting the absorbed energy in the IR. This is a strongly nonlinear process that makes independent studies of the UV-optical and IR susceptible to large uncertainties and degeneracies. Over the years, UV to IR spectral energy distribution (SED) fitting utilizing varying approximations has revealed important results on the stellar and dust properties of galaxies. Yet the approximations limit the fidelity of the derived properties. There is sufficient computer power now available that it is now possible to remove these approximations and map out of landscape of galaxy SEDs using full dust radiative transfer. This improves upon previous work by directly connecting the UV, optical, and IR through dust grain physics. We present the DirtyGrid, a grid of radiative transfer models of SEDs of dusty stellar populations in galactic environments designed to span the full range of physical parameters of galaxies. Using the stellar and gas radiation input from the stellar population synthesis model PEGASE, our radiative transfer model DIRTY self-consistently computes the UV to far-IR/sub-mm SEDs for each set of parameters in our grid. DIRTY computes the dust absorption, scattering, and emission from the local radiation field and a dust grain model, thereby physically connecting the UV-optical to the IR. We describe the computational method and explain the choices of parameters in DirtyGrid. The computation took millions of CPU hours on supercomputers, and the SEDs produced are an invaluable tool for fitting multi-wavelength data sets. We provide the complete set of SEDs in an online table.

Re-evaluating old stellar populations

Determining the properties of old stellar populations (those with age >1 Gyr) has long involved the comparison of their integrated light, either in the form of photometry or spectroscopic indexes, with empirical or synthetic templates. Here we reevaluate the properties of old stellar populations using a new set of stellar population synthesis models, designed to incorporate the effects of binary stellar evolution pathways as a function of stellar mass and age. We find that single-aged stellar population models incorporating binary stars, as well as new stellar evolution and atmosphere models, can reproduce the colours and spectral indices observed in both globular clusters and quiescent galaxies. The best fitting model populations are often younger than those derived from older spectral synthesis models, and may also lie at slightly higher metallicities.

Multiple populations within globular clusters in early-type galaxies exploring their effect on stellar initial mass function estimates

It is now well-established that most (if not all) ancient globular clusters host multiple popula- tions, that are characterised by distinct chemical features such as helium abundance variations along with N-C and Na-O anti-correlations, at fixed [Fe/H]. These very distinct chemical fea- tures are similar to what is found in the centres of the massive early-type galaxies and may influence measurements of the global properties of the galaxies. Additionally, recent results have suggested that M/L variations found in the centres of massive early-type galaxies might be due to a bottom-heavy stellar initial mass function. We present an analysis of the effects of globular cluster-like multiple populations on the integrated properties of early-type galaxies. In particular, we focus on spectral features in the integrated optical spectrum and the global mass-to-light ratio that have been used to infer variations in the stellar initial mass function. To achieve this we develop appropriate stellar population synthesis models and take into account, for the first time, an initial-final mass relation which takes into consideration a varying He abundance. We conclude that while the multiple populations may be present in massive early-type galaxies, they are likely not responsible for the observed variations in the mass- to-light ratio and IMF sensitive line strengths. Finally, we estimate the fraction of stars with multiple populations chemistry that come from disrupted globular clusters within massive ellipticals and find that they may explain some of the observed chemical patterns in the centres of these galaxies.

Star-forming Galaxies as AGN Imposters? A Theoretical Investigation of the Mid-infrared Colors of AGNs and Extreme Starbursts

Abstract We conduct for the first time a theoretical investigation of the mid-infrared spectral energy distribution (SED) produced by dust heated by an active galactic nucleus (AGN) and an extreme starburst. These models employ an integrated modeling approach using photoionization and stellar population synthesis models in which both the line and emergent continuum is predicted from gas exposed to the ionizing radiation from a young starburst and an AGN. In this work, we focus on the infrared colors from the Wide-field Infrared Survey Explorer, predicting the dependence of the colors on the input radiation field, the interstellar medium conditions, the obscuring column, and the metallicity. We find that an extreme starburst can mimic an AGN in two band mid-infrared color cuts employed in the literature. However, the three-band color cuts employed in the literature require starbursts with extremely high ionization parameters or gas densities. We show that the extreme mid-infrared colors seen in some blue compact dwarf galaxies are not due to metallicity but rather a combination of high ionization parameters and high column densities. Based on our theoretical calculations, we present a theoretical mid-infrared color cut that will exclude even the most extreme starburst that we have modeled in this work. The theoretical AGN demarcation region presented here can be used to identify elusive AGN candidates for future follow-up studies with the James Webb Space Telescope. The full suite of simulated SEDs are available online.