Predicting observable infrared signatures of nanosilicates in the diffuse interstellar medium

Gas-Phase Production of Hydroxylated Silicon Oxide Cluster Cations: Structure, Infrared Spectroscopy, and Astronomical Relevance

The interaction of free cationic silicon oxide clusters, Si x O y + (x = 2-5, y ≥ x), with dilute water vapor, was investigated in a flow tube reactor. Product mass distributions indicate cluster size-dependent dissociative water adsorption. To probe the structure and vibrational spectra of the resulting Si x O y H2 + (x = 2-4) clusters, we employed infrared multiple photon dissociation spectroscopy and density functional theory calculations. The planar rhombic cluster core of the disilicon oxides (x = 2) appears to be retained upon dissociative adsorption of one H2O unit, whereas a significant structural transformation of the tri- and tetra-silicon oxides (x = 3 and 4) is induced, resulting in an increased coordination of the Si atoms and more 3D cluster structures. In an astronomical context, we discuss the potential relevance of Si x O y H z + clusters as seeds for dust nucleation and catalysts for carbon-based chemistry in diffuse or translucent interstellar clouds, where all the necessary conditions for producing these species are found. In the produced clusters, the frequency of the isolated silanol Si-OH stretching vibrational mode is considerably blue-shifted compared to that in hydroxylated bulk silica and small inorganic compounds. This mode has a characteristic frequency range between 1200 cm-1 (8.3 μm) and 1090 cm-1 (9.2 μm) and is associated with the anomalously small Si-OH bond lengths in these ionised species. In infrared observations such high frequency Si-O stretching modes are usually associated with a pure bulk silica component of silicate cosmic dust. The presence of Si x O y H2 + clusters in low silica astrophysical environments could thus potentially be detected via their signature Si-O band using the James Webb space telescope.

Search for Chirality in Hydrogenated Magnesium Nanosilicates: A DFT and TD-DFT Investigation

The formation of silicate grains in the interstellar medium (ISM), especially those containing chiral surfaces such as clinopyroxenes, is poorly understood. Moreover, silicate interactions with various forms of hydrogen-proton (H+), neutral H (HI), and molecular hydrogen (H2) are of high importance as hydrogen comprises >90% of the ISM gas budget, and these species play important roles in the formation of new molecules in space. Furthermore, silicate surfaces catalyze the formation of H2 in the interstellar medium formed on dust grain surfaces by H-H association. The technical difficulty of in situ laboratory investigations of nanosilicate nucleation using astrophysically relevant environmental conditions makes computational chemistry a useful tool for studying potential nanosilicate structures. Furthermore, chiral surfaces interacting with chiral organic molecules could serve as templates that lead to the enantiomeric excess of l-amino acids and d-polyols detected in carbonaceous meteorites. However, in order for this effect to take place, an excess of one chiral form of a mineral is required to break the symmetry. This symmetry-breaking event could have been due to the asymmetric absorption of circularly polarized light by the nanosilicate as it traverses star-forming regions. We investigate this possibility using a metastable chiral form of an enstatite dimer as an input for density functional theory (DFT) and time-dependent (TD)-DFT calculations to obtain various properties and circular dichroism spectra. All in all, twenty-six magnesium nanosilicate structures were studied using varying degrees of hydrogenation: none, with HI, with H+, and with H2. The HSE06/aug-cc-pVQZ level of theory was used for the DFT calculations. TD-DFT calculations utilized the CAM-B3LYP/cc-pVQZ and ωB97X-D3/cc-pVQZ functional and basic set pairings. Results show that (1) all twenty-six structures have absorption bands that fall within the 0.6-28.3 μm range available with the newly launched James Webb Space Telescope and (2) there is a small enantioselective effect by UV-CPL on the eight chiral enstatite dimers (predicted g-values of up to 0.007). While the observed effect is small, it opens up the possibility that it is the inorganic material that becomes enantiomerically biased by UV-CPL, driving chiral enhancements in meteoric organic molecules.

Quantum Chemistry and Astrochemistry: A Match Made in the Heavens

Quantum chemistry can uniquely answer astrochemical questions that no other technique can provide. Computations can be parallelized, automated, and left to run continuously providing exceptional molecular throughput that cannot be done through experimentation. Additionally, the granularity of the individual computations that are required of potential energy surfaces, reaction mechanism pathways, or other quantum chemically derived observables produces a unique mosaic that make up the larger whole. These pieces can be dissected for their individual contributions or evaluated in an ad hoc fashion for each of their roles in generating the larger whole. No other scientific approach is capable of reporting such fine-grained insights. Quantum chemistry also works from a bottom-up approach in providing properties directly from the desired molecule instead of a top-down perspective as required of experiment where molecules have to be linked to observed phenomena. Furthermore, modern quantum chemistry is well within the range of "chemical accuracy" and is approaching "spectroscopic accuracy." As such, the seemingly difficult questions asked by astrochemistry that would not be asked initially for any other application require quantum chemical reference data. While the results of quantum chemical computations are needed to interpret astrochemical observation, modeling, or laboratory experimentation, such hard questions, regardless of the original need to answer them, produce unique solutions. While questions in astrochemistry often require novel developments in and implementations of quantum chemistry as outlined herein, the applications of these solutions will stretch beyond astrochemistry and may yet impact fields much closer to Earth.

A machine learning potential for simulating infrared spectra of nanosilicate clusters

The use of machine learning (ML) in chemical physics has enabled the construction of interatomic potentials having the accuracy of ab initio methods and a computational cost comparable to that of classical force fields. Training an ML model requires an efficient method for the generation of training data. Here, we apply an accurate and efficient protocol to collect training data for constructing a neural network-based ML interatomic potential for nanosilicate clusters. Initial training data are taken from normal modes and farthest point sampling. Later on, the set of training data is extended via an active learning strategy in which new data are identified by the disagreement between an ensemble of ML models. The whole process is further accelerated by parallel sampling over structures. We use the ML model to run molecular dynamics simulations of nanosilicate clusters with various sizes, from which infrared spectra with anharmonicity included can be extracted. Such spectroscopic data are needed for understanding the properties of silicate dust grains in the interstellar medium and in circumstellar environments.