Far-infrared (FIR) emission lines are a powerful tool to investigate the properties of the interstellar medium, especially in high-redshift galaxies, where ALMA observations have provided unprecedented information. Interpreting such data with state-of-the-art cosmological simulations post-processed with CLOUDY, has provided insights on the internal structure and gas dynamics of these systems. However, no detailed investigation of the consistency and uncertainties of this kind of analysis has been performed to date. Here, we compare different approaches to estimate FIR line emission from state-of-the-art cosmological simulations, either with CLOUDY or with on-the-fly non-equilibrium chemistry. We find that [CII]$_{158\mu}$ predictions are robust. [OI] emission lines are instead model-dependent, as these lines are strongly affected by the thermodynamic state of the gas and non-equilibrium photoionisation effects. For the same reasons, [OI] lines represent an excellent ISM diagnostic. Future observations targeting these lines will be also crucial to constrain models.

Cosmic rays are a global source of ionisation, and the ionisation fraction represents a fundamental parameter in the interstellar medium. Ions couple to magnetic fields, affect the chemistry, and the dynamics of star-forming regions as well as planetary atmospheres. However, the cosmic-ray ionisation rate represents one of the bottlenecks for astrochemical models, and its determination is one of the most puzzling problems in astrophysics. While for diffuse clouds reasonable values have been provided from H$_3^+$ observations, for dense clouds, due to the lack of rotational transitions, this is not possible, and estimates are strongly biased by the employed model. We present here an analytical expression, obtained from first principles, to estimate the cosmic-ray ionisation rate from observational quantities. The theoretical predictions are validated with high-resolution three-dimensional numerical simulations and applied to the well known core L1544; we obtained an estimate of $\zeta_2 \sim 2-3 \times 10^{-17}$ s$^{-1}$. Our results and the analytical formulae provided represent the first model-independent, robust tool to probe the cosmic-ray ionisation rate in the densest part of star-forming regions (on spatial scales of $R \leq 0.05$ pc). An error analysis is presented to give statistical relevance to our study.

Current galaxy observations suggest that a roughly linear correlation exists between the [CII] emission and the star formation rate, either as spatially-resolved or integrated quantities. Observationally, this correlation seems to be independent of metallicity, but the very large scatter does not allow to properly assess whether this is true. On the other hand, theoretical models tend to suggest a metallicity dependence of the correlation. In this study, we investigate the metallicity evolution of the correlation via a high-resolution zoom-in cosmological simulation of a dwarf galaxy employing state-of-the-art sub-grid modelling for gas cooling, star formation, and stellar feedback, and that self-consistently evolves the abundances of metal elements out of equilibrium. Our results suggest that the correlation should evolve with metallicity, in agreement with theoretical predictions, but also that this evolution can be hardly detected in observations, because of the large scatter. We also find that most of the [CII] emission is associated with neutral gas at low-intermediate densities, whereas the highest emissivity is produced by the densest regions around star-forming regions.

Disentangling the different stages of the star-formation process, in particular in the high-mass regime, is a challenge in astrophysics. Chemical clocks could help alleviating this problem, but their evolution strongly depends on many parameters, leading to degeneracy in the interpretation of the observational data. One of these uncertainties is the degree of CO depletion. We present here the first self-consistent magneto-hydrodynamic simulations of high-mass star-forming regions at different scales, fully coupled with a non-equilibrium chemical network, which includes C-N-O bearing molecules. Depletion and desorption processes are treated time-dependently. The results show that full CO-depletion (i.e. all gas-phase CO frozen-out on the surface of dust grains), can be reached very quickly, in one third or even smaller fractions of the free-fall time, whether the collapse proceeds on slow or fast timescales. This leads to a high level of deuteration in a short time both for typical tracers like N$_2$H$^+$, as well as for the main ion H$_3^+$, the latter being in general larger and more extended. N$_2$ depletion is slightly less efficient, and no direct effects on N-bearing molecules and deuterium fractionation are observed. We show that CO depletion is not the only driver of deuteration, and that there is a strong impact on $D_{frac}$ when changing the grain-size. We finally apply a two-dimensional gaussian Point Spread Function to our results to mimic observations with single-dish and interferometers. Our findings suggest that the low-values observed in high-mass star-forming clumps are in reality masking a full-depletion stage in the inner 0.1 pc region.

An estimate of the degree of CO-depletion ($f_D$) provides information on the physical conditions occurring in the innermost and densest regions of molecular clouds. A key parameter in these studies is the size of the depletion radius, i.e. the radius within which the C-bearing species, and in particular CO, are largely frozen onto dust grains. A strong depletion state (i.e. $f_D>10$, as assumed in our models) is highly favoured in the innermost regions of dark clouds, where the temperature is $<20$ K and the number density of molecular hydrogen exceeds a few $\times$10$^{4}$ cm$^{-3}$. In this work, we estimate the size of the depleted region by studying the Infrared Dark Cloud (IRDC) G351.77-0.51. Continuum observations performed with the $Herschel$ $Space$ $Observatory$ and the $LArge$ $APEX$ $BOlometer$ $CAmera$, together with APEX C$^{18}$O and C$^{17}$O J=2$\rightarrow$1 line observations, allowed us to recover the large-scale beam- and line-of-sight-averaged depletion map of the cloud. We built a simple model to investigate the depletion in the inner regions of the clumps in the filament and the filament itself. The model suggests that the depletion radius ranges from 0.02 to 0.15 pc, comparable with the typical filament width (i.e.$\sim$0.1 pc). At these radii, the number density of H$_2$ reaches values between 0.2 and 5.5$\times$10$^{5}$ cm$^{-3}$. These results provide information on the approximate spatial scales on which different chemical processes operate in high-mass star-forming regions and also suggest caution when using CO for kinematical studies in IRDCs.

Observations of high-redshift quasars provide information on the massive black holes (MBHs) powering them and the galaxies hosting them. Current observations of $z \gtrsim 6$ hosts, at sub-mm wavelengths, trace the properties of cold gas, and these are used to compare with the correlations between MBHs and galaxies characterising the $z=0$ population. The relations at $z=0$, however, rely on stellar-based tracers of the galaxy properties. We perform a very-high resolution cosmological zoom-in simulation of a $z=7$ quasar including state-of-the-art non-equilibrium chemistry, MBH formation, growth and feedback, to assess the evolution of the galaxy host and the central MBH, and compare the results with recent ALMA observations of high-redshift quasars. We measure both the stellar-based quantities used to establish the $z=0$ correlations, as well as the gas-based quantities available in $z \gtrsim 6$ observations, adopting the same assumptions and techniques used in observational studies. The high-redshift studies argued that MBHs at high redshift deviate from the local MBH-galaxy correlations. In our analysis of the single galaxy we evolve, we find that the high-redshift population sits on the same correlations as the local one, when using the same tracers used at $z=0$. When using the gas-based tracers, however, MBHs appear to be over-massive. The discrepancy between local and high-redshift MBHs seems caused by the different tracers employed, and necessary assumptions, and not by an intrinsic difference. Better calibration of the tracers, higher resolution data and availability of facilities that can probe the stellar population will be crucial to assess precisely and accurately high-redshift quasar hosts.

Young massive stars are usually found embedded in dense and massive molecular clumps and are known for being highly obscured and distant. During their formation process, deuteration is regarded as a potentially good indicator of the formation stage. Therefore, proper observations of such deuterated molecules are crucial, but still, hard to perform. In this work, we test the observability of the transition o-H$_2$D$^+(1_{10}$-$1_{11})$, using a synthetic source, to understand how the physical characteristics are reflected in observations through interferometers and single-dish telescopes. In order to perform such tests, we post-processed a magneto-hydrodynamic simulation of a collapsing magnetized core using the radiative transfer code POLARIS. Using the resulting intensity distributions as input, we performed single-dish (APEX) and interferometric (ALMA) synthetic observations at different evolutionary times, always mimicking realistic configurations. Finally, column densities were derived to compare our simulations with real observations previously performed. Our derivations for o-H$_2$D$^+$ are in agreement with values reported in the literature, in the range of 10$^{\!10-11}$cm$^{\!-2}$ and 10$^{\!12-13}$cm$^{\!-2}$ for single-dish and interferometric measurements, respectively.

Context: In cold and dense gas prior to the formation of young stellar objects, heavy molecular species (including CO) are accreted onto dust grains. Under these conditions H$_3^+$ and its deuterated isotopologues become more abundant, enhancing the deuterium fraction of molecules such as N$_2$H$^+$ that are formed via ion-neutral reactions. Because this process is extremely temperature sensitive, the abundance of these species is likely linked to the evolutionary stage of the source. Aims: We investigate how the abundances of o-H$_2$D$^+$ and N$_2$D$^+$ vary with evolution in high-mass clumps. Methods: We observed with APEX the ground-state transitions of o-H$_2$D$^+$ near 372 GHz, and N$_2$D$^+$(3-2) near 231 GHz for three massive clumps in different evolutionary stages. The sources were selected within the G351.77-0.51 complex to minimise the variation of initial chemical conditions, and to remove distance effects. We modelled their dust continuum emission to estimate their physical properties, and also modelled their spectra under the assumption of local thermodynamic equilibrium to calculate beam-averaged abundances. Results: We find an anticorrelation between the abundance of o-H$_2$D$^+$ and that of N$_2$D$^+$, with the former decreasing and the latter increasing with evolution. With the new observations we are also able to provide a qualitative upper limit to the age of the youngest clump of about 10$^5$ yr, comparable to its current free-fall time. Conclusions: We can explain the evolution of the two tracers with simple considerations on the chemical formation paths, depletion of heavy elements, and evaporation from the grains. We therefore propose that the joint observation and the relative abundance of o-H$_2$D$^+$ and N$_2$D$^+$ can act as an efficient tracer of the evolutionary stages of the star-formation process.

Black holes with masses of $\rm 10^6-10^9~M_{\odot}$ dwell in the centers of most galaxies, but their formation mechanisms are not well known. A subdominant dissipative component of dark matter with similar properties to the ordinary baryons, known as mirror dark matter, may collapse to form massive black holes during the epoch of first galaxies formation. In this study, we explore the possibility of massive black hole formation via this alternative scenario. We perform three-dimensional cosmological simulations for four distinct halos and compare their thermal, chemical and dynamical evolution in both the ordinary and the mirror sectors. We find that the collapse of halos is significantly delayed in the mirror sector due to the lack of $\rm H_2$ cooling and only halos with masses above $ \rm \geq 10^7~ M_{\odot}$ are formed. Overall, the mass inflow rates are $\rm \geq 10^{-2}~M_{\odot}/yr$ and there is less fragmentation. This suggests that the conditions for the formation of massive objects, including black holes, are more favorable in the mirror sector.

To explain the observed population of supermassive black holes at z~7, very massive seed black holes or, alternatively, super-Eddington scenarios are needed to reach final masses of the order of 10^9 solar masses. A popular explanation for massive seeds has been the direct collapse model, which predicts the formation of a single massive object due to the direct collapse of a massive gas cloud. Simulations over the last years have however shown that such a scenario is very difficult to achieve. A realistic model of black hole formation should therefore take fragmentation into account, and consider the interaction between stellar-dynamical and gas-dynamical processes. We present here numerical simulations pursued with the AMUSE code, employing an approximate treatment of the gas. Based on these simulations, we show that very massive black holes of 10^4-10^5 solar masses may form depending on the gas supply and the accretion onto the protostars.

The presence of supermassive black holes (SMBHs) of a few billion solar masses at very high redshift has motivated us to study how these massive objects formed during the first billion years after the Big Bang. The most promising model that has been proposed to explain this is the direct collapse of protogalactic gas clouds. In this scenario, very high accretion rates are needed to form massive objects early on and the suppression of $\mathrm H_2$ cooling is important in regulating the fragmentation. Recent studies have shown that if we use a strong radiation background, the hydrogen molecules are destroyed, favoring the high accretion rates and therefore producing objects of very high mass. In this work we study the impact of UV radiation fields in a primordial gas cloud using the recently coupled code GRADSPH-KROME for the modeling of gravitational collapse including primordial chemistry to explore the fragmentation in AU scales and hence the formation of first SMBHs. We found that to suppress the formation of $\mathrm H_2$ a very high value of $J_{21}$ is required, because of that we conclude that the direct collapse black holes (DCBHs) are very unlikely to be an explanation for the formation of the first SMBHs.

We present here a three-dimesional hydrodynamical simulation for star formation. Our aim is to explore the effect of the metal-line cooling on the thermodynamics of the star-formation process. We explore the effect of changing the metallicty of the gas from $Z/Z_{\odot}=10^{-4}$ to $Z/Z_{\odot}=10^{-2}$. Furthermore, we explore the implications of using the observational abundance pattern of a CEMP-no star, which have been considered to be the missing second-generation stars, the so-called Pop. III.2 stars. In order to pursue our aim, we modelled the microphysics by employing the public astrochemistry package KROME, using a chemical network which includes sixteen chemical species (H, H$^{+}$, H$^{-}$, He, He$^{+}$, He$^{++}$, e$^{-}$, H$_{2}$, H$_{2}^{+}$, C, C$^{+}$, O, O$^{+}$, Si, Si$^{+}$, and Si$^{++}$). We couple KROME with the fully three-dimensional Smoothed-particle hydrodynamics (SPH) code GRADSPH. With this framework we investigate the collapse of a metal-enhanced cloud, exploring the fragmentation process and the formation of stars. We found that the metallicity has a clear impact on the thermodynamics of the collapse, allowing the cloud to reach the CMB temperature floor for a metallicity $Z/Z_{\odot}=10^{-2}$, which is in agreement with previous work. Moreover, we found that adopting the abundance pattern given by the star SMSS J031300.36-670839.3 the thermodynamics behavior is very similar to simulations with a metallicity of $Z/Z_{\odot}=10^{-2}$, due to the high carbon abundance. As long as only metal line cooling is considered, our results support the metallicity threshold proposed by previous works, which will very likely regulate the first episode of fragmentation and potentially determine the masses of the resulting star clusters.

The interplay between chemistry and thermodynamics determines the final outcome of the process of gravitational collapse and sets the conditions for the formation of the first cosmological objects, including primordial supermassive black holes. In this chapter, we will review the main chemical reactions and the most important heating/cooling processes taking place in a gas of primordial composition, including the effects of local and cosmological radiation backgrounds.

The first stars are known to form in primordial gas, either in minihalos with about $10^6$~M$_\odot$ or so-called atomic cooling halos of about $10^8$~M$_\odot$. Simulations have shown that gravitational collapse and disk formation in primordial gas yield dense stellar clusters. In this paper, we focus particularly on the formation of protostellar binary systems, and aim to quantify their properties during the early stage of their evolution. For this purpose, we combine the smoothed particle hydrodynamics code GRADSPH with the astrochemistry package KROME. The GRADSPH-KROME framework is employed to investigate the collapse of primordial clouds in the high-density regime, exploring the fragmentation process and the formation of binary systems. We observe a strong dependence of fragmentation on the strength of the turbulent Mach number $\mathcal{M}$ and the rotational support parameter $\beta{}$. Rotating clouds show significant fragmentation, and have produced several Pop.~III proto-binary systems. We report maximum and minimum mass accretion rates of $2.31 \times 10^{-1}$~M$_{\odot}$ yr$^{-1}$ and $2.18\times 10^{-4}$~M$_{\odot}$ yr$^{-1}$. The mass spectrum of the individual Pop III proto-binary components ranges from $0.88$~M$_{\odot}$ to $31.96$~M$_{\odot}$ and has a sensitive dependence on the Mach number $\mathcal{M}$ as well as on the rotational parameter $\beta{}$. We also report a range from $\sim0.01$ to $\sim1$ for the mass ratio of our proto-binary systems.

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.

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.

In this study, we present a suite of high-resolution numerical simulations of an isolated galaxy to test a sub-grid framework to consistently follow the formation and dissociation of H$_2$ with non-equilibrium chemistry. The latter is solved via the package KROME, coupled to the mesh-less hydrodynamic code GIZMO. We include the effect of star formation (SF), modelled with a physically motivated prescription independent of H$_2$, supernova feedback and mass losses from low-mass stars, extragalactic and local stellar radiation, and dust and H$_2$ shielding, to investigate the emergence of the observed correlation between H$_2$ and SF rate surface densities. We present two different sub-grid models and compare them with on-the-fly radiative transfer (RT) calculations, to assess the main differences and limits of the different approaches. We also discuss a sub-grid clumping factor model to enhance the H$_2$ formation, consistent with our SF prescription, which is crucial, at the achieved resolution, to reproduce the correlation with H$_2$. We find that both sub-grid models perform very well relative to the RT simulation, giving comparable results, with moderate differences, but at much lower computational cost. We also find that, while the Kennicutt-Schmidt relation for the total gas is not strongly affected by the different ingredients included in the simulations, the H$_2$-based counterpart is much more sensitive, because of the crucial role played by the dissociating radiative flux and the gas shielding.

We show that a subdominant component of dissipative dark matter resembling the Standard Model can form many intermediate-mass black hole seeds during the first structure formation epoch. We also observe that, in the presence of this matter sector, the black holes will grow at a much faster rate with respect to the ordinary case. These facts can explain the observed abundance of supermassive black holes feeding high-redshift quasars. The scenario will have interesting observational consequences for dark substructures and gravitational wave production.

Deuterium fractionation is considered as an important process to infer the chemical ages of prestellar cores in filaments. We present here the first magneto-hydrodynamical simulations including a chemical network to study deuterium fractionation in magnetized and turbulent filaments and their substructures. The filaments typically show widespread deuterium fractionation with average values $\gtrsim0.01$. For individual cores of similar age, we observe the deuteration fraction to increase with time, but also to be independent of their average properties such as density, virial or mass-to-magnetic flux ratio. We further find a correlation of the deuteration fraction with core mass, average H$_2$ density and virial parameter only at late evolutionary stages of the filament and attribute this to the lifetime of the individual cores. Specifically, chemically old cores reveal higher deuteration fractions. Within the radial profiles of selected cores, we notice differences in the structure of the deuteration fraction or surface density, which we can attribute to their different turbulent properties. High deuteration fractions of the order $0.01-0.1$ may be reached within approximately $200$~kyrs, corresponding to two free-fall times, as defined for cylindrical systems, of the filaments

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.

To improve our understanding of high-z galaxies we study the impact of H$_{2}$ chemistry on their evolution, morphology and observed properties. We compare two zoom-in high-resolution (30 pc) simulations of prototypical $M_{\star}\sim 10^{10} {\rm M}_{\odot}$ galaxies at $z=6$. The first, "Dahlia", adopts an equilibrium model for H$_{2}$ formation, while the second, "Alth{\ae}a", features an improved non-equilibrium chemistry network. The star formation rate (SFR) of the two galaxies is similar (within 50\%), and increases with time reaching values close to 100 ${\rm M}_{\odot}/\rm yr$ at $z=6$. They both have SFR-stellar mass relation consistent with observations, and a specific SFR of $\simeq 5\, {\rm Gyr}^{-1}$. The main differences arise in the gas properties. The non-equilibrium chemistry determines the H$\rightarrow$ H$_{2}$~transition to occur at densities $> 300\,{cm}^{-3}$, i.e. about 10 times larger than predicted by the equilibrium model used for Dahlia. As a result, Alth{\ae}a features a more clumpy and fragmented morphology, in turn making SN feedback more effective. Also, because of the lower density and weaker feedback, Dahlia sits $3\sigma$ away from the Schmidt-Kennicutt relation; Alth{\ae}a, instead nicely agrees with observations. The different gas properties result in widely different observables. Alth{\ae}a outshines Dahlia by a factor of 7 (15) in [CII]~$157.74\,\mu{\rm m}$ (H$_{2}$~$17.03\,\mu{\rm m}$) line emission. Yet, Alth{\ae}a is under-luminous with respect to the locally observed [CII]-SFR relation. Whether this relation does not apply at high-z or the line luminosity is reduced by CMB and metallicity effects remains as an open question.

High-mass stars are expected to form from dense prestellar cores. Their precise formation conditions are widely discussed, including their virial condition, which results in slow collapse for super-virial cores with strong support by turbulence or magnetic fields, or fast collapse for sub-virial sources. To disentangle their formation processes, measurements of the deuterium fractions are frequently employed to approximately estimate the ages of these cores and to obtain constraints on their dynamical evolution. We here present 3D magneto-hydrodynamical simulations including for the first time an accurate non-equilibrium chemical network with 21 gas-phase species plus dust grains and 213 reactions. With this network we model the deuteration process in fully depleted prestellar cores in great detail and determine its response to variations in the initial conditions. We explore the dependence on the initial gas column density, the turbulent Mach number, the mass-to-magnetic flux ratio and the distribution of the magnetic field, as well as the initial ortho-to-para ratio of H2. We find excellent agreement with recent observations of deuterium fractions in quiescent sources. Our results show that deuteration is rather efficient, even when assuming a conservative ortho-to-para ratio of 3 and highly sub-virial initial conditions, leading to large deuterium fractions already within roughly a free-fall time. We discuss the implications of our results and give an outlook to relevant future investigations.

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.

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.

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 $\mu$m emission both for high-redshift as well as for local galaxies.

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.

The origin of strong magnetic fields in the Universe can be explained by amplifying weak seed fields via turbulent motions on small spatial scales and subsequently transporting the magnetic energy to larger scales. This process is known as the turbulent dynamo and depends on the properties of turbulence, i.e. on the hydrodynamical Reynolds number and the compressibility of the gas, and on the magnetic diffusivity. While we know the growth rate the magnetic energy in the linear regime, the saturation level, i.e. the ratio of magnetic energy to turbulent kinetic energy that can be reached, is not known from analytical calculations. In this paper we present the first scale-dependent saturation model based on an effective turbulent resistivity which is determined by the turnover timescale of turbulent eddies and the magnetic energy density. The magnetic resistivity increases compared to the Spitzer value and the effective scale on which the magnetic energy spectrum is at its maximum moves to larger spatial scales. This process ends when the peak reaches a characteristic wavenumber k* which is determined by the critical magnetic Reynolds number. The saturation level of the dynamo also depends on the type of turbulence and differs for the limits of large and small magnetic Prandtl numbers Pm. With our model we find saturation levels between 43.8% and 1.3% for Pm>>1 and between 2.43% and 0.135% for Pm<<1, where the higher values refer to incompressible turbulence and the lower ones to highly compressible turbulence.

The turbulent dynamo may explain the origin of cosmic magnetism. While the exponential amplification of magnetic fields has been studied for incompressible gases, little is known about dynamo action in highly-compressible, supersonic plasmas, such as the interstellar medium of galaxies and the early Universe. Here we perform the first quantitative comparison of theoretical models of the dynamo growth rate and saturation level with three-dimensional magnetohydrodynamical simulations of supersonic turbulence with grid resolutions of up to 1024^3 cells. We obtain numerical convergence and find that dynamo action occurs for both low and high magnetic Prandtl numbers Pm = nu/eta = 0.1-10 (the ratio of viscous to magnetic dissipation), which had so far only been seen for Pm >= 1 in supersonic turbulence. We measure the critical magnetic Reynolds number, Rm_crit = 129 (+43, -31), showing that the compressible dynamo is almost as efficient as in incompressible gas. Considering the physical conditions of the present and early Universe, we conclude that magnetic fields need to be taken into account during structure formation from the early to the present cosmic ages, because they suppress gas fragmentation and drive powerful jets and outflows, both greatly affecting the initial mass function of stars.

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.

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.

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.

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.

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.

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.

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

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.

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.

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.

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.

(Abridged) The small-scale dynamo may play a substantial role in magnetizing the Universe under a large range of conditions, including subsonic turbulence at low Mach numbers, highly supersonic turbulence at high Mach numbers and a large range of magnetic Prandtl numbers Pm, i.e. the ratio of kinetic viscosity to magnetic resistivity. Low Mach numbers may in particular lead to the well-known, incompressible Kolmogorov turbulence, while for high Mach numbers, we are in the highly compressible regime, thus close to Burgers turbulence. In this study, we explore whether in this large range of conditions, a universal behavior can be expected. Our starting point are previous investigations in the kinematic regime. Here, analytic studies based on the Kazantsev model have shown that the behavior of the dynamo depends significantly on Pm and the type of turbulence, and numerical simulations indicate a strong dependence of the growth rate on the Mach number of the flow. Once the magnetic field saturates on the current amplification scale, backreactions occur and the growth is shifted to the next-larger scale. We employ a Fokker-Planck model to calculate the magnetic field amplification during the non-linear regime, and find a resulting power-law growth that depends on the type of turbulence invoked. For Kolmogorov turbulence, we confirm previous results suggesting a linear growth of magnetic energy. For more general turbulent spectra, where the turbulent velocity v_t scales with the characteristic length scale as u_\ell\propto \ell^{\vartheta}, we find that the magnetic energy grows as (t/T_{ed})^{2\vartheta/(1-\vartheta)}, with t the time-coordinate and T_{ed} the eddy-turnover time on the forcing scale of turbulence. For Burgers turbulence, \vartheta=1/2, a quadratic rather than linear growth may thus be expected, and a larger timescale until saturation is reached.

The present-day Universe is highly magnetized, even though the first magnetic seed fields were most probably extremely weak. To explain the growth of the magnetic field strength over many orders of magnitude fast amplification processes need to operate. The most efficient mechanism known today is the small-scale dynamo, which converts turbulent kinetic energy into magnetic energy leading to an exponential growth of the magnetic field. The efficiency of the dynamo depends on the type of turbulence indicated by the slope of the turbulence spectrum v(l) \propto l^{theta}, where v(l) is the eddy velocity at a scale l. We explore turbulent spectra ranging from incompressible Kolmogorov turbulence with theta = 1/3 to highly compressible Burgers turbulence with theta = 1/2. In this work we analyze the properties of the small-scale dynamo for low magnetic Prandtl numbers Pm, which denotes the ratio of the magnetic Reynolds number, Rm, to the hydrodynamical one, Re. We solve the Kazantsev equation, which describes the evolution of the small-scale magnetic field, using the WKB approximation. In the limit of low magnetic Prandtl numbers the growth rate is proportional to Rm^{(1-theta)/(1+theta)}. We furthermore discuss the critical magnetic Reynolds number Rm_crit, which is required for small-scale dynamo action. The value of Rm_crit is roughly 100 for Kolmogorov turbulence and 2700 for Burgers. Furthermore, we discuss that Rm_crit provides a stronger constraint in the limit of low Pm than it does for large Pm. We conclude that the small-scale dynamo can operate in the regime of low magnetic Prandtl numbers, if the magnetic Reynolds number is large enough. Thus, the magnetic field amplification on small scales can take place in a broad range of physical environments and amplify week magnetic seed fields on short timescales.

The small-scale dynamo provides a highly efficient mechanism for the conversion of turbulent into magnetic energy. In astrophysical environments, such turbulence often occurs at high Mach numbers, implying steep slopes in the turbulent spectra. It is thus a central question whether the small-scale dynamo can amplify magnetic fields in the interstellar or intergalactic media, where such Mach numbers occur. To address this long-standing issue, we employ the Kazantsev model for turbulent magnetic field amplification, systematically exploring the effect of different turbulent slopes, as expected for Kolmogorov, Burgers, the Larson laws and results derived from numerical simulations. With the framework employed here, we give the first solution encompassing the complete range of magnetic Prandtl numbers, including Pm << 1, Pm ~ 1 and Pm >> 1. We derive scaling laws of the growth rate as a function of hydrodynamic and magnetic Reynolds number for Pm << 1 and Pm >> 1 for all types of turbulence. A central result concerns the regime of Pm ~ 1, where the magnetic field amplification rate increases rapidly as a function of Pm. This phenomenon occurs for all types of turbulence we explored. We further find that the dynamo growth rate can be decreased by a few orders of magnitude for turbulence spectra steeper than Kolmogorov. We calculate the critical magnetic Reynolds number Rm_c for magnetic field amplification, which is highest for the Burgers case. As expected, our calculation shows a linear behaviour of the amplification rate close to the threshold proportional to Rm-Rm_c. Based on the Kazantsev model, we therefore expect the existence of the small-scale dynamo for any given value of Pm, as long as the magnetic Reynolds number is above the critical threshold.

We explore the amplification of magnetic fields in the high-redshift Universe. For this purpose, we perform high-resolution cosmological simulations following the formation of primordial halos with \sim10^7 M_solar, revealing the presence of turbulent structures and complex morphologies at resolutions of at least 32 cells per Jeans length. Employing a turbulence subgrid-scale model, we quantify the amount of unresolved turbulence and show that the resulting turbulent viscosity has a significant impact on the gas morphology, suppressing the formation of low-mass clumps. We further demonstrate that such turbulence implies the efficient amplification of magnetic fields via the small-scale dynamo. We discuss the properties of the dynamo in the kinematic and non-linear regime, and explore the resulting magnetic field amplification during primordial star formation. We show that field strengths of \sim10^{-5} G can be expected at number densities of \sim5 cm^{-3}.

The first galaxies form due to gravitational collapse of primordial halos. During this collapse, weak magnetic seed fields get amplified exponentially by the small-scale dynamo - a process converting kinetic energy from turbulence into magnetic energy. We use the Kazantsev theory, which describes the small-scale dynamo analytically, to study magnetic field amplification for different turbulent velocity correlation functions. For incompressible turbulence (Kolmogorov turbulence), we find that the growth rate is proportional to the square root of the hydrodynamic Reynolds number, Re^(1/2). In the case of highly compressible turbulence (Burgers turbulence) the growth rate increases proportional to Re^(1/3). With a detailed chemical network we are able to follow the chemical evolution and determine the kinetic and magnetic viscosities (due to Ohmic and ambipolar diffusion) during the collapse of the halo. This way, we can calculate the growth rate of the small-scale dynamo quantitatively and predict the evolution of the small-scale magnetic field. As the magnetic energy is transported to larger scales on the local eddy-timescale, we obtain an estimate for the magnetic field on the Jeans scale. Even there, we find that equipartition with the kinetic energy is reached on small timescales. Dynamically relevant field structures can thus be expected already during the formation of the first objects in the Universe.

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.

We present the results of quantum calculations based on entirely ab initio methods for a variety of molecular processes and chemical reactions involving the LiHe$^+$ ionic polar molecule. With the aid of these calculations we derive accurate reaction rates and fitting expressions valid over a range of gas temperatures representative of the typical conditions of the pregalactic gas. With the help of a full chemical network, we then compute the evolution of the abundance of LiHe$^+$ as function of redshift in the early Universe. Finally, we compare the relative abundance of LiHe$^+$ with that of other polar cations formed in the same redshift interval.

The role of HeH+ has been newly assessed with the aid of newly calculated rates which use entirely ab initio methods, thereby allowing us to compute more accurately the relevant abundances within the global chemical network of the early universe. A comparison with the similar role of the ionic molecule LiH+ is also presented. Quantum calculations have been carried out for the gas-phase reaction of HeH+ with H atoms with our new in-house code, based on the negative imaginary potential method. Integral cross sections and reactive rate coefficients obtained under the general conditions of early universe chemistry are presented and discussed. With the new reaction rate, the abundance of HeH+ in the early universe is more than one order of magnitude larger than in previous studies. Our more accurate findings further buttress the possibility to detect cosmological signatures of HeH+.

The relative efficiencies of the chemical pathways that can lead to the destruction of LiH and LiH+ molecules, conjectured to be present in the primordial gas and to control molecular cooling processes in the gravitational collapse of the post-recombination era, are revisited by using accurate quantum calculations for the several reactions involved. The new rates are employed to survey the behavior of the relative abundance of these molecules at redshifts of interest for early universe conditions. We find significant differences with respect to previous calculations, the present ones yielding LIH abundances higher than LiH+ at all redshifts.

We present a quantum-mechanical study of the exothermic 7LiH reaction with H. Accurate reactive probabilities and rate coefficients are obtained by solving the Schrodinger equation for the motion of the three nuclei on a single Born-Oppenheimer potential energy surface (PES) and using a coupled-channel hyperspherical coordinate method. Our new rates indeed confirm earlier, qualitative predictions and some previous theoretical calculations, as discussed in the main text. In the astrophysical domain we find that the depletion process largely dominates for redshift (z) between 400 and 100, a range significant for early Universe models. This new result from first-principle calculations leads us to definitively surmise that LiH should be already destroyed when the survival processes become important. Because of this very rapid depletion reaction, the fractional abundance of LiH is found to be drastically reduced, so that it should be very difficult to manage to observe it as an imprinted species in the cosmic background radiation (CBR). The present findings appear to settle the question of LiH observability in the early Universe. We further report several state-to-state computed reaction rates in the same range of temperatures of interest for the present problem.