Analytic solution of Chemical Evolution Models with Type Ia SNe

P. A. Palicio,E. Spitoni,A. Recio-Blanco,F. Matteucci,S. Peirani,L. Greggio
DOI: https://doi.org/10.1051/0004-6361/202346567
2023-07-23
Abstract:Context: In the last years, a significant number of works have focused 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. Aims: 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. Methods: We integrate the chemical evolution model by extending the instantaneous recycling approximation with the contribution of Type Ia SNe. 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. Results: We obtain the exact solution for a chemical model with Type Ia SNe widely used in previous works. 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. Conclusions: 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.
Astrophysics of Galaxies
What problem does this paper attempt to address?
### Problems the Paper Attempts to Solve This paper aims to address the issue in the chemical evolution models of the Milky Way, where the presence of the delay time distribution (DTD) of Type Ia supernovae (Type Ia SNe) prevents obtaining an analytical solution. Specifically, many existing chemical evolution models can analytically describe the changes in metallicity in the Milky Way, but these models often ignore the contribution of Type Ia supernovae because the DTD of Type Ia supernovae complicates the analytical solution. Therefore, this paper proposes a new analytical solution method that can handle different forms of DTD and is applicable to both single and multiple gas inflow scenarios. ### Main Contributions 1. **Proposal of Analytical Solution**: The paper proposes a new analytical solution for a chemical evolution model that includes the contribution of Type Ia supernovae. This analytical solution can handle different forms of DTD, including single-degenerate and double-degenerate scenarios. 2. **Model Applicability**: The analytical solution can be used not only to describe the chemical evolution of the Milky Way but also to apply to other Milky Way-like galaxies. 3. **Computational Efficiency**: Compared to numerical integration methods, the analytical solution can accurately predict the chemical evolution of α-elements and iron-peak elements in a shorter computation time. 4. **Application Examples**: The paper validates the effectiveness of the new analytical solution by comparing the APOGEE DR17 observational data with the model predictions. Additionally, it demonstrates how to use this analytical solution to study the bimodal distribution of α-element abundance in the Galactic disk. ### Method Overview - **Model Extension**: The paper extends the instantaneous recycling approximation (IRA) by introducing the contribution of Type Ia supernovae, which requires adding an extra term dependent on DTD to the model. - **DTD Handling**: For DTDs that lead to non-analytical integrals, the paper describes them using a combination of Gaussian, exponential, and 1/t functions and fits them using a constrained least squares fitting method. - **Derivation of Analytical Solution**: The paper derives the evolution equations of the surface mass density of various elements over time in the chemical model through analytical methods and decomposes them into IRA contributions and Type Ia supernova contributions. ### Results and Conclusions - **Exact Solution**: The paper obtains an exact solution for a widely used chemical model that includes Type Ia supernovae, avoiding numerical integration errors. - **Computational Efficiency**: The analytical solution is superior in computation time compared to numerical integration methods, allowing for rapid prediction of the chemical evolution of α-elements and iron-peak elements. - **Model Validation**: The model's effectiveness is validated by comparing it with APOGEE DR17 observational data. In particular, the low-α sequence can be explained by delayed gas inflow. - **Application Prospects**: The analytical solution provides a powerful tool for Galactic archaeology, capable of simulating the relationship between α-element abundance and [Fe/H] in the solar neighborhood. ### Summary This paper proposes a new analytical solution method that can handle chemical evolution models including the contribution of Type Ia supernovae, providing a more accurate and efficient tool for studying the chemical evolution of the Milky Way and other similar galaxies.