-
Chemical evolution models: the role of type Ia supernovae in the $α$-elements over Iron relative abundances and their variations in time and space
Authors:
O. Cavichia,
M. Mollá,
J. J. Bazán,
A. Castrillo,
L. Galbany,
I. Millán-Irigoyen,
Y. Ascasibar,
A. I Díaz,
H. Monteiro
Abstract:
The role of type Ia supernovae (SN Ia), mainly the Delay Time Distributions (DTDs) determined by the binary systems, and the yields of elements created by different explosion mechanisms, are studied by using the MulChem chemical evolution model applied to our Galaxy. We explored 15 DTDs and 12 tables of elemental yields produced by different SN Ia explosion mechanisms, doing a total of 180 models.…
▽ More
The role of type Ia supernovae (SN Ia), mainly the Delay Time Distributions (DTDs) determined by the binary systems, and the yields of elements created by different explosion mechanisms, are studied by using the MulChem chemical evolution model applied to our Galaxy. We explored 15 DTDs and 12 tables of elemental yields produced by different SN Ia explosion mechanisms, doing a total of 180 models. Chemical abundances for $α$-elements (O, Mg, Si, Ca) and Fe derived from these models are compared with recent solar region observational data of $α$-elements over Fe relative abundances, [X/Fe], as a function of [Fe/H] and age. A multidimensional maximum likelihood analysis shows that 52 models are able to fit all these datasets simultaneously, considering the 1$σ$ confidence level. The combination of STROLG1 DTD and LN20181 SN Ia yields provides the best fit. The exponential model with very prompt events is a possible DTD, but a combination of several channels is more probable. The SN Ia yields that include MCh or Near MCh correspond to 39 (75%) of the 52 best models. Regarding the DTD, 31 (60%) of the 52 most probable models correspond to the SD scenario, while the remaining 21 (40%) are based on the DD scenario. Our results also show that the relatively large dispersion of the observational data may be explained by the stellar migration from other radial regions, and/or perhaps a combination of DTDs and explosion channels.
△ Less
Submitted 30 June, 2024; v1 submitted 31 October, 2023;
originally announced November 2023.
-
A resolved study of the inner regions of nearby galaxies with an excess of young massive stars: missing link in the AGN-starburst connection?
Authors:
Guinevere Kauffmann,
Iker Millán-Irigoyen,
Paul Crowther,
Claudia Maraston
Abstract:
We have selected galaxies with very high levels of H$α$ emission (EQW(H$α$)$>$700 Å.) in their central regions from the final data release of the MaNGA survey . Our study focuses on 14 very well-resolved nearby galaxies with stellar masses in the range $9.5 < \log M_*/(M_{\odot}) < 11.5$. We investigate a variety of procedures for selecting galaxy regions that are likely to harbour excess populati…
▽ More
We have selected galaxies with very high levels of H$α$ emission (EQW(H$α$)$>$700 Å.) in their central regions from the final data release of the MaNGA survey . Our study focuses on 14 very well-resolved nearby galaxies with stellar masses in the range $9.5 < \log M_*/(M_{\odot}) < 11.5$. We investigate a variety of procedures for selecting galaxy regions that are likely to harbour excess populations of young massive stars, finding that selection in the 2-dimensional space of extinction-corrected H$α$ EQW and [SIII]/[SII] line ratio produces the best results. By comparing stacked spectra covering these regions with stacked spectra covering normal starburst regions with 100Å$<$EQW(H$α$)$<$200Å, we obtain the following main results: 1) Clear signatures of excess Wolf Rayet stars are found in half of the H$α$ excess regions, 2) Galaxy regions containing excess Wolf Rayet stars are more often associated with the presence of high-ionization emission lines characteristic of accreting black holes. Excess [NeIII] is detected in 4 out of 8 of the WR regions and there are tentative [FeX] detections in 2 galaxies. 3) Regions of the galaxy with excess Wolf Rayet stars are located where the interstellar medium has larger ionized gas turbulent velocities and higher neutral gas overdensities. We make a first attempt to constrain changes in the high mass end of the stellar initial mass function (IMF) using the HR-pyPopStar evolutionary population synthesis models that include high wavelength-resolution theoretical atmosphere libraries for Wolf Rayet stars.
△ Less
Submitted 9 October, 2023;
originally announced October 2023.
-
MEGASTAR (III). Stellar parameters and data products for DR1 late-type stars
Authors:
M. Mollá,
M. L. García-Vargas,
I. Millán-Irigoyen,
N. Cardiel,
E. Carrasco,
A. Gil de Paz,
S. R. Berlanas,
P. Gómez-Álvarez
Abstract:
MEGARA is the optical integral field and multi-object spectrograph at the Gran Telescopio Canarias. We have created MEGASTAR, an empirical library of stellar spectra obtained using MEGARA at high resolution $R=20\,000$ (FWHM), available in two wavelength ranges: one centered in H$α$, from 6420 to 6790\,Å and the other centered in the \ion{Ca}{ii} triplet, from 8370 to 8885\,Å (\mbox{HR-R} and \mbo…
▽ More
MEGARA is the optical integral field and multi-object spectrograph at the Gran Telescopio Canarias. We have created MEGASTAR, an empirical library of stellar spectra obtained using MEGARA at high resolution $R=20\,000$ (FWHM), available in two wavelength ranges: one centered in H$α$, from 6420 to 6790\,Å and the other centered in the \ion{Ca}{ii} triplet, from 8370 to 8885\,Å (\mbox{HR-R} and \mbox{HR-I} VPH-grating configurations). In this work, we use MEGASTAR spectra, combination of these two short wavelength intervals, to estimate the stellar parameters namely effective temperature, surface gravity and metallicity (and their associated errors) for a sample of 351 MEGASTAR members with spectral types earlier than B2. We have applied a $χ^2$ technique by comparing MEGASTAR data to theoretical stellar models. For those stars with stellar parameters derived in the literature, we have obtained a good agreement between those published parameters and ours. Besides the stellar parameters, we also provide several products like the rectified spectra, radial velocities and stellar indices for this sample of stars. In a near future, we will use MEGASTAR spectra and their derived stellar parameters to compute stellar population evolutionary synthesis models, which will contribute to a better interpretation of star clusters and galaxies spectra obtained with MEGARA.
△ Less
Submitted 21 December, 2022;
originally announced December 2022.
-
Stellar Populations in type Ia supernova host galaxies at intermediate-high redshift: Star formation and metallicity enrichment histories
Authors:
I. Millan-Irigoyen,
M. G. del Valle-Espinosa,
R. Fernández-Aranda,
L. Galbany,
J. M. Gomes,
M. Moreno-Raya,
Á. R. López-Sánchez,
M. Mollá
Abstract:
We present a summary of our project that studies galaxies hosting type Ia supernova (SN Ia) at different redshifts. We present Gran Telescopio de Canarias (GTC) optical spectroscopy of six SN Ia host galaxies at redshift $z\sim 0.4-0.5$. They are joined to a set of SN Ia host galaxies at intermediate-high redshift, which include galaxies from surveys SDSS and COSMOS. The final sample, after a sele…
▽ More
We present a summary of our project that studies galaxies hosting type Ia supernova (SN Ia) at different redshifts. We present Gran Telescopio de Canarias (GTC) optical spectroscopy of six SN Ia host galaxies at redshift $z\sim 0.4-0.5$. They are joined to a set of SN Ia host galaxies at intermediate-high redshift, which include galaxies from surveys SDSS and COSMOS. The final sample, after a selection of galaxy spectra in terms of signal-to-noise and other characteristics, consists of 680 galaxies with redshift in the range $0.04 < z < 1$. We perform an inverse stellar population synthesis with the code {\sc fado} to estimate the star formation and enrichment histories of this set of galaxies, simultaneously obtaining their mean stellar age and metallicity and stellar mass. After analysing the correlations among these characteristics, we look for possible dependencies of the Hubble diagram residuals and supernova features (luminosity, color and strength parameter) on these stellar parameters. We find that the Hubble residuals show a clear dependence on the stellar metallicity weighted by mass with a slope of -0.061\,mag\,dex$^{-1}$, when represented in logarithmic scale, $\log{ \langle Z_{M}/Z_{\odot}\rangle }$. This result supports our previous findings obtained from gas oxygen abundances for local and SDSS-survey galaxies. Comparing with other works from the literature that also use the stellar metallicity, we find a similar value, but with more precision and a better significance (2.08 vs $\sim$ 1.1), due to the higher number of objects and wider range of redshift of our sample.
△ Less
Submitted 21 September, 2022;
originally announced September 2022.
-
Predicting interstellar radiation fields from chemical evolution models
Authors:
Mario Romero,
Pablo Corcho-Caballero,
Iker Millán-Irigoyen,
Mercedes Mollá,
Yago Ascasibar
Abstract:
We present a self-consistent prediction of the interstellar radiation field (ISRF), from the extreme ultraviolet (EUV) to sub-mm range, based on two chemical evolution models of a Milky Way-like galaxy (MWG). To this end, we develop a new tool called Mixclask to include gas emission, absorption and scattering from the photoionization code Cloudy into the Monte Carlo radiative transfer code Skirt.…
▽ More
We present a self-consistent prediction of the interstellar radiation field (ISRF), from the extreme ultraviolet (EUV) to sub-mm range, based on two chemical evolution models of a Milky Way-like galaxy (MWG). To this end, we develop a new tool called Mixclask to include gas emission, absorption and scattering from the photoionization code Cloudy into the Monte Carlo radiative transfer code Skirt. Both algorithms are invoked iteratively, until the physical properties of the ISM converge. We have designed a first test, reminiscent of a HII region, and we find that the results of Mixclask are in good agreement with a spherically symmetric Cloudy simulation. Both MWG models based on chemical evolution codes give results broadly consistent with previous empirical models reported in the literature for the ISRF of our Galaxy, albeit they systematically underestimate the mid-infrared emission. We also find significant differences between our two models in the whole ultraviolet range, not fully explored in previous ISRF models. These results show the feasibility of our method of combining radiative transfer with chemical evolution models: there is increased predictive power and the interstellar radiation field obtained provides further constraints on the model parameters. Python source code to implement our method is publicly available at https://github.com/MarioRomeroC/Mixclask.
△ Less
Submitted 8 February, 2023; v1 submitted 9 March, 2022;
originally announced March 2022.
-
HR-pyPopStar: high wavelength-resolution stellar populations evolutionary synthesis model
Authors:
I. Millán-Irigoyen,
M. Mollá,
M. Cerviño,
Y. Ascasibar,
M. L. García-Vargas,
P. R. T. Coelho
Abstract:
We present the HR-pyPopStar model, which provides a complete set (in ages) of high resolution (HR) Spectral Energy Distributions of Single Stellar Populations. The model uses the most recent high wavelength-resolution theoretical atmosphere libraries for main sequence, post-AGB/planetary nebulae and Wolf-Rayet stars. The Spectral Energy Distributions are given for more than a hundred ages ranging…
▽ More
We present the HR-pyPopStar model, which provides a complete set (in ages) of high resolution (HR) Spectral Energy Distributions of Single Stellar Populations. The model uses the most recent high wavelength-resolution theoretical atmosphere libraries for main sequence, post-AGB/planetary nebulae and Wolf-Rayet stars. The Spectral Energy Distributions are given for more than a hundred ages ranging from 0.1 Myr to 13.8 Gyr, at four different values of the metallicity (Z = 0.004, 0.008, 0.019 and 0.05), considering four different IMFs. The wavelength range goes from 91 to 24 000 Å in linear steps δλ = 0.1 Å, giving a theoretical resolving power R_{th,5000} ~ 50 000 at 5000 Å. This is the main novelty of these spectra, unique for their age and wavelength ranges. The models include the ionising stellar populations that are relevant both at young (massive hot stars) as well as old (planetary nebulae) ages. We have tested the results with some examples of HR spectra recently observed with MEGARA at GTC. We highlight the importance of wavelength-resolution in reproducing and interpreting the observational data from the last and forthcoming generations of astronomical instruments operating at 8-10m class telescopes, with higher spectral resolution than their predecessors.
△ Less
Submitted 27 July, 2021;
originally announced July 2021.
-
Chemical evolution of galaxies: emerging dust and the different gas phases in a new multiphase code
Authors:
I. Millán-Irigoyen,
M. Mollá,
Y. Ascasibar
Abstract:
Dust plays an important role in the evolution of a galaxy, since it is one of the main ingredients for efficient star formation. Dust grains are also a sink/source of metals when they are created/destroyed, and, therefore, a self-consistent treatment is key in order to correctly model chemical evolution. In this work, we discuss the implementation of dust physics into our current multiphase model,…
▽ More
Dust plays an important role in the evolution of a galaxy, since it is one of the main ingredients for efficient star formation. Dust grains are also a sink/source of metals when they are created/destroyed, and, therefore, a self-consistent treatment is key in order to correctly model chemical evolution. In this work, we discuss the implementation of dust physics into our current multiphase model, which also follows the evolution of atomic, ionised and molecular gas. Our goal is to model the conversion rates among the different phases of the interstellar medium, including the creation, growth and destruction of dust, based on physical principles rather than phenomenological recipes inasmuch as possible. We first present the updated set of differential equations and then discuss the results. We calibrate our model against observations of the Milky Way Galaxy and compare its predictions with extant data. Our results are broadly consistent with the observed data for intermediate and high metallicities, but the models tend to produce more dust than observed in the low metallicity regime.
△ Less
Submitted 5 March, 2020; v1 submitted 25 April, 2019;
originally announced April 2019.