Showing 1–50 of 53 results for author: Grassi, T

Planet Formation and Disk Chemistry: Dust and Gas Evolution during Planet Formation

Authors: G. Perotti, L. Cacciapuoti, N. -D. Tung, T. Grassi, E. Schisano, L. Testi

Abstract: Over the past decade, progress in observational capabilities, combined with theoretical advancements, have transformed our comprehension of the physics and chemistry during planet formation. Despite these important steps forward, open questions persist on the chemical and physical evolution of solids in their journey from the collapsing molecular cores to disks and planetary bodies. This chapter i… ▽ More Over the past decade, progress in observational capabilities, combined with theoretical advancements, have transformed our comprehension of the physics and chemistry during planet formation. Despite these important steps forward, open questions persist on the chemical and physical evolution of solids in their journey from the collapsing molecular cores to disks and planetary bodies. This chapter is a repository of such burning questions. It has the ambition to identify the most promising avenues for future research based on current observational and modeling opportunities. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: 9 pages, 5 figures, Chapter IV of the National Congress of (proto-)Planetary Astrochemistry (CNAP II) Proceedings

Testing analytical methods to derive the cosmic-ray ionisation rate in cold regions via synthetic observations

Authors: E. Redaelli, S. Bovino, A. Lupi, T. Grassi, D. Gaete-Espinoza, G. Sabatini, P. Caselli

Abstract: Cosmic rays (CRs) heavily impact the chemistry and physics of cold and dense star-forming regions. However, characterising their ionisation rate is still challenging from an observational point of view. In the past, a few analytical formulas have been proposed to infer the cosmic-ray ionization rate $ζ_2$ from molecular line observations. These have been derived from the chemical kinetics of the i… ▽ More Cosmic rays (CRs) heavily impact the chemistry and physics of cold and dense star-forming regions. However, characterising their ionisation rate is still challenging from an observational point of view. In the past, a few analytical formulas have been proposed to infer the cosmic-ray ionization rate $ζ_2$ from molecular line observations. These have been derived from the chemical kinetics of the involved species, but they have not been validated using synthetic data processed with a standard observative pipeline. We aim to bridge this gap. We perform the radiative transfer on a set of three-dimensional magneto-hydrodynamical simulations of prestellar cores, exploring different initial $ζ_2$, evolutionary stages, types of radiative transfer (e.g. assuming local-thermodynamic-equilibrium conditions), and telescope responses. We then compute the column densities of the involved tracers to determine $ζ_2$, using, in particular, the equation proposed by Bovino et. al (2020) and by Caselli et al. (1998) both used nowadays. Our results confirm that the method of Bovino et al. (2020) accurately retrieves the actual $ζ_2$ within a factor of $2-3$, in the physical conditions explored in our tests. Since we also explore a non-local thermodynamic equilibrium radiative transfer, this work indirectly offers insights into the excitation temperatures of common transitions at moderate volume densities ($n\approx 10^5 \, \rm cm^{-3}$). We have also performed a few tests using the formula proposed by Caselli et al. (1998), which overestimates the actual $ζ_2$ by at least two orders of magnitudes. We also consider a new derivation of this method, which, however, still leads to large overestimates. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Comments: Accepted for publication in A&A (15th Feb.)

Ultraviolet H$_2$ luminescence in molecular clouds induced by cosmic rays

Authors: Marco Padovani, Daniele Galli, Liam H. Scarlett, Tommaso Grassi, Una S. Rehill, Mark C. Zammit, Igor Bray, Dmitry V. Fursa

Abstract: Galactic cosmic rays (CRs) play a crucial role in ionisation, dissociation, and excitation processes within dense cloud regions where UV radiation is absorbed by dust grains and gas species. CRs regulate the abundance of ions and radicals, leading to the formation of more and more complex molecular species, and determine the charge distribution on dust grains. A quantitative analysis of these effe… ▽ More Galactic cosmic rays (CRs) play a crucial role in ionisation, dissociation, and excitation processes within dense cloud regions where UV radiation is absorbed by dust grains and gas species. CRs regulate the abundance of ions and radicals, leading to the formation of more and more complex molecular species, and determine the charge distribution on dust grains. A quantitative analysis of these effects is essential for understanding the dynamical and chemical evolution of star-forming regions. The CR-induced photon flux has a significant impact on the evolution of the dense molecular medium in its gas and dust components. This study is intended to evaluate the flux of UV photons generated by CRs to calculate the photon-induced dissociation and ionisation rates of a vast number of atomic and molecular species, as well as the integrated UV photon flux. Our study takes advantage of recent developments in the determination of the spectra of secondary electrons, in the calculation of state-resolved excitation cross sections of H$_2$ by electron impact, and of photodissociation and photoionisation cross sections. We calculate the H$_2$ level population of each rovibrational level of the $X$, $B$, $C$, $B'$, $D$, $B''$, $D'$ and $a$ states. We then compute the UV photon spectrum of H$_2$ in its line and continuum components between 72 and 700 nm, with unprecedented accuracy as a function of the CR spectrum incident on a molecular cloud, the H$_2$ column density, the isomeric H$_2$ composition, and the dust properties. The resulting photodissociation and photoionisation rates are, on average, smaller than previous determinations by a factor of about 2, with deviations up to a factor of 5. A special focus is given to the photoionisation rates of H$_2$, HF, and H$_2$, as well as to the photodissociation of H$_2$, which we find to be orders of magnitude higher than previous estimates. △ Less

Submitted 14 December, 2023; v1 submitted 4 December, 2023; originally announced December 2023.

Comments: 22 pages, 14 figures, accepted by Astronomy & Astrophysics

The Astrochemistry Low-energy Electron Cross-Section (ALeCS) database I. Semi-empirical electron-impact ionization cross-section calculations and ionization rates

Authors: Brandt A. L. Gaches, Tommaso Grassi, Stefan Vogt-Geisse, Giulia M. Bovolenta, Claire Vallance, David Heathcote, Marco Padovani, Stefano Bovino, Prasanta Gorai

Abstract: (Abridged) Electron-molecule interaction is a fundamental process in radiation-driven chemistry in space, from the interstellar medium to comets. Therefore, knowledge of interaction cross-sections is key. While there has been a plethora of studies of total ionization cross-sections, data is often spread over many sources, or not public or readily available. We introduce the Astrochemistry Low-ener… ▽ More (Abridged) Electron-molecule interaction is a fundamental process in radiation-driven chemistry in space, from the interstellar medium to comets. Therefore, knowledge of interaction cross-sections is key. While there has been a plethora of studies of total ionization cross-sections, data is often spread over many sources, or not public or readily available. We introduce the Astrochemistry Low-energy Electron Cross-Section (ALeCS) database, a public database for electron interaction cross-sections and ionization rates for molecules of astrochemical interest. In this work, we present the first data release comprising total ionization cross-sections and ionization rates for over 200 neutral molecules. We include optimized geometries and molecular orbital energies at various levels of theory, and for a subset of the molecules, the ionization potentials. We compute total ionization cross-sections using the binary-encounter Bethe model and screening-corrected additivity rule, and ionization rates and reaction network coefficients for molecular cloud environments for $>$200 neutral molecules ranging from diatomics to complex organics. We demonstrate that our binary-encounter Bethe cross-sections agree well with experimental data. We show that the ionization rates scale roughly linearly with the number of constituent atoms in the molecule. We introduce and describe the public ALeCS database. For the initial release, we include total ionization cross-sections for $>$200 neutral molecules and several cations and anions calculated with different levels of quantum chemistry theory, the chemical reaction rates for the ionization, and network files in the formats of the two most popular astrochemical networks, the KIDA and UMIST. The database will be continuously updated for more molecules and interactions. △ Less

Submitted 13 February, 2024; v1 submitted 16 October, 2023; originally announced October 2023.

Comments: Accepted to A&A. Public database online at https://github.com/AstroBrandt/ALeCS

A 3D physico-chemical model of a pre-stellar core. I. Environmental and structural impact on the distribution of CH$_3$OH and $c$-C$_3$H$_2$

Authors: S. S. Jensen, S. Spezzano, P. Caselli, T. Grassi, T. Haugbølle

Abstract: Pre-stellar cores represent the earliest stage of the star- and planet-formation process. By characterizing the physical and chemical structure of these cores we can establish the initial conditions for star and planet formation and determine to what degree the chemical composition of pre-stellar cores is inherited to the later stages. A 3D MHD model of a pre-stellar core embedded in a dynamic sta… ▽ More Pre-stellar cores represent the earliest stage of the star- and planet-formation process. By characterizing the physical and chemical structure of these cores we can establish the initial conditions for star and planet formation and determine to what degree the chemical composition of pre-stellar cores is inherited to the later stages. A 3D MHD model of a pre-stellar core embedded in a dynamic star-forming cloud is post-processed using sequentially continuum radiative transfer, a gas-grain chemical model, and a line-radiative transfer model. Results are analyzed and compared to observations of CH$_3$OH and $c$-C$_3$H$_2$ in L1544. Nine different chemical models are compared to the observations to determine which initial conditions are compatible with the observed chemical segregation in the prototypical pre-stellar core L1544. The model is able to reproduce several aspects of the observed chemical differentiation in L1544. Extended methanol emission is shifted towards colder and more shielded regions of the core envelope while $c$-C$_3$H$_2$ emission overlaps with the dust continuum, consistent with the observed chemical structure. Increasing the strength of the interstellar radiation field or the cosmic-ray ionization rate with respect to the typical values assumed in nearby star-forming regions leads to synthetic maps that are inconsistent with the observed chemical structure. Our model shows that the observed chemical dichotomy in L1544 can arise as a result of uneven illumination due to the asymmetrical structure of the 3D core and the environment within which the core has formed. This highlights the importance of the 3D structure at the core-cloud transition on the chemistry of pre-stellar cores. △ Less

Submitted 10 May, 2023; originally announced May 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 675, A34 (2023)

Implementation of chemistry in the Athena++ code

Authors: Munan Gong, Ka-Wai Ho, James M. Stone, Eve C. Ostriker, Paola Caselli, Tommaso Grassi, Chang-Goo Kim, Jeong-Gyu Kim, Goni Halevi

Abstract: Chemistry plays a key role in many aspects of astrophysical fluids. Atoms and molecules are agents for heating and cooling, determine the ionization fraction, serve as observational tracers, and build the molecular foundation of life. We present the implementation of a chemistry module in the publicly available magneto-hydrodynamic code Athena++. We implement several chemical networks and heating… ▽ More Chemistry plays a key role in many aspects of astrophysical fluids. Atoms and molecules are agents for heating and cooling, determine the ionization fraction, serve as observational tracers, and build the molecular foundation of life. We present the implementation of a chemistry module in the publicly available magneto-hydrodynamic code Athena++. We implement several chemical networks and heating and cooling processes suitable for simulating the interstellar medium (ISM). A general chemical network framework in the KIDA format is also included, allowing the user to easily implement their own chemistry. Radiation transfer and cosmic-ray ionization are coupled with chemistry and solved with the simple six-ray approximation. The chemical and thermal processes are evolved as a system of coupled ODEs with an implicit solver from the CVODE library. We perform and present a series of tests to ensure the numerical accuracy and convergence of the code. Many tests combine chemistry with gas dynamics, including comparisons with analytic solutions, 1D problems of the photo-dissociation regions and shocks, and realistic 3D simulations of the turbulent ISM. We release the code with the new public version of Athena++, aiming to provide a robust and flexible code for the astrochemical simulation community. △ Less

Submitted 8 May, 2023; originally announced May 2023.

Presence of liquid water during the evolution of exomoons orbiting ejected free-floating planets

Authors: Giulia Roccetti, Tommaso Grassi, Barbara Ercolano, Karan Molaverdikhani, Aurélien Crida, Dieter Braun, Andrea Chiavassa

Abstract: Free-floating planets (FFPs) can result from dynamical scattering processes happening in the first few million years of a planetary system's life. Several models predict the possibility, for these isolated planetary-mass objects, to retain exomoons after their ejection. The tidal heating mechanism and the presence of an atmosphere with a relatively high optical thickness may support the formation… ▽ More Free-floating planets (FFPs) can result from dynamical scattering processes happening in the first few million years of a planetary system's life. Several models predict the possibility, for these isolated planetary-mass objects, to retain exomoons after their ejection. The tidal heating mechanism and the presence of an atmosphere with a relatively high optical thickness may support the formation and maintenance of oceans of liquid water on the surface of these satellites. In order to study the timescales over which liquid water can be maintained, we perform dynamical simulations of the ejection process and infer the resulting statistics of the population of surviving exomoons around free-floating planets. The subsequent tidal evolution of the moons' orbital parameters is a pivotal step to determine when the orbits will circularize, with a consequential decay of the tidal heating. We find that close-in ($a \lesssim 25 $R$_{\rm J}$) Earth-mass moons with CO$_2$-dominated atmospheres could retain liquid water on their surfaces for long timescales, depending on the mass of the atmospheric envelope and the surface pressure assumed. Massive atmospheres are needed to trap the heat produced by tidal friction that makes these moons habitable. For Earth-like pressure conditions ($p_0$ = 1 bar), satellites could sustain liquid water on their surfaces up to 52 Myr. For higher surface pressures (10 and 100 bar), moons could be habitable up to 276 Myr and 1.6 Gyr, respectively. Close-in satellites experience habitable conditions for long timescales, and during the ejection of the FFP remain bound with the escaping planet, being less affected by the close encounter. △ Less

Submitted 9 February, 2023; originally announced February 2023.

Comments: 21 pages, 12 figures, accepted for publication on International Journal of Astrobiology

Non-ideal magnetohydrodynamics of self-gravitating filaments

Authors: Nicol Gutiérrez-Vera, Tommaso Grassi, Stefano Bovino, Alessandro Lupi, Daniele Galli, Dominik R. G. Schleicher

Abstract: Filaments have been studied in detail through observations and simulations. A range of numerical works have separately investigated how chemistry and diffusion effects, as well as magnetic fields and their structure impact the gas dynamics of the filament. However, non-ideal effects have hardly been explored thus far. We investigate how non-ideal magnetohydrodynamic (MHD) effects, combined with a… ▽ More Filaments have been studied in detail through observations and simulations. A range of numerical works have separately investigated how chemistry and diffusion effects, as well as magnetic fields and their structure impact the gas dynamics of the filament. However, non-ideal effects have hardly been explored thus far. We investigate how non-ideal magnetohydrodynamic (MHD) effects, combined with a simplified chemical model affect the evolution and accretion of a star-forming filament. We modeled an accreting self-gravitating turbulent filament using lemongrab, a one-dimensional (1D) non-ideal MHD code that includes chemistry. We explore the influence of non-ideal MHD, the orientation and strength of the magnetic field, and the cosmic ray ionization rate, on the evolution of the filament, with particular focus on the width and accretion rate. We find that the filament width and the accretion rate are determined by the magnetic field properties, including the initial strength, the coupling with the gas controlled by the cosmic ray ionization rate, and the orientation of the magnetic field with respect to the accretion flow direction. Increasing the cosmic-ray ionization rate leads to a behavior closer to that of ideal MHD, reducing the magnetic pressure support and, hence, damping the accretion efficiency with a consequent broadening of the filament width. For the same reason, we obtained a narrower width and a larger accretion rate when we reduced the initial magnetic field strength. Overall, while these factors affect the final results by approximately a factor of~2, removing the non-ideal MHD effects results in a much greater variation (up to a factor of~7). The inclusion of non-ideal MHD effects and the cosmic-ray ionization is crucial for the study of self-gravitating filaments and in determining critical observable quantities, such as the filament width and accretion rate. △ Less

Submitted 28 November, 2022; originally announced November 2022.

Comments: 11 pages, 9 figures, to be published in Astronomy & Astrophysics

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

Interpreting molecular hydrogen and atomic oxygen line emission of T Tauri disks with photoevaporative disk-wind models

Authors: Ch. Rab, M. Weber, T. Grassi, B. Ercolano, G. Picogna, P. Caselli, W. -F. Thi, I. Kamp, P. Woitke

Abstract: Winds in protoplanetary disks play an important role in their evolution and dispersal. However, what physical process is driving the winds is still unclear (i.e. magnetically vs thermally driven), and can only be understood by directly confronting theoretical models with observational data. We use hydrodynamic photoevaporative disk-wind models and post-process them with a thermo-chemical model to… ▽ More Winds in protoplanetary disks play an important role in their evolution and dispersal. However, what physical process is driving the winds is still unclear (i.e. magnetically vs thermally driven), and can only be understood by directly confronting theoretical models with observational data. We use hydrodynamic photoevaporative disk-wind models and post-process them with a thermo-chemical model to produce synthetic observables for the o-H$_2$ at 2.12 micron and [OI] at 0.63 micron spectral lines and directly compare the results to a sample of observations. Our photoevaporative disk-wind model is consistent with the observed signatures of the blueshifted narrow low-velocity component (NLVC), which is usually associated with slow disk winds, for both tracers. Only for one out of seven targets that show blueshifted NLVCs does the photoevaporative model fail to explain the observed line kinematics. Our results also indicate that interpreting spectral line profiles by simple methods, such as the thin-disk approximation, to determine the line emitting region can yield misleading conclusions. The photoevaporative disk-wind models are largely consistent with the studied observational data set, but it is not possible to clearly discriminate between different wind-driving mechanisms. Further improvements to the models, such as consistent modelling of the dynamics and chemistry and detailed modelling of individual targets would be beneficial. Furthermore, a direct comparison of magnetically driven disk-wind models to the observational data set is necessary in order to determine whether or not spatially unresolved observations of multiple wind tracers are sufficient to discriminate between theoretical models. △ Less

Submitted 18 November, 2022; v1 submitted 27 October, 2022; originally announced October 2022.

Comments: 14 pages, 13 figures, accepted for publication in A&A

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

Binding Energy Evaluation Platform: A database of quantum chemical binding energy distributions for the astrochemical community

Authors: Giulia M. Bovolenta, Stefan Vogt-Geisse, Stefano Bovino, Tommaso Grassi

Abstract: The quality of astrochemical models is highly dependent on reliable binding energy (BE) values that consider the morphological and energetic variety of binding sites on the surface of ice-grain mantles. Here, we present the Binding Energy Evaluation Platform (BEEP) and database that, using quantum chemical methods, produces full BE distributions of molecules bound to an amorphous solid water (ASW)… ▽ More The quality of astrochemical models is highly dependent on reliable binding energy (BE) values that consider the morphological and energetic variety of binding sites on the surface of ice-grain mantles. Here, we present the Binding Energy Evaluation Platform (BEEP) and database that, using quantum chemical methods, produces full BE distributions of molecules bound to an amorphous solid water (ASW) surface model. BEEP is highly automatized and allows to sample binding sites on set of water clusters and to compute accurate BEs. Using our protocol, we computed 21 BE distributions of interstellar molecules and radicals on an amorphized set of 15-18 water clusters of 22 molecules each. The distributions contain between 225 and 250 unique binding sites. We apply a Gaussian fit and report the mean and standard deviation for each distribution. We compare with existing experimental results and find that the low and high coverage experimental BEs coincide well with the high BE tail and mean value of our distributions, respectively. Previously reported single BE theoretical values are broadly in line with ours, even though in some cases significant differences can be appreciated. We show how the use of different BE values impact a typical problem in astrophysics, such as the computation of snow lines in protoplanetary discs. BEEP will be publicly released so that the database can be expanded to other molecules or ice-models in a community effort. △ Less

Submitted 26 July, 2022; originally announced July 2022.

Comments: Accepted on ApJS

Broadband spectroscopy of astrophysical ice analogues: II. Optical constants of CO and CO$_2$ ices in the terahertz and infrared ranges

Authors: A. A. Gavdush, F. Kruczkiewicz, B. M. Giuliano, B. Müller, G. A. Komandin, T. Grassi, P. Theulé, K. I. Zaytsev, A. V. Ivlev, P. Caselli

Abstract: Context: Broadband optical constants of astrophysical ice analogues in the infrared (IR) and terahertz (THz) ranges are required for modeling the dust continuum emission and radiative transfer in dense and cold regions, where thick icy mantles are formed on the surface of dust grains. Aims: In this paper, the THz time-domain spectroscopy (TDS) and the Fourier-transform IR spectroscopy (FTIR) are c… ▽ More Context: Broadband optical constants of astrophysical ice analogues in the infrared (IR) and terahertz (THz) ranges are required for modeling the dust continuum emission and radiative transfer in dense and cold regions, where thick icy mantles are formed on the surface of dust grains. Aims: In this paper, the THz time-domain spectroscopy (TDS) and the Fourier-transform IR spectroscopy (FTIR) are combined to study optical constants of CO and CO$_2$ ices in the broad THz-IR spectral range. Methods: The measured ices are grown at cryogenic temperatures by gas deposition on a cold Si window. A method to quantify the broadband THz-IR optical constants of ices is developed based on the direct reconstruction of the complex refractive index of ices in the THz range from the TDS data, and the use of the Kramers-Kronig relation in the IR range for the reconstruction from the FTIR data. Uncertainties of the Kramers-Kronig relation are eliminated by merging the THz and IR spectra. The reconstructed THz-IR response is then analyzed using classical models of complex dielectric permittivity. Results: The complex refractive index of CO and CO$_2$ ices deposited at the temperature of $28$ K is obtained in the range of 0.3--12.0 THz. Based on the measured dielectric constants, opacities of the astrophysical dust with CO and CO$_2$ icy mantles are computed. Conclusions: The developed method can be used for a model-independent reconstruction of optical constants of various astrophysical ice analogs in a broad THz-IR range. Such data can provide important benchmarks to interpret the broadband observations from the existing and future ground-based facilities and space telescopes. The reported results will be useful to model sources that show a drastic molecular freeze-out, such as central regions of prestellar cores and mid-planes of protoplanetary disks, as well as CO and CO$_2$ snow lines in disks. △ Less

Submitted 22 July, 2022; originally announced July 2022.

Comments: Accepted for publication in A&A, 9 pages, 7 figures

Journal ref: A&A 667, A49 (2022)

Cosmic rays in molecular clouds probed by H$_{2}$ rovibrational lines -- Perspectives for the James Webb Space Telescope

Authors: Marco Padovani, Shmuel Bialy, Daniele Galli, Alexei V. Ivlev, Tommaso Grassi, Liam H. Scarlett, Una S. Rehill, Mark C. Zammit, Dmitry V. Fursa, Igor Bray

Abstract: Cosmic rays (CRs) at sub-TeV energies play a fundamental role in the chemical and dynamical evolution of molecular clouds, as they control the ionisation, dissociation, and excitation of H$_{2}$. Their characterisation is important both for the interpretation of observations and for the development of theoretical models. The methods used so far for estimating the CR ionisation rate ($ζ$) in molecu… ▽ More Cosmic rays (CRs) at sub-TeV energies play a fundamental role in the chemical and dynamical evolution of molecular clouds, as they control the ionisation, dissociation, and excitation of H$_{2}$. Their characterisation is important both for the interpretation of observations and for the development of theoretical models. The methods used so far for estimating the CR ionisation rate ($ζ$) in molecular clouds have several limitations due to uncertainties in the adopted chemical networks. We refine and extend the method proposed by Bialy (2020) to estimate $ζ$ by observing rovibrational transitions of H$_{2}$ at near-infrared wavelengths, which are mainly excited by secondary CR electrons. Combining models of interstellar CR propagation and attenuation with the calculation of the expected secondary electron spectrum and updated H$_{2}$ excitation cross sections by electron collisions, we derive the intensity of the four H$_{2}$ rovibrational transitions observable in dense, cold gas: (1-0)O(2), (1-0)Q(2), (1-0)S(0), and (1-0)O(4). The proposed method allows the estimation of $ζ$ for a given observed line intensity and H$_{2}$ column density. We are also able to deduce the shape of the low-energy CR proton spectrum impinging upon the molecular cloud. We present a look-up plot and a web-based application that can be used to constrain the low-energy spectral slope of the interstellar CR proton spectrum. We comment on the capability of the James Webb Space Telescope to detect these near-infrared H$_{2}$ lines, making it possible to derive for the first time spatial variation of $ζ$ in dense gas. Besides the implications for the interpretation of the chemical-dynamic evolution of a molecular cloud, it will be possible to test competing models of CR propagation and attenuation in the interstellar medium, as well as compare CR spectra in different Galactic regions. △ Less

Submitted 20 January, 2022; originally announced January 2022.

Comments: 14 pages, 11 figures, accepted by A&A

Journal ref: A&A 658, A189 (2022)

On the low ortho-to-para H2 ratio in star-forming filaments

Authors: Alessandro Lupi, Stefano Bovino, Tommaso Grassi

Abstract: The formation of stars and planetary systems is a complex phenomenon, which relies on the interplay of multiple physical processes. Nonetheless, it represents a crucial stage for our understanding of the Universe, and in particular of the conditions leading to the formation of key molecules (e.g. water) on comets and planets. Herschel observations demonstrated that stars form out of gaseous filame… ▽ More The formation of stars and planetary systems is a complex phenomenon, which relies on the interplay of multiple physical processes. Nonetheless, it represents a crucial stage for our understanding of the Universe, and in particular of the conditions leading to the formation of key molecules (e.g. water) on comets and planets. Herschel observations demonstrated that stars form out of gaseous filamentary structures in which the main constituent is molecular hydrogen (H$_2$). Depending on its nuclear spin H$_2$ can be found in two forms: `ortho' with parallel spins and `para' where the spins are anti-parallel. The relative ratio among these isomers, i.e. the ortho-to-para ratio (OPR), plays a crucial role in a variety of processes related to the thermodynamics of star-forming gas and to the fundamental chemistry affecting the deuteration of water in molecular clouds, commonly used to determine the origin of water in Solar System's bodies. Here, for the first time, we assess the evolution of the OPR starting from the warm neutral medium, by means of state-of-the-art three-dimensional magneto-hydrodynamic simulations of turbulent molecular clouds. Our results show that star-forming clouds exhibit a low OPR ($\ll 0.1$) already at moderate densities ($\sim$1000 cm$^{-3}$). We also constrain the cosmic rays ionisation rate, finding that $10^{-16}\,\rm s^{-1}$ is the lower limit required to explain the observations of diffuse clouds. Our results represent a step forward in the understanding of the star and planet formation process providing a robust determination of the chemical initial conditions for both theoretical and observational studies. △ Less

Submitted 4 January, 2022; v1 submitted 6 September, 2021; originally announced September 2021.

Comments: 12 pages, 9 figures, 2 tables, matches the published version

Presence of water on exomoons orbiting free-floating planets: a case study

Authors: Patricio Javier Ávila, Tommaso Grassi, Stefano Bovino, Andrea Chiavassa, Barbara Ercolano, Sebastian Oscar Danielache, Eugenio Simoncini

Abstract: A free-floating planet is a planetary-mass object that orbits around a non-stellar massive object (e.g. a brown dwarf) or around the Galactic Center. The presence of exomoons orbiting free-floating planets has been theoretically predicted by several models. Under specific conditions, these moons are able to retain an atmosphere capable of ensuring the long-term thermal stability of liquid water on… ▽ More A free-floating planet is a planetary-mass object that orbits around a non-stellar massive object (e.g. a brown dwarf) or around the Galactic Center. The presence of exomoons orbiting free-floating planets has been theoretically predicted by several models. Under specific conditions, these moons are able to retain an atmosphere capable of ensuring the long-term thermal stability of liquid water on their surface. We model this environment with a one-dimensional radiative-convective code coupled to a gas-phase chemical network including cosmic rays and ion-neutral reactions. We find that, under specific conditions and assuming stable orbital parameters over time, liquid water can be formed on the surface of the exomoon. The final amount of water for an Earth-mass exomonoon is smaller than the amount of water in Earth oceans, but enough to host the potential development of primordial life. The chemical equilibrium time-scale is controlled by cosmic rays, the main ionization driver in our model of the exomoon atmosphere. △ Less

Submitted 15 June, 2021; originally announced June 2021.

Comments: 12 pages, 6 figures

Journal ref: International Journal of Astrobiology (2021) 1-12

Establishing the evolutionary timescales of the massive star formation process through chemistry

Authors: G. Sabatini, S. Bovino, A. Giannetti, T. Grassi, J. Brand, E. Schisano, F. Wyrowski, S. Leurini, K. M. Menten

Abstract: (Abridged) Understanding the details of the formation process of massive (i.e. M<8-10M$_\odot$) stars is a long-standing problem in astrophysics. [...] We present a method to derive accurate timescales of the different evolutionary phases of the high-mass star formation process. We model a representative number of massive clumps of the ATLASGAL-TOP100 sample which cover all the evolutionary stages… ▽ More (Abridged) Understanding the details of the formation process of massive (i.e. M<8-10M$_\odot$) stars is a long-standing problem in astrophysics. [...] We present a method to derive accurate timescales of the different evolutionary phases of the high-mass star formation process. We model a representative number of massive clumps of the ATLASGAL-TOP100 sample which cover all the evolutionary stages. The models describe an isothermal collapse and the subsequent warm-up phase, for which we follow their chemical evolution. The timescale of each phase is derived by comparing the results of the models with the properties of the sources of the ATLASGAL-TOP100 sample, taking into account the mass and luminosity of the clumps, and the column densities of methyl acetylene (CH$_3$CCH), acetonitrile (CH$_3$CN), formaldehyde (H$_2$CO) and methanol (CH$_3$OH). We find that the chosen molecular tracers are affected by the thermal evolution of the clumps, showing steep ice evaporation gradients from 10$^3$ to 10$^5$ AU during the warm-up phase. We succeed in reproducing the observed column densities of CH$_3$CCH and CH$_3$CN, while H$_2$CO and CH$_3$OH show a poorer agreement with the observed values. The total (massive) star formation time is found to be $\sim5.2\times10^5$ yr, which is defined by the timescales of the individual evolutionary phases of the ATLASGAL-TOP100 sample: $\sim5\times10^4$ yr for 70-$μ$m weak, $\sim1.2\times10^5$ yr for mid-IR weak, $\sim2.4\times10^5$ yr for mid-IR bright and $\sim1.1\times10^5$ yr for HII-regions phases. Our models, with an appropriate selection of molecular tracers that can act as chemical clocks, allow to get robust estimates of the duration of the individual phases of the high-mass star formation process, with the advantage of being capable to include additional tracers aimed at increasing the accuracy of the estimated timescales. △ Less

Submitted 13 October, 2021; v1 submitted 1 June, 2021; originally announced June 2021.

Comments: Published on A&A (19 pages, 9 figures, 7 tables)

Journal ref: A&A 652, A71 (2021)

Reducing the complexity of chemical networks via interpretable autoencoders

Authors: T. Grassi, F. Nauman, J. P. Ramsey, S. Bovino, G. Picogna, B. Ercolano

Abstract: In many astrophysical applications, the cost of solving a chemical network represented by a system of ordinary differential equations (ODEs) grows significantly with the size of the network, and can often represent a significant computational bottleneck, particularly in coupled chemo-dynamical models. Although standard numerical techniques and complex solutions tailored to thermochemistry can some… ▽ More In many astrophysical applications, the cost of solving a chemical network represented by a system of ordinary differential equations (ODEs) grows significantly with the size of the network, and can often represent a significant computational bottleneck, particularly in coupled chemo-dynamical models. Although standard numerical techniques and complex solutions tailored to thermochemistry can somewhat reduce the cost, more recently, machine learning algorithms have begun to attack this challenge via data-driven dimensional reduction techniques. In this work, we present a new class of methods that take advantage of machine learning techniques to reduce complex data sets (autoencoders), the optimization of multi-parameter systems (standard backpropagation), and the robustness of well-established ODE solvers to to explicitly incorporate time-dependence. This new method allows us to find a compressed and simplified version of a large chemical network in a semi-automated fashion that can be solved with a standard ODE solver, while also enabling interpretability of the compressed, latent network. As a proof of concept, we tested the method on an astrophysically-relevant chemical network with 29 species and 224 reactions, obtaining a reduced but representative network with only 5 species and 12 reactions, and a x65 speed-up. △ Less

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

Comments: accepted A&A, code available at https://bitbucket.org/tgrassi/latent_ode_paper/

Modelling H$_2$ and its effects on star formation using a joint implementation of GADGET-3 and KROME

Authors: Emanuel Sillero, Patricia B. Tissera, Diego G. Lambas, Stefano Bovino, Dominik R. Schleicher, Tommaso Grassi, Gustavo Bruzual, Stéphane Charlot

Abstract: We present P-GADGET3-K, an updated version of GADGET3, that incorporates the chemistry package KROME. P-GADGET3-K follows the hydrodynamical and chemical evolution of cosmic structures, incorporating the chemistry and cooling of H$_2$ and metal cooling in non-equilibrium. We performed different runs of the same ICs to assess the impact of various physical parameters and prescriptions, namely gas m… ▽ More We present P-GADGET3-K, an updated version of GADGET3, that incorporates the chemistry package KROME. P-GADGET3-K follows the hydrodynamical and chemical evolution of cosmic structures, incorporating the chemistry and cooling of H$_2$ and metal cooling in non-equilibrium. We performed different runs of the same ICs to assess the impact of various physical parameters and prescriptions, namely gas metallicity, molecular hydrogen formation on dust, star formation recipes including or not H$_2$ dependence, and the effects of numerical resolution. We find that the characteristics of the simulated systems, both globally and at kpc-scales, are in good agreement with several observable properties of molecular gas in star-forming galaxies. The surface density profiles of SFR and H$_2$ are found to vary with the clumping factor and resolution. In agreement with previous results, the chemical enrichment of the gas component is found to be a key ingredient to model the formation and distribution of H$_2$ as a function of gas density and temperature. A SF algorithm that takes into account the H$_2$ fraction together with a treatment for the local stellar radiation field improves the agreement with observed H$_2$ abundances over a wide range of gas densities and with the molecular Kennicutt-Schmidt law, implying a more realistic modelling of the star formation process. △ Less

Submitted 9 April, 2021; originally announced April 2021.

Comments: 21 pages, 16 figures. Accepted for publication in MNRAS

High level ab initio binding energy distribution of molecules on interstellar ices: Hydrogen fluoride

Authors: Giulia Bovolenta, Stefano Bovino, Esteban Vöhringer-Martinez, David A. Saez, Tommaso Grassi, Stefan Vogt-Geisse

Abstract: The knowledge of the binding energy of molecules on astrophysically relevant ices can help to obtain an estimate of the desorption rate, i.e. the molecules residence time on the surface. This represents an important parameter for astrochemical models, crucial to determine the chemical fate of complex organic molecules formed on dust grains and observed in the densest regions of the interstellar me… ▽ More The knowledge of the binding energy of molecules on astrophysically relevant ices can help to obtain an estimate of the desorption rate, i.e. the molecules residence time on the surface. This represents an important parameter for astrochemical models, crucial to determine the chemical fate of complex organic molecules formed on dust grains and observed in the densest regions of the interstellar medium. In this work, we propose a new robust procedure to study the interaction of atoms and molecules with interstellar ices, based on \textit{ab initio} molecular dynamics and density functional theory, validated by high-level \textit{ab initio} methods at a CCSD(T)/CBS level. We have applied this procedure to hydrogen fluoride (HF), a promising tracer of the molecular content of galaxies. In total we found 13 unique equilibrium structures of HF binding to small water clusters of up to 4 molecules, with binding energies ranging from 1208 to 7162 K. We computed a 22-molecules model of amorphous solid water (ASW) surface using \textit{ab initio} molecular dynamics simulations and carried out a systematic analysis of the binding sites of HF, in terms of binding modes and binding energies. Considering 10 different water clusters, we found a binding energy distribution with an average value of $5313\pm74$ K, and a dispersion of $921\pm115$ K. Finally, the effect of the electrostatic field of the 22 water molecules on the binding energies was investigated incrementally by symmetry adapted perturbation theory, in order to gauge the effect of the water environment. The results indicate that the extent of the electrostatic interaction of HF with ASW depends strongly on the properties of the binding site. We expect that this work will provide a solid foundation for a systematic development of a binding energy distribution database of molecules on interstellar surfaces. △ Less

Submitted 18 October, 2020; originally announced October 2020.

Comments: 12 pages, 9 figures, 6 tables, accepted to Mol. Astrophys

A novel framework to study the impact of binding energy distributions on the chemistry of dust grains

Authors: T. Grassi, S. Bovino, P. Caselli, G. Bovolenta, S. Vogt-Geisse, B. Ercolano

Abstract: The evaporation of molecules from dust grains is crucial to understand some key aspects of the star- and the planet-formation processes. During the warm-up phase the presence of young protostellar objects induces molecules to evaporate from the dust surface into the gas phase, enhancing its chemical complexity. Similarly, in circumstellar disks, the position of the so-called snow-lines is determin… ▽ More The evaporation of molecules from dust grains is crucial to understand some key aspects of the star- and the planet-formation processes. During the warm-up phase the presence of young protostellar objects induces molecules to evaporate from the dust surface into the gas phase, enhancing its chemical complexity. Similarly, in circumstellar disks, the position of the so-called snow-lines is determined by evaporation, with important consequences for the formation of planets. The amount of molecules that are desorbed depends on the interaction between the species and the grain surface, which is controlled by the binding energy. Recent theoretical and experimental works point towards a distribution of values for this parameter instead of the single value often employed in astrochemical models.We present here a new "multi-binding energy" framework, to assess the effects that a distribution of binding energies has on the amount of species bound to the grains. We find that the efficiency of the surface chemistry is significantly influenced by this process with crucial consequences on the theoretical estimates of the desorbed species. △ Less

Submitted 8 September, 2020; originally announced September 2020.

Comments: Accepted A&A

Modelling thermochemical processes in protoplanetary disks I: numerical methods

Authors: T. Grassi, B. Ercolano, L. Szűcs, J. Jennings, G. Picogna

Abstract: The dispersal phase of planet-forming disks via winds driven by irradiation from the central star and/or magnetic fields in the disk itself is likely to play an important role in the formation and evolution of planetary systems. Current theoretical models lack predictive power to adequately constrain observations. We present PRIZMO, a code for evolving thermochemistry in protoplanetary disks capab… ▽ More The dispersal phase of planet-forming disks via winds driven by irradiation from the central star and/or magnetic fields in the disk itself is likely to play an important role in the formation and evolution of planetary systems. Current theoretical models lack predictive power to adequately constrain observations. We present PRIZMO, a code for evolving thermochemistry in protoplanetary disks capable of being coupled with hydrodynamical and multi-frequency radiative transfer codes. We describe the main features of the code, including gas and surface chemistry, photochemistry, microphysics, and the main cooling and heating processes. The results of a suite of benchmarks, which include photon-dominated regions, slabs illuminated by radiation spectra that include X-ray, and well-established cooling functions evaluated at different temperatures show good agreement both in terms of chemical and thermal structures. The development of this code is an important step to perform quantitative spectroscopy of disk winds, and ultimately the calculation of line profiles, which is urgently needed to shed light on the nature of observed disk winds. △ Less

Submitted 9 April, 2020; originally announced April 2020.

Comments: MNRAS accepted

Broadband spectroscopy of astrophysical ice analogs. I. Direct measurement of complex refractive index of CO ice using terahertz time-domain spectroscopy

Authors: B. M. Giuliano, A. A. Gavdush, B. Müller, K. I. Zaytsev, T. Grassi, A. V. Ivlev, M. E. Palumbo, G. A. Baratta, C. Scirè, G. A. Komandin, S. O. Yurchenko, P. Caselli

Abstract: Context: Reliable, directly measured optical properties of astrophysical ice analogs in the infrared (IR) and terahertz (THz) range are missing. These parameters are of great importance to model the dust continuum radiative transfer in dense and cold regions, here thick ice mantles are present, and are necessary for the interpretation of future observations planned in the far-IR region. Aims: Cohe… ▽ More Context: Reliable, directly measured optical properties of astrophysical ice analogs in the infrared (IR) and terahertz (THz) range are missing. These parameters are of great importance to model the dust continuum radiative transfer in dense and cold regions, here thick ice mantles are present, and are necessary for the interpretation of future observations planned in the far-IR region. Aims: Coherent THz radiation allows direct measurement of the complex dielectric function (refractive index) of astrophysically relevant ice species in the THz range. Methods: The time-domain waveforms and the frequency-domain spectra of reference samples of CO ice, deposited at a temperature of 28.5 K and annealed to 33 K at different thicknesses, have been recorded. A new algorithm is developed to reconstruct the real and imaginary parts of the refractive index from the time-domain THz data. Results: The complex refractive index in the wavelength range of 1 mm - 150 $μ$m (0.3 - 2.0 THz) has been determined for the studied ice samples, and compared with available data found in the literature. Conclusions: The developed algorithm of reconstructing the real and imaginary parts of the refractive index from the time-domain THz data enables, for the first time, the determination of optical properties of astrophysical ice analogs without using the Kramers-Kronig relations. The obtained data provide a benchmark to interpret the observational data from current ground based facilities as well as future space telescope missions, and have been used to estimate the opacities of the dust grains in presence of CO ice mantles. △ Less

Submitted 26 July, 2019; originally announced July 2019.

Journal ref: A&A 629, A112 (2019)

The challenges of modelling microphysics: ambipolar diffusion, chemistry, and cosmic rays in MHD shocks

Authors: T. Grassi, M. Padovani, J. P. Ramsey, D. Galli, N. Vaytet, B. Ercolano, T. Haugboelle

Abstract: From molecular clouds to protoplanetary disks, non-ideal magnetic effects are important in many astrophysical environments. Indeed, in star and disk formation processes, it has become clear that these effects are critical to the evolution of the system. The efficacy of non-ideal effects are, however, determined by the complex interplay between magnetic fields, ionising radiation, cosmic rays, micr… ▽ More From molecular clouds to protoplanetary disks, non-ideal magnetic effects are important in many astrophysical environments. Indeed, in star and disk formation processes, it has become clear that these effects are critical to the evolution of the system. The efficacy of non-ideal effects are, however, determined by the complex interplay between magnetic fields, ionising radiation, cosmic rays, microphysics, and chemistry. In order to understand these key microphysical parameters, we present a one-dimensional non-ideal magnetohydrodynamics code and apply it to a model of a time-dependent, oblique, magnetic shock wave. By varying the microphysical ingredients of the model, we find that cosmic rays and dust play a major role, and that, despite the uncertainties, the inclusion of microphysics is essential to obtain a realistic outcome in magnetic astrophysical simulations. △ Less

Submitted 2 January, 2019; originally announced January 2019.

Comments: MNRAS accepted

Lampray: Multi-group long characteristics ray tracing for adaptive mesh radiation hydrodynamics

Authors: Troels Frostholm, Troels Haugbølle, Tommaso Grassi

Abstract: We present Lampray: a multi-group long characteristics ray tracing method for adaptive mesh radiation hydrodynamics in the Ramses code. It avoids diffusion, captures shadows, and treats colliding beams correctly, and therefore complements existing moment-based ray tracing in Ramses. Lampray includes different options for interpolation between ray and cell domain, and use either integral, Fourier,… ▽ More We present Lampray: a multi-group long characteristics ray tracing method for adaptive mesh radiation hydrodynamics in the Ramses code. It avoids diffusion, captures shadows, and treats colliding beams correctly, and therefore complements existing moment-based ray tracing in Ramses. Lampray includes different options for interpolation between ray and cell domain, and use either integral, Fourier, or an implicit method for hydrogen ionization to solve the radiative transfer. The opacity can either be tabulated or computed through a coupling to the general non-equilibrium astro-chemistry framework Krome. We use an H-He-C-O network with 36 species and 240 reactions to track the photo-chemistry in the interstellar medium across 6 and 10 orders of magnitude in temperature and density. Self-shielding prescriptions for H$_2$ and CO are used together with a new model for the diffuse interstellar UV-field. We also track the dust temperature, formation of H$_2$ on grains, and H$_2$O and CO ices in detail. Lampray is tested against standard benchmarks for molecular cloud and star formation simulations, including the formation of a Strömgren sphere, the expansion of an ionization front, the photo-evaporation of a dense clump, and the H-He-C-O chemistry in a static photo-dissociation front. Efficient parallelisation is achieved with a separate domain decomposition for rays where points along a ray reside in the same memory space, and data movement from cell- to ray-domain is done with a direct hash-table lookup algorithm. Point sources are treated without splitting rays, and therefore the method currently only scales to a few point sources, while diffuse radiation has excellent scaling. △ Less

Submitted 14 September, 2018; originally announced September 2018.

Comments: 31 pages, 38 figures

The effect of non-equilibrium metal cooling on the interstellar medium

Authors: Pedro R. Capelo, Stefano Bovino, Alessandro Lupi, Dominik R. G. Schleicher, Tommaso Grassi

Abstract: By using a novel interface between the modern smoothed particle hydrodynamics code GASOLINE2 and the chemistry package KROME, we follow the hydrodynamical and chemical evolution of an isolated galaxy. In order to assess the relevance of different physical parameters and prescriptions, we constructed a suite of ten simulations, in which we vary the chemical network (primordial and metal species), h… ▽ More By using a novel interface between the modern smoothed particle hydrodynamics code GASOLINE2 and the chemistry package KROME, we follow the hydrodynamical and chemical evolution of an isolated galaxy. In order to assess the relevance of different physical parameters and prescriptions, we constructed a suite of ten simulations, in which we vary the chemical network (primordial and metal species), how metal cooling is modelled (non-equilibrium versus equilibrium; optically thin versus thick approximation), the initial gas metallicity (from ten to hundred per cent solar), and how molecular hydrogen forms on dust. This is the first work in which metal injection from supernovae, turbulent metal diffusion, and a metal network with non-equilibrium metal cooling are self-consistently included in a galaxy simulation. We find that properly modelling the chemical evolution of several metal species and the corresponding non-equilibrium metal cooling has important effects on the thermodynamics of the gas, the chemical abundances, and the appearance of the galaxy: the gas is typically warmer, has a larger molecular gas mass fraction, and has a smoother disc. We also conclude that, at relatively high metallicity, the choice of molecular-hydrogen formation rates on dust is not crucial. Moreover, we confirm that a higher initial metallicity produces a colder gas and a larger fraction of molecular gas, with the low-metallicity simulation best matching the observed molecular Kennicutt-Schmidt relation. Finally, our simulations agree quite well with observations which link star formation rate to metal emission lines. △ Less

Submitted 4 January, 2018; v1 submitted 3 October, 2017; originally announced October 2017.

Comments: Accepted for publication in MNRAS

Astrochemistry: From primordial gas to present-day clouds

Authors: Dominik R. G. Schleicher, Stefano Bovino, Bastian Körtgen, Tommaso Grassi, Robi Banerjee

Abstract: Astrochemistry plays a central role during the process of star formation, both in the primordial regime as well as in the present-day Universe. We revisit here the chemistry in both regimes, focusing first on the chemistry under close to primordial conditions, as observed in the so-called Caffau star SDSS J102915+172927, and subsequently discuss deuteration processes in present-day star-forming co… ▽ More Astrochemistry plays a central role during the process of star formation, both in the primordial regime as well as in the present-day Universe. We revisit here the chemistry in both regimes, focusing first on the chemistry under close to primordial conditions, as observed in the so-called Caffau star SDSS J102915+172927, and subsequently discuss deuteration processes in present-day star-forming cores. In models of the high-redshift Universe, the chemistry is particularly relevant to determine the cooling, while it also serves as an important diagnostic in the case of present-day star formation. △ Less

Submitted 23 August, 2017; originally announced August 2017.

Comments: 3 pages, 2 figures. Proceedings article for "Francesco's Legacy: Star Formation in Space and Time", to be published in "Memorie della Società Astronomica Italiana"

H$_2$ ortho-to-para conversion on grains: A route to fast deuterium fractionation in dense cloud cores?

Authors: Stefano Bovino, Tommaso Grassi, Dominik R. Schleicher, Paola Caselli

Abstract: Deuterium fractionation, i.e. the enhancement of deuterated species with respect to the non-deuterated ones, is considered to be a reliable chemical clock of star-forming regions. This process is strongly affected by the ortho-to-para (o-p) H$_2$ ratio. In this letter we explore the effect of the o-p H$_2$ conversion on grains on the deuteration timescale in fully depleted dense cores, including t… ▽ More Deuterium fractionation, i.e. the enhancement of deuterated species with respect to the non-deuterated ones, is considered to be a reliable chemical clock of star-forming regions. This process is strongly affected by the ortho-to-para (o-p) H$_2$ ratio. In this letter we explore the effect of the o-p H$_2$ conversion on grains on the deuteration timescale in fully depleted dense cores, including the most relevant uncertainties that affect this complex process. We show that (i) the o-p H$_2$ conversion on grains is not strongly influenced by the uncertainties on the conversion time and the sticking coefficient and (ii) that the process is controlled by the temperature and the residence time of ortho-H$_2$ on the surface, i.e. by the binding energy. We find that for binding energies in between 330-550 K, depending on the temperature, the o-p H$_2$ conversion on grains can shorten the deuterium fractionation timescale by orders of magnitude, opening a new route to explain the large observed deuteration fraction $D_\mathrm{frac}$ in dense molecular cloud cores. Our results suggest that the star formation timescale, when estimated through the timescale to reach the observed deuteration fractions, might be shorter than previously proposed. However, more accurate measurements of the binding energy are needed to better assess the overall role of this process. △ Less

Submitted 23 October, 2017; v1 submitted 7 August, 2017; originally announced August 2017.

Comments: Accepted for publication in ApJ Letters

Modeling the role of electron attachment rates on column density ratios for CnH-/CnH (n=4,6,8) in dense molecular clouds

Authors: F. A. Gianturco, T. Grassi, R. Wester

Abstract: (abridged) The fairly recent detection of a variety of anions in the Interstellar Molecular Clouds have underlined the importance of realistically modeling the processes governing their abundance. To this aim, our earlier calculations for the radiative electron attachment (REA) rates for C4H-, C6H-, and C8H- are employed to generate the corresponding column density ratios of anion/neutral (A/N) re… ▽ More (abridged) The fairly recent detection of a variety of anions in the Interstellar Molecular Clouds have underlined the importance of realistically modeling the processes governing their abundance. To this aim, our earlier calculations for the radiative electron attachment (REA) rates for C4H-, C6H-, and C8H- are employed to generate the corresponding column density ratios of anion/neutral (A/N) relative abundances. The latter are then compared with those obtained from observational measurements. The calculations involved the time-dependent solutions of a large network of chemical processes over an extended time interval and included a series of runs in which the values of REA rates were repeatedly scaled. Macroscopic parameters for the clouds' modeling were also varied to cover a broad range of physical environments. It was found that, within the range and quality of the processes included in the present network,and selected from state-of-the-art astrophysical databases, the REA values required to match the observed A/N ratios needed to be reduced by orders of magnitude for C4H- case, while the same rates for C6H- and C8H- only needed to be scaled by much smaller factors. The results suggest that the generally proposed formation of interstellar anions by REA mechanism is overestimated by current models for the C4H- case, for which is likely to be an inefficient path to formation. This path is thus providing a rather marginal contribution to the observed abundances of C4H-, the latter being more likely to originate from other chemical processes in the network, as we discuss in some detail in the present work.Possible physical reasons for the much smaller differences against observations found instead for the values of the (A/N) ratios in two other, longer members of the series are put forward and analyzed within the evolutionary modeling discussed in the present work. △ Less

Submitted 19 July, 2016; originally announced July 2016.

Comments: Journal of Physics B, accepted

A detailed framework to incorporate dust in hydrodynamical simulations

Authors: T. Grassi, S. Bovino, T. Haugboelle, D. R. G. Schleicher

Abstract: Dust plays a key role in the evolution of the ISM and its correct modelling in numerical simulations is therefore fundamental. We present a new and self-consistent model that treats grain thermal coupling with the gas, radiation balance, and surface chemistry for molecular hydrogen. This method can be applied to any dust distribution with an arbitrary number of grain types without affecting the ov… ▽ More Dust plays a key role in the evolution of the ISM and its correct modelling in numerical simulations is therefore fundamental. We present a new and self-consistent model that treats grain thermal coupling with the gas, radiation balance, and surface chemistry for molecular hydrogen. This method can be applied to any dust distribution with an arbitrary number of grain types without affecting the overall computational cost. In this paper we describe in detail the physics and the algorithm behind our approach, and in order to test the methodology, we present some examples of astrophysical interest, namely (i) a one-zone collapse with complete gas chemistry and thermochemical processes, (ii) a 3D model of a low-metallicity collapse of a minihalo starting from cosmological initial conditions, and (iii) a turbulent molecular cloud with H-C-O chemistry (277 reactions), together with self-consistent cooling and heating solved on the fly. Although these examples employ the publicly available code KROME, our approach can be easily integrated into any computational framework. △ Less

Submitted 4 November, 2016; v1 submitted 3 June, 2016; originally announced June 2016.

Comments: 17 pages, 12 figures, MNRAS accepted

The formation of the primitive star SDSS J102915+172927: effect of the dust mass and the grain-size distribution

Authors: S. Bovino, T. Grassi, D. R. G. Schleicher, R. Banerjee

Abstract: Understanding the formation of the extremely metal poor star SDSS-J102915+172927 is of fundamental importance to improve our knowledge on the transition between the first and second generation of stars in the Universe. In this paper, we perform three-dimensional cosmological hydrodynamical simulations of dust-enriched halos during the early stages of the collapse process including a detailed treat… ▽ More Understanding the formation of the extremely metal poor star SDSS-J102915+172927 is of fundamental importance to improve our knowledge on the transition between the first and second generation of stars in the Universe. In this paper, we perform three-dimensional cosmological hydrodynamical simulations of dust-enriched halos during the early stages of the collapse process including a detailed treatment of the dust physics. We employ the astrochemistry package \krome coupled with the hydrodynamical code \textsc{enzo} assuming grain size distributions produced by the explosion of core-collapse supernovae of 20 and 35 M$_\odot$ primordial stars which are suitable to reproduce the chemical pattern of the SDSS-J102915+172927 star. We find that the dust mass yield produced from Population III supernovae explosions is the most important factor which drives the thermal evolution and the dynamical properties of the halos. Hence, for the specific distributions relevant in this context, the composition, the dust optical properties, and the size-range have only minor effects on the results due to similar cooling functions. We also show that the critical dust mass to enable fragmentation provided by semi-analytical models should be revised, as we obtain values one order of magnitude larger. This determines the transition from disk fragmentation to a more filamentary fragmentation mode, and suggests that likely more than one single supernova event or efficient dust growth should be invoked to get such a high dust content. △ Less

Submitted 7 September, 2016; v1 submitted 18 January, 2016; originally announced January 2016.

Comments: Accepted on ApJ

A chemical model for the interstellar medium in galaxies

Authors: Stefano Bovino, Tommaso Grassi, Pedro R. Capelo, Dominik R. G. Schleicher, R. Banerjee

Abstract: We present and test chemical models for three-dimensional hydrodynamical simulations of galaxies. We explore the effect of changing key parameters such as metallicity, radiation and non-equilibrium versus equilibrium metal cooling approximations on the transition between the gas phases in the interstellar medium. The microphysics is modelled by employing the public chemistry package KROME and the… ▽ More We present and test chemical models for three-dimensional hydrodynamical simulations of galaxies. We explore the effect of changing key parameters such as metallicity, radiation and non-equilibrium versus equilibrium metal cooling approximations on the transition between the gas phases in the interstellar medium. The microphysics is modelled by employing the public chemistry package KROME and the chemical networks have been tested to work in a wide range of densities and temperatures. We describe a simple H/He network following the formation of H$_2$, and a more sophisticated network which includes metals. Photochemistry, thermal processes, and different prescriptions for the H$_2$ catalysis on dust are presented and tested within a one-zone framework. The resulting network is made publicly available on the KROME webpage. We find that employing an accurate treatment of the dust-related processes induces a faster HI--H$_2$ transition. In addition, we show when the equilibrium assumption for metal cooling holds, and how a non-equilibrium approach affects the thermal evolution of the gas and the HII--HI transition. These models can be employed in any hydrodynamical code via an interface to KROME and can be applied to different problems including isolated galaxies, cosmological simulations of galaxy formation and evolution, supernova explosions in molecular clouds, and the modelling of star-forming regions. The metal network can be used for a comparison with observational data of CII 158 $μ$m emission both for high-redshift as well as for local galaxies. △ Less

Submitted 15 February, 2016; v1 submitted 23 October, 2015; originally announced October 2015.

Comments: A&A accepted

Journal ref: A&A 590, A15 (2016)

The chemical evolution of self-gravitating primordial disks

Authors: Dominik R. G. Schleicher, Stefano Bovino, Muhammad A. Latif, Andrea Ferrara, Tommaso Grassi

Abstract: Numerical simulations show the formation of self-gravitating primordial disks during the assembly of the first structures in the Universe, in particular during the formation of Pop.~III and supermassive stars. Their subsequent evolution is expected to be crucial to determine the mass scale of the first cosmological objects, which depends on the temperature of the gas and the dominant cooling mecha… ▽ More Numerical simulations show the formation of self-gravitating primordial disks during the assembly of the first structures in the Universe, in particular during the formation of Pop.~III and supermassive stars. Their subsequent evolution is expected to be crucial to determine the mass scale of the first cosmological objects, which depends on the temperature of the gas and the dominant cooling mechanism. Here, we derive a one-zone framework to explore the chemical evolution of such disks and show that viscous heating leads to the collisional dissociation of an initially molecular gas. The effect is relevant on scales of 10 AU (1000 AU) for a central mass of 10 M_sun (10^4 M_sun) at an accretion rate of 10^{-1} M_sun yr^{-1}, and provides a substantial heat input to stabilize the disk. If the gas is initially atomic, it remains atomic during the further evolution, and the effect of viscous heating is less significant. The additional thermal support is particularly relevant for the formation of very massive objects, such as the progenitors of the first supermassive black holes. The stabilizing impact of viscous heating thus alleviates the need for a strong radiation background as a means of keeping the gas atomic. △ Less

Submitted 8 October, 2015; v1 submitted 23 April, 2015; originally announced April 2015.

Comments: 13 pages, 5 figures, 6 tables, accepted at A&A

Journal ref: A&A 585, A11 (2016)

How realistic UV spectra and X-rays suppress the abundance of direct collapse black holes

Authors: M. A. Latif, S. Bovino, T. Grassi, D. R. G. Schleicher, M. Spaans

Abstract: Observations of high redshift quasars at $z>6$ indicate that they harbor supermassive black holes (SMBHs) of a billion solar masses. The direct collapse scenario has emerged as the most plausible way to assemble SMBHs. The nurseries for the direct collapse black holes are massive primordial halos illuminated with an intense UV flux emitted by population II (Pop II) stars. In this study, we compute… ▽ More Observations of high redshift quasars at $z>6$ indicate that they harbor supermassive black holes (SMBHs) of a billion solar masses. The direct collapse scenario has emerged as the most plausible way to assemble SMBHs. The nurseries for the direct collapse black holes are massive primordial halos illuminated with an intense UV flux emitted by population II (Pop II) stars. In this study, we compute the critical value of such a flux ($J_{21}^{\rm crit}$) for realistic spectra of Pop II stars through three-dimensional cosmological simulations. We derive the dependence of $J_{21}^{\rm crit}$ on the radiation spectra, on variations from halo to halo, and on the impact of X-ray ionization. Our findings show that the value of $J_{21}^{\rm crit}$ is a few times $\rm 10^4$ and only weakly depends on the adopted radiation spectra in the range between $T_{\rm rad}=2 \times 10^4-10^5$ K. For three simulated halos of a few times $\rm 10^{7}$~M$_{\odot}$, $J_{21}^{\rm crit}$ varies from $\rm 2 \times 10^4 - 5 \times 10^4$. The impact of X-ray ionization is almost negligible and within the expected scatter of $J_{21}^{\rm crit}$ for background fluxes of $J_{\rm X,21} \leq 0.1$. The computed estimates of $J_{21}^{\rm crit}$ have profound implications for the quasar abundance at $z=10$ as it lowers the number density of black holes forming through an isothermal direct collapse by a few orders of magnitude below the observed black holes density. However, the sites with moderate amounts of $\rm H_2$ cooling may still form massive objects sufficient to be compatible with observations. △ Less

Submitted 31 October, 2014; v1 submitted 13 August, 2014; originally announced August 2014.

Comments: Accepted for publication in MNRAS, comments are welcome

Effects of turbulence and rotation on protostar formation as a precursor to seed black holes

Authors: C. Van Borm, S. Bovino, M. A. Latif, D. R. G. Schleicher, M. Spaans, T. Grassi

Abstract: Context. The seeds of the first supermassive black holes may have resulted from the direct collapse of hot primordial gas in $\gtrsim 10^4$ K haloes, forming a supermassive or quasistar as an intermediate stage. Aims. We explore the formation of a protostar resulting from the collapse of primordial gas in the presence of a strong Lyman-Werner radiation background. Particularly, we investigate th… ▽ More Context. The seeds of the first supermassive black holes may have resulted from the direct collapse of hot primordial gas in $\gtrsim 10^4$ K haloes, forming a supermassive or quasistar as an intermediate stage. Aims. We explore the formation of a protostar resulting from the collapse of primordial gas in the presence of a strong Lyman-Werner radiation background. Particularly, we investigate the impact of turbulence and rotation on the fragmentation behaviour of the gas cloud. We accomplish this goal by varying the initial turbulent and rotational velocities. Methods. We performed 3D adaptive mesh refinement simulations with a resolution of 64 cells per Jeans length using the ENZO code, simulating the formation of a protostar up to unprecedentedly high central densities of $10^{21}$ cm$^{-3}$, and spatial scales of a few solar radii. To achieve this goal, we employed the KROME package to improve modelling of the chemical and thermal processes. Results. We find that the physical properties of the simulated gas clouds become similar on small scales, irrespective of the initial amount of turbulence and rotation. After the highest level of refinement was reached, the simulations have been evolved for an additional ~5 freefall times. A single bound clump with a radius of $2 \times 10^{-2}$ AU and a mass of ~$7 \times 10^{-2}$ M$_{\odot}$ is formed at the end of each simulation, marking the onset of protostar formation. No strong fragmentation is observed by the end of the simulations, regardless of the initial amount of turbulence or rotation, and high accretion rates of a few solar masses per year are found. Conclusions. Given such high accretion rates, a quasistar of $10^5$ M$_{\odot}$ is expected to form within $10^5$ years. △ Less

Submitted 29 October, 2014; v1 submitted 6 August, 2014; originally announced August 2014.

Comments: 18 pages, 7 figures, fixed typos, added references and clarified some details; accepted for publication in A&A

Journal ref: A&A 572, A22 (2014)

Formation of carbon-enhanced metal-poor stars in the presence of far ultraviolet radiation

Authors: S. Bovino, T. Grassi, D. R. G. Schleicher, M. A. Latif

Abstract: Recent discoveries of carbon-enhanced metal-poor stars like SMSS J031300.36-670839.3 provide increasing observational insights into the formation conditions of the first second-generation stars in the Universe, reflecting the chemical conditions after the first supernova explosion. Here, we present the first cosmological simulations with a detailed chemical network including primordial species as… ▽ More Recent discoveries of carbon-enhanced metal-poor stars like SMSS J031300.36-670839.3 provide increasing observational insights into the formation conditions of the first second-generation stars in the Universe, reflecting the chemical conditions after the first supernova explosion. Here, we present the first cosmological simulations with a detailed chemical network including primordial species as well as C, C$^+$, O, O$^+$, Si, Si$^+$, and Si$^{2+}$ following the formation of carbon-enhanced metal poor stars. The presence of background UV flux delays the collapse from $z=21$ to $z=15$ and cool the gas down to the CMB temperature for a metallicity of Z/Z$_\odot$=10$^{-3}$. This can potentially lead to the formation of lower mass stars. Overall, we find that the metals have a stronger effect on the collapse than the radiation, yielding a comparable thermal structure for large variations in the radiative background. We further find that radiative backgrounds are not able to delay the collapse for Z/Z$_\odot$=10$^{-2}$ or a carbon abundance as in SMSS J031300.36-670839.3. △ Less

Submitted 12 June, 2014; originally announced June 2014.

Comments: submitted to ApJL

The formation of massive primordial stars in the presence of moderate UV backgrounds

Authors: M. A. Latif, D. R. G. Schleicher, S. Bovino, T. Grassi, M. Spaans

Abstract: Radiative feedback from populations II stars played a vital role in early structure formation. Particularly, photons below the Lyman limit can escape the star forming regions and produce a background ultraviolet (UV) flux which consequently may influence the pristine halos far away from the radiation sources. These photons can quench the formation of molecular hydrogen by photo-detachment of… ▽ More Radiative feedback from populations II stars played a vital role in early structure formation. Particularly, photons below the Lyman limit can escape the star forming regions and produce a background ultraviolet (UV) flux which consequently may influence the pristine halos far away from the radiation sources. These photons can quench the formation of molecular hydrogen by photo-detachment of $\rm H^{-}$. In this study, we explore the impact of such UV radiation on fragmentation in massive primordial halos of a few times $\rm 10^{7}$~M${_\odot}$. To accomplish this goal, we perform high resolution cosmological simulations for two distinct halos and vary the strength of the impinging background UV field in units of $\rm J_{21}$. We further make use of sink particles to follow the evolution for 10,000 years after reaching the maximum refinement level. No vigorous fragmentation is observed in UV illuminated halos while the accretion rate changes according to the thermal properties. Our findings show that a few 100-10, 000 solar mass protostars are formed when halos are irradiated by $\rm J_{21}=10-500$ at $\rm z>10$ and suggest a strong relation between the strength of UV flux and mass of a protostar. This mode of star formation is quite different from minihalos, as higher accretion rates of about $\rm 0.01-0.1$ M$_{\odot}$/yr are observed by the end of our simulations. The resulting massive stars are the potential cradles for the formation of intermediate mass black holes at earlier cosmic times and contribute to the formation of a global X-ray background. △ Less

Submitted 5 June, 2014; originally announced June 2014.

Comments: Submitted to APJ, comments are welcome. High resolution copy is available at http://www.astro.physik.uni-goettingen.de/~mlatif/IMBHs_apj.pdf

A UV flux constraint on the formation of direct collapse black holes

Authors: M. A. Latif, S. Bovino, C. Van Borm, T. Grassi, D. R. G. Schleicher, M. Spaans

Abstract: The ability of metal free gas to cool by molecular hydrogen in primordial halos is strongly associated with the strength of ultraviolet (UV) flux produced by the stellar populations in the first galaxies. Depending on the stellar spectrum, these UV photons can either dissociate $\rm H_{2}$ molecules directly or indirectly by photo-detachment of $\rm H^{-}$ as the latter provides the main pathway f… ▽ More The ability of metal free gas to cool by molecular hydrogen in primordial halos is strongly associated with the strength of ultraviolet (UV) flux produced by the stellar populations in the first galaxies. Depending on the stellar spectrum, these UV photons can either dissociate $\rm H_{2}$ molecules directly or indirectly by photo-detachment of $\rm H^{-}$ as the latter provides the main pathway for $\rm H_{2}$ formation in the early universe. In this study, we aim to determine the critical strength of the UV flux above which the formation of molecular hydrogen remains suppressed for a sample of five distinct halos at $z>10$ by employing a higher order chemical solver and a Jeans resolution of 32 cells. We presume that such flux is emitted by PopII stars implying atmospheric temperatures of $\rm 10^{4}$~K. We performed three-dimensional cosmological simulations and varied the strength of the UV flux below the Lyman limit in units of $\rm J_{21}$. Our findings show that the value of $\rm J_{21}^{crit}$ varies from halo to halo and is sensitive to the local thermal conditions of the gas. For the simulated halos it varies from 400-700 with the exception of one halo where $\rm J_{21}^{crit} \geq 1500$. This has important implications for the formation of direct collapse black holes and their estimated population at z > 6. It reduces the number density of direct collapse black holes by almost three orders of magnitude compared to the previous estimates. △ Less

Submitted 19 June, 2014; v1 submitted 23 April, 2014; originally announced April 2014.

Comments: 10 pages, 6 figures, matches the accepted version to ber published in MNRAS, higher resolution version is available at http://www.astro.physik.uni-goettingen.de/~mlatif/Jcrit.pdf

Dark-matter halo mergers as a fertile environment for low-mass Population III star formation

Authors: S. Bovino, M. A. Latif, T. Grassi, D. R. G. Schleicher

Abstract: While Population III stars are typically thought to be massive, pathways towards lower-mass Pop III stars may exist when the cooling of the gas is particularly enhanced. A possible route is enhanced HD cooling during the merging of dark-matter halos. The mergers can lead to a high ionization degree catalysing the formation of HD molecules and may cool the gas down to the cosmic microwave backgroun… ▽ More While Population III stars are typically thought to be massive, pathways towards lower-mass Pop III stars may exist when the cooling of the gas is particularly enhanced. A possible route is enhanced HD cooling during the merging of dark-matter halos. The mergers can lead to a high ionization degree catalysing the formation of HD molecules and may cool the gas down to the cosmic microwave background (CMB) temperature. In this paper, we investigate the merging of mini-halos with masses of a few 10$^5$ M$_\odot$ and explore the feasibility of this scenario. We have performed three-dimensional cosmological hydrodynamics calculations with the ENZO code, solving the thermal and chemical evolution of the gas by employing the astrochemistry package KROME. Our results show that the HD abundance is increased by two orders of magnitude compared to the no-merging case and the halo cools down to $\sim$60 K triggering fragmentation. Based on Jeans estimates the expected stellar masses are about 10 M$_\odot$. Our findings show that the merging scenario is a potential pathway for the formation of low-mass stars. △ Less

Submitted 8 April, 2014; v1 submitted 18 February, 2014; originally announced February 2014.

Comments: Submitted to MNRAS

KROME - a package to embed chemistry in astrophysical simulations

Authors: T. Grassi, S. Bovino, D. R. G. Schleicher, J. Prieto, D. Seifried, E. Simoncini, F. A. Gianturco

Abstract: Chemistry plays a key role in many astrophysical situations regulating the cooling and the thermal properties of the gas, which are relevant during gravitational collapse, the evolution of disks and the fragmentation process. In order to simplify the usage of chemical networks in large numerical simulations, we present the chemistry package KROME, consisting of a Python pre-processor which gener… ▽ More Chemistry plays a key role in many astrophysical situations regulating the cooling and the thermal properties of the gas, which are relevant during gravitational collapse, the evolution of disks and the fragmentation process. In order to simplify the usage of chemical networks in large numerical simulations, we present the chemistry package KROME, consisting of a Python pre-processor which generates a subroutine for the solution of chemical networks which can be embedded in any numerical code. For the solution of the rate equations, we make use of the high-order solver DLSODES, which was shown to be both accurate and efficient for sparse networks, which are typical in astrophysical applications. KROME also provides a large set of physical processes connected to chemistry, including photochemistry, cooling, heating, dust treatment, and reverse kinetics. The package presented here already contains a network for primordial chemistry, a small metal network appropriate for the modelling of low metallicities environments, a detailed network for the modelling of molecular clouds, a network for planetary atmospheres, as well as a framework for the modelling of the dust grain population. In this paper, we present an extended test suite ranging from one-zone and 1D-models to first applications including cosmological simulations with ENZO and RAMSES and 3D collapse simulations with the FLASH code. The package presented here is publicly available at http://kromepackage.org/ and https://bitbucket.org/krome/krome_stable △ Less

Submitted 15 January, 2014; v1 submitted 5 November, 2013; originally announced November 2013.

Comments: accepted for publication in MNRAS

Primordial star formation: relative impact of H2 three-body rates and initial conditions

Authors: S. Bovino, D. R. G. Schleicher, T. Grassi

Abstract: Population III stars are the first stars in the Universe to form at z=20-30 out of a pure hydrogen and helium gas in minihalos of 10^5-10^6 M$_\odot$ . Cooling and fragmentation is thus regulated via molecular hydrogen. At densities above 10^8 cm$^{-3}$, the three-body H2 formation rates are particularly important for making the gas fully molecular. These rates were considered to be uncertain by a… ▽ More Population III stars are the first stars in the Universe to form at z=20-30 out of a pure hydrogen and helium gas in minihalos of 10^5-10^6 M$_\odot$ . Cooling and fragmentation is thus regulated via molecular hydrogen. At densities above 10^8 cm$^{-3}$, the three-body H2 formation rates are particularly important for making the gas fully molecular. These rates were considered to be uncertain by at least a few orders of magnitude. We explore the impact of new accurate three-body H2 formation rates derived by Forrey (2013) for three different minihalos, and compare to the results obtained with three-body rates employed in previous studies. The calculations are performed with the cosmological hydrodynamics code ENZO (release 2.2) coupled with the chemistry package KROME (including a network for primordial chemistry), which was previously shown to be accurate in high resolution simulations. While the new rates can shift the point where the gas becomes fully molecular, leading to a different thermal evolution, there is no trivial trend in how this occurs. While one might naively expect the results to be inbetween the calculations based on Palla et al. (1983) and Abel et al. (2002), the behavior can be close to the former or the latter depending on the dark matter halo that is explored. We conclude that employing the correct three-body rates is about as equally important as the use of appropriate initial conditions, and that the resulting thermal evolution needs to be calculated for every halo individually. △ Less

Submitted 17 December, 2013; v1 submitted 29 July, 2013; originally announced July 2013.

Comments: 10 pages, 9 figures, A&A, 561, A13 (2014)

CH+ depletion by atomic hydrogen: accuracy of new rates in photo-dominated and self-shielded environments

Authors: S. Bovino, T. Grassi, M. Tacconi, F. A. Gianturco

Abstract: A detailed quantum analysis of a ionic reaction with a crucial role in the ISM is carried out to generate ab initio reactive cross sections with a quantum method. From them we obtain the corresponding CH+ depletion rates over a broad range of temperatures. The new rates are further linked to a complex chemical network that shows the evolution in time of the CH+ abundance in photodissociation regio… ▽ More A detailed quantum analysis of a ionic reaction with a crucial role in the ISM is carried out to generate ab initio reactive cross sections with a quantum method. From them we obtain the corresponding CH+ depletion rates over a broad range of temperatures. The new rates are further linked to a complex chemical network that shows the evolution in time of the CH+ abundance in photodissociation region (PDR) and molecular cloud (MC) environments. The evolutionary abundances of CH+ are given by numerical solutions of a large set of coupled, first-order kinetics equations by employing the new chemical package KROME. The differences found between all existing calculations from low-T experiments are explained via a simple numerical model that links the low-T cross section reductions to collinear approaches where nonadiabatic crossings dominate. The analysis of evolutionary abundance of CH+ reveals that the important region for the depletion reaction of this study is that above 100 K, hence showing that, at least for this reaction, the differences with the existing low-temperature experiments are of essentially no importance within the astrochemical environments. A detailed analysis of the chemical network involving CH+ also shows that a slight decrease in the initial oxygen abundance might lead to higher CH+ abundance since the main chemical carbon ion depletion channel is reduced in efficiency. This simplified observation might provide an alternative starting point to understand the problem of astrochemical models in matching the observed CH+ abundances. △ Less

Submitted 20 December, 2013; v1 submitted 20 June, 2013; originally announced June 2013.

Comments: 9 pags, 6 figs, submitted to MNRAS

Impact of an accurate modeling of primordial chemistry in high resolution studies

Authors: S. Bovino, T. Grassi, M. A. Latif, D. R. G. Schleicher

Abstract: The formation of the first stars in the Universe is regulated by a sensitive interplay of chemistry and cooling with the dynamics of a self-gravitating system. As the outcome of the collapse and the final stellar masses depend sensitively on the thermal evolution, it is necessary to accurately model the thermal evolution in high resolution simulations. As previous investigations raised doubts rega… ▽ More The formation of the first stars in the Universe is regulated by a sensitive interplay of chemistry and cooling with the dynamics of a self-gravitating system. As the outcome of the collapse and the final stellar masses depend sensitively on the thermal evolution, it is necessary to accurately model the thermal evolution in high resolution simulations. As previous investigations raised doubts regarding the convergence of the temperature at high resolution, we investigate the role of the numerical method employed to model the chemistry and the thermodynamics. Here we compare the standard implementation in the adaptive-mesh refinement code \verb|ENZO|, employing a first order backward differentiation formula (BDF), with the 5th order accurate BDF solver \verb|DLSODES|. While the standard implementation in \verb|ENZO| shows a strong dependence on the employed resolution, the results obtained with \verb|DLSODES| are considerably more robust, both with respect to the chemistry and thermodynamics, but also for dynamical quantities such as density, total energy or the accretion rate. We conclude that an accurate modeling of the chemistry and thermodynamics is central for primordial star formation. △ Less

Submitted 23 May, 2013; v1 submitted 7 May, 2013; originally announced May 2013.

Comments: accepted by MNRAS Letters

Chemical complexity in astrophysical simulations: optimization and reduction techniques

Authors: T. Grassi, S. Bovino, D. Schleicher, F. A. Gianturco

Abstract: Chemistry has a key role in the evolution of the interstellar medium (ISM), so it is highly desirable to follow its evolution in numerical simulations. However, it may easily dominate the computational cost when applied to large systems. In this paper we discuss two approaches to reduce these costs: (i) based on computational strategies, and (ii) based on the properties and on the topology of the… ▽ More Chemistry has a key role in the evolution of the interstellar medium (ISM), so it is highly desirable to follow its evolution in numerical simulations. However, it may easily dominate the computational cost when applied to large systems. In this paper we discuss two approaches to reduce these costs: (i) based on computational strategies, and (ii) based on the properties and on the topology of the chemical network. The first methods are more robust, while the second are meant to be giving important information on the structure of large, complex networks. To this aim we first discuss the numerical solvers for integrating the system of ordinary differential equations (ODE) associated with the chemical network. We then propose a buffer method that decreases the computational time spent in solving the ODE system. We further discuss a flux-based method that allows one to determine and then cut on the fly the less active reactions. In addition we also present a topological approach for selecting the most probable species that will be active during the chemical evolution, thus gaining information on the chemical network that otherwise would be difficult to retrieve. This topological technique can also be used as an a priori reduction method for any size network. We implemented these methods into a 1D Lagrangian hydrodynamical code to test their effects: both classes lead to large computational speed-ups, ranging from x2 to x5. We have also tested some hybrid approaches finding that coupling the flux method with a buffer strategy gives the best trade-off between robustness and speed-up of calculations. △ Less

Submitted 13 February, 2013; v1 submitted 21 November, 2012; originally announced November 2012.

Comments: accepted for publication in MNRAS

Carbon Cations and Silicon Atoms in the ISM: modeling their charge exchange reaction

Authors: M. Satta, T. Grassi, F. A. Gianturco

Abstract: The time-dependent rate coefficients for the charge exchange reaction C+ + Si -> C + Si+ for doublet and quartet states have been determined with ab initio quantum calculations coupled with a non-adiabatic transition model based on a simple Landau-Zener picture. This reaction plays a key role in determining the abundances of C, Si, and their ions, in the ISM since these abundances affect the fine… ▽ More The time-dependent rate coefficients for the charge exchange reaction C+ + Si -> C + Si+ for doublet and quartet states have been determined with ab initio quantum calculations coupled with a non-adiabatic transition model based on a simple Landau-Zener picture. This reaction plays a key role in determining the abundances of C, Si, and their ions, in the ISM since these abundances affect the fine structure cooling and hence the star formation rates. We also provide additional calculations to evaluate the differences between the gas evolution as obtained by using the empirical rate estimates found in the current literature and the calculations presented in this work which are based on our more realistic evaluation of such rates from ab initio transition probabilities . We shall thus show here that the new rates yield important differences for metal-rich environments where $T<10^4$ K and the UV flux is almost negligible, while becoming less important at higher T values and higher photon fluxes. △ Less

Submitted 5 November, 2012; v1 submitted 4 October, 2012; originally announced October 2012.

Comments: Accepted to be published on MNRAS

Electron attachment rates for PAH anions in the ISM and dark molecular clouds: dependence on their chemical properties

Authors: F. Carelli, T. Grassi, F. A. Gianturco

Abstract: CONTEXT: The attachment of free electrons to polycondensed aromatic ring molecules (PAHs) is studied for the variety of these molecules with different numbers of condensed rings and over a broad range of electron temperatures, using a multichannel quantum scattering approach. The calculations of the relevant cross sections are used in turn to model the corresponding attachment rates for each of th… ▽ More CONTEXT: The attachment of free electrons to polycondensed aromatic ring molecules (PAHs) is studied for the variety of these molecules with different numbers of condensed rings and over a broad range of electron temperatures, using a multichannel quantum scattering approach. The calculations of the relevant cross sections are used in turn to model the corresponding attachment rates for each of the systems under study, and these rates are parametrized as a function of temperature using a commonly employed expression for two-body processes in the interstellar medium (ISM). AIM: The scope of this work is to use first principles to establish the influence of chemical properties on the efficiency of the electron-attachment process for PAHs. METHODS: Quantum multichannel scattering methods are employed to generate the relevant cross sections, hence the attachment rates, using integral elastic cross sections computed over a broad range of relevant energies, from threshold up to 1000 K and linking the attachment to low-energy resonant collisions. RESULTS: The rates obtained for the present molecules are found to markedly vary within the test ensemble of the present work and to be lower than the earlier values used for the entire class of PAHs anions, when modelling their evolutions in ISM environments. The effects of such differences on the evolutions of chemical networks that include both PAH and PAH- species are analysed in some detail and related to previous calculations. △ Less

Submitted 4 December, 2012; v1 submitted 11 July, 2012; originally announced July 2012.

Comments: accepted to be published on A&A

Electron-attachment rates for carbon-rich molecules in protoplanetary atmospheres: the role of chemical differences

Authors: F. Carelli, T. Grassi, F. Sebastianelli, F. A. Gianturco

Abstract: The formation of anionic species in the interstellar medium from interaction of linear molecules containing carbon, nitrogen and hydrogen as atomic components (polyynes) with free electrons in the environment is modelled via a quantum treatment of the collision dynamics. The ensuing integral cross sections are employed to obtain the corresponding attachment rates over a broad range of temperatures… ▽ More The formation of anionic species in the interstellar medium from interaction of linear molecules containing carbon, nitrogen and hydrogen as atomic components (polyynes) with free electrons in the environment is modelled via a quantum treatment of the collision dynamics. The ensuing integral cross sections are employed to obtain the corresponding attachment rates over a broad range of temperatures for the electrons. The calculations unequivocally show that a parametrization form often employed for such rates yields a broad range of values that turn out to be specific for each molecular species considered, thus excluding using a unique set for the whole class of polyynes. △ Less

Submitted 1 October, 2012; v1 submitted 11 July, 2012; originally announced July 2012.

Comments: accepted to be published on MNRAS

Complexity reduction of astrochemical networks

Authors: T. Grassi, S. Bovino, F. A. Gianturco, P. Baiocchi, E. Merlin

Abstract: We present a new computational scheme aimed at reducing the complexity of the chemical networks in astrophysical models, one which is shown to markedly improve their computational efficiency. It contains a flux-reduction scheme that permits to deal with both large and small systems. This procedure is shown to yield a large speed-up of the corresponding numerical codes and provides good accord with… ▽ More We present a new computational scheme aimed at reducing the complexity of the chemical networks in astrophysical models, one which is shown to markedly improve their computational efficiency. It contains a flux-reduction scheme that permits to deal with both large and small systems. This procedure is shown to yield a large speed-up of the corresponding numerical codes and provides good accord with the full network results. We analyse and discuss two examples involving chemistry networks of the interstellar medium and show that the results from the present reduction technique reproduce very well the results from fuller calculations. △ Less

Submitted 18 June, 2012; v1 submitted 30 April, 2012; originally announced April 2012.

Comments: 9 pages, 7 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

Formation and Evolution of Early-Type Galaxies. III Star formation history as a function of mass and over-density

Authors: Emiliano Merlin, Cesare Chiosi, Lorenzo Piovan, Tommaso Grassi, Umberto Buonomo, Francesco La Barbera

Abstract: We investigate the influence of the initial proto-galaxies over-densities and masses on their evolution, to understand whether the internal properties of the proto-galactic haloes are sufficient to account for the varied properties of the galactic populations. By means of fully hydrodynamical N-body simulations performed with the code EvoL we produce twelve self-similar models of early-type galaxi… ▽ More We investigate the influence of the initial proto-galaxies over-densities and masses on their evolution, to understand whether the internal properties of the proto-galactic haloes are sufficient to account for the varied properties of the galactic populations. By means of fully hydrodynamical N-body simulations performed with the code EvoL we produce twelve self-similar models of early-type galaxies of different initial masses and over-densities, following their evolution from z \geq 20 down to z \leq 1. The simulations include radiative cooling, star formation, stellar energy feedback, a reionizing photoheating background, and chemical enrichment of the ISM. We find a strong correlation between the initial properties of the proto-haloes and their star formation histories. Massive (10^13M\odot) haloes experience a single, intense burst of star formation (with rates \geq 10^3M\odot/yr) at early epochs, consistently with observations, with a less pronounced dependence on the initial over-density; intermediate mass (10^11M\odot) haloes histories strongly depend on their initial over-density, whereas small (10^9M\odot) haloes always have fragmented histories, resulting in multiple stellar populations, due to the "galactic breathing" phenomenon. The galaxy models have morphological, structural and photometric properties comparable to real galaxies, often closely matching the observed data; even though some disagreement is still there, likely a consequence of some numerical choices. We conclude that internal properties are essentially sufficient to explain many of the observed features of early type galaxies, particularly the complicated and different star formation histories shown by haloes of very different mass. In this picture, nature seems to play the dominant role, whereas nurture has a secondary importance. △ Less

Submitted 17 August, 2012; v1 submitted 23 April, 2012; originally announced April 2012.

Comments: 27 pages, 25 figures, 5 tables

Formation and Evolution of the Dust in Galaxies. III. The Disk of the Milky Way

Authors: L. Piovan, C. Chiosi, E. Merlin, T. Grassi, R. Tantalo, U. Buonomo, L. P. Cassarà

Abstract: Models of chemical evolution of galaxies including the dust are nowadays required to decipher the high-z universe. In a series of three papers we have tackled the problem and set a modern chemical evolution model. In the first paper (Piovan et al., 2011a) we revised the condensation coefficients for the elements that typically are present in the dust. In the second paper (Piovan et al., 2011b) we… ▽ More Models of chemical evolution of galaxies including the dust are nowadays required to decipher the high-z universe. In a series of three papers we have tackled the problem and set a modern chemical evolution model. In the first paper (Piovan et al., 2011a) we revised the condensation coefficients for the elements that typically are present in the dust. In the second paper (Piovan et al., 2011b) we have implemented the dust into the Padova chemical model and tested it against the observational data for the Solar Neighbourhood. In this paper we extend it to the whole Disk of the Milky Way (MW). The Disk is used as a laboratory to analyze the spatial and temporal behaviour of (i) several dust grain families with the aid of which we can describe the ISM, (ii) the abundances in the gas, dust, and total ISM of the elements present in the dust and (iii) the depletion of the same elements. The temporal evolution of the dust and gas across the Disk is calculated under the effect of radial flows and a central Bar. The gradients of the abundances of C, N, O, Mg, Si, S, Ca and Fe in gas and dust across the Disk are derived as a function of time. The theoretical gradients nicely reproduce those derived from Cepheids, OB stars, Red Giants and HII regions. This provides the backbone for the companion processes of dust formation and evolution across the Disk. We examine in detail the contributions to dust by AGB stars, SNae and grain accretion in the ISM at different galacto-centric distances. Furthermore, we examine the variation of the ratio silicates/carbonaceous grains with time and position in the Disk. Finally, some hints about the depletion of the elements in regions of high and low SFR (inner and outer Disk) are presented. The results obtained make it possible to extend the model to other astrophysical situations or different theoretical models like the chemo-dynamical N-Body simulations. △ Less

Submitted 22 July, 2011; originally announced July 2011.

Formation and Evolution of the Dust in Galaxies. II. The Solar Neighbourhood

Authors: L. Piovan, C. Chiosi, E. Merlin, T. Grassi, R. Tantalo, U. Buonomo, L. P. Cassarà

Abstract: Over the past decade a new generation of chemical models have included the dust in the treatment of the ISM. This major accomplishment has been spurred by the growing amounts of data on the highly obscured high-z Universe and the intriguing local properties of the Solar Neighbourhood (SoNE). We present here a new model able to simulate the formation and evolution of dust in the ISM. The model foll… ▽ More Over the past decade a new generation of chemical models have included the dust in the treatment of the ISM. This major accomplishment has been spurred by the growing amounts of data on the highly obscured high-z Universe and the intriguing local properties of the Solar Neighbourhood (SoNE). We present here a new model able to simulate the formation and evolution of dust in the ISM. The model follows the evolution of 16 elemental species, with particular attention to those that are simultaneously present in form of gas and dust, e.g. C, N, O, Mg, Si, S, Ca and Fe. In this study we focus on the SoNe and the MW Disk as a whole which are considered as laboratories to test the physical ingredients governing the dust evolution. Infall of primordial gas, birth and death of stars, radial flows of matter between contiguous shells, presence of a central bar, star-dust emission by SNae and AGB stars, dust destruction and accretion are taken into account. The model reproduces the local depletion of the elements in the gas, and simultaneously satisfies other constraints obtained from the observations. The evolution of the element abundances in the gas and dust has been well reproduced for plausible choices of the parameters. The Mg/Si ratio, in particular, drives the formation of silicates. We show that for most of the evolution of the MW, the main process for dust enrichment is the accretion in the cold regions of the ISM. SNae dominate in the early phases of the evolution. We have examined the main factors controlling the temporal window in which SNae govern the dust budget both in low and high star forming environments. The role played by AGB stars is discussed. We find that IMFs with regular slope in the range of massive stars better reproduce the observed depletions. The results obtained for the SoNe lead us to safely extend the model. △ Less

Submitted 22 July, 2011; originally announced July 2011.

Comments: 37 pages, 18 figures and 5 tables

Formation and Evolution of the Dust in Galaxies. I. The Condensation Efficiencies

Authors: L. Piovan, C. Chiosi, E. Merlin, T. Grassi, R. Tantalo, U. Buonomo, L. P. Cassarà

Abstract: The growing interest in the high-z universe, where strongly obscured objects are present, has determined an effort to improve the simulations of dust formation and evolution in galaxies. Three main basic ingredients enter the problem influencing the total dust budget and the mixture of the dust grains: the types and amounts of dust injected by AGB stars and SNae and the accretion/destruction proce… ▽ More The growing interest in the high-z universe, where strongly obscured objects are present, has determined an effort to improve the simulations of dust formation and evolution in galaxies. Three main basic ingredients enter the problem influencing the total dust budget and the mixture of the dust grains: the types and amounts of dust injected by AGB stars and SNae and the accretion/destruction processes of dust in the ISM. They govern the relative abundances of the gas and dust in the ISM. In this study, we focus on star-dust and present a database of condensation efficiencies for the refractory elements C, O, Mg, Si, S, Ca and Fe in AGB stars and SNae that can be easily applied to the gaseous ejecta, in order to determine the amount and kind of refractory elements locally embedded into dust and injected into the ISM. The best theoretical recipes available nowadays in literature to estimate the amount of dust produced by SNae and AGB stars have been discussed and for SNae compared to the observations to get clues on the problem. The condensation efficiencies have been analyzed in the context of a classical chemical model of dust formation and evolution in the Solar Neighbourhood and Galactic Disk. Tables of coefficients are presented for (i) AGB stars at varying the metallicity and (ii) SNae at varying the density n_H of the ISM where the SNa explosions took place. In particular, we show how the controversial CNT approximation widely adopted to form dust in SNae, still gives good results and agrees with some clues coming from the observations. A new generation of dust formation models in SNae is however required to solve some contradictions recently emerged. A simple database of condensation efficiencies is set up to be used in chemical models including the effect of dust and meant to simulate galaxies going from primordial proto-galaxies to those currently seen in the local universe. △ Less

Submitted 22 July, 2011; originally announced July 2011.

Comments: 28 pages, 15 figures and 3+11 tables