Authors: Garcia-Arroyo, Esther; Campos-Martinez, Jose; Bartolomei, Massimiliano; Hernandez, Marta I.; Pirani, Fernando; Halberstadt, Nadine


Publication date: 2023/10/29

DOI: 10.1002/cphc.202300424

Abstract: Interactions between molecular hydrogen and ions are of interest in cluster science, astrochemistry and hydrogen storage. In dynamical simulations, H-2 molecules are usually modelled as point particles, an approximation that can fail for anisotropic interactions. Here, we apply an adiabatic separation of the H-2 rotational motion to build effective pseudoatom-ion potentials and in turn study the properties of (H-2)(n)Na+/Cl- clusters. These interaction potentials are based on high-level ab initio calcu- lations and Improved Lennard-Jones parametrizations, while the subsequent dynamics has been performed by quantum Monte Carlo calculations. By comparisons with simulations explicitly describing the molecular rotations, it is concluded that the present adiabatic model is very adequate. Interestingly, we find differences in the cluster stabilities and coordination shells depending on the spin isomer considered (para- or ortho-H-2), especially for the anionic clusters.