J/A+A/678/A61 Analytic solution of the CEM with Type Ia SNe (Palicio+, 2023)
Analytic solution of chemical evolution models with Type Ia supernovae.
I. Disc bimodality in the [alpha/Fe] versus [Fe/H] plane and other applications.
Palicio P.A., Spitoni E., Recio-Blanco A., Matteucci F., Peirani S.,
Greggio L.
<Astron. Astrophys. 678, A61 (2023)>
=2023A&A...678A..61P 2023A&A...678A..61P (SIMBAD/NED BibCode)
ADC_Keywords: Supernovae ; Models
Keywords: Galaxy: abundances - Galaxy: evolution - solar neighborhood -
ISM: general - evolution
Abstract:
In the last years, a significant number of works have focussed on
finding analytic solutions for the chemical enrichment models of
galactic systems, including the Milky Way. Some of these solutions,
however, cannot account for the enrichment produced by Type Ia SNe due
to the presence of the delay time distributions (DTDs) in the models.
We present a new analytic solution for the chemical evolution model of
the Galaxy. This solution can be used with different prescriptions of
the DTD, including the single and double degenerate scenarios, and
allows the inclusion of an arbitrary number of pristine gas infalls.
We integrate the chemical evolution model by extending the
instantaneous recycling approximation with the contribution of Type Ia
SNe. This implies an extra term in the modelling that depends on the
DTD. For those DTDs that lead to non-analytic integrals, we describe
them as a superposition of Gaussian, exponential and 1/t functions
using a restricted least-squares fitting method.
We obtain the exact solution for a chemical model with Type Ia SNe
widely used in previous works, avoiding numerical integration errors.
This solution can reproduce the expected chemical evolution of the
alpha and iron-peak elements in less computing time than numerical
integration methods. We compare the pattern in the [Si/Fe] vs. [Fe/H]
plane observed by APOGEE DR17 with that predicted by the model. We
find the low alpha sequence can be explained by a delayed gas infall.
We exploit the applicability of our solution by modelling the chemical
evolution of a simulated Milky Way-like galaxy from its star formation
history. The implementation of our solution is released as a Python
package.
Our solution constitutes a promising tool for the Galactic Archaeology
and is able to model the observed trends in alpha element abundances
versus [Fe/H] in the solar neighbourhood. We infer the chemical
information of a simulated galaxy modelled without Chemistry.
Description:
Values of the parameters introduced in the modelling of the delay time
distributions (DTDs) considered in the paper. These quantities were
computed by performing a restricted least-squares fitting algorithm
that accounts for the shape of each DTD (see Appendix A of the paper.)
File Summary:
--------------------------------------------------------------------------------
FileName Lrecl Records Explanations
--------------------------------------------------------------------------------
ReadMe 80 . This file
table1.dat 130 36 Values of the parameters of the DTDs considered
in this work
--------------------------------------------------------------------------------
Byte-by-byte Description of file: table1.dat
--------------------------------------------------------------------------------
Bytes Format Units Label Explanations
--------------------------------------------------------------------------------
1- 8 A8 --- DTDname Identifier of the DTD (1)
10- 21 F12.9 1/Gyr AG ? Amplitude of the Gaussian term
23- 33 F11.9 Gyr taup ? Offset of the Gaussian term
35- 45 F11.9 Gyr sigmap ? Width of the Gaussian term
47- 51 F5.3 Gyr tau1G ? Lower time limit of the Gaussian term
53- 58 F6.3 Gyr tau2G ? Upper time limit of the Gaussian term
60- 72 F13.9 1/Gyr AE ? Amplitude of the exponential term
74- 85 F12.9 Gyr tauD ? Characteristic timescale of the
exponential term
87- 91 F5.3 Gyr tau1E ? Lower time limit of the exponential term
93- 98 F6.3 Gyr tau2E ? Upper time limit of the exponential term
100-111 F12.9 1/Gyr AI ? Amplitude of the inverse-of-time term
113-117 F5.3 Gyr tau0 ? Offset of the inverse-of-time term
119-123 F5.3 Gyr tau1I ? Lower time limit of the inverse-of-time term
125-130 F6.3 Gyr tau2I ? Upper time limit of the inverse-of-time term
--------------------------------------------------------------------------------
Note (1): Identifier of the DTD as follows:
MR01 = Single Degenerate (Matteucci & Recchi, 2001ApJ...558..351M 2001ApJ...558..351M)
G05WIDE = WIDE Double Degenerate (Greggio, 2005A&A...441.1055G 2005A&A...441.1055G)
G05CLOSE = CLOSE Double Degenerate (Greggio, 2005A&A...441.1055G 2005A&A...441.1055G)
MVP06 = Empirical bimodal distribution (Mannucci et al.,
2006MNRAS.372..265M 2006MNRAS.372..265M)
T08 = Empirical ∝ t-1 (Totani et al., 2008PASJ...60.1327T 2008PASJ...60.1327T)
P08 = Empirical ∝ t-1/2 (Pritchet et al., 2008ApJ...683L..25P 2008ApJ...683L..25P)
S05 = Empirical Gaussian (Strolger et al., 2004ApJ...613..200S 2004ApJ...613..200S,
Cat. J/ApJ/613/200, 2005ApJ...635.1370S 2005ApJ...635.1370S)
--------------------------------------------------------------------------------
Acknowledgements:
Pedro Alonso Palicio, pedro.alonso-palicio(at)oca.eu
(End) Patricia Vannier [CDS] 28-Jul-2023