Search

One of the most fascinating discoveries in recent years, in the cold and low pressure regions of the universe, was the detection of ArH and HeH species. The identification of such noble gas-containing molecules in space is the key to understanding noble gas chemistry. In the present work, we discuss the possibility of [ArH] existence as a potentially detectable molecule in the interstellar medium, providing new data on possible astronomical pathways and energetics of this compound. As a first step, a data-driven approach is proposed to construct a full 3D machine-learning potential energy surface (ML-PES) via the reproducing kernel Hilbert space (RKHS) method. The training and testing data sets are generated from CCSD(T)/CBS[56] computations, while a validation protocol is introduced to ensure the quality of the potential. In turn, the resulting ML-PES is employed to compute vibrational levels and molecular spectroscopic constants for the cation. In this way, the most common isotopologue in ISM, [ArH], was characterized for the first time, while simultaneously, comparisons with previously reported values available for [ArH] are discussed. Our present data could serve as a benchmark for future studies on this system, as well as on higher-order cationic Ar-hydrides of astrophysical interest.

Computational quantum chemistry has become a powerful tool with a wide range of possibilities to solve chemical-physical problems. As a result of this, the interest in the applications of computational quantum chemistry has expanded considerably, and has opened up novel research opportunities. In particular, those related to the characterization of heavy-atoms complexes, as most electronic structure calculations for such systems struggle with the problem posed by the large number of electrons present in them, and consequently, the introduction of relativistic effects. The present study performed an exhaustive assess to characterized the uncommon NgH+ (Ng = Kr, Xe, and Rn) hydride cations in order to provide accurate rovibrational data of their isotopes to assist in the laboratory characterization or even their astronomical detection. Scalar relativistic effects were included, and the ground and first electronically exited states potential curves were obtained from benchmark ab initio CCSD(T)/CBS and MRCI+Q electronic structure calculations. Next, such interaction potentials, correctly extended to long-range asymptotic regions, were employed in quantum bound state calculations and molecular spectroscopic constants were determined for the most abundance 84Kr, 132Xe, and 222Rn isotopes. Our results were discussed in comparison with available experimental and previous theoretical estimates, aiming to treat accuracy issues. The new sets provide reference data that could serve for spectroscopic characterization of such low abundance and high radioactive species.

We report new results on the translational-rotational (T-R) states of the CO2 molecule inside the sI clathrate-hydrate cages. We adopted the multiconfiguration time-dependent Hartree methodology to solve the nuclear molecular Hamiltonian, and to address issues on the T-R couplings. Motivated by experimental X-ray observations on the CO2 orientation in the D and T sI cages, we aim to evaluate the effect of the CO2 -water interaction on quantum dynamics. Thus, we first compared semiempirical and ab initio-based pair interaction model potentials against first-principles DFT-D calculations for ascertaining the importance of nonadditive many-body effects on such guest-host interactions. Our results reveal that the rotational and translational excited states quantum dynamics is remarkably different, with the pattern and density of states clearly affected by the underlying potential model. By analyzing the corresponding the probability density distributions of the calculated T-R eigenstates on both semiempirical and ab initio pair CO2 -water nanocage potentials, we have extracted information on the altered CO2 guest local structure, and we discussed it in connection with experimental data on the orientation of the CO2 molecule in the D and T sI clathrate cages available from neutron diffraction and 13C solid-state NMR studies, as well as in comparison with previous molecular dynamics simulations. Our calculations provide a very sensitive test of the potential quality by predicting the low-lying T-R states and corresponding transitions for the encapsulated CO2 molecule. As such spectroscopic observables have not been measured so far, our results could trigger further detailed experimental and theoretical investigations leading to a quantitative description of the present guest-host interactions.

The structures and energetics of Li 2+ -doped He clusters have been determined by means of evolutionary programming optimizations and classical molecular dynamics simulations. The underlying interactions in the He N Li 2+ complexes are described by sum-of-potentials ab initio-based models. The classical picture of the He atoms surrounding the cationic dimer shows a selective growth of the clusters. By analyzing spatial distribution probabilities, and single-atom evaporative energies from molecular dynamics calculations, we found pronounced drops in the computed energy for N= 2 , 4 and 6, and smaller ones for N= 10 , 13, 15, 18, 20, 22, 24, 27 and 29. The most energetic favored structure (compared to its neighbors) is found when six He atoms are attached in the Li 2+ -cation, forming He 3 -motifs at each side and leading to the formation of the first solvation shell. In turn, as extra He atoms are added, a new shell is started to form at N= 10 , with the He 3 triangles being the building blocks in all classical solid-like structural arrangements in He N Li 2+ clusters. By combining the outcome of our work on the existence of such local solvation microstructures in Li 2+ -doped He clusters together with those from previous studies on alkali ions, we may speculate that such microsolvation effects could influence the short-time solute–solvent dynamics and thus contribute to the observed slow ions’ mobility in He droplets. Graphic Abstract: [Figure not available: see fulltext.].

A water molecule confined inside the C70 fullerene was quantum-mechanically described using a computational approach within the MCTDH framework. Such procedure involves the development of a full-dimensional coupled hamiltonian, with an exact kinetic energy operator, including all rotational, translational and vibrational degrees of freedom of the endofullerene system. In turn, through an effective pairwise potential model, the ground and rotationally excited states of the encapsulated H2 O inside the C70 cage were calculated, and traced back to the isotropic case of the H2 O@C60 endofullerene in order to understand the nature and physical origin of the symmetry breaking observed experimentally in the latter system. Moreover, the computational scheme used here allows to study the quantization of the translational movement of the encapsulated water molecule inside the C70 fullerene, and to investigate the confinement effects in the vibrational energy levels of the H2 O@C70 system.

In the realm of solid water hydrostructures, helium atoms have a tendency to occupy the interstitial spaces formed within the crystal lattice of ice structures. The primary objective of this study is to examine the stability of various ice crystals when influenced by the presence of He atoms. Presenting a first attempt at a detailed computational description of the whole energy components (guest-water, water-water, guest-guest) in the complete crystal unit cells contributes to enhancing the knowledge available about these relatively unexplored helium-water systems, which could potentially benefit future experiments. For this purpose, two different ice structures were considered: the previously established He@ice II system, and the predicted (but currently nonexistent) He@ice XVII system. One of the main features of these He-filled structures is the stability conferred by the weak van der Waals dispersion forces that occur between the host lattice and the guest atoms, in addition to the hydrogen bonds established among the water molecules. Hence, it is crucial to accurately describe these interactions. Therefore, the first part of this research is devoted examining the performance and accuracy of various semi-local and non-local DFT/DFT-D functionals, in comparison with previous experimental and/or high-level computational data. Once the best-performing DFT functional has been identified, the stability of these empty and He-filled structures, including different number of He atoms within the lattices, is analysed in terms of their structural (lattice deformation), mechanical (pressure compression effects) and energetic properties (binding and saturation energies). In this manner, the potential formation of these structures under zero temperature and pressure conditions can be evaluated, while their maximum storage capacity is also determined. The obtained results reveal that, despite the weak underlying interactions, the He encapsulation has a rather notable effect o

The reactive collision between 36 Ar and the 36 ArH+ species has been investigated by means of quantum mechanical (QM), quasiclassical trajectories (QCT) and statistical quantum mechanical (SQM) approaches. Reaction probabilities, cross sections as a function of the energy and rate constants in terms of the temperature have been obtained. Cumulative distributions as a function of the collision time and the inspection of selected QCT corresponding to specific dynamical mechanisms have been analysed. Predictions by means of the SQM method are in good agreement with the QM results, thus supporting the complex-forming nature of the process.

Clathrate hydrates have the ability to encapsulate atoms and molecules within their cavities, and thus they could be potentially large storage capacity materials. The present work studies the multiple cage occupancy effects in the recently discovered He@sII crystal. On the basis of previous theoretical and experimental findings, the stability of He(1/1)@sII, He(1/4)@sII and He(2/4)@sII crystals was analysed in terms of structural, mechanical and energetic properties. For this purpose, first-principles DFT/DFT-D computations were performed by using both semi-local and non-local functionals, not only to elucidate which configuration is the most energetically favoured, but also to scrutinize the relevance of the long-range dispersion interactions in these kinds of compounds. We have encountered that dispersion interactions play a fundamental role in describing the underlying interactions, and different tendencies were observed depending on the choice of the functional. We found that PW86PBE-XDM shows the best performance, while the non-local functionals tested here were not able to correctly account for them. The present results reveal that the most stable configuration is the one presenting singly occupied D cages and tetrahedrally occupied H cages (He(1/4)@sII) in line with the experimental observation.

We performed first-principles computations to investigate the complex interplay of molecular interaction energies in determining the lattice structure and stability of CO2@sH clathrate hydrates. Density functional theory computations using periodic boundary conditions were employed to characterize energetics and the key structural properties of the sH clathrate crystal under pressure, such as equilibrium lattice volume and bulk modulus. The performance of exchange−correlation functionals together with recently developed dispersion-corrected schemes was evaluated in describing interactions in both short-range and long-range regions of the potential. Structural relaxations of the fully CO2-filled and empty sH unit cells yield crystal structure and lattice energies, while their compressibility parameters were derived by including the pressure dependencies. The present quantum chemistry computations suggest anisotropy in the compressibility of the sH clathrate hydrates, with the crystal being less compressible along the a-axis direction than along the c-axis one, in distinction from nearly isotropic sI and sII structures. The detailed results presented here give insight into the complex nature of the underlying guest−host interactions, checking earlier assumptions, providing critical tests, and improving estimates. Such entries may eventually lead to better predictions of thermodynamic properties and formation conditions, with a direct impact on emerging hydrate-based technologies.

Dissolved ions in aqueous media are ubiquitous in many physicochemical processes, with a direct impact on research fields, such as chemistry, climate, biology, and industry. Ions play a crucial role in the structure of the surrounding network of water molecules as they can either weaken or strengthen it. Gaining a thorough understanding of the underlying forces from small clusters to bulk solutions is still challenging, which motivates further investigations. Through a systematic analysis of the interaction energies obtained from high-level electronic structure methodologies, we assessed various dispersion-corrected density functional approaches, as well as ab initio-based data-driven potential models for halide ion-water clusters. We introduced an active learning scheme to automate the generation of optimally weighted datasets, required for the development of efficient bottom-up anion-water models. Using an evolutionary programming procedure, we determined optimized and reference configurations for such polarizable and first-principles-based representation of the potentials, and we analyzed their structural characteristics and energetics in comparison with estimates from DF-MP2 and DFT+D quantum chemistry computations. Moreover, we presented new benchmark datasets, considering both equilibrium and non-equilibrium configurations of higher-order species with an increasing number of water molecules up to 54 for each F, Cl, Br, and I anions, and we proposed a validation protocol to cross-check methods and approaches. In this way, we aim to improve the predictive ability of future molecular computer simulations for determining the ongoing conflicting distribution of different ions in aqueous environments, as well as the transition from nanoscale clusters to macroscopic condensed phases.

We report the results of a detailed and accurate investigation focused on structures and energetics of poly-hydrated halides employing first-principles polarizable halide-water potentials to describe the underlying forces. Following a bottom-up data-driven potential approach, we initially looked into the classical behavior of higher-order X(HO) clusters. We have located several low-lying energies, such as global and local minima, structures for each cluster, with various water molecules (up to N = 8) surrounding the halide anion (X = F, Cl, Br, I), employing an evolutionary programming method. It is found that the F-water clusters exhibit different structural configurations than the heavier halides, however independently of the halide anion, all clusters show in general a selective growth with the anion preferring to be connected to the outer shell of the water molecule arrangements. In turn, path-integral molecular dynamics simulations are performed to incorporate explicitly nuclear quantum and thermal effects in describing the nature of halide ion microsolvation in such prototypical model systems. Our data reveal that at low finite temperatures, nuclear quantum effects affect certain structural properties, such as weakening hydrogen bonding between the halide anion and water molecules, with minor distortions in the water network beyond the first hydration shell, indicating local structure rearrangements. Such structural characteristics and the promising cluster size trends observed in the single-ion solvation energies motivated us to draw connections of small size cluster data to the limits of continuum bulk values, toward the investigation of the challenging computational modeling of bulk single ion hydration.

We explore the origin of the anomalous splitting of the 1 levels reported experimentally for the HO@C endofullerene, in order to give some insight about the physical interpretations of the symmetry breaking observed. We performed fully-coupled quantum computations within the multiconfiguration time-dependent Hartree approach employing a rigorous procedure to handle such computationally challenging problems. We introduce two competing physical models, and discuss the observed unconventional quantum patterns in terms of anisotropy in the interfullerene interactions, caused by the change in the off-center position of the encapsulated water molecules inside the cage or the uniaxial C-cage distortion, arising from noncovalent bonding upon water's encapsulation, or exohedral fullerene perturbations. Our results show that both scenarios could reproduce the experimentally observed rotational degeneracy pattern, although quantitative agreement with the available experimental rotational levels splitting value has been achieved by the model that considers an uniaxial elongation of the C-cage. Such finding supports that the observed symmetry breaking could be mainly caused by the distortion of the fullerene cage. However, as nuclear quantum treatments rely on the underlying interactions, a decisive conclusion hinges on the availability of their improved description, taken into account both endofullerene and exohedral environments, from forthcoming highly demanding electronic structure many-body interaction studies.

The dynamics of the formation of a specific clathrate hydrate as well as its thermodynamic transitions depend on the interactions between the trapped molecules and the host water lattice. The molecular-level understanding of the different underlying processes benefits not only the description of the properties of the system, but also allows the development of multiple technological applications such as gas storage, gas separation, energy transport, etc. In this work we investigate the stability of periodic crystalline structures, such as He@sI and He@sII clathrate hydrates by first-principles computations. We consider such host water networks interacting with a guest He atom using selected density functional theory approaches, in order to explore the effects on the encapsulation of a light atom in the sI/sII crystals, by deriving all energy components (guest-water, water-water, guest-guest). Structural properties and energies were first computed by structural relaxations of the He-filled and empty sI/sII unit cells, yielding lattice and compressibility parameters comparable to experimental and theoretical values available for those hydrates. According to the results obtained, the He enclathration in the sI/sII unit cells is a stabilizing process, and both He@sI and He@sII clathrates, considering single cage occupancy, are predicted to be stable whatever the XDM or D4 dispersion correction applied. Our results further reveal that despite the weak underlying interactions the He encapsulation has a rather notable effect on both lattice parameters and energetics, with the He@sII being the most energetically favorable in accord with recent experimental observations.

Dissolved ions in aqueous media are ubiquitous in many physicochemical processes, with a direct impact on research fields, such as chemistry, climate, biology, and industry. Ions play a crucial role in the structure of the surrounding network of water molecules as they can either weaken or strengthen it. Gaining a thorough understanding of the underlying forces from small clusters to bulk solutions is still challenging, which motivates further investigations. Through a systematic analysis of the interaction energies obtained from high-level electronic structure methodologies, we assessed various dispersion-corrected density functional approaches, as well as ab initio-based data-driven potential models for halide ion-water clusters. We introduced an active learning scheme to automate the generation of optimally weighted datasets, required for the development of efficient bottom-up anion-water models. Using an evolutionary programming procedure, we determined optimized and reference configurations for such polarizable and first-principles-based representation of the potentials, and we analyzed their structural characteristics and energetics in comparison with estimates from DF-MP2 and DFT+D quantum chemistry computations. Moreover, we presented new benchmark datasets, considering both equilibrium and non-equilibrium configurations of higher-order species with an increasing number of water molecules up to 54 for each F, Cl, Br, and I anions, and we proposed a validation protocol to cross-check methods and approaches. In this way, we aim to improve the predictive ability of future molecular computer simulations for determining the ongoing conflicting distribution of different ions in aqueous environments, as well as the transition from nanoscale clusters to macroscopic condensed phases.