Showing 1–50 of 78 results for author: Spitoni, E

A Census of Sun's Ancestors and their Contributions to the Solar System Chemical Composition

Authors: F. Fiore, F. Matteucci, E. Spitoni, M. Molero, P. Salucci, D. Romano, A. Vasini

Abstract: In this work we compute the rates and numbers of different types of stars and phenomena (SNe, novae, white dwarfs, merging neutron stars, black holes) that contributed to the chemical composition of the Solar System. Stars die and restore the newly formed elements into the interstellar gas. This process is called "chemical evolution". In particular, we analyse the death rates of stars of all masse… ▽ More In this work we compute the rates and numbers of different types of stars and phenomena (SNe, novae, white dwarfs, merging neutron stars, black holes) that contributed to the chemical composition of the Solar System. Stars die and restore the newly formed elements into the interstellar gas. This process is called "chemical evolution". In particular, we analyse the death rates of stars of all masses, dying either quiescently or explosively. These rates and total star numbers are computed in the context of a revised version of the two-infall model for the chemical evolution of the Milky Way, which reproduces fairly well the observed abundance patterns of several chemical species, as well as the global solar metallicity. We compute also the total number of stars ever born and still alive as well as the number of stars born up to the formation of the Solar System and with a mass and metallicity like the Sun. This latter number will account for all the possible existing Solar Systems which can host life in the solar vicinity. Among all the stars (from 0.8 to 100 M$_{\odot}$) born and died from the beginning up to the Solar System formation epoch, which contributed to its chemical composition, 93.00$\%$ are represented by stars dying as single white dwarfs (without interacting significantly with a companion star) and originating in the mass range 0.8-8 M$_{\odot}$, while 5.24$\%$ are neutron stars and 0.73$\%$ are black holes, both originating from SNe core-collapse (M>8 M$_{\odot}$); 0.64$\%$ are Type Ia SNe and 0.40$\%$ are nova systems, both originating from the same mass range as the white dwarfs. The number of stars similar to the Sun born from the beginning up to the Solar System formation, with metallicity in the range 12+log(Fe/H)= 7.50 $\pm$ 0.04 dex is 3.1732$\cdot$ 10$^{7}$, and in particular our Sun is the 2.6092$\cdot$ 10$^7$-th star of this kind, born in the solar vicinity. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: Submitted to A&A, 11 pages, 9 figures

Chemical evolution of the Galactic bulge with different stellar populations

Authors: M. Molero, F. Matteucci, E. Spitoni, A. Rojas-Arriagada, R. M. Rich

Abstract: The metallicity distribution function (MDF) of the Galactic bulge features a multi-peak shape, with a metal-poor peak at [Fe/H]=-0.3 dex and a metal-rich peak at [Fe/H]=+0.3 dex. This bimodality is also seen in [alpha/Fe] versus [Fe/H] ratios, indicating different stellar populations in the bulge. We aim to replicate the observed MDF by proposing a scenario where the metal-poor bulge stars formed… ▽ More The metallicity distribution function (MDF) of the Galactic bulge features a multi-peak shape, with a metal-poor peak at [Fe/H]=-0.3 dex and a metal-rich peak at [Fe/H]=+0.3 dex. This bimodality is also seen in [alpha/Fe] versus [Fe/H] ratios, indicating different stellar populations in the bulge. We aim to replicate the observed MDF by proposing a scenario where the metal-poor bulge stars formed in situ during an intense star formation burst, while the metal-rich stars formed during a second burst and/or were accreted from the inner Galactic disk due to a growing bar. We used a chemical evolution model that tracks various chemical species with detailed nucleosynthesis, focusing on Fe production from both Type Ia supernovae and massive stars, including rotating massive stars with varying velocities. Our model also accounts for gas infall, outflow, and the effect of stellar migration. Results are compared to 13,000 stars from the SDSS/APOGEE survey within 3.5 kpc of the Galactic center. Our model successfully reproduces the double-peak shape of the bulge MDF and the alpha-element abundance trends relative to Fe by assuming (i) a multi-burst star formation history with a 250 Myr quenching of the first burst and (ii) stellar migration from the inner disk due to a growing bar. We estimate that about 40% of the bulge-bar's stellar mass originates from the inner disk. Nucleosynthesis models that assume either no rotation for massive stars or a rotational velocity distribution favoring slow rotation at high metallicities best match the observed MDF and [alpha/Fe] and [Ce/Fe] versus [Fe/H] abundance patterns. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: Accepted for publication in A&A, 12 pages, 12 figures

(Re)mind the gap: a hiatus in star formation history unveiled by APOGEE DR17

Authors: E. Spitoni, F. Matteucci, R. Gratton, B. Ratcliffe, I. Minchev, G. Cescutti

Abstract: The analysis of several spectroscopic surveys indicates the presence of a bimodality between the disc stars in the abundance ratio space of [$α$/Fe] versus [Fe/H]. The two stellar groups are commonly referred to as the high-$α$ and low-$α$ sequences. Some models capable of reproducing such a bimodality, invoke the presence of a hiatus in the star formation history in our Galaxy, whereas other mode… ▽ More The analysis of several spectroscopic surveys indicates the presence of a bimodality between the disc stars in the abundance ratio space of [$α$/Fe] versus [Fe/H]. The two stellar groups are commonly referred to as the high-$α$ and low-$α$ sequences. Some models capable of reproducing such a bimodality, invoke the presence of a hiatus in the star formation history in our Galaxy, whereas other models explain the two sequences by means of stellar migration. Our aim is to show that the existence of the gap in the star formation rate between high-$α$ and low-$α$ is evident in the stars of APOGEE DR17, if one plots [Fe/$α$] versus [$α$/H], thus confirming previous suggestions by Gratton et al. (1996) and Fuhrmann (1998). Then we try to interpret the data by means of detailed chemical models. We compare the APOGEE DR17 red giant stars with the predictions of a detailed chemical evolution model based on the two-infall paradigm, taking also into account possible accretion of dwarf satellites. The APOGEE DR17 abundance ratios [Fe/$α$] versus [$α$/H] exhibit a sharp increase of [Fe/$α$] at a nearly constant [$α$/H] (where $α$ elements considered are Mg, Si, O) during the transition between the two disc phases. This observation strongly supports the hypothesis that a hiatus in star formation occurred during this evolutionary phase. Notably, the most pronounced growth in the [Fe/$α$] versus [$α$/H] relation is observed for oxygen, as this element is exclusively synthesised in core-collapse supernovae. A chemical model predicting a stop in the star formation of a duration of roughly 3.5 Gyr, and where the high-$α$ disc starts forming from pre-enriched gas by a previous encounter with a dwarf galaxy can well explain the observations. △ Less

Submitted 17 May, 2024; originally announced May 2024.

Comments: Submitted to A&A, 14 pages, 10 figures

Discovery of a dormant 33 solar-mass black hole in pre-release Gaia astrometry

Authors: Gaia Collaboration, P. Panuzzo, T. Mazeh, F. Arenou, B. Holl, E. Caffau, A. Jorissen, C. Babusiaux, P. Gavras, J. Sahlmann, U. Bastian, Ł. Wyrzykowski, L. Eyer, N. Leclerc, N. Bauchet, A. Bombrun, N. Mowlavi, G. M. Seabroke, D. Teyssier, E. Balbinot, A. Helmi, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne , et al. (390 additional authors not shown)

Abstract: Gravitational waves from black-hole merging events have revealed a population of extra-galactic BHs residing in short-period binaries with masses that are higher than expected based on most stellar evolution models - and also higher than known stellar-origin black holes in our Galaxy. It has been proposed that those high-mass BHs are the remnants of massive metal-poor stars. Gaia astrometry is exp… ▽ More Gravitational waves from black-hole merging events have revealed a population of extra-galactic BHs residing in short-period binaries with masses that are higher than expected based on most stellar evolution models - and also higher than known stellar-origin black holes in our Galaxy. It has been proposed that those high-mass BHs are the remnants of massive metal-poor stars. Gaia astrometry is expected to uncover many Galactic wide-binary systems containing dormant BHs, which may not have been detected before. The study of this population will provide new information on the BH-mass distribution in binaries and shed light on their formation mechanisms and progenitors. As part of the validation efforts in preparation for the fourth Gaia data release (DR4), we analysed the preliminary astrometric binary solutions, obtained by the Gaia Non-Single Star pipeline, to verify their significance and to minimise false-detection rates in high-mass-function orbital solutions. The astrometric binary solution of one source, Gaia BH3, implies the presence of a 32.70 \pm 0.82 M\odot BH in a binary system with a period of 11.6 yr. Gaia radial velocities independently validate the astrometric orbit. Broad-band photometric and spectroscopic data show that the visible component is an old, very metal-poor giant of the Galactic halo, at a distance of 590 pc. The BH in the Gaia BH3 system is more massive than any other Galactic stellar-origin BH known thus far. The low metallicity of the star companion supports the scenario that metal-poor massive stars are progenitors of the high-mass BHs detected by gravitational-wave telescopes. The Galactic orbit of the system and its metallicity indicate that it might belong to the Sequoia halo substructure. Alternatively, and more plausibly, it could belong to the ED-2 stream, which likely originated from a globular cluster that had been disrupted by the Milky Way. △ Less

Submitted 19 April, 2024; v1 submitted 16 April, 2024; originally announced April 2024.

Comments: 23 pages, accepted fro publication in A&A Letters. New version with small fixes

MINCE II. Neutron capture elements

Authors: P. François, G. Cescutti, P. Bonifacio, E. Caffau, L. Monaco, M. Steffen, J. Puschnig, F. Calura, S. Cristallo, P. Di Marcantonio, V. Dobrovolskas, M. Franchini, A. J. Gallagher, C. J. Hansen, A. Korn, A. Kuvinskas, R. Lallement, L. Lombardo, F. Lucertini, L. Magrini, A. M. Matas Pinto, F. Matteucci, A. Mucciarelli, L. Sbordone, M. Spite , et al. (2 additional authors not shown)

Abstract: The MINCE (Measuring at Intermediate metallicity Neutron-Capture Elements) project aims to gather the abundances of neutron-capture elements but also of light elements and iron peak elements in a large sample of giant stars in this metallicity range. T The aim of this work is to study the chemical evolution of galactic sub-components recently identified (i.e. Gaia Sausage Enceladus (GSE), Sequoia)… ▽ More The MINCE (Measuring at Intermediate metallicity Neutron-Capture Elements) project aims to gather the abundances of neutron-capture elements but also of light elements and iron peak elements in a large sample of giant stars in this metallicity range. T The aim of this work is to study the chemical evolution of galactic sub-components recently identified (i.e. Gaia Sausage Enceladus (GSE), Sequoia). We used high signal-to-noise ratios, high-resolution spectra and standard 1D LTE spectrum synthesis to determine the detailed abundances. We could determine the abundances for up to 10 neutron-capture elements (Sr, Y, Zr, Ba, La, Ce, Pr, Nd, Sm and Eu) in 33 stars. The general trends of abundance ratios [n-capture element/Fe] versus [Fe/H] are in agreement with the results found in the literature. When our sample is divided in sub-groups depending on their kinematics, we found that the run of [Sr/Ba] vs [Ba/H] for the stars belonging to the GSE accretion event shows a tight anti-correlation. The results for the Sequoia stars, although based on a very limited sample, shows a [Sr/Ba] systematically higher than the [Sr/Ba] found in the GSE stars at a given [Ba/H] hinting at a different nucleosynthetic history. Stochastic chemical evolution models have been computed to understand the evolution of the GSE chemical composition of Sr and Ba. The first conclusions are that the GSE chemical evolution is similar to the evolution of a dwarf galaxy with galactic winds and inefficient star formation. Detailed abundances of neutron-capture elements have been measured in high-resolution, high signal-to-noise spectra of intermediate metal-poor stars, the metallicity range covered by the MINCE project. These abundances have been compared to detailed stochastic models of galactic chemical evolution. △ Less

Submitted 12 April, 2024; originally announced April 2024.

Comments: 18 pages, 15 figures

Cosmic Type Ia SN rate and constraints on SN Ia progenitors

Authors: P. A. Palicio, F. Matteucci, M. Della Valle, E. Spitoni

Abstract: Type Ia supernovae play a key role in the evolution of galaxies by polluting the interstellar medium with a fraction of iron peak elements larger than that released in the core collapse supernova events. Their light-curve, moreover, is widely used in cosmological studies as it constitutes a reliable distance indicator at extra-galactic scales. Among the mechanisms proposed to explain the Type Ia S… ▽ More Type Ia supernovae play a key role in the evolution of galaxies by polluting the interstellar medium with a fraction of iron peak elements larger than that released in the core collapse supernova events. Their light-curve, moreover, is widely used in cosmological studies as it constitutes a reliable distance indicator at extra-galactic scales. Among the mechanisms proposed to explain the Type Ia SNe, the single and double degenerate channels are thought to be the dominant ones, which imply a different distribution of time delays between the progenitor formation and the explosion. In this paper, we aim at determining the dominant mechanism by comparing a compilation of Type Ia SN rates with those computed from various cosmic star formation histories coupled with different delay time distribution functions, and evaluating the relative contributions of both channels. By using a least-squares fitting procedure, we model the observations of Type Ia SN rates assuming different combinations of three recent cosmic star formation rates and seven delay time distributions. The goodness of these fits are statistically quantified by the chi-squared test. For two of the three cosmic star formation rates, the single degenerate scenario provides the most accurate explanation for the observations, while a combination of 34% single degenerate and 66% double degenerate delay time distributions is more plausible for the remaining tested cosmic star formation rates. The wide double degenerate scenario mechanism slightly under-predicts the observations at redshift z>1, unless the cosmic SFR flattens in that regime. On the contrary, although the purely close double degenerate scenario can be ruled out, we cannot rule out a mixed scenario with single and double degenerate progenitors. △ Less

Submitted 19 June, 2024; v1 submitted 26 February, 2024; originally announced February 2024.

Comments: Accepted for publication in A&A

Stellar Population Astrophysics (SPA) with TNG, Fluorine abundances in seven open clusters

Authors: Shilpa Bijavara Seshashayana, Henrik Jönsson, Valentina D'Orazi, Govind Nandakumar, Ernesto Oliva, Angela Bragaglia, Nicoletta Sanna, Donatella Romano, Emanuele Spitoni, Amanda Karakas, Maria Lugaro, Livia Origlia

Abstract: The age, evolution, and chemical properties of the Galactic disk can be effectively ascertained using open clusters. Within the large program Stellar Populations Astrophysics at the Telescopio Nazionale Galileo, we specifically focused on stars in open clusters, to investigate various astrophysical topics, from the chemical content of very young systems to the abundance patterns of lesser studied… ▽ More The age, evolution, and chemical properties of the Galactic disk can be effectively ascertained using open clusters. Within the large program Stellar Populations Astrophysics at the Telescopio Nazionale Galileo, we specifically focused on stars in open clusters, to investigate various astrophysical topics, from the chemical content of very young systems to the abundance patterns of lesser studied intermediate-age and old open clusters. We investigate the astrophysically interesting element fluorine (F), which has an uncertain and intriguing cosmic origin. We also determine the abundance of cerium (Ce), as F abundance is expected to correlate with the s-process elements. High-resolution near-infrared spectra were obtained using the GIANO-B spectrograph. The Python version of Spectroscopy Made Easy (PySME), was used to derive atmospheric parameters and abundances. The stellar parameters were determined using OH, CN, and CO molecular lines along with Fe I lines. This paper presents the first F Galactic radial abundance gradient. Our results are also compared with literature estimates and with Galactic chemical evolution models that have been generated using different F production channels. Our results indicate a constant, solar pattern in the [F/Fe] ratios across clusters of different ages, supporting the latest findings that fluorine levels do not exhibit any secondary behavior for stars with solar or above-solar metallicity. By comparing our sample stars with the predictions of Galactic chemical evolution models, we came to the conclusion that both asymptotic giant branch stars and massive stars, including a fraction of fast rotators that increase with decreasing metallicity, are needed to explain the cosmic origin of F. △ Less

Submitted 14 February, 2024; originally announced February 2024.

Comments: 15 pages, 6 figures, 4 tables, Accepted by A&A

2D chemical evolution models II. Effects of multiple spiral arm patterns on O, Eu, Fe and Ba abundance gradients

Authors: E. Spitoni, G. Cescutti, A. Recio-Blanco, I. Minchev, E. Poggio, P. A. Palicio, F. Matteucci, S. Peirani, M. Barbillon, A. Vasini

Abstract: According to observations and numerical simulations, the Milky Way could exhibit several spiral arm modes with multiple pattern speeds, wherein the slower patterns are located at larger Galactocentric distances. Our aim is to quantify the effects of the spiral arms on the azimuthal variations of the chemical abundances for oxygen, iron and for the first time for neutron-capture elements (europium… ▽ More According to observations and numerical simulations, the Milky Way could exhibit several spiral arm modes with multiple pattern speeds, wherein the slower patterns are located at larger Galactocentric distances. Our aim is to quantify the effects of the spiral arms on the azimuthal variations of the chemical abundances for oxygen, iron and for the first time for neutron-capture elements (europium and barium) in the Galactic disc. We assume a model based on multiple spiral arm modes with different pattern speeds. The resulting model represents an updated version of previous 2D chemical evolution models. We apply new analytical prescriptions for the spiral arms in a 2D Galactic disc chemical evolution model, exploring the possibility that the spiral structure is formed by the overlap of chunks with different pattern speeds and spatial extent. The predicted azimuthal variations in abundance gradients are dependent on the considered chemical element. Elements synthesised on short time scales (i.e., oxygen and europium in this study) exhibit larger abundance fluctuations. In fact, for progenitors with short lifetimes, the chemical elements restored into the ISM perfectly trace the star formation perturbed by the passage of the spiral arms. The map of the star formation rate predicted by our chemical evolution model with multiple patterns of spiral arms presents arcs and arms compatible with those revealed by multiple tracers (young upper main sequence stars, Cepheids, and distribution of stars with low radial actions). Finally, our model predictions are in good agreement with the azimuthal variations that emerged from the analysis of Gaia DR3 GSP-Spec [M/H] abundance ratios, if at most recent times the pattern speeds match the Galactic rotational curve at all radii. △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: Accepted for publication in Astronomy and Astrophysics (A&A), 18 pages, 17 figures

Journal ref: A&A 680, A85 (2023)

Gaia Focused Product Release: Sources from Service Interface Function image analysis -- Half a million new sources in omega Centauri

Authors: Gaia Collaboration, K. Weingrill, A. Mints, J. Castañeda, Z. Kostrzewa-Rutkowska, M. Davidson, F. De Angeli, J. Hernández, F. Torra, M. Ramos-Lerate, C. Babusiaux, M. Biermann, C. Crowley, D. W. Evans, L. Lindegren, J. M. Martín-Fleitas, L. Palaversa, D. Ruz Mieres, K. Tisanić, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, F. Arenou, A. Barbier , et al. (378 additional authors not shown)

Abstract: Gaia's readout window strategy is challenged by very dense fields in the sky. Therefore, in addition to standard Gaia observations, full Sky Mapper (SM) images were recorded for nine selected regions in the sky. A new software pipeline exploits these Service Interface Function (SIF) images of crowded fields (CFs), making use of the availability of the full two-dimensional (2D) information. This ne… ▽ More Gaia's readout window strategy is challenged by very dense fields in the sky. Therefore, in addition to standard Gaia observations, full Sky Mapper (SM) images were recorded for nine selected regions in the sky. A new software pipeline exploits these Service Interface Function (SIF) images of crowded fields (CFs), making use of the availability of the full two-dimensional (2D) information. This new pipeline produced half a million additional Gaia sources in the region of the omega Centauri ($ω$ Cen) cluster, which are published with this Focused Product Release. We discuss the dedicated SIF CF data reduction pipeline, validate its data products, and introduce their Gaia archive table. Our aim is to improve the completeness of the {\it Gaia} source inventory in a very dense region in the sky, $ω$ Cen. An adapted version of {\it Gaia}'s Source Detection and Image Parameter Determination software located sources in the 2D SIF CF images. We validated the results by comparing them to the public {\it Gaia} DR3 catalogue and external Hubble Space Telescope data. With this Focused Product Release, 526\,587 new sources have been added to the {\it Gaia} catalogue in $ω$ Cen. Apart from positions and brightnesses, the additional catalogue contains parallaxes and proper motions, but no meaningful colour information. While SIF CF source parameters generally have a lower precision than nominal {\it Gaia} sources, in the cluster centre they increase the depth of the combined catalogue by three magnitudes and improve the source density by a factor of ten. This first SIF CF data publication already adds great value to the {\it Gaia} catalogue. It demonstrates what to expect for the fourth {\it Gaia} catalogue, which will contain additional sources for all nine SIF CF regions. △ Less

Submitted 8 November, 2023; v1 submitted 10 October, 2023; originally announced October 2023.

Journal ref: A&A 680, A35 (2023)

Gaia Focused Product Release: A catalogue of sources around quasars to search for strongly lensed quasars

Authors: Gaia Collaboration, A. Krone-Martins, C. Ducourant, L. Galluccio, L. Delchambre, I. Oreshina-Slezak, R. Teixeira, J. Braine, J. -F. Le Campion, F. Mignard, W. Roux, A. Blazere, L. Pegoraro, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, F. Arenou, C. Babusiaux, A. Barbier, M. Biermann, O. L. Creevey, D. W. Evans, L. Eyer, R. Guerra , et al. (376 additional authors not shown)

Abstract: Context. Strongly lensed quasars are fundamental sources for cosmology. The Gaia space mission covers the entire sky with the unprecedented resolution of $0.18$" in the optical, making it an ideal instrument to search for gravitational lenses down to the limiting magnitude of 21. Nevertheless, the previous Gaia Data Releases are known to be incomplete for small angular separations such as those ex… ▽ More Context. Strongly lensed quasars are fundamental sources for cosmology. The Gaia space mission covers the entire sky with the unprecedented resolution of $0.18$" in the optical, making it an ideal instrument to search for gravitational lenses down to the limiting magnitude of 21. Nevertheless, the previous Gaia Data Releases are known to be incomplete for small angular separations such as those expected for most lenses. Aims. We present the Data Processing and Analysis Consortium GravLens pipeline, which was built to analyse all Gaia detections around quasars and to cluster them into sources, thus producing a catalogue of secondary sources around each quasar. We analysed the resulting catalogue to produce scores that indicate source configurations that are compatible with strongly lensed quasars. Methods. GravLens uses the DBSCAN unsupervised clustering algorithm to detect sources around quasars. The resulting catalogue of multiplets is then analysed with several methods to identify potential gravitational lenses. We developed and applied an outlier scoring method, a comparison between the average BP and RP spectra of the components, and we also used an extremely randomised tree algorithm. These methods produce scores to identify the most probable configurations and to establish a list of lens candidates. Results. We analysed the environment of 3 760 032 quasars. A total of 4 760 920 sources, including the quasars, were found within 6" of the quasar positions. This list is given in the Gaia archive. In 87\% of cases, the quasar remains a single source, and in 501 385 cases neighbouring sources were detected. We propose a list of 381 lensed candidates, of which we identified 49 as the most promising. Beyond these candidates, the associate tables in this Focused Product Release allow the entire community to explore the unique Gaia data for strong lensing studies further. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Comments: 35 pages, 60 figures, accepted for publication by Astronomy and Astrophysics

Journal ref: A&A 685, A130 (2024)

Gaia Focused Product Release: Radial velocity time series of long-period variables

Authors: Gaia Collaboration, Gaia Collaboration, M. Trabucchi, N. Mowlavi, T. Lebzelter, I. Lecoeur-Taibi, M. Audard, L. Eyer, P. García-Lario, P. Gavras, B. Holl, G. Jevardat de Fombelle, K. Nienartowicz, L. Rimoldini, P. Sartoretti, R. Blomme, Y. Frémat, O. Marchal, Y. Damerdji, A. G. A. Brown, A. Guerrier, P. Panuzzo, D. Katz, G. M. Seabroke, K. Benson , et al. (382 additional authors not shown)

Abstract: The third Gaia Data Release (DR3) provided photometric time series of more than 2 million long-period variable (LPV) candidates. Anticipating the publication of full radial-velocity (RV) in DR4, this Focused Product Release (FPR) provides RV time series for a selection of LPVs with high-quality observations. We describe the production and content of the Gaia catalog of LPV RV time series, and the… ▽ More The third Gaia Data Release (DR3) provided photometric time series of more than 2 million long-period variable (LPV) candidates. Anticipating the publication of full radial-velocity (RV) in DR4, this Focused Product Release (FPR) provides RV time series for a selection of LPVs with high-quality observations. We describe the production and content of the Gaia catalog of LPV RV time series, and the methods used to compute variability parameters published in the Gaia FPR. Starting from the DR3 LPVs catalog, we applied filters to construct a sample of sources with high-quality RV measurements. We modeled their RV and photometric time series to derive their periods and amplitudes, and further refined the sample by requiring compatibility between the RV period and at least one of the $G$, $G_{\rm BP}$, or $G_{\rm RP}$ photometric periods. The catalog includes RV time series and variability parameters for 9\,614 sources in the magnitude range $6\lesssim G/{\rm mag}\lesssim 14$, including a flagged top-quality subsample of 6\,093 stars whose RV periods are fully compatible with the values derived from the $G$, $G_{\rm BP}$, and $G_{\rm RP}$ photometric time series. The RV time series contain a mean of 24 measurements per source taken unevenly over a duration of about three years. We identify the great most sources (88%) as genuine LPVs, with about half of them showing a pulsation period and the other half displaying a long secondary period. The remaining 12% consists of candidate ellipsoidal binaries. Quality checks against RVs available in the literature show excellent agreement. We provide illustrative examples and cautionary remarks. The publication of RV time series for almost 10\,000 LPVs constitutes, by far, the largest such database available to date in the literature. The availability of simultaneous photometric measurements gives a unique added value to the Gaia catalog (abridged) △ Less

Submitted 9 October, 2023; originally announced October 2023.

Comments: 36 pages, 38 figures

Galactic Archaeology with [Mg/Mn] versus [Al/Fe] abundance ratios -- Uncertainties and caveats

Authors: Arianna Vasini, Emanuele Spitoni, Francesca Matteucci

Abstract: The diagram depicting the abundance ratios [Mg/Mn] vs. [Al/Fe] has gained attention in recent literature as a valuable tool for exploring fundamental aspects of the evolution of the Milky Way and the Local Group. In particular, this combination of elements is supposed to be highly sensitive to the star formation history (SFH), unveiled by the imprints left on those abundances. Unfortunately, a com… ▽ More The diagram depicting the abundance ratios [Mg/Mn] vs. [Al/Fe] has gained attention in recent literature as a valuable tool for exploring fundamental aspects of the evolution of the Milky Way and the Local Group. In particular, this combination of elements is supposed to be highly sensitive to the star formation history (SFH), unveiled by the imprints left on those abundances. Unfortunately, a complete discussion on the uncertainties associated is still missing, making it difficult to know how reliable the associated results are. In this paper we analyze, by means of detailed chemical evolution models, the nuclear uncertainties of Mg, Al, Mn and Fe to show how different yields can affect the trends in the [Mg/Mn] vs. [Al/Fe] plane. In fact, if different yield assumptions produce conflicting results, then the [Mg/Mn] vs. [Al/Fe] diagram does not represent a strong diagnostic for the SFH of a galaxy. We discuss the results on the [Mg/Mn] vs. [Al/Fe] diagram, as predicted by several Milky Way (MW) and Large Magellanic Cloud (LMC) chemical evolution models adopting different nucleosynthesis prescriptions. The results show that the literature yields require some corrective factors to reproduce the APOGEE DR17 abundances of Mg, Al and Mn in the MW and that the same factors can also improve the results for the LMC. In particular, we show that by modifying the massive stars yields of Mg and Al the behaviour of the [Mg/Mn] vs. [Al/Fe] plot changes substantially. In conclusion, by changing the yields within their error bars, one obtains trends which differ strongly, making it difficult to draw any reliable conclusion on the SFH of galaxies. The proposed diagram is therefore uncertain from a theoretical point of view and it could represent a good diagnostic for SFH if the uncertainties on the nucleosynthesis of these elements (Mg, Mn, Al and Fe) could be reduced by future stellar calculations. △ Less

Submitted 6 October, 2023; originally announced October 2023.

Comments: 10 pages, 8 figures, Submitted to A&A

Oxygen, sulfur, and iron radial abundance gradients of classical Cepheids across the Galactic thin disk

Authors: R. da Silva, V. D'Orazi, M. Palla, G. Bono, V. F. Braga, M. Fabrizio, B. Lemasle, E. Spitoni, F. Matteucci, H. Jonsson, V. Kovtyukh, L. Magrini, M. Bergemann, M. Dall'Ora, I. Ferraro, G. Fiorentino, P. Francois, G. Iannicola, L. Inno, R. -P. Kudritzki, N. Matsunaga, M. Monelli, M. Nonino, C. Sneden, J. Storm , et al. (3 additional authors not shown)

Abstract: Classical Cepheids (CCs) are solid distance indicators and tracers of young stellar populations. Our aim is to provide iron, oxygen, and sulfur abundances for the largest and most homogeneous sample of Galactic CCs ever analyzed. The current sample covers a wide range in Galactocentric distances (RG), pulsation modes and periods. High-resolution and high S/N spectra collected with different spectr… ▽ More Classical Cepheids (CCs) are solid distance indicators and tracers of young stellar populations. Our aim is to provide iron, oxygen, and sulfur abundances for the largest and most homogeneous sample of Galactic CCs ever analyzed. The current sample covers a wide range in Galactocentric distances (RG), pulsation modes and periods. High-resolution and high S/N spectra collected with different spectrographs were adopted to estimate the atmospheric parameters. Individual distances are based on Gaia trigonometric parallaxes or on near-infrared Period-Luminosity relations. We found that Fe and alpha-element radial gradients based on CCs display a well-defined change in the slope for RG larger than 12 kpc. Radial gradients based on open clusters, covering a wide range in age, display similar trends, meaning that the flattening in the outer disk is an intrinsic feature of the radial gradients since it is independent of age. Empirical evidence indicates that the radial gradient for S is steeper than for Fe. The difference in the slope is a factor of two in the linear fit. We also found that S is, on average, under-abundant compared with O. We performed a detailed comparison with Galactic chemical evolution models and we found that a constant Star Formation Efficiency for RG larger than 12 kpc takes account for the flattening in both Fe and alpha-elements. To further constrain the impact that predicted S yields for massive stars have on radial gradients, we adopted a "toy model" and we found that the flattening in the outermost regions requires a decrease of a factor of four in the current S predictions. Sulfur photospheric abundances, compared with other alpha-elements, have the key advantage of being a volatile element. Therefore, stellar S abundances can be directly compared with nebular S abundances in external galaxies. △ Less

Submitted 2 August, 2023; originally announced August 2023.

Comments: 21 pages, 13 figures, 3 tables; to be published in the A&A journal

Journal ref: A&A 678, A195 (2023)

Chemical clocks and their time zones: understanding the [s/Mg]--age relation with birth radii

Authors: Bridget Ratcliffe, Ivan Minchev, Gabriele Cescutti, Emanuele Spitoni, Henrik Jönsson, Friedrich Anders, Anna Queiroz, Matthias Steinmetz

Abstract: The relative enrichment of s-process to $α$-elements ([s/$α$]) has been linked with age, providing a potentially useful avenue in exploring the Milky Way's chemical evolution. However, the age--[s/$α$] relationship is non-universal, with dependencies on metallicity and current location in the Galaxy. In this work, we examine these chemical clock tracers across birth radii (… ▽ More The relative enrichment of s-process to $α$-elements ([s/$α$]) has been linked with age, providing a potentially useful avenue in exploring the Milky Way's chemical evolution. However, the age--[s/$α$] relationship is non-universal, with dependencies on metallicity and current location in the Galaxy. In this work, we examine these chemical clock tracers across birth radii ($\rm \text{R}_\text{birth}$), recovering the inherent trends between the variables. We derive $\rm \text{R}_\text{birth}$ and explore the [s/$α$]--age--$\rm \text{R}_\text{birth}$ relationship for 36,652 APOGEE DR17 red giant and 24,467 GALAH DR3 main sequence turnoff and subgiant branch disk stars using [Ce/Mg], [Ba/Mg], and [Y/Mg]. We discover that the age--[s/Mg] relation is strongly dependent on birth location in the Milky Way, with stars born in the inner disk having the weakest correlation. This is congruent with the Galaxy's initially weak, negative [s/Mg] radial gradient, which becomes positive and steep with time. We show that the non-universal relations of chemical clocks is caused by their fundamental trends with $\rm \text{R}_\text{birth}$ over time, and suggest that the tight age--[s/Mg] relation obtained with solar-like stars is due to similar $\rm \text{R}_\text{birth}$ for a given age. Our results are put into context with a Galactic chemical evolution model, where we demonstrate the need for data-driven nucleosynthetic yields. △ Less

Submitted 20 July, 2023; originally announced July 2023.

Comments: submitted to MNRAS

A Bayesian chemical evolution model of the DustPedia Galaxy M74

Authors: Francesco Calura, Marco Palla, Laura Morselli, Emanuele Spitoni, Viviana Casasola, Kuldeep Verma, Andrea Enia, Massimo Meneghetti, Simone Bianchi, Francesca Pozzi, Carlotta Gruppioni

Abstract: We introduce a new, multi-zone chemical evolution model of the DustPedia galaxy M74, calibrated by means of MCMC methods. We take into account the observed stellar and gas density profiles and use Bayesian analysis to constrain two fundamental parameters characterising the gas accretion and star formation timescale, i.e. the infall timescale tau and the SF efficiency nu, respectively, as a funct… ▽ More We introduce a new, multi-zone chemical evolution model of the DustPedia galaxy M74, calibrated by means of MCMC methods. We take into account the observed stellar and gas density profiles and use Bayesian analysis to constrain two fundamental parameters characterising the gas accretion and star formation timescale, i.e. the infall timescale tau and the SF efficiency nu, respectively, as a function of galactocentric radius R. Our analysis supports an infall timescale increasing with R and a star formation efficiency decreasing with R, thus supporting an 'Inside-Out' formation for M74. For both tau and nu, we find a weaker radial dependence than in the Milky Way. We also investigate the dust content of M74, comparing the observed dust density profile with the results of our chemical evolution models. Various prescriptions have been considered for two key parameters, i.e. the typical dust accretion timescale and the mass of gas cleared out of the dust by a supernova remnant, regulating the dust growth and destruction rate, respectively. Two models with a different current balance between destruction and accretion, i.e. with equilibrium and dominion of accretion over destruction, can equally reproduce the observed dust profile of M74. This outlines the degeneracy between these parameters in shaping the interstellar dust content in galaxies. Our methods will be extended to more DustPedia galaxies to shed more light on the relative roles of dust production and destruction. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Comments: MNRAS, accepted for publication, 19 pages, 14 figures

Evolution of radioactive elements in the LMC: predictions for future gamma-ray surveys

Authors: Arianna Vasini, Francesca Matteucci, Emanuele Spitoni, Thomas Siegert

Abstract: Short-lived radionuclides, such as 26Al and 60Fe, are tracers of star formation. Therefore, their abundances can unravel the recent star formation history of the host galaxy. In view of future gamma-ray surveys, we predict the masses and fluxes of these two elements in the Large Magellanic Cloud (LMC) using new chemical evolution models. Our best model reproduces the abundance patterns of alpha/Fe… ▽ More Short-lived radionuclides, such as 26Al and 60Fe, are tracers of star formation. Therefore, their abundances can unravel the recent star formation history of the host galaxy. In view of future gamma-ray surveys, we predict the masses and fluxes of these two elements in the Large Magellanic Cloud (LMC) using new chemical evolution models. Our best model reproduces the abundance patterns of alpha/Fe ratios, the gas mass, the average metallicity, the present time supernova and nova rates observed in LMC. We show three main results: i) the best model for the LMC suggests a star formation rate very mild at the beginning with a recent burst, and a Salpeter-like initial mass function. ii) The predicted mass of 26Al is 0.33 M_Sun, 2/3 produced by massive stars and 1/3 by novae. iii) The predicted mass of 60Fe is 0.44 M_Sun, entirely produced by massive stars. This result suggests a larger fraction of 60Fe, at variance with the Milky Way. The explanation for this lies in the adopted initial mass function, that for the LMC contains more massive stars than for the Milky Way. These predictions can be useful for the COSI-SMEX mission planned for launch in 2027. The expected gamma-ray line fluxes for the 1.809 MeV line of 26Al and the 1.173 and 1.332 MeV lines of 60Fe are in the range of (0.2-2.7)x10^-6 ph cm^-2 s^-1 and (0.7-2.8)x10^-7 ph cm^-2 s^-1, respectively. This new instrument could have the sensitivity to detect the upper end of the predicted 1.8 MeV flux within its nominal two-year mission. △ Less

Submitted 8 May, 2023; originally announced May 2023.

Comments: 11 pages, 2 figures, 4 tables, accepted for publication by MNRAS

Origin of neutron capture elements with the Gaia-ESO survey: the evolution of s- and r-process elements across the Milky Way

Authors: Marta Molero, Laura Magrini, Francesca Matteucci, Donatella Romano, Marco Palla, Gabriele Cescutti, Carlos Viscasillas Vázquez, Emanuele Spitoni

Abstract: We study the abundance patterns and the radial gradients of s-process elements (Y, Zr, Ba, La and Ce), r-process elements (Eu) and mixed-process elements (Mo, Nd and Pr) in the Galactic thin disc by means of a detailed two-infall chemical evolution model for the Milky Way with state-of-the-art nucleosynthesis prescriptions. We consider r-process nucleosynthesis from merging neutron stars (MNS), ma… ▽ More We study the abundance patterns and the radial gradients of s-process elements (Y, Zr, Ba, La and Ce), r-process elements (Eu) and mixed-process elements (Mo, Nd and Pr) in the Galactic thin disc by means of a detailed two-infall chemical evolution model for the Milky Way with state-of-the-art nucleosynthesis prescriptions. We consider r-process nucleosynthesis from merging neutron stars (MNS), magneto-rotational supernovae (MR-SNe) and s-process synthesis from low- and intermediate- mass stars (LIMS) and rotating massive stars. The predictions of our model are compared with data from the sixth data release of the Gaia-ESO survey, from which we consider 62 open clusters with age > 0.1 Gyr and 1300 Milky Way disc field stars. We conclude that: i) the [Eu/Fe] vs. [Fe/H] is reproduced by both a prompt and a delayed source, but the quick source completely dominates the Eu production; ii) rotation in massive stars contribute substantially to the s-process elements of the first peak, but MNS and MR-SNe are necessary in order to reproduce the observations; iii) due to the adopted yields, our model overpredicts Pr and underpredicts Nd, while the [Mo/Fe] vs. [Fe/H] is nicely reproduced. For the radial gradients, we conclude that: i) our predicted slope of the [Fe/H] gradient is in agreement with the one observed in open clusters by Gaia-ESO and other high-resolution spectroscopic surveys. ii) The predicted slope of the [Eu/H] radial gradient is steeper than the observed one, independently on how quick the production of Eu is. We discuss the possible causes of this discrepancy in terms of both different Galaxy formation scenarios and stellar radial migration effects. iii) For all the elements belonging to the second s-process peak (Ba, La, Ce) as well as for Pr, we predict a plateau at low Galactocentric distances, which is probably due to the enhanced enrichment from LIMS in the inner regions. △ Less

Submitted 28 August, 2023; v1 submitted 13 April, 2023; originally announced April 2023.

Comments: 17 pages, 11 figures. Published on MNRAS

Journal ref: MNRAS, 523, 2974 (2023)

Analytic solution of Chemical Evolution Models with Type Ia SNe

Authors: P. A. Palicio, E. Spitoni, A. Recio-Blanco, F. Matteucci, S. Peirani, L. Greggio

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 c… ▽ More 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. △ Less

Submitted 22 July, 2023; v1 submitted 31 March, 2023; originally announced April 2023.

Comments: Accepted for publication in A&A. 22 pages, 16 figures. ChEAP code available at https://bitbucket.org/pedroap/cheap/src/master/

Journal ref: A&A 678, A61 (2023)

Detailed $α$ abundance trends in the inner Galactic bulge

Authors: N. Nieuwmunster, G. Nandakumar, E. Spitoni, N. Ryde, M. Schultheis, R. M. Rich, P. S. Barklem, O. Agertz, F. Renaud, F. Matteucci

Abstract: In this paper, we aim to derive high-precision alpha-element abundances using CRIRES high-resolution IR spectra of 72 cool M giants of the inner Galactic bulge. Silicon, magnesium, and calcium abundances were determined by fitting a synthetic spectrum for each star. We also incorporated recent theoretical data into our spectroscopic analysis (i.e. updated K-band line list, better broadening parame… ▽ More In this paper, we aim to derive high-precision alpha-element abundances using CRIRES high-resolution IR spectra of 72 cool M giants of the inner Galactic bulge. Silicon, magnesium, and calcium abundances were determined by fitting a synthetic spectrum for each star. We also incorporated recent theoretical data into our spectroscopic analysis (i.e. updated K-band line list, better broadening parameter estimation, non-local thermodynamic equilibrium (NLTE) corrections). We compare these inner bulge alpha abundance trends with those of solar neighbourhood stars observed with IGRINS using the same line list and analysis technique; we also compare our sample to APOGEE DR17 abundances for inner bulge stars. We investigate bulge membership using spectro-photometric distances and orbital simulations. We construct a chemical-evolution model that fits our metallicity distribution function (MDF) and our alpha-element trends. Among our 72 stars, we find four that are not bulge members. [Si/Fe] and [Mg/Fe] versus [Fe/H] trends show a typical thick disc alpha-element behaviour, except that we do not see any plateau at supersolar metallicities as seen in other works. The NLTE analysis lowers [Mg/Fe] typically by $\sim$0.1 dex, resulting in a noticeably lower trend of [Mg/Fe] versus [Fe/H]. The derived [Ca/Fe] versus [Fe/H] trend has a larger scatter than those for Si and Mg, but is in excellent agreement with local thin and thick disc trends. With our updated analysis, we constructed one of the most detailed studies of the alpha abundance trends of cool M giants in the inner Galactic bulge. We modelled these abundances by adopting a two-infall chemical-evolution model with two distinct gas-infall episodes with timescales of 0.4 Gyr and 2 Gyr, respectively. Based on a very meticulous spectral analysis, we have constructed detailed and precise chemical abundances of Mg, Si, and Ca for cool M giants. △ Less

Submitted 24 January, 2023; originally announced January 2023.

Comments: 19 pages, 11 figures

Journal ref: A&A 671, A94 (2023)

GalCEM I -- An Open-Source Detailed Isotopic Chemical Evolution Code

Authors: Eda Gjergo, Aleksei G. Sorokin, Anthony Ruth, Emanuele Spitoni, Francesca Matteucci, Xilong Fan, Jinning Liang, Marco Limongi, Yuta Yamazaki, Motohiko Kusakabe, Toshitaka Kajino

Abstract: This is the first of a series of papers that will introduce a user-friendly, detailed, and modular GALactic Chemical Evolution Model, GalCEM, that tracks isotope masses as a function of time in a given galaxy. The list of tracked isotopes automatically adapts to the complete set provided by the input yields. The present iteration of GalCEM tracks 86 elements broken down in 451 isotopes. The prescr… ▽ More This is the first of a series of papers that will introduce a user-friendly, detailed, and modular GALactic Chemical Evolution Model, GalCEM, that tracks isotope masses as a function of time in a given galaxy. The list of tracked isotopes automatically adapts to the complete set provided by the input yields. The present iteration of GalCEM tracks 86 elements broken down in 451 isotopes. The prescription includes massive stars, low-to-intermediate mass stars, and Type Ia supernovae as enrichment channels. We have developed a preprocessing tool that extracts multi-dimensional interpolation curves from the input yield tables. These interpolation curves improve the computation speeds of the full convolution integrals, which are computed for each isotope and for each enrichment channel. We map the integrand quantities onto consistent array grids in order to perform the numerical integration at each time step. The differential equation is solved with a fourth-order Runge-Kutta method. We constrain our analysis to the evolution of all the light and intermediate elements from carbon to zinc, and lithium. Our results are consistent up to the extremely metal poor regime with Galactic abundances. We provide tools to track the mass rate change of individual isotopes on a typical spiral galaxy with a final baryonic mass of $5\times 10^{10} M_{\odot}$. Future iterations of the work will extend to the full periodic table by including the enrichment from neutron-capture channels as well as spatially-dependent treatments of galaxy properties. GalCEM is publicly available at https://github.com/egjergo/GalCEM/ △ Less

Submitted 5 January, 2023; originally announced January 2023.

Comments: Accepted for publication by The Astrophysical Journal Supplement Series. 29 pages, 13 figures, 2 tables, 2 listings

MINCE I. Presentation of the project and of the first year sample

Authors: G. Cescutti, P. Bonifacio, E. Caffau, L. Monaco, M. Franchini, L. Lombardo, A. M. Matas Pinto, F. Lucertini, P. François, E. Spitoni, R. Lallement, L. Sbordone, A. Mucciarelli, M. Spite, C. J. Hansen, P. Di Marcantonio, A. Kučinskas, V. Dobrovolskas, A. J. Korn, M. Valentini, L. Magrini, S. Cristallo, F. Matteucci

Abstract: In recent years, Galactic archaeology has become a particularly vibrant field of astronomy, with its main focus set on the oldest stars of our Galaxy. In most cases, these stars have been identified as the most metal-poor. However, the struggle to find these ancient fossils has produced an important bias in the observations - in particular, the intermediate metal-poor stars (-2.5<[Fe/H]< -1.5) hav… ▽ More In recent years, Galactic archaeology has become a particularly vibrant field of astronomy, with its main focus set on the oldest stars of our Galaxy. In most cases, these stars have been identified as the most metal-poor. However, the struggle to find these ancient fossils has produced an important bias in the observations - in particular, the intermediate metal-poor stars (-2.5<[Fe/H]< -1.5) have been frequently overlooked. The missing information has consequences for the precise study of the chemical enrichment of our Galaxy, in particular for what concerns neutron-capture elements and it will be only partially covered by future multi-object spectroscopic surveys such as WEAVE and 4MOST. Measuring at Intermediate Metallicity Neutron Capture Elements (MINCE) is gathering the first high-quality spectra (high S/N ratio and high resolution) for several hundreds of bright and metal-poor stars, mainly located in our Galactic halo. We compiled our selection mainly on the basis of Gaia data and determined the stellar atmospheres of our sample and the chemical abundances of each star. In this paper, we present the first sample of 59 spectra of 46 stars. We measured the radial velocities and computed the Galactic orbits for all stars. We found that 8 stars belong to the thin disc, 15 to disrupted satellites, and the remaining cannot be associated to the mentioned structures, and we call them halo stars. For 33 of these stars, we provide abundances for the elements up to zinc. We also show the chemical evolution results for eleven chemical elements, based on recent models. Our observational strategy of using multiple telescopes and spectrographs to acquire high S/N and high-resolution spectra has proven to be very efficient since the present sample was acquired over only about one year of observations. Finally, our target selection strategy proved satisfactory for our purposes. △ Less

Submitted 11 November, 2022; originally announced November 2022.

Comments: 21 pages, 23 figures, online material. Accepted for publication in A&A

Journal ref: A&A 668, A168 (2022)

The Gaia-ESO survey: mapping the shape and evolution of the radial abundance gradients with open clusters

Authors: L. Magrini, C. Viscasillas Vazquez, L. Spina, S. Randich, D. Romano, E. Franciosini, A. Recio-Blanco, T. Nordlander, V. D'Orazi, M. Baratella, R. Smiljanic, M. L. L. Dantas, L. Pasquini, E. Spitoni, G. Casali, M. Van der Swaelmen, T. Bensby, E. Stonkute, S. Feltzing. G. G. Sacco, A. Bragaglia, E. Pancino, U. Heiter, K. Biazzo, G. Gilmore, M. Bergemann , et al. (5 additional authors not shown)

Abstract: The spatial distribution of elemental abundances and their time evolution are among the major constraints to disentangle the scenarios of formation and evolution of the Galaxy. We used the sample of open clusters available in the final release of the Gaia-ESO survey to trace the Galactic radial abundance and abundance to iron ratio gradients, and their time evolution. We selected member stars in 6… ▽ More The spatial distribution of elemental abundances and their time evolution are among the major constraints to disentangle the scenarios of formation and evolution of the Galaxy. We used the sample of open clusters available in the final release of the Gaia-ESO survey to trace the Galactic radial abundance and abundance to iron ratio gradients, and their time evolution. We selected member stars in 62 open clusters, with ages from 0.1 to about 7~Gyr, located in the Galactic thin disc at Galactocentric radii from about 6 to 21~kpc. We analysed the shape of the resulting [Fe/H] gradient, the average gradients [El/H] and [El/Fe] combining elements belonging to four different nucleosynthesis channels, and their individual abundance and abundance ratio gradients. We also investigated the time evolution of the gradients dividing open clusters in three age bins. The[Fe/H] gradient has a slope of -0.054 dex~kpc-1. We saw different behaviours for elements belonging to different channels. We found that the youngest clusters in the inner disc have lower metallicity than their older counterpart and they outline a flatter gradient. We considered some possible explanations, including the effects of gas inflow and migration. We suggested that it might be a bias introduced by the standard spectroscopic analysis producing lower metallicities in low gravity stars. To delineate the shape of the `true' gradient, we should limit our analysis to stars with low surface gravity logg>2.5 and xi<1.8 km~s-1. Based on this reduced sample, we can conclude that the gradient has minimally evolved over the time-frame outlined by the open clusters, indicating a slow and stationary formation of the thin disc in the latest Gyr. We found a secondary role of clusters' migration in shaping the gradient, with a more prominent role of migration for the oldest clusters. △ Less

Submitted 27 October, 2022; originally announced October 2022.

Comments: 25 pages, 14 figures and 4 tables in the main text, 3 figures and 7 tables in the Appendix. Accepted for publication in A&A

Journal ref: A&A 669, A119 (2023)

Spiral-like features in the disc revealed by Gaia DR3 radial actions

Authors: P. A. Palicio, A. Recio-Blanco, E. Poggio, T. Antoja, P. J. McMillan, E. Spitoni

Abstract: We aim to explore the actions of the new Gaia DR3 astrometry to find structures in the Galactic disc. We compute the actions and the orbital parameters of the Gaia DR3 stars with full astrometry and velocities assuming an axisymmetric model for the Milky Way. Using Gaia DR3 photometric data, we select a subset of giants stars with better astrometry as control sample. The maps of the percentiles of… ▽ More We aim to explore the actions of the new Gaia DR3 astrometry to find structures in the Galactic disc. We compute the actions and the orbital parameters of the Gaia DR3 stars with full astrometry and velocities assuming an axisymmetric model for the Milky Way. Using Gaia DR3 photometric data, we select a subset of giants stars with better astrometry as control sample. The maps of the percentiles of the radial action J_R reveal spiral-like shape structures. We find a high J_R region centered at R~10.5 kpc of 1 kpc width, as well as three arc-shape regions dominated by circular orbits at inner radii. We also identify the spiral arms in the overdensities of the giant population. We find a good agreement with the literature in the innermost region for the Scutum-Sagittarius spiral arms. At larger radii, the low J_R structure tracks the Local arm at negative X, while for the Perseus arm the agreement is restricted to the X<2 kpc region, with a displacement with respect to the literature at more negative longitudes. We detect a high J_R area at a Galactocentric radii of ~10.5 kpc, consistent with some estimations of the Outer Lindblad Resonance location. We conclude that the pattern in the dynamics of the old stars is consistent in several places with spatial distribution of the spiral arms traced by young populations, with small potential contributions from the moving groups. △ Less

Submitted 27 January, 2023; v1 submitted 20 September, 2022; originally announced September 2022.

Comments: 11 pages, 9 figures. Accepted for publication in A&A. Acceptance date: 04 December 2022

Journal ref: A&A 670, L7 (2023)

The cerium content of the Milky Way as revealed by Gaia DR3 GSP-Spec abundances

Authors: G. Contursi, P. de Laverny, A. Recio-Blanco, E. Spitoni, P. A. Palicio, E. Poggio, V. Grisoni, G. Cescutti, F. Matteucci, L. Spina, M. A. Alvarez, G. Kordopatis, C. Ordenovic, I. Oreshina-Slezak, H. Zhao

Abstract: The recent Gaia Third Data Release contains a homogeneous analysis of millions of high-quality Radial Velocity Spectrometer (RVS) stellar spectra by the GSP-Spec module. This led to the estimation of millions of individual chemical abundances and allows us to chemically map the Milky Way. Among the published GSP-Spec abundances, three heavy-elements produced by neutron-captures in stellar interior… ▽ More The recent Gaia Third Data Release contains a homogeneous analysis of millions of high-quality Radial Velocity Spectrometer (RVS) stellar spectra by the GSP-Spec module. This led to the estimation of millions of individual chemical abundances and allows us to chemically map the Milky Way. Among the published GSP-Spec abundances, three heavy-elements produced by neutron-captures in stellar interiors can be found: Ce, Zr and Nd. We use a sample of about 30,000 LTE Ce abundances, selected after applying different combinations of GSP-Spec flags. Thanks to the Gaia DR3 astrometric data and radial velocities, we explore the cerium content in the Milky Way and, in particular, in its halo and disc components. The high quality of the Ce GSP-Spec abundances is quantified thanks to literature comparisons. We found a rather flat [Ce/Fe] versus [M/H] trend. We also found a flat radial gradient in the disc derived from field stars and, independently, from about 50 open clusters, in agreement with previous studies. The [Ce/Fe] vertical gradient has also been estimated. We also report an increasing [Ce/Ca] vs [Ca/H] in the disc, illustrating the late contribution of AGB with respect to SN II. Our cerium abundances in the disc, including the young massive population, are well reproduced by a new three-infall chemical evolution model. Among the halo population, the M 4 globular cluster is found to be enriched in cerium. Moreover, eleven stars with cerium abundances belonging to the Thamnos, Helmi Stream and Gaia-Sausage-Enceladus accreted systems were identified from chemo-dynamical diagnostics. We found that the Helmi Stream could be slightly underabundant in cerium, compared to the two other systems. This work illustrates the high quality of the GSP-Spec chemical abundances, that significantly contributes to unveil the heavy elements evolution history of the Milky Way. △ Less

Submitted 15 December, 2022; v1 submitted 12 July, 2022; originally announced July 2022.

Comments: 17 pages, 12 figures, accepted by A&A

The chemical signature of the Galactic spiral arms revealed by Gaia DR3

Authors: E. Poggio, A. Recio-Blanco, P. A. Palicio, P. Re Fiorentin, P. de Laverny, R. Drimmel, G. Kordopatis, M. G. Lattanzi, M. Schultheis, A. Spagna, E. Spitoni

Abstract: Taking advantage of the recent Gaia Data Release 3 (DR3), we map chemical inhomogeneities in the Milky Way's disc out to a distance of $\sim$ 4 kpc from the Sun, using different samples of bright giant stars. The samples are selected using effective temperatures and surface gravities from the GSP-Spec module, and are expected to trace stellar populations of different typical age. The cool (old) gi… ▽ More Taking advantage of the recent Gaia Data Release 3 (DR3), we map chemical inhomogeneities in the Milky Way's disc out to a distance of $\sim$ 4 kpc from the Sun, using different samples of bright giant stars. The samples are selected using effective temperatures and surface gravities from the GSP-Spec module, and are expected to trace stellar populations of different typical age. The cool (old) giants exhibit a relatively smooth radial metallicity gradient with an azimuthal dependence. Binning in Galactic azimuth $φ$, the slope gradually varies from $d$[M/H]$/dR \sim -0.054$ dex kpc$^{-1}$ at $φ\sim -20^{\circ}$ to $\sim -0.035$ dex kpc$^{-1}$ at $φ\sim 20^{\circ}$. On the other hand, the relatively hotter (and younger) stars present remarkable inhomogeneities, apparent as three (possibly four) metal-rich elongated features in correspondence of the spiral arms' locations in the Galactic disc. When projected onto Galactic radius, those features manifest themselves as statistically significant bumps on top of the observed radial metallicity gradients with amplitudes up to $ \sim 0.05-0.1$ dex, making the assumption of a linear radial decrease not applicable to this sample. The strong correlation between the spiral structure of the Galaxy and the observed chemical pattern in the young sample indicates that the spiral arms might be at the origin for the detected chemical inhomogeneities. In this scenario, the spiral arms would leave in the younger stars a strong signature, which progressively disappears when cooler (and older) giants are considered. △ Less

Submitted 1 September, 2022; v1 submitted 29 June, 2022; originally announced June 2022.

Comments: 9 pages, 6 figures. Accepted for publication in A&A Letters

Journal ref: A&A 666, L4 (2022)

Beyond the two-infall model I. Indications for a recent gas infall with Gaia DR3 chemical abundances

Authors: E. Spitoni, A. Recio-Blanco, P. de Laverny, P. A. Palicio, G. Kordopatis, M. Schultheis, G. Contursi, E. Poggio, D. Romano, F. Matteucci

Abstract: The recent Gaia Data Release 3 (DR3) represents an unparalleled revolution in Galactic Archaeology, providing us with numerous radial velocities chemical abundances for millions of stars, with all-sky coverage. We present a new chemical evolution model for the Galactic disc components (high- and low- $α$ sequence stars) designed to reproduce the new abundance ratios provided by the GSP-spec module… ▽ More The recent Gaia Data Release 3 (DR3) represents an unparalleled revolution in Galactic Archaeology, providing us with numerous radial velocities chemical abundances for millions of stars, with all-sky coverage. We present a new chemical evolution model for the Galactic disc components (high- and low- $α$ sequence stars) designed to reproduce the new abundance ratios provided by the GSP-spec module for the Gaia DR3 and also constrained by the detailed star formation histories for both the thick and thin disc stars inferred from previous Gaia releases. Sophisticated modeling based on previous Gaia releases has found evidence for narrow episodes of enhanced SF inferred in recent time. Additionally, Gaia DR3 highlighted the presence of young (massive) low-$α$ disc stars which show evidence of a recent chemical impoverishment in several elements. Hence, in this study, we compare Gaia DR3 chemical abundances with the predictions of a three-infall chemical evolution model for the high- and low-$α$ components. The proposed three-infall chemical evolution model nicely reproduces the main features of the abundance ratio [X/Fe] versus [M/H] (X=Mg, Si, Ca, Ti, $α$) of Gaia DR3 stars in different age bins for the considered $α$ elements. Moreover, the most recent gas infall - which started $\sim$ 2.7 Gyr ago - allows us to predict well the Gaia DR3 young population which has experienced a recent chemical impoverishment. △ Less

Submitted 13 December, 2022; v1 submitted 24 June, 2022; originally announced June 2022.

Comments: Accepted for publication in Astronomy and Astrophysics (A&A), 17 pages, 16 figures

Journal ref: A&A 670, A109 (2023)

Stellar ages, masses, extinctions and orbital parameters based on spectroscopic parameters of Gaia DR3

Authors: G. Kordopatis, M. Schultheis, P. J. McMillan, P. A. Palicio, P. de Laverny, A. Recio-Blanco, O. Creevey, M. A. Álvarez, R. Andrae, E. Poggio, E. Spitoni, G. Contursi, H. Zhao, I. Oreshina-Slezak, C. Ordenovic, A. Bijaoui

Abstract: Gaia DR3 provides radial velocities for 33 million stars and spectroscopically derived atmospheric parameters for more than five million targets. When combined with the astrometric data, these allow us to derive orbital and stellar parameters that are key in order to understand the stellar populations of the Milky Way and perform galactic archaeology. We use the calibrated atmospheric parameters,… ▽ More Gaia DR3 provides radial velocities for 33 million stars and spectroscopically derived atmospheric parameters for more than five million targets. When combined with the astrometric data, these allow us to derive orbital and stellar parameters that are key in order to understand the stellar populations of the Milky Way and perform galactic archaeology. We use the calibrated atmospheric parameters, 2MASS and Gaia-EDR3 photometry, and parallax-based distances to compute, via an isochrone fitting method, the ages, initial stellar masses and reddenings for the stars with spectroscopic parameters. We also derive the orbits (actions, eccentricities, apocentre, pericentre and Zmax) for all of the stars with measured radial velocities and astrometry, adopting two sets of line-of-sight distances from the literature and an axisymmetric potential of the Galaxy. Comparisons with reference catalogues of field and cluster stars suggest that reliable ages are obtained for stars younger than 9-10Gyr when the estimated relative age uncertainty is <50%. For older stars, ages tend to be under-estimated. The most reliable stellar type for age determination are turn-off stars, even when the input atmospheric parameters have large uncertainties. Ages for giants and main-sequence stars are retrieved with uncertainties of ~2Gyr when extinction towards the star's line-of sight is smaller than A_V<2.5mag. The full catalogue is made publicly available to be downloaded. With it, the full chemo-dynamical properties of the extended Solar neighbourhood unfold, and allow us to better identify the properties of the spiral arms, to parameterise the dynamical heating of the disc, or to thoroughly study the chemical enrichment of the Milky Way. △ Less

Submitted 16 June, 2022; originally announced June 2022.

Comments: The Table (in .hdf format) and the readme file can be downloaded here: https://ftp.oca.eu/pub/gkordo/GDR3/ The parameters computed in this paper have been prepared in the context of Gaia's Performance verification paper concerning the Chemical cartography of the Milky Way (Gaia collaboration, Recio-Blanco et al, 2022). 22 pages, 20 figures, 3 appendixes. Submitted to A&A

Journal ref: A&A 669, A104 (2023)

Gaia Data Release 3: Analysis of RVS spectra using the General Stellar Parametriser from spectroscopy

Authors: A. Recio-Blanco, P. de Laverny, P. A. Palicio, G. Kordopatis, M. A. Álvarez, M. Schultheis, G. Contursi, H. Zhao, G. Torralba Elipe, C. Ordenovic, M. Manteiga, C. Dafonte, I. Oreshina-Slezak, A. Bijaoui, Y. Fremat, G. Seabroke, F. Pailler, E. Spitoni, E. Poggio, O. L. Creevey, A. Abreu Aramburu, S. Accart, R. Andrae, C. A. L. Bailer-Jones, I. Bellas-Velidis , et al. (55 additional authors not shown)

Abstract: The chemo-physical parametrisation of stellar spectra is essential for understanding the nature and evolution of stars and of Galactic stellar populations. Gaia DR3 contains the parametrisation of RVS data performed by the General Stellar Parametriser-spectroscopy, module. Here we describe the parametrisation of the first 34 months of RVS observations. GSP-spec estimates the chemo-physical paramet… ▽ More The chemo-physical parametrisation of stellar spectra is essential for understanding the nature and evolution of stars and of Galactic stellar populations. Gaia DR3 contains the parametrisation of RVS data performed by the General Stellar Parametriser-spectroscopy, module. Here we describe the parametrisation of the first 34 months of RVS observations. GSP-spec estimates the chemo-physical parameters from combined RVS spectra of single stars. The main analysis workflow described here, MatisseGauguin, is based on projection and optimisation methods and provides the stellar atmospheric parameters; the individual chemical abundances of N, Mg, Si, S, Ca, Ti, Cr, FeI, FeII, Ni, Zr, Ce and Nd; the differential equivalent width of a cyanogen line; and the parameters of a DIB feature. Another workflow, based on an artificial neural network, provides a second set of atmospheric parameters that are useful for classification control. We implement a detailed quality flag chain considering different error sources. With about 5.6 million stars, the Gaia DR3 GSP-spec all-sky catalogue is the largest compilation of stellar chemo-physical parameters ever published and the first one from space data. Internal and external biases have been studied taking into account the implemented flags. In some cases, simple calibrations with low degree polynomials are suggested. The homogeneity and quality of the estimated parameters enables chemo-dynamical studies of Galactic stellar populations, interstellar extinction studies from individual spectra, and clear constraints on stellar evolution models. We highly recommend that users adopt the provided quality flags for scientific exploitation . The Gaia DR3 GSP-spec catalogue is a major step in the scientific exploration of Milky Way stellar populations, confirming the Gaia promise of a new Galactic vision (abridged). △ Less

Submitted 11 June, 2022; originally announced June 2022.

Comments: Astronomy and Astrophysics (accepted, in press)

Journal ref: A&A 674, A29 (2023)

Disc dichotomy signature in the vertical distribution of [Mg/Fe] and the delayed gas infall scenario

Authors: E. Spitoni, V. Aguirre Børsen-Koch, K. Verma, A. Stokholm

Abstract: The analysis of the APOGEE data suggests the existence of a clear distinction between two sequences of disc stars in the [$α$/Fe] vs. [Fe/H] abundance ratio space. We aim to test if the two-infall chemical evolution models designed to reproduce these two sequences in the solar neighbourhood are also capable to predict the disc bimodality observed in the vertical distribution of [Mg/Fe] in APOGEE D… ▽ More The analysis of the APOGEE data suggests the existence of a clear distinction between two sequences of disc stars in the [$α$/Fe] vs. [Fe/H] abundance ratio space. We aim to test if the two-infall chemical evolution models designed to reproduce these two sequences in the solar neighbourhood are also capable to predict the disc bimodality observed in the vertical distribution of [Mg/Fe] in APOGEE DR16 data. Along with the predicted chemical composition of SSPs born at different Galactic times in the solar vicinity, we provide their maximum vertical height |zmax| above the Galactic plane computed assuming the relation between the vertical action and stellar age in thin disc stars. The predicted vertical distribution of the [Mg/Fe] abundance ratio is in agreement with the one observed combining the APOGEE DR16 data and the astroNN catalogue (stellar ages, orbital parameters) for stars younger than 8 Gyr (only low-$α$ sequence stars). Including the high-$α$ disc component, the dichotomy in the vertical [Mg/Fe] abundance distribution is reproduced considering the observational cut in the Galactic height of |z| < 2 kpc. However, our model predicts a too flat growth of the |zmax| as a function of [Mg/Fe] for high-$α$ objects in contrast with the median values from APOGEE data. Possible explanations for such a tension are: i) the data sample with |z| < 2 kpc is more likely contaminated by halo stars, causing the median values to be kinematically hotter, ii) external perturbations such as minor mergers could have heated up the disc, and the heating of the orbits cannot be modelled by only scattering processes. Assuming for the data a disc dissection based on chemistry, the observed |zmax| distributions for high-$α$ and low-$α$ sequences are in good agreement with our model predictions if we consider in the calculation the errors in the vertical action estimates. △ Less

Submitted 15 April, 2022; originally announced April 2022.

Comments: Accepted for publication in Astronomy and Astrophysics (A&A), 17 pages, 15 figures

Journal ref: A&A 663, A174 (2022)

Chemical evolution of $^{26}$Al and $^{60}$Fe in the Milky Way

Authors: A. Vasini, F. Matteucci, E. Spitoni

Abstract: We present theoretical mass estimates of $^{26}$Al and $^{60}$Fe throughout the Galaxy, performed with a numerical chemical evolution model including detailed nucleosynthesis prescriptions for stable and radioactive nuclides. We compared the results for several sets of stellar yields taken from the literature, for massive, low and intermediate mass stars, nova systems (only for $^{26}$Al) and supe… ▽ More We present theoretical mass estimates of $^{26}$Al and $^{60}$Fe throughout the Galaxy, performed with a numerical chemical evolution model including detailed nucleosynthesis prescriptions for stable and radioactive nuclides. We compared the results for several sets of stellar yields taken from the literature, for massive, low and intermediate mass stars, nova systems (only for $^{26}$Al) and supernovae Type Ia.We then computed the total masses of $^{26}$Al and $^{60}$Fe in the Galaxy. We studied the bulge and the disc of the Galaxy in a galactocentric radius range 0-22 kpc. We assumed that the bulge region (within 2 kpc) evolved quickly suffering a strong star formation burst, while the disc formed more slowly and inside-out. We compared our results with the $^{26}$Al mass observed by the $γ$-ray surveys COMPTEL and INTEGRAL to select the best model. Concerning $^{60}$Fe, we do not have any observed mass value so we just performed a theoretical prediction for future observations. In conclusion, low, intermediate mass stars and Type Ia supernovae contribute negligibly to the two isotopes, while massive stars are the dominant source. The nova contribution is, however, necessary to reproduce the observations of $^{26}$Al. Our best model predicts $2.12$ M$_{\odot}$ of $^{26}$Al, in agreement with observations, while for $^{60}$Fe our best mass estimate is $\sim 1.05$ M$_{\odot}$. We also predicted the present injection rate of $^{26}$Al and $^{60}$Fe in the Galaxy and compared it with previous results, and we found a larger present time injection rate along the disc. △ Less

Submitted 23 May, 2022; v1 submitted 1 April, 2022; originally announced April 2022.

Comments: accepted,10 pages, 6 figures

The G-dwarf distribution in star-forming galaxies: a tug-of-war between infall and outflow

Authors: E. Spitoni, F. Calura, V. Silva Aguirre, R. Gilli

Abstract: In the past, the cumulative metallicity distribution function (CMDF) turned out as a useful tool to constrain the accretion history of various components of the Milky Way. In this Letter, by means of analytical, leaky-box chemical evolution models (i.e. including both infall and galactic outflows) we study the CMDF of local star-forming galaxies that follow two fundamental empirical scaling relati… ▽ More In the past, the cumulative metallicity distribution function (CMDF) turned out as a useful tool to constrain the accretion history of various components of the Milky Way. In this Letter, by means of analytical, leaky-box chemical evolution models (i.e. including both infall and galactic outflows) we study the CMDF of local star-forming galaxies that follow two fundamental empirical scaling relations, namely the mass-metallicity and main sequence relations. Our analysis shows that galactic winds, which are dominant mostly in low-mass systems, play a fundamental role in shaping this function and, in particular, in determining its steepness and curvature. We show that the CMDF of low-mass (M$_{\star}$/M$_{\odot} \le 10^{9.5}$) and high-mass (M$_{\star}$/M$_{\odot}$>10$^{10.5}$) galaxies deviate substantially from the results of a 'closed-box' model, as the evolution of the former (latter) systems is mostly dominated by outflows (infall). In the context of galactic downsizing, we show that downward-concave CMDFs (associated with systems with extremely small infall timescales and with very strong winds) are more frequent in low-mass galaxies, which include larger fractions of young systems and present more substantial deviations from equilibrium between gas accretion and reprocessing (either via star formation or winds). △ Less

Submitted 6 April, 2021; v1 submitted 19 March, 2021; originally announced March 2021.

Comments: Accepted for publication in Astronomy and Astrophysics Letters (A&A Letters), 6 pages, 4 figures

Journal ref: A&A 648, L5 (2021)

APOGEE DR16: a multi-zone chemical evolution model for the Galactic disc based on MCMC methods

Authors: E. Spitoni, K. Verma, V. Silva Aguirre, F. Vincenzo, F. Matteucci, B. Vaičekauskaitė, M. Palla, V. Grisoni, F. Calura

Abstract: The analysis of the APOGEE DR16 data suggests the existence of a clear distinction between two sequences of disc stars at different Galactocentric distances in the [$α$/Fe] vs. [Fe/H] abundance ratio space: the so-called high-$α$ sequence, classically associated to an old population of stars in the thick disc, and the low-$α$ sequence, which mostly comprises relatively young stars in the thin disc… ▽ More The analysis of the APOGEE DR16 data suggests the existence of a clear distinction between two sequences of disc stars at different Galactocentric distances in the [$α$/Fe] vs. [Fe/H] abundance ratio space: the so-called high-$α$ sequence, classically associated to an old population of stars in the thick disc, and the low-$α$ sequence, which mostly comprises relatively young stars in the thin disc. We perform a Bayesian analysis based on a Markov Chain Monte Carlo method to constrain a multi-zone two-infall chemical evolution model designed for regions at different Galactocentric distances using measured chemical abundances from the APOGEE DR16 sample. An inside-out formation of the Galaxy disc naturally emerges from the best fit of our two-infall chemical-evolution model to APOGEE-DR16: inner Galactic regions are assembled on shorter time-scales compared to the external ones. In the outer disc (with radii $R>6$ kpc), the chemical dilution due to a late accretion event of gas with primordial chemical composition is the main driver of the [Mg/Fe] vs. [Fe/H] abundance pattern in the low-$α$ sequence. In the inner disc, in the framework of the two-infall model, we confirm the presence of an enriched gas infall in the low-$α$ phase as suggested by chemo-dynamical models. Our Bayesian analysis of the recent APOGEE DR16 data suggests a significant delay time, ranging from $\sim$3.0 to 4.7 Gyr, between the first and second gas infall events for all the analyzed Galactocentric regions. Our results propose a clear interpretation of the [Mg/Fe] vs. [Fe/H] relations along the Galactic discs. The signatures of a delayed gas-rich merger which gives rise to a hiatus in the star formation history of the Galaxy are impressed in the [Mg/Fe] vs. [Fe/H] relation, determining how the low-$α$ stars are distributed in the abundance space at different Galactocentric distances. △ Less

Submitted 21 January, 2021; originally announced January 2021.

Comments: Accepted for publication in Astronomy and Astrophysics (A&A), 16 pages, 19 figures

Journal ref: A&A 647, A73 (2021)

Chemical evolution of the Milky Way: constraints on the formation of the thick and thin discs

Authors: Marco Palla, Francesca Matteucci, Emanuele Spitoni, Fiorenzo Vincenzo, Valeria Grisoni

Abstract: We study the evolution of Milky Way thick and thin discs in the light of the most recent observational data. In particular, we analyze abundance gradients of O, N, Fe and Mg along the thin disc as well as the [Mg/Fe] vs. [Fe/H] relations and the metallicity distribution functions at different Galactocentric distances. We run several models starting from the two-infall paradigm, assuming that the t… ▽ More We study the evolution of Milky Way thick and thin discs in the light of the most recent observational data. In particular, we analyze abundance gradients of O, N, Fe and Mg along the thin disc as well as the [Mg/Fe] vs. [Fe/H] relations and the metallicity distribution functions at different Galactocentric distances. We run several models starting from the two-infall paradigm, assuming that the thick and thin discs formed by means of two different infall episodes, and we explore several physical parameters, such as radial gas flows, variable efficiency of star formation, different times for the maximum infall onto the disc, different distributions of the total surface mass density of the thick disc and enriched gas infall. Our best model suggests that radial gas flows and variable efficiency of star formation should be acting together with the inside-out mechanism for the thin disc formation. The timescale for maximum infall onto the thin disc, which determines the gap between the formation of the two discs, should be $t_{max}\simeq 3.25$ Gyr. The thick disc should have an exponential, small scale length density profile and gas infall on the inner thin disc should be enriched. We compute also the evolution of Gaia-Enceladus system and study the effects of possible interactions with the thick and thin discs. We conclude that the gas lost by Enceladus or even part of it could have been responsible for the formation of the thick disc but not the thin disc. △ Less

Submitted 17 August, 2020; originally announced August 2020.

Comments: 18 pages, 12 figures, 5 tables. Accepted for publication in MNRAS

Fluorine in the solar neighborhood: modelling the Galactic thick and thin discs

Authors: Valeria Grisoni, Donatella Romano, Emanuele Spitoni, Francesca Matteucci, Nils Ryde, Henrik Jönsson

Abstract: We investigate the evolution of the abundance of fluorine in the Milky Way thick and thin discs by means of detailed chemical evolution models compared with recent observational data. The chemical evolution models adopted here have already been shown to fit the observed abundance patterns of CNO and $α$-elements as well as the metallicity distribution functions for the Galactic thick and thin disc… ▽ More We investigate the evolution of the abundance of fluorine in the Milky Way thick and thin discs by means of detailed chemical evolution models compared with recent observational data. The chemical evolution models adopted here have already been shown to fit the observed abundance patterns of CNO and $α$-elements as well as the metallicity distribution functions for the Galactic thick and thin disc stars. We apply them here to the study of the origin and evolution of fluorine, which is still a matter of debate. First, we study the importance of the various sites proposed for the production of fluorine. Then, we apply the reference models to follow the evolution of the two different Galactic components. We conclude that rotating massive stars are important producers of F and they can set a plateau in F abundance below [Fe/H]=-0.5 dex, though its existence for [Fe/H]<-1 has yet to be confimed by extensive observations of halo stars. In order to reproduce the F abundance increase in the discs at late times, instead, a contribution from lower mass stars - single asymptotic giant branch stars and/or novae - is required. The dichotomy between the thick and thin discs is more evident in the [F/O] vs. [O/H] plot than in the [F/Fe] vs. [Fe/H] one, and we confirm that the thick disc has evolved much faster than the thin disc, in agreement with findings from the abundance patterns of other chemical elements. △ Less

Submitted 3 August, 2020; originally announced August 2020.

Comments: 8 pages, 2 figures. Accepted for publication in MNRAS

High-precision abundances of elements in solar-type stars. Evidence of two distinct sequences in abundance-age relations

Authors: P. E. Nissen, J. Christensen-Dalsgaard, J. R. Mosumgaard, V. Silva Aguirre, E. Spitoni, K. Verma.

Abstract: HARPS spectra with signal-to-noise ratios S/N > 600 at 6000 A were analysed with MARCS model atmospheres to obtain 1D LTE abundances of C, O, Na, Mg, Al, Si, Ca, Ti, Cr, Fe, Ni, Sr, and Y for 72 nearby solar-type stars with metallicities in the range -0.3 < [Fe/H] < +0.3 and ASTEC stellar models were used to determine stellar ages from effective temperatures, luminosities obtained via Gaia DR2 par… ▽ More HARPS spectra with signal-to-noise ratios S/N > 600 at 6000 A were analysed with MARCS model atmospheres to obtain 1D LTE abundances of C, O, Na, Mg, Al, Si, Ca, Ti, Cr, Fe, Ni, Sr, and Y for 72 nearby solar-type stars with metallicities in the range -0.3 < [Fe/H] < +0.3 and ASTEC stellar models were used to determine stellar ages from effective temperatures, luminosities obtained via Gaia DR2 parallaxes, and heavy element abundances. The resulting age-metallicity distribution appears to consist of two distinct populations: a sequence of old stars with a steep rise of [Fe/H] to about +0.3 dex at an age of ~7 Gyr and a younger sequence with [Fe/H] increasing from -0.3 dex to about +0.2 dex over the last six Gyr. Furthermore, the trends of several abundance ratios, [O/Fe], [Na/Fe], [Ca/Fe], and [Ni/Fe], as a function of stellar age split into two corresponding sequences. The [Y/Mg]-age relation, on the other hand, shows no offset between the two age sequences and has no significant dependence on [Fe/H], but the components of a visual binary star, zeta Reticuli, have a large and puzzling deviation. The split of the age-metallicity distribution into two sequences may be interpreted as evidence of two episodes of accretion of gas onto the Galactic disk with a quenching of star formation in between. Some of the [X/Fe]-age relations support this scenario but other relations are not so easy to explain, which calls for a deeper study of systematic errors in the derived abundances as a function of [Fe/H], in particular 3D non-LTE effects. △ Less

Submitted 10 June, 2020; originally announced June 2020.

Comments: 15 pages, 19 figures. Accepted for publication in A&A

Journal ref: A&A 640, A81 (2020)

On the variation of carbon abundance in galaxies and its implications

Authors: Donatella Romano, Mariagrazia Franchini, Valeria Grisoni, Emanuele Spitoni, Francesca Matteucci, Carlo Morossi

Abstract: The trends of chemical abundances and abundance ratios observed in stars of different ages, kinematics, and metallicities bear the imprints of several physical processes that concur to shape the host galaxy properties. By inspecting these trends, we get precious information on stellar nucleosynthesis, the stellar mass spectrum, the timescale of structure formation, the efficiency of star formation… ▽ More The trends of chemical abundances and abundance ratios observed in stars of different ages, kinematics, and metallicities bear the imprints of several physical processes that concur to shape the host galaxy properties. By inspecting these trends, we get precious information on stellar nucleosynthesis, the stellar mass spectrum, the timescale of structure formation, the efficiency of star formation, as well as any inward or outward flows of gas. In this paper, we analyse recent determinations of carbon-to-iron and carbon-to-oxygen abundance ratios in different environments (the Milky Way and elliptical galaxies) using our latest chemical evolution models that implement up-to-date stellar yields and rely on the tight constraints provided by asteroseismic stellar ages (whenever available). A scenario where most carbon is produced by rotating massive stars, with yields largely dependent on the metallicity of the parent proto-star clouds, allows us to fit simultaneously the high-quality data available for the local Galactic components (thick and thin discs) and for microlensed dwarf stars in the Galactic bulge, as well as the abundance ratios inferred for massive elliptical galaxies. Yet, more efforts are needed from both observers and theoreticians in order to base these conclusions on firmer grounds. △ Less

Submitted 12 May, 2020; originally announced May 2020.

Comments: 10 pages, 8 figures, 1 table, accepted for publication in A&A

Journal ref: A&A 639, A37 (2020)

Detailed abundances in the Galactic center: Evidence of a metal-rich alpha-enhanced stellar population

Authors: B. Thorsbro, N. Ryde, R. M. Rich, M. Schultheis, F. Renaud, E. Spitoni, T. K. Fritz, A. Mastrobuono-Battisti, L. Origlia, F. Matteucci, R. Schödel

Abstract: We present a detailed study of the composition of 20 M giants in the Galactic center with 15 of them confirmed to be in the Nuclear Star Cluster. As a control sample we have also observed 7 M giants in the Milky Way Disk with similar stellar parameters. All 27 stars are observed using the NIRSPEC spectograph on the KECK II telescope in the K-band at a resolving power of R=23,000. We report the fir… ▽ More We present a detailed study of the composition of 20 M giants in the Galactic center with 15 of them confirmed to be in the Nuclear Star Cluster. As a control sample we have also observed 7 M giants in the Milky Way Disk with similar stellar parameters. All 27 stars are observed using the NIRSPEC spectograph on the KECK II telescope in the K-band at a resolving power of R=23,000. We report the first silicon abundance trends versus [Fe/H] for stars in the Galactic center. While finding a disk/bulge like trend at subsolar metallicities, we find that [Si/Fe] is enhanced at supersolar metallicities. We speculate on possible enrichment scenarios to explain such a trend. However, the sample size is modest and the result needs to be confirmed by additional measurements of silicon and other \textalpha-elements. We also derive a new distribution of [Fe/H] and find the most metal rich stars at [Fe/H]=+0.5 dex, confirming our earlier conclusions that the Galactic center hosts no stars with extreme chemical composition. △ Less

Submitted 26 March, 2020; v1 submitted 24 March, 2020; originally announced March 2020.

Comments: 13 pages, 10 figures, accepted for publication in ApJ, Figure 7 corrected

Fluorine in the solar neighbourhood: the need for several cosmic sources

Authors: Nils Ryde, Henrik Jönsson, Gregory Mace, Katia Cunha, Emanuele Spitoni, Melike Afsar, Daniel Jaffe, Rebecca Forsberg, Kyle F. Kaplan, Benjamin T. Kidder, Jae-Joon Lee, Heeyoung Oh, Verne V. Smith, Christopher Sneden, Kimberly R. Sokal, Emily Strickland, Brian Thorsbro

Abstract: The cosmic origin of fluorine is still not well constrained. Several nucleosynthetic channels at different phases of stellar evolution have been suggested, but these must be constrained by observations. For this, the fluorine abundance trend with metallicity spanning a wide range is required. Our aim is to determine stellar abundances of fluorine for -1.1<[Fe/H]<+0.4. We determine the abundances f… ▽ More The cosmic origin of fluorine is still not well constrained. Several nucleosynthetic channels at different phases of stellar evolution have been suggested, but these must be constrained by observations. For this, the fluorine abundance trend with metallicity spanning a wide range is required. Our aim is to determine stellar abundances of fluorine for -1.1<[Fe/H]<+0.4. We determine the abundances from HF lines in infrared K-band spectra (at 2.3 microns) of cool giants, observed with the IGRINS and Phoenix high-resolution spectrographs. We derive accurate stellar parameters for all our observed K giants, which is important since the HF lines are very temperature sensitive. We find that [F/Fe] is flat as a function of metallicity at [F/Fe]~0, but increases as the metallicity increases. The fluorine slope shows a clear secondary behavior in this metallicity range. We also find that the [F/Ce] ratio is relatively flat for -0.6<[Fe/H]<0, and that for two metal-poor ([Fe/H]<-0.8), s-process element enhanced giants, we do not detect an elevated fluorine abundance. We interpret all these observational constraints to indicate that several major processes are at play for the cosmic budget of fluorine over time; from those in massive stars at low metallicities, through the asymptotic giant branch-star contribution at -0.6<[Fe/H]<0, to processes with increasing yields with metallicity at super-solar metallicities. The origins of the latter, and whether or not Wolf-Rayet stars and/or novae could contribute at super-solar metallicities, is currently not known. To quantify these observational results, theoretical modelling is required. More observations in the metal-poor region are required to clarify the processes there. △ Less

Submitted 10 March, 2020; originally announced March 2020.

Comments: Accepted by ApJ, 17 pages, 10 figures

On the connection between galactic downsizing and the most fundamental galactic scaling relations

Authors: E. Spitoni, F. Calura, M. Mignoli, R. Gilli, V. Silva Aguirre, A. Gallazzi

Abstract: In their evolution, star-forming galaxies are known to follow scaling relations between some fundamental physical quantities, such as the mass-metallicity and the main sequence relations. We aim at studying the evolution of galaxies that, at a given redshift, lie simultaneously on the mass-metallicity and main sequence relations (MZR, MSR). To this aim, we use the analytical, 'leaky-box' chemical… ▽ More In their evolution, star-forming galaxies are known to follow scaling relations between some fundamental physical quantities, such as the mass-metallicity and the main sequence relations. We aim at studying the evolution of galaxies that, at a given redshift, lie simultaneously on the mass-metallicity and main sequence relations (MZR, MSR). To this aim, we use the analytical, 'leaky-box' chemical evolution model of Spitoni et al. (2017), in which galaxy evolution is described by an infall timescale $τ$ and a wind efficiency $λ$. We provide a detailed analysis of the temporal evolution of galactic metallicity, stellar mass, mass-weighted age and gas fraction. The evolution of the galaxies lying on the MZR and MSR at $z\sim0.1$ suggests that the average infall time-scale in two different bins of stellar masses ($M_{\star}<10^{10} M_{\odot}$ and $M_{\star}>10^{10} M_{\odot}$) decreases with decreasing redshift. This means that at each redshift, only the youngest galaxies can be assembled on the shortest timescales and still belong to the star-forming MSR. In the lowest mass bin, a decrease of the median $τ$ is accompanied by an increase of the median $λ$ value. This implies that systems which have formed at more recent times will need to eject a larger amount of mass to keep their metallicity at low values. Another important result is that galactic downsizing, as traced by the age-mass relation, is naturally recovered by imposing that local galaxies lie on both the MZR and MSR. Finally, we study the evolution of the hosts of C$_{\rm IV}$ -selected AGN, which at $z\sim 2$ follow a flat MZR, as found by Mignoli et al. (2019). If we impose that these systems lie on the MSR, at lower redshifts we find an 'inverted' MZR, meaning that some additional processes must be at play in their evolution. △ Less

Submitted 1 August, 2020; v1 submitted 4 March, 2020; originally announced March 2020.

Comments: Accepted for publication in Astronomy and Astrophysics (A&A), 20 pages, 26 figures

Journal ref: A&A 642, A113 (2020)

Detection and characterisation of oscillating red giants: first results from the TESS satellite

Authors: Víctor Silva Aguirre, Dennis Stello, Amalie Stokholm, Jakob R. Mosumgaard, Warrick Ball, Sarbani Basu, Diego Bossini, Lisa Bugnet, Derek Buzasi, Tiago L. Campante, Lindsey Carboneau, William J. Chaplin, Enrico Corsaro, Guy R. Davies, Yvonne Elsworth, Rafael A. García, Patrick Gaulme, Oliver J. Hall, Rasmus Handberg, Marc Hon, Thomas Kallinger, Liu Kang, Mikkel N. Lund, Savita Mathur, Alexey Mints , et al. (56 additional authors not shown)

Abstract: Since the onset of the `space revolution' of high-precision high-cadence photometry, asteroseismology has been demonstrated as a powerful tool for informing Galactic archaeology investigations. The launch of the NASA TESS mission has enabled seismic-based inferences to go full sky -- providing a clear advantage for large ensemble studies of the different Milky Way components. Here we demonstrate i… ▽ More Since the onset of the `space revolution' of high-precision high-cadence photometry, asteroseismology has been demonstrated as a powerful tool for informing Galactic archaeology investigations. The launch of the NASA TESS mission has enabled seismic-based inferences to go full sky -- providing a clear advantage for large ensemble studies of the different Milky Way components. Here we demonstrate its potential for investigating the Galaxy by carrying out the first asteroseismic ensemble study of red giant stars observed by TESS. We use a sample of 25 stars for which we measure their global asteroseimic observables and estimate their fundamental stellar properties, such as radius, mass, and age. Significant improvements are seen in the uncertainties of our estimates when combining seismic observables from TESS with astrometric measurements from the Gaia mission compared to when the seismology and astrometry are applied separately. Specifically, when combined we show that stellar radii can be determined to a precision of a few percent, masses to 5-10% and ages to the 20% level. This is comparable to the precision typically obtained using end-of-mission Kepler data △ Less

Submitted 5 February, 2020; v1 submitted 16 December, 2019; originally announced December 2019.

Comments: Accepted for publication in ApJ Letters

Galactic Archaeology with asteroseismic ages. II. Confirmation of a delayed gas infall using Bayesian analysis based on MCMC methods

Authors: E. Spitoni, K. Verma, V. Silva Aguirre, F. Calura

Abstract: With the wealth of information from large surveys and observational campaigns in the contemporary era, it is critical to properly exploit the data to constrain the parameters of Galactic chemical evolution models and quantify the associated uncertainties. We aim at constraining the two-infall chemical evolution models for the solar annulus using the measured chemical abundance ratios and seismical… ▽ More With the wealth of information from large surveys and observational campaigns in the contemporary era, it is critical to properly exploit the data to constrain the parameters of Galactic chemical evolution models and quantify the associated uncertainties. We aim at constraining the two-infall chemical evolution models for the solar annulus using the measured chemical abundance ratios and seismically inferred age of stars in the APOKASC sample. In the revised two-infall chemical evolution models by Spitoni et al. (2019), a significant delay of $\sim 4.3$ Gyr has been invoked between the two episodes of gas accretion. In this work, we wish to test its robustness and statistically confirm/quantify the delay. For the first time, a Bayesian framework based on Markov Chain Monte Carlo methods has been used for fitting the two-infall chemical evolution models to the data. In addition to fitting the data for stars in the APOKASC sample, our best fit models also reproduce other important observational constraints of the chemical evolution of the disk: i) present day stellar surface density; ii) present day supernova and star formation rates; iii) the metallicity distribution function; and iv) solar abundance values. We found a significant delay between the two gas accretion episodes for various models explored with different values for the star formation efficiencies. The values for the delay lie in the range $4.5-5.5$ Gyr. The results suggest that the APOKASC sample carries the signature of delayed gas-rich merger, with the dilution as main process determining the shape of low-$α$ stars in the abundance ratios space. △ Less

Submitted 28 January, 2020; v1 submitted 9 December, 2019; originally announced December 2019.

Comments: Accepted for publication in Astronomy and Astrophysics (A&A), 12 pages, 13 figures

Journal ref: A&A 635, A58 (2020)

The origin of stellar populations in the Galactic bulge from chemical abundances

Authors: Francesca Matteucci, Valeria Grisoni, Emanuele Spitoni, Anna Zulianello, Alvaro Rojas-Arriagada, Mathias Schultheis, Nils Ryde

Abstract: In this work, we study the formation and chemical evolution of the Galactic bulge with particular focus on the abundance pattern ([Mg/Fe] vs. [Fe/H]), metallicity and age distribution functions. We consider detailed chemical evolution models for the Galactic bulge and inner disc, with the aim of shedding light on the connection between these components and the origin of bulge stars. In particular,… ▽ More In this work, we study the formation and chemical evolution of the Galactic bulge with particular focus on the abundance pattern ([Mg/Fe] vs. [Fe/H]), metallicity and age distribution functions. We consider detailed chemical evolution models for the Galactic bulge and inner disc, with the aim of shedding light on the connection between these components and the origin of bulge stars. In particular, we first present a model assuming a fast and intense star formation, with the majority of bulge stars forming on a timescale less than 1 Gyr. Then we analyze the possibility of two distinct stellar populations in the bulge, as suggested by Gaia-ESO and APOGEE data. These two populations, one metal poor and the other metal rich, can have had two different origins: i) the metal rich formed after a stop of roughly 250 Myr in the star formation rate of the bulge, or ii) the metal rich population is made of stars formed in the inner disc and brought into the bulge by the early secular evolution of the bar. We also examine the case of multiple star bursts in the bulge with consequent formation of multiple populations, as suggested by studies of microlensed stars. After comparing model results and observations, we suggest that the most likely scenario is that there are two main stellar populations, both made mainly by old stars (> 10 Gyr), with the metal rich and younger one formed from inner thin disc stars, in agreement with kinematical arguments. However, on the basis of dynamical simulations, we cannot completely exclude that the second population formed after a stop in the star formation during the bulge evolution, so that all the stars formed in situ. △ Less

Submitted 21 June, 2019; v1 submitted 29 March, 2019; originally announced March 2019.

Comments: 11 pages, 9 figures, accepted by MNRAS

The Fall of a Giant. Chemical evolution of Enceladus, alias the Gaia Sausage

Authors: Fiorenzo Vincenzo, Emanuele Spitoni, Francesco Calura, Francesca Matteucci, Victor Silva Aguirre, Andrea Miglio, Gabriele Cescutti

Abstract: We present the first chemical evolution model for Enceladus, alias the Gaia Sausage, to investigate the star formation history of one of the most massive satellites accreted by the Milky Way during a major merger event. Our best chemical evolution model for Enceladus nicely fits the observed stellar [$α$/Fe]-[Fe/H] chemical abundance trends, and reproduces the observed stellar metallicity distribu… ▽ More We present the first chemical evolution model for Enceladus, alias the Gaia Sausage, to investigate the star formation history of one of the most massive satellites accreted by the Milky Way during a major merger event. Our best chemical evolution model for Enceladus nicely fits the observed stellar [$α$/Fe]-[Fe/H] chemical abundance trends, and reproduces the observed stellar metallicity distribution function, by assuming low star formation efficiency, fast infall time scale, and mild outflow intensity. We predict a median age for Enceladus stars $12.33^{+0.92}_{-1.36}$ Gyr, and - at the time of the merger with our Galaxy ($\approx10$ Gyr ago from Helmi et al.) - we predict for Enceladus a total stellar mass $M_{\star} \approx 5 \times 10^{9}\,\text{M}_{\odot}$. By looking at the predictions of our best model, we discuss that merger events between the Galaxy and systems like Enceladus may have inhibited the gas accretion onto the Galaxy disc at high redshifts, heating up the gas in the halo. This scenario could explain the extended period of quenching in the star formation activity of our Galaxy about 10 Gyr ago, which is predicted by Milky Way chemical evolution models, in order to reproduce the observed bimodality in [$α$/Fe]-[Fe/H] between thick- and thin-disc stars. △ Less

Submitted 12 May, 2019; v1 submitted 8 March, 2019; originally announced March 2019.

Comments: Accepted for publication in MNRAS Letters

2D chemical evolution model: the impact of galactic disc asymmetries on azimuthal chemical abundance variations

Authors: E. Spitoni, G. Cescutti, I. Minchev, F. Matteucci, V. Silva Aguirre, M. Martig, G. Bono, C. Chiappini

Abstract: Galactic disc chemical evolution models generally ignore azimuthal surface density variation that can introduce chemical abundance azimuthal gradients. Recent observations, however, have revealed chemical abundance changes with azimuth in the gas and stellar components of both the Milky Way and external galaxies. To quantify the effects of spiral arm density fluctuations on the azimuthal variation… ▽ More Galactic disc chemical evolution models generally ignore azimuthal surface density variation that can introduce chemical abundance azimuthal gradients. Recent observations, however, have revealed chemical abundance changes with azimuth in the gas and stellar components of both the Milky Way and external galaxies. To quantify the effects of spiral arm density fluctuations on the azimuthal variations of the oxygen and iron abundances in disc galaxies. We develop a new 2D galactic disc chemical evolution model, capable of following not just radial but also azimuthal inhomogeneities. The density fluctuations resulting from a Milky Way-like N-body disc formation simulation produce azimuthal variations in the oxygen abundance gradients of the order of 0.1 dex. Moreover, in agreement with the most recent observations in external galaxies, the azimuthal variations are more evident in the outer galactic regions. Using a simple analytical model, we show that the largest fluctuations with azimuth result near the spiral structure corotation resonance, where the relative speed between spiral and gaseous disc is the slowest. In conclusion we provided a new 2D chemical evolution model capable of following azimuthal density variations. Density fluctuations extracted from a Milky Way-like dynamical model lead to a scatter in the azimuthal variations of the oxygen abundance gradient in agreement with observations in external galaxies. We interpret the presence of azimuthal scatter at all radii by the presence of multiple spiral modes moving at different pattern speeds, as found in both observations and numerical simulations. △ Less

Submitted 24 June, 2019; v1 submitted 27 November, 2018; originally announced November 2018.

Comments: Accepted for publication in Astronomy and Astrophysics (A&A), 17 pages, 18 Figures

Journal ref: A&A 628, A38 (2019)

Galactic Archaeology with asteroseismic ages: evidence for delayed gas infall in the formation of the Milky Way disc

Authors: E. Spitoni, V. Silva Aguirre, F. Matteucci, F. Calura, V. Grisoni

Abstract: Precise stellar ages from asteroseismology have become available and can help setting stronger constraints on the evolution of the Galactic disc components. Recently, asteroseismology has confirmed a clear age difference in the solar annulus between two distinct sequences in the [$α$/Fe] versus [Fe/H] abundance ratios relation: the high-$α$ and low-$α$ stellar populations. We aim at reproducing th… ▽ More Precise stellar ages from asteroseismology have become available and can help setting stronger constraints on the evolution of the Galactic disc components. Recently, asteroseismology has confirmed a clear age difference in the solar annulus between two distinct sequences in the [$α$/Fe] versus [Fe/H] abundance ratios relation: the high-$α$ and low-$α$ stellar populations. We aim at reproducing these new data with chemical evolution models including different assumptions for the history and number of accretion events. We tested two different approaches: a revised version of the `two-infall' model where the high-$α$ phase forms by a fast gas accretion episode and the low-$α$ sequence follows later from a slower gas infall rate, and the parallel formation scenario where the two disc sequences form coevally and independently. The revised `two-infall' model including uncertainties in age and metallicity is capable of reproducing: i) the [$α$/Fe] vs. [Fe/H] abundance relation at different Galactic epochs, ii) the age$-$metallicity relation and the time evolution [$α$/Fe]; iii) the age distribution of the high-$α$ and low-$α$ stellar populations, iv) the metallicity distribution function. The parallel approach is not capable of properly reproduce the stellar age distribution, in particular at old ages. In conclusion, the best chemical evolution model is the revised `two-infall' one, where a consistent delay of $\sim$4.3 Gyr in the beginning of the second gas accretion episode is a crucial assumption to reproduce stellar abundances and ages. △ Less

Submitted 29 January, 2019; v1 submitted 4 September, 2018; originally announced September 2018.

Comments: Accepted for publication in A&A; 20 pages, 20 figures

Journal ref: A&A 623, A60 (2019)

Abundance gradients along the Galactic disc from chemical evolution models

Authors: V. Grisoni, E. Spitoni, F. Matteucci

Abstract: In this paper, we study the formation and chemical evolution of the Milky Way disc with particular focus on the abundance patterns ([$α$/Fe] vs. [Fe/H]) at different Galactocentric distances, the present-time abundance gradients along the disc and the time evolution of abundance gradients. We consider the chemical evolution models for the Galactic disc developed by Grisoni et al. (2017) for the so… ▽ More In this paper, we study the formation and chemical evolution of the Milky Way disc with particular focus on the abundance patterns ([$α$/Fe] vs. [Fe/H]) at different Galactocentric distances, the present-time abundance gradients along the disc and the time evolution of abundance gradients. We consider the chemical evolution models for the Galactic disc developed by Grisoni et al. (2017) for the solar neighborhood, both the two-infall and the one-infall ones, and we extend our analysis to the other Galactocentric distances. In particular, we examine the processes which mainly influence the formation of the abundance gradients: the inside-out scenario, a variable star formation efficiency, and radial gas flows. We compare our model results with recent abundance patterns obtained along the Galactic disc from the APOGEE survey and with abundance gradients observed from Cepheids, open clusters, HII regions and PNe. We conclude that the inside-out scenario is a key ingredient, but cannot be the only one to explain abundance patterns at different Galactocentric distances and abundance gradients. Further ingredients, such as radial gas flows and variable star formation efficiency, are needed to reproduce the observed features in the thin disc. The evolution of abundance gradients with time is also shown, although firm conclusions cannot still be drawn. △ Less

Submitted 10 September, 2018; v1 submitted 28 May, 2018; originally announced May 2018.

Comments: 12 pages, 7 figures. Accepted for publication in MNRAS

Is the IMF in ellipticals bottom-heavy? Clues from their chemical abundances

Authors: C. De Masi, F. Vincenzo, F. Matteucci, G. Rosani, La Barbera, A. Pasquali, E. Spitoni

Abstract: We tested the implementation of different IMFs in our model for the chemical evolution of ellipticals, with the aim of reproducing the observed relations of [Fe/H] and [Mg/Fe] abundances with galaxy mass in a sample of early-type galaxies selected from the SPIDER-SDSS catalog. Abundances in the catalog were derived from averaged spectra, obtained by stacking individual spectra according to central… ▽ More We tested the implementation of different IMFs in our model for the chemical evolution of ellipticals, with the aim of reproducing the observed relations of [Fe/H] and [Mg/Fe] abundances with galaxy mass in a sample of early-type galaxies selected from the SPIDER-SDSS catalog. Abundances in the catalog were derived from averaged spectra, obtained by stacking individual spectra according to central velocity dispersion, as a proxy of galaxy mass. We tested initial mass functions already used in a previous work, as well as two new models, based on low-mass tapered ("bimodal") IMFs, where the IMF becomes either (1) bottom-heavy in more massive galaxies, or (2) is time-dependent, switching from top-heavy to bottom-heavy in the course of galactic evolution. We found that observations could only be reproduced by models assuming either a constant, Salpeter IMF, or a time-dependent distribution, as other IMFs failed. We further tested the models by calculating their M/L ratios. We conclude that a constant, time-independent bottom-heavy IMF does not reproduce the data, especially the increase of the $[α/Fe]$ ratio with galactic stellar mass, whereas a variable IMF, switching from top to bottom-heavy, can match observations. For the latter models, the IMF switch always occurs at the earliest possible considered time, i.e. $t_{\text{switch}}= 0.1$ Gyr. △ Less

Submitted 15 November, 2018; v1 submitted 17 May, 2018; originally announced May 2018.

Comments: 24 pages, 25 figures, accepted for publication on MNRAS

Fluorine in the Solar Neighborhood: Chemical Evolution Models

Authors: E. Spitoni, F. Matteucci, H. Jönsson, N. Ryde, D. Romano

Abstract: In the light of the new observational data related to fluorine abundances in the solar neighborhood stars, we present here chemical evolution models testing different fluorine nucleosynthesis prescriptions with the aim to best fit those new data related to the abundance ratios [F/O] vs. [O/H] and [F/Fe] vs. [Fe/H]. The adopted chemical evolution models are: i) the classical "two-infall" model whic… ▽ More In the light of the new observational data related to fluorine abundances in the solar neighborhood stars, we present here chemical evolution models testing different fluorine nucleosynthesis prescriptions with the aim to best fit those new data related to the abundance ratios [F/O] vs. [O/H] and [F/Fe] vs. [Fe/H]. The adopted chemical evolution models are: i) the classical "two-infall" model which follows the chemical evolution of halo-thick disk and thin disk phases, ii) and the "one-infall" model designed only for the thin disk evolution. We tested the effects on the predicted fluorine abundance ratios of different nucleosynthesis yield sources: AGB stars, Wolf-Rayet stars, Type II and Type Ia supernovae, and novae. We find that the fluorine production is dominated by AGB stars but the Wolf-Rayet stars are required to reproduce the trend of the observed data in the solar neighborhood by Jönsson et al. (2017a) with our chemical evolution models. In particular, the best model both for the "two-infall" and "one-infall" cases requires an increase by a factor of two of the Wolf-Rayet yields given by Meynet & Arnould (2000). We also show that the novae, even if their yields are still uncertain, could help to better reproduce the secondary behavior of F in the [F/O] vs. [O/H] relation. The inclusion of the fluorine production by Wolf-Rayet stars seems to be essential to reproduce the observed ratio [F/O] vs [O/H] in the solar neighborhood by Jönsson et al. (2017a). Moreover, the inclusion of novae helps substantially to reproduce the observed fluorine secondary behavior. △ Less

Submitted 27 November, 2017; originally announced November 2017.

Comments: 10 pages, 10 figure. Accepted for publication in A&A

Journal ref: A&A 612, A16 (2018)

Origin of the Galactic Halo: accretion vs. in situ formation

Authors: E. Spitoni, F. Vincenzo, F. Matteucci, D. Romano

Abstract: We test the hypothesis that the classical and ultra-faint dwarf spheroidal satellites of the our Galaxy have been the building blocks of the Galactic halo by comparing their [O/Fe] and [Ba/Fe] vs. [Fe/H] patterns with the ones observed in Galactic halo stars. The [O/Fe] ratio deviates substantially from the observed abundance ratios in the Galactic halo stars for [Fe/H] > -2 dex, while they overla… ▽ More We test the hypothesis that the classical and ultra-faint dwarf spheroidal satellites of the our Galaxy have been the building blocks of the Galactic halo by comparing their [O/Fe] and [Ba/Fe] vs. [Fe/H] patterns with the ones observed in Galactic halo stars. The [O/Fe] ratio deviates substantially from the observed abundance ratios in the Galactic halo stars for [Fe/H] > -2 dex, while they overlap for lower metallicities. On the other hand, for the neutron capture elements, the discrepancy is extended at all the metallicities, suggesting that the majority of stars in the halo are likely to have been formed in situ. We present the results for a model considering the effects of an enriched gas stripped from dwarf satellites on the chemical evolution of the Galactic halo. We find that the resulting chemical abundances of the halo stars depend on the adopted infall time-scale, and the presence of a threshold in the gas for star formation. △ Less

Submitted 25 July, 2017; v1 submitted 24 July, 2017; originally announced July 2017.

Comments: 4 pages, 2 figures, to appear in Proceedings of IAU Symposium 334 "Rediscovering our Galaxy" in Potsdam (Germany), July 10-14, 2017

The AMBRE Project: chemical evolution models for the Milky Way thick and thin discs

Authors: V. Grisoni, E. Spitoni, F. Matteucci, A. Recio-Blanco, P. de Laverny, M. Hayden, Š. Mikolaitis, C. C. Worley

Abstract: We study the chemical evolution of the thick and thin discs of the Galaxy by comparing detailed chemical evolution models with recent data from the AMBRE Project. The data suggest that the stars in the thick and thin discs form two distinct sequences with the thick disc stars showing higher [α/Fe] ratios. We adopt two different approaches to model the evolution of thick and thin discs. In particul… ▽ More We study the chemical evolution of the thick and thin discs of the Galaxy by comparing detailed chemical evolution models with recent data from the AMBRE Project. The data suggest that the stars in the thick and thin discs form two distinct sequences with the thick disc stars showing higher [α/Fe] ratios. We adopt two different approaches to model the evolution of thick and thin discs. In particular, we adopt: i) a two-infall approach where the thick disc forms fast and before the thin disc and by means of a fast gas accretion episode, whereas the thin disc forms by means of a second accretion episode on a longer timescale; ii) a parallel approach, where the two discs form in parallel but at different rates. By comparing our model results with the observed [Mg/Fe] vs. [Fe/H] and the metallicity distribution functions in the two Galactic components, we conclude that the parallel approach can account for a group of α-enhanced metal rich stars present in the data, whereas the two-infall approach cannot explain these stars unless they are the result of stellar migration. In both approaches, the thick disc has formed on a timescale of accretion of 0.1 Gyr, whereas the thin disc formed on a timescale of 7 Gyr in the solar region. In the two-infall approach a gap in star formation between the thick and thin disc formation of several hundreds of Myr should be present, at variance with the parallel approach where no gap is present. △ Less

Submitted 29 August, 2017; v1 submitted 8 June, 2017; originally announced June 2017.

Comments: 12 pages, 10 figures. Accepted for publication in MNRAS