Structural and quantum effects from anionic centers in rare gas clusters: The (Ne)nH− and (Ne)n+1 systems

Solvation of Isoelectronic Halide and Alkali Metal Ions by Argon Atoms

This work systematically examines the interactions of alkali metal cations and their isoelectronic halide counterparts with up to six solvating Ar atoms (M⁺Ar_n and X^-Ar_n, where M = Li, Na, K, and Rb; X = H, F, Cl, and Br; and n = 1-6) via full geometry optimizations with the MP2 method and robust, correlation-consistent quadruple-ζ (QZ) basis sets. 116 unique M⁺Ar_n and X^-Ar_n stationary points have been characterized on the MP2/QZ potential energy surface. To the best of our knowledge, approximately two dozen of these stationary points have been reported here for the first time. Some of these new structures are either the lowest-energy stationary point for a particular cluster or energetically competitive with it. The CCSD(T) method was employed to perform additional single-point energy computations upon all MP2/QZ-optimized structures using the same basis set. CCSD(T)/QZ results indicate that internally solvated structures with the ion at/near the geometric center of the cluster have appreciably higher energies than those placing the ion on the periphery. While this study extends the prior investigations of M⁺Ar_n clusters found within the literature, it notably provides one of the first thorough characterizations of and comparisons to the corresponding negatively charged X^-Ar_n clusters.

Predicted M(H2)12n+ (M = Ac, Th, Pa, U, La and n = 3, 4) complexes with twenty-four hydrogen atoms bound to the metal ion

Herein, we have shown that La(iii), Ac(iii), Th(iii), Th(iv), Pa(iv) and U(iv) can directly bind with a maximum of 24 hydrogen atoms in M(H2)12 in the first sphere of coordination, which would be a new record in any metal-hydrogen complex investigated at the molecular level, where all the hydrogen atoms are directly connected to the central metal ion through M-η2(H2) bonds. Moreover, Ac(H2)n3+ (n = 9-12) systems satisfy the 18-electron rule.

Competitive bond breaking in floppy molecular trimers: HeNeH and HeNeH- calculations

Calculations for the J=0 bound states associated with the two title molecular systems are carried out using distributed Gaussian function approach and using the sum of two-body potential approximation to describe the overall interactions. The results yield one bound state for the neutral trimer and three bound states for the anionic triatom. The relative values of dissociation energies place the first dissociation threshold at the H emission for the neutral complex and at the He emission for the corresponding anion. The general spatial properties of the various bound states are analyzed on both systems.

Rovibrational Structures in Floppy Triatomics: Distributed Gaussian Functions Treatment for the Ne2H- System

The full sequence of the bound states for a very floppy triatomic complex, Ne2H- in its ground electronic state, are initially computed for the rotationless situation and employing a variational approach that expands the total nuclear wave function over a large set of symmetry-adapted, distributed Gaussian functions and employs accurate atom-atom potential energy data. The results are tested for numerical convergence, compared with the behavior of both its diatomic fragments, Ne2 and NeH-, and further compared with the results for the Ne3 case. The computational analysis is extended to the production of the rotational constants for the very nonclassical ground state vibrational configuration by making use of the previous findings. The method is shown to provide us with several illuminating details on the nanoscopic internal dynamics of this very weakly bound quantum aggregate.

Predicting atomic dopant solvation in helium clusters: The MgHen case

We present a quantum Monte Carlo study of the solvation and spectroscopic properties of the Mg-doped helium clusters MgHe(n) with n=2-50. Three high-level [MP4, CCSD(T), and CCSDT] MgHe interaction potentials have been used to study the sensitivity of the dopant location on the shape of the pair interaction. Despite the similar MgHe well depth, the pair distribution functions obtained in the diffusion Monte Carlo simulations markedly differ for the three pair potentials, therefore indicating different solubility properties for Mg in He(n). Moreover, we found interesting size effects for the behavior of the Mg impurity. As a sensitive probe of the solvation properties, the Mg excitation spectra have been simulated for various cluster sizes and compared with the available experimental results. The interaction between the excited 1P Mg atom and the He moiety has been approximated using the diatomics-in-molecules method and the two excited 1pi and 1sigma MgHe potentials. The shape of the simulated MgHe50 spectra shows a substantial dependency on the location of the Mg impurity, and hence on the MgHe pair interaction employed. To unravel the dependency of the solvation behavior on the shape of the computed potentials, exact density-functional theory has been adapted to the case of doped He(n) and various energy distributions have been computed. The results indicate the shape of the repulsive part of the MgHe potential as an important cause of the different behaviors.

Microsolvation of LiH+ in helium clusters: many-body effects and additivity models for the interaction forces

The ab initio calculation of the interaction forces between the LiH+ molecular ion, at its equilibrium geometry, and several He atoms is carried out in order to isolate and assess the importance of many-body contributions in the search for realistic energy and geometry data. The full potential energy surface (PES) with a single helium partner is obtained first by using an aug-cc-pVQZ basis set for He and higher quality ones for Li and H. The calculations were performed at the CAS-SCF plus MRCI level for the lowest potential energy surface over a total of 480 grid points of the two intermolecular Jacobi coordinates, whereas the excited state surface has also been examined in order to exclude the presence of any significant nonadiabatic interaction between the two PESs. A numerical fit of the lower surface is presented and the general physical changes of the ionic interaction when going from the lower to the upper of the two potentials are described and discussed. The fairly limited importance of many-body effects for such systems is seen from further ab initio calculations including several He atoms: our results suggest that, at least in the present case, no strong charge migration occurs after He attachment, and therefore, one could realistically model larger clusters by implementing a sum-of-potentials approach via the presently computed PES.

A complete configurational study for the bound states of Ne trimers

The structural properties and the energetics of the ground and the excited bound states of Ne(3) for zero total angular momentum are examined using different modelings for the two-body interactions. We employ a method consisting of a variational approach with a distributed Gaussian functions (DGF) basis set expansion. We discuss at length the advantages and possible limitations of such an approach, comparing it to other methods which have been applied in the literature to the same system. The DGF method turns out to be very accurate in giving us the bound states energetics and also provides in a natural way a convincing pictorial description of all the states, including those with dominant linear configurations. Additional bound states are found for the Ne(3) system with respect to those indicated in previous works and we suggest a "stabilization" procedure that can be used to assess the truly bound nature of a state. Some considerations on the relative reliability of the examined two-body interactions are also reported.

Replacement equivalence of H- and argon in small (Ar)nH- clusters from optimized structure calculations

The structural properties of some of the smaller ionic clusters of argon atoms containing the atomic impurity H-, ArnH- with n from 2 up to 7, are examined using different modeling for the interactions within each cluster and by employing different theoretical treatments, both classical and quantum, for the energetics. The same calculations are also carried out for the corresponding neutral homogeneous clusters Ar(n+1). The results of the calculations, the physical reliability of the interactions modeling, and the similarities and the difference between the anionic and the neutral complexes are discussed in some detail. The emerging picture shows that, due to specific features of the employed atom-atom potentials, the ArnH- and Ar(n+1) clusters present very similar structures, where the H- dopant substitutes for one of the outer Ar atoms but does not undergo as yet solvation within such small clusters.