Source author record

T. Grassi

T. Grassi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2026arXiv

New Estimate for the Cosmic Ray-Induced $\rm H_2$ Photodissociation Rate in the Interstellar Medium

In the interstellar medium, cosmic rays (CRs) generate a field of ultraviolet (UV) photons via the excitation and subsequent radiative decay of $\rm H_2$ molecules. This UV field is a major agent of ionization and dissociation in the inner regions of molecular clouds that are shielded from the effects of the interstellar radiation field. In particular, the dissociation of $\rm H_2$, by far the most abundant molecule in interstellar clouds, leads to the production of atomic hydrogen which then takes part in the production of a multitude of molecules, in particular complex organics on the surfaces of interstellar dust grains. Precise knowledge of the rates of CR-induced dissociation processes is thus crucial for constructing reliable chemical models. For the present paper, we have derived a new value of $k_{\rm diss, CR}(\mbox{$\rm H_2$})=0.831ζ$ for the rate of $\rm H_2$ dissociation, where $ζ$ is the CR ionization rate of $\rm H_2$. This prediction contrasts a previous value from the Leiden database which overestimated the rate due to an inconsistent treatment of the $\rm H_2$ abundances and photodissociation cross sections. By running a series of chemical models, we show that the overestimated dissociation rate has a large effect on the results of chemical simulations, with the abundance of methanol being overestimated by over one order of magnitude. Hence, we strongly recommend the adoption of our new estimate $k_{\rm diss, CR}(\mbox{$\rm H_2$})=0.831ζ$ in all chemical models that include this process. Our newly derived value corresponds to $\rm H_2$ being purely in the para form ($J^{\prime\prime} = 0$). However, in the interiors of molecular clouds the $\rm H_2$ ortho-to-para ratio is low and using the rate for para-$\rm H_2$ is an adequate approximation.

preprint2020arXiv

Modelling thermochemical processes in protoplanetary disks I: numerical methods

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.

preprint2016arXiv

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

(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.

preprint2016arXiv

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

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.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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.

preprint2014arXiv

KROME - a package to embed chemistry in astrophysical simulations

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

preprint2014arXiv

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

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.

preprint2013arXiv

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

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.

preprint2013arXiv

Chemical complexity in astrophysical simulations: optimization and reduction techniques

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.

preprint2013arXiv

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

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.

preprint2013arXiv

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

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.

preprint2012arXiv

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

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.

preprint2012arXiv

Complexity reduction of astrochemical networks

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.

preprint2012arXiv

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

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.

preprint2012arXiv

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

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.

preprint2011arXiv

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

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.

preprint2011arXiv

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

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.

preprint2011arXiv

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

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.

preprint2011arXiv

MaNN: Multiple Artificial Neural Networks for modelling the Interstellar Medium

Modelling the complex physics of the Interstellar Medium (ISM) in the context of large-scale numerical simulations is a challenging task. A number of methods have been proposed to embed a description of the ISM into different codes. We propose a new way to achieve this task: Artificial Neural Networks (ANNs). The ANN has been trained on a pre-compiled model database, and its predictions have been compared to the expected theoretical ones, finding good agreement both in static and in dynamical tests run using the Padova Tree-SPH code \textsc{EvoL}. A neural network can reproduce the details of the interstellar gas evolution, requiring limited computational resources. We suggest that such an algorithm can replace a real-time calculation of mass elements chemical evolution in hydrodynamical codes.

preprint2011arXiv

ROBO: a Model and a Code for the Study of the Interstellar Medium

We present ROBO, a model and its companion code for the study of the interstellar medium (ISM). The aim is to provide an accurate description of the physical evolution of the ISM and to set the ground for an ancillary tool to be inserted in NBody-Tree-SPH (NB-TSPH) simulations of large scale structures in cosmological context or of the formation and evolution of individual galaxies. The ISM model consists of gas and dust. The gas chemical composition is regulated by a network of reactions that includes a large number of species (hydrogen and deuterium based molecules, helium, and metals). New reaction rates for the charge transfer in $\mathrm H^+$ and $\mathrm H_2$ collisions are presented. The dust contains the standard mixture of carbonaceous grains (graphite grains and PAHs) and silicates of which the model follows the formation and destruction by several processes. The model takes into account an accurate treatment of the cooling process, based on several physical mechanisms, and cooling functions recently reported in the literature. The model is applied to a wide range of the input parameters and the results for important quantities describing the physical state of the gas and dust are presented. The results are organized in a database suited to the artificial neural networks (ANNs). Once trained, the ANNs yield the same results obtained by ROBO, with great accuracy. We plan to develop ANNs suitably tailored for applications to NB-TSPH simulations of cosmological structures and/or galaxies.