Computational Modeling: Up-to-Date Approaches and Cutting-Edge Applications from Clusters, Nanostructures to Bulk Systems

The contributions in this special theme collection, in honor to Prof. P. Villarreal, cover a broad variety of computational methodologies and experimental techniques, containing studies on gas phase, clusters and condensed phase systems.

Trihydrogen Cation Helium Clusters: A New Potential Energy Surface

We present a new analytical potential energy surface (PES) for the interaction between the trihydrogen cation and a He atom,H 3 + - H e ${{H}_{3}^{+}-He}$ , in its electronic ground state. The proposed PES has been built as a sum of two contributions: a polarization energy term due to the electric field generated by the molecular cation at the position of the polarizable He atom, and an exchange-repulsion and dispersion interactions represented by a sum of "atom-bond" potentials between the three bonds ofH 3 + ${{H}_{3}^{+}}$ and the He atom. All parameters of this new PES have been chosen and fitted from data obtained from high-level ab-initio calculations. Using this new PES plus the Aziz-Slaman potential for the interaction between Helium atoms and assuming pair-wise interactions, we carry out classical Basin-Hopping (BH) global optimization, semiclassical BH with Zero Point Energy corrections, and quantum Diffusion Monte Carlo simulations. We have found the minimum energy configurations of small He clusters doped withH 3 + ${{H}_{3}^{+}}$ ,H 3 + H e N ${{H}_{3}^{+}{\left(He\right)}_{N}}$ , with N=1-16. The study of the energies of these clusters allows us to find a pronounced anomaly for N=12, in perfect agreement with previous experimental findings, which we relate to a greater relative stability of this aggregate.

Helium nanodroplets as an efficient tool to investigate hydrogen attachment to alkali cations

We report a novel method to reversibly attach and detach hydrogen molecules to positively charged sodium clusters formed inside a helium nanodroplet host matrix. It is based on the controlled production of multiply charged helium droplets which, after picking up sodium atoms and exposure to H2 vapor, lead to the formation of Nam+(H2)n clusters, whose population was accurately measured using a time-of-flight mass spectrometer. The mass spectra reveal particularly favorable Na+(H2)n and Na2+(H2)n clusters for specific "magic" numbers of attached hydrogen molecules. The energies and structures of these clusters have been investigated by means of quantum-mechanical calculations employing analytical interaction potentials based on ab initio electronic structure calculations. A good agreement is found between the experimental and the theoretical magic numbers.

From Molecular Aggregation to a One-Dimensional Quantum Crystal of Deuterium Inside a Carbon Nanotube of 1 nm Diameter

The quantum motion of clusters of up to four deuterium molecules under confinement in a single-wall (1 nm diameter) carbon nanotube is investigated by applying a highly accurate full quantum treatment of the most relevant nuclear degrees of freedom and an ab initio-derived potential model of the underlying dispersion-dominated intermolecular interactions. The wave functions and energies are calculated using an ad hoc-developed discrete variable representation (DVR) numerical approach in internal coordinates, with the space grid approaching a few billion grid points. We unambiguously demonstrate the formation of a solid-like pyramidal one-dimensional chain structure of molecules under the cylindrical nanotube confinement. The onset of solid-like packing is explained by analyzing the potential minima landscape. The stabilization of collective rotational motion through "rigid rotations" of four deuterium molecules provides conclusive evidence for the onset of a quantum solid-like behavior resembling that of quantum rings featuring persistent current (charged particles) or persistent flow (neutral particles).

A combined experimental and theoretical investigation of Cs+ ions solvated in HeN clusters

Solvation of Cs⁺ ions inside helium droplets has been investigated both experimentally and theoretically. On the one hand, mass spectra of doped helium clusters ionized with a crossed electron beam, He_NCs⁺, have been recorded for sizes up to N = 60. The analysis of the ratio between the observed peaks for each size N reveals evidences of the closure of the first solvation shell when 17 He atoms surround the alkali ion. On the other hand, we have obtained energies and geometrical structures of the title clusters by means of basin-hopping, diffusion Monte Carlo (DMC), and path integral Monte Carlo (PIMC) methods. The analytical He-Cs⁺ interaction potential employed in our calculations is represented by the improved Lennard-Jones expression optimized on high level ab initio energies. The weakness of the existing interaction between helium and Cs⁺ in comparison with some other alkali ions such as Li⁺ is found to play a crucial role. Our theoretical findings confirm that the first solvation layer is completed at N = 17 and both evaporation and second difference energies obtained with the PIMC calculation seem to reproduce a feature observed at N = 12 for the experimental ion abundance. The analysis of the DMC probability distributions reveals the important contribution from the icosahedral structure to the overall configuration for He₁₂Cs⁺.

Spectroscopy of a rotating hydrogen molecule in carbon nanotubes

Computing the energy levels of molecular hydrogen rotating in carbon nanotubes of increasing size.

Snowball formation for Cs+ solvation in molecular hydrogen and deuterium

Interactions of atomic cations with molecular hydrogen are of interest for a wide range of applications in hydrogen technologies. These interactions are fairly strong despite being non-covalent, hence one can ask whether hydrogen molecules would form dense, solid-like, solvation shells around the ion (snowballs) or rather a more weakly bound compound. In this work, the interactions between Cs+ and H2 are studied both experimentally and computationally. Isotopic substitution of H2 by D2 is also investigated. On the one hand, helium nanodroplets doped with cesium and hydrogen or deuterium are ionized by electron impact and the (H2/D2)nCs+ (up to n = 30) clusters formed are identified via mass spectrometry. On the other hand, a new analytical potential energy surface, based on ab initio calculations, is developed and used to study cluster energies and structures by means of classical and quantum-mechanical Monte Carlo methods. The most salient features of the measured ion abundances are remarkably mimicked by the computed evaporation energies, particularly for the clusters composed of deuterium. This result supports the reliability of the present potential energy surface and allows us to recommend its use in related systems. Clusters with either twelve H2 or D2 molecules stand out for their stability and quasi-rigid icosahedral structures. However, the first solvation shell involves thirteen or fourteen molecules for hydrogenated or deuterated clusters, respectively. This shell retains its internal structure when extra molecules are added to the second shell and is nearly solid-like, especially for the deuterated clusters. The role played by three-body induction interactions as well as the rotational degrees of freedom is analyzed and they are found to be significant (up to 15% and 18%, respectively) for the molecules belonging to the first solvation shell.

Lithium ions solvated in helium

We report on a combined experimental and theoretical study of Li+ ions solvated by up to 50 He atoms. The experiments show clear enhanced abundances associated with HenLi+ clusters where n = 2, 6, 8, and 14. We find that classical methods, e.g. basin-hopping (BH), give results that qualitatively agree with quantum mechanical methods such as path integral Monte Carlo, diffusion Monte Carlo and quantum free energy, regarding both energies and the solvation structures that are formed. The theory identifies particularly stable structures for n = 4, 6 and 8 which line up with some of the most abundant features in the experiments.

DFT Studies of Single Lithium Adsorption on Coronene

Density functional theory (DFT) calculations were performed to study the adsorption of neutral and cationic forms of single lithium (Li) on representative original and each of oxygen/sulfur-terminated coronene monolayer surfaces. First, the monolayers of coronene structures were prepared. Next, Li/Li+ adsorptions were investigated on the surfaces of the already optimised coronene models. The results indicate that the singular coronene models can be considered as appropriate surfaces for Li/Li+ adsorption, with stronger Li+ adsorption. Localisations of LI/Li+ species were carefully examined at the central carbon zone of the monolayer surface; however, only one model showed discrepancy by getting localised at the monolayer edge. Energy levels and distribution patters for the molecular orbitals indicate the effects of atomic terminations and Li/Li+ adsorptions, in which the singular and Li+-adsorbed models reveal identical results. Atomic-scale nuclear quadrupole resonance (NQR) properties were also evaluated, with the results indicating that the atomic properties can determine the major electronic properties for applying the coronene structure for specific applications.

Quantum confinement of molecular deuterium clusters in carbon nanotubes: ab initio evidence for hexagonal close packing

This study shows ab initio evidence for hexagonal close packing of D₂ molecules in carbon nanotubes, with a = 3.6 Å and .