Computational development of models and tools for the kinetic study of astrochemical gas-phase reactions

This PhD thesis focuses on the application and development of computational tools and methodologies for the modeling of the kinetics of gas-phase reactions of astrophysical interest in the interstellar medium (ISM). The complexity related to the investigation of chemical reactivity in space is mostly due to the extreme physical conditions of temperature, pressure and exposure to high-energy radiation, which in turn also lead to the formation of exotic species, like radicals and ions. Nevertheless, there is still much to be understood about the formation of molecules, the major issue being the lack of sufficient laboratory (experimental and computational) studies. A more detailed and accurate study of all the chemical processes occurring in the ISM will allow us to obtain the data necessary to simulate the chemical evolution of an interstellar cloud over time using kinetic models including thousands of reactions that involve hundreds of species. The collection of the kinetic parameters required for the relevant reactions has led to the growth of different astrochemical databases, such as KIDA and UMIST. However, the data gathered in these catalogues are incomplete, and rely extensively on crude estimations and extrapolations. These rates are of paramount importance to get a better comprehension of the relative abundances of the chemical compounds extrapolated by the astronomers from the spectral data recorded through the radio telescopes and the in-orbit devices, like the satellites. Accurate state-of-the-art computational approaches play a fundamental role in analyzing feasible reaction mechanisms and in accurately predicting the associated kinetics. Such approaches usually rely on chemical intuition where a by-hand search of the most likely pathways is performed. Unfortunately, thisprocedure can lead to overlook significant mechanisms, especially when large molecular systems are investigated. Increasing the size of a molecule can also increase the number of its possible conformers which can show a different chemical reactivity with respect to the same chemical partner. This brings to get very complex chemical reaction networks in which hundreds of chemical species are involved and thousands of chemical reactions can occur.During the last decades, a lot of effort has been done to develop computational techniques able to perform extensive and thorough investigations of complex reaction mechanisms. Such approaches rely on automated computational protocols which drastically decrease the risk of making blunders during the search for significant reaction pathways.Furthermore, the accurate characterization of the potential energy surfaces (PESs) critical points, like reactants, intermediates, transition states and products involved in the reaction mechanism, is crucial in order to carry out a reliable kinetic investigation. The kinetic analysis of an erroneous potential energy surface, would lead to gross errors in the estimation of the rate constants of the chemical species involved in the reaction.In order to avoid such errors, the combination of high-level electronic structure calculations via composite scheme can be helpful to get a more precise estimation of the energy barriers involved in the reaction mechanism. It has been proven that "cheap"[1] composite schemes can achieve subchemical accuracy without any empirical parameters and with convenient computation times, making them perfect for the purpose of this thesis.In recent decades, many efforts have been made to develop theoretical and computational methodologies to perform accurate numerical simulations of the kinetics of such complex reaction mechanisms in a wide range of thermodynamic conditions that mimic extreme reaction environmentsas for combustion systems, the atmosphere and the ISM. Such methodologies are based on the ab initio-transition-state-theory-based master equation approach, which allows the determination of rate coefficients and branching ratios of chemical species involved in complex chemical reactions. This methodology allows to make accurate predictions of the relative abundances of the reaction products for complex reactions even under conditions of temperature and pressure not experimentally accessible, such as those that characterize the ISM. Based on these premises, this dissertation has been focused on the application of a computational protocol for the ab initio-based computational modeling and kinetic investigation of gas-phase reactions which can occur in the ISM.This protocol is based on the application of validated methodologies for the automated discovery of complex reaction mechanisms by means of the AutoMeKin[2] program, the accurate calculation of the energetic of the potential energy surfaces (PESs) through the junChS and junChS-F12a "cheap" composite schemes and the kinetic investigation using the StarRate computer program specifically designed to study gas-phase reactions of astrochemical interest in conjunction with the MESS program. Furthermore, this dissertation has been also focused on the development and implementation of StarRate, a computer program for the accurate calculation of kinetics through a chemical master equation approach of multi-step chemical reactions. StarRate is an object-based program written in the so-called F language. It is structured in three main modules, namely molecules, steps and reactions, which extract the properties needed to calculate the kinetics for the single-step reactions partecipating in the overall reaction. Another module, in_out, handles program’s input and output operations. The main program,starrate, controls the sequences of the calling of the procedures contained in each of the three main modules.Through these modular structure, StarRate[3] can compute canonical and microcanonical rate coefficients taking into account for the tunneling effect and the energy-dependent and time-dependent evolution of the species concentrations involved in the reaction mechanism. Such protocol has been applied to investigate the formation reaction mechanisms of some complex interstellar polyatomic molecules, named interstellar complex organic molecules (iCOMs). More specifically, the formation of prebiotic iCOMs in space has raised considerable interest in the scientific community, because they are considered as precursors of more complex biological systems involved in the origin of life in the Universe. Debate on the origins of these biomolecular building blocks has been further stimulated by the discovery of nucleobases and amino acids in meteorites and other extraterrestrial sources. However, few insights on the chemistry which brings to the formation of such compounds is known.  References: [1] Jacopo Lupi,Silvia Alessandrini,Cristina Puzzarini,and Vincenzo Barone.junchs and junchs-F12 models:Parameter-free efficient yet accurate compositeschemes for energies and structures of noncovalent complexes. Journal of Chem-ical Theory and Computation, 17(11):6974–6992, 2021. PMID: 34677974.[2] Emilio Martínez-Núñez, George L. Barnes, David R. Glowacki, Sabine Kopec,Daniel Peláez, Aurelio Rodríguez, Roberto Rodríguez-Fernández, Robin J. Shan-non, James J. P. Stewart, Pablo G. Tahoces, and Saulo A. Vazquez.Au-tomekin2021: An open-source program for automated reaction discovery. Journalof Computational Chemistry, 42(28):2036–2048, 2021.[3] Surajit Nandi, Bernardo Ballotta, Sergio Rampino, and Vincenzo Barone.Ageneral user-friendly tool for kinetic calculations of multi-step reactions withinthe virtual multifrequency spectrometer project. Applied Sciences, 10(5), 2020.