This paper deals with the possible formation of two molecular anions often considered likely components in the physical environments of the Interstellar Medium ( ISM) : $\mathrm{C_2H^-}$ and $\mathrm{NCO^-}$. They are both discussed here by computationally following the radiative association (RA) mechanism starting from $\mathrm{C_2^-}$, $\mathrm{H}$, $\mathrm{NC^-}$ and $\mathrm{O}$ as partners. The corresponding RA total cross sections produced by the calculations are in turn employed to generate the overall association rates over the relevant range of temperatures. The latter are found to be in line with other molecular ions formed by RA but not large enough to uniquivocally suggest this path as the main route to the anions formation in the ISM. Other possible paths of formation are also analysed and discussed. The presence of resonant structures during the association dynamics for both systems is found by the calculations and their consequences are discussed in some detail in the present study.

We investigate the relative efficiencies of low-temperature chemical reactions in the Interstellar medium (ISM) with H$^-$ anion reacting in the gas phase with cyanopolyyne neutral molecules, leading to the formation of anionic C$_x$N$^-$ linear chains of different length and of H$_2$. All the reactions turn out to be without barriers, highly exothermic reactions which provide a chemical route to the formation of anionic chains of the same length . Some of the anions have been observed in the dark molecular clouds and in the diffuse interstellar envelopes.Quantum calculations are carried for the corresponding reactive potential energy surfaces (RPESs) for all the odd-numbered members of the series (x=1, 3, 5, 7). We employ the Minimum Energy paths (MEPs) to obtain the relevant Transition State (TS) configurations and use the latter within the Variational Transition State ( VTS) model to obtain the chemical rates. The present results indicate that, at typical temperatures around 100 K, a set of significantly larger rate values exists for x=3 and x=5, while are smaller for CN$^-$ and C$_7$N$^-$. At those temperatures, however, all the rates turn out to be larger than the estimates in the current literature for the Radiative Electron Attachment (REA) rates, thus indicating the greater importance of the present chemical path with respect to REA processes at those temperatures. The physical reasons for our findings are discussed in detail and linked with the existing observational findings.

The direct dynamical paths leading to Associative Detachment (AD) in the gas-phase, and specifically in the low-temperature regions of the Dark Molecular Clouds (DMC) in the ISM, or in cold trap laboratory experiments, are investigated with quantum chemical methods by using a high-level multi-reference Configuration Interaction (CI) approach that employs single and double excitations plus Davidson perturbative correction [MRSDCI(Q)] and the d-aug-cc-pV5Z basis set. The potential energy curves for H + CN- are constructed for different directions of the H partner approaching the CN- anion within the framework of the Born-Oppenheimer approximation. The present calculations found that the AD energetics at low temperature becomes favorable only along a selected range of approaching directions, thus showing that there is a preferred path of forming HCN at low temperatures, while that of forming its HNC isomer is found to be energetically forbidden. Given the existence in the ISM of different HCN/HNC ratios in different environments, we discuss the implications of our findings for selective formation of either isomer in the low-temperature conditions of the Molecular Cloud Cores.

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

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.

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.

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.

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.

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.

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.