Quantum Features of Anionic Species He*⁻ and He2*⁻ in Small He(N) Clusters.

We present variational calculations on systems containing a few boson helium atoms attached to electronically excited atomic and molecular helium anions He*⁻ and He2*⁻ and characterize their structures and energetics. Previously reported high-level ab initio results [Huber, S. E.; Mauracher, A. Mol. Phys. 2014, 112, 794] to describe the interactions between excited (metastable) anions and a neutral He atom have been employed. For the case of the atomic species He*⁻, the corresponding interaction with He suggests large anharmonicity effects due to the presence of a deep well of ∼17,500 cm⁻¹ at short distances, together with a more external shallow secondary well of ∼4 cm⁻¹, both supporting bound levels. Moreover, when a sum of pairwise interactions is assumed to describe the full PES corresponding to the presence of several neutral He atoms, geometrical constraints already predict the complete solvation of the anionic impurity by six helium atoms, giving rise to a bipyramidal structure. In turn, for the anisotropic weak interaction He-He2*⁻, where the anionic dimer is considered as a rigid rotor, the obtained structures show the tendency of the helium atoms to pack themselves together and largely far away from the dopant, thereby confirming the heliophobic character of He2*⁻.

Reactive scattering calculations for (87)Rb+(87)RbHe→Rb2((3)Σ(u)(+),v)+He from ultralow to intermediate energies.

We investigate atom-diatom reactive collisions, as a preliminary step,in order to assess the possibility of forming Rb(2) molecules in their lowest triplet electronic state by cold collisions of rubidium atoms on the surface of helium nanodroplets [corrected]. A simple model related to the well-known Rosen treatment of linear triatomic molecules [N. Rosen, J. Chem. Phys. 1, 319 (1933)] in relative coordinates is used, allowing to estimate reactive probabilities for different values of the total angular momentum. The best available full dimensional potential energy surface [Guillon et al., J. Chem. Phys. 136, 174307 (2012)] is employed through the calculations. Noticeable values of the probabilities in the ultracold regime, which numerically fulfill the Wigner threshold law, support the feasibility of the process. The rubidium dimer is mainly produced at high vibrational states, and the reactivity is more efficient for a bosonic helium partner than when the fermion species is considered.

Computational study of interactions and nuclear magnetic shielding constants in linear chains of formamide clusters.

We investigated the energetic, structural, dielectric, and nuclear magnetic shielding properties of linear n-formamide clusters, with n up to 6, to quantitatively characterize cooperative effects in model biological systems. The geometries of the complexes were optimized at the MP2 and DFT/B3LYP levels by using the pc-2 and pc-3 basis sets, while the nuclear magnetic shielding constants were calculated by employing pcS-n type basis sets, which have been optimized specifically for density functional calculations of these properties. The interaction energies show the cooperative effect, which favors the successive addition of monomers. In addition, by analyzing structural changes in the intermolecular C=O, C-N and hydrogen O⋯H bonds, as well as in the average dipole moments as cluster size increases, we found that the cooperative interaction far exceeds that expected for electrostatic interactions. Such non-pairwise-additive effects are also reflected in the changes of the nuclear magnetic shielding constants. In particular, the negativity of O shielding decreases around 23% from the monomer to the 6-formamide chain. It is possible to note the decrease in the shielding of H and in the deshielding of O as a result of their hydrogen bonding. However, the results obtained show that these variations in the extremes of formamide chains tend to zero, and the respective shielding values tend to stabilize as the number of monomers increases in the chain. Also, the cooperative effect increases in the middle of the chains, by decreasing the shielding for all atoms except that of O, which decreases its deshielding. These results could serve to guide improvements in current conventional models for simulating hydrogen bonded systems.

Simulating liquid water for determining its structural and transport properties.

Molecular dynamics simulations are carried out for calculating structural and transport properties of pure liquid water, such as radial distribution functions and self-diffusion and viscosity coefficients, respectively. We employed reparameterized versions of the ab initio water potential by Niesar, Clementi and Corongiu (NCC). In order to investigate the role of the electrostatic contribution, the partial charges of the NCC model are adjusted so that to reproduce the dipole moment values of the SPC/E, SPC/Fw and TIP4P/2005 water models. The single and collective transport coefficients are obtained by employing the Green-Kubo relations at various temperatures. Additionally, in order to overcome convergence difficulties arising from the long correlation times of the stress-tensor autocorrelation functions, a previously reported fitting scheme was employed. The present results indicate that there is a significant relationship between the dipole moment value of the model, and the calculated transport coefficients. We found that by adjusting the molecular dipole moment of the NCC to the value of the TIP4P/2005, the obtained values for the self-diffusion and viscosity coefficients are in better agreement with experiment, compared to the values obtained with the original NCC model. Even though the predictions of the present model exhibits an overall correct behavior, we conclude that further improvements are still required. In order to achieve that, a careful reparameterization of the repulsion-dispersion terms of the potential model is proposed. Also, the effect of the inclusion of many-body effects such as polarizability, should also be investigated.

Including nuclear quantum effects into highly correlated electronic structure calculations of weakly bound systems.

An interface between the APMO code and the electronic structure package MOLPRO is presented. The any particle molecular orbital APMO code [González et al., Int. J. Quantum Chem. 108, 1742 (2008)] implements the model where electrons and light nuclei are treated simultaneously at Hartree-Fock or second-order Möller-Plesset levels of theory. The APMO-MOLPRO interface allows to include high-level electronic correlation as implemented in the MOLPRO package and to describe nuclear quantum effects at Hartree-Fock level of theory with the APMO code. Different model systems illustrate the implementation: (4)He2 dimer as a protype of a weakly bound van der Waals system; isotopomers of [He-H-He](+) molecule as an example of a hydrogen bonded system; and molecular hydrogen to compare with very accurate non-Born-Oppenheimer calculations. The possible improvements and future developments are outlined.

Solvent states and spectroscopy of doped helium clusters as a quantum-chemistry-like problem.

The Full-Configuration-Interaction Nuclear-Orbital (FCI-NO) approach [J. Chem. Phys., 2009, 131, 19401], as the implementation of the quantum-chemistry ansatz, is overviewed and applied to (He)N-Cl2(X) clusters (N≤ 4). The ground and excited states of both fermionic (3)He and bosonic (4)He [see also, J. Phys. Chem. Lett., 2012, 2, 2145] clusters are studied. It is shown that the FCI-NO approach allows us to overcome three main difficulties: (1) the Fermi-Dirac (Bose-Einstein) nuclear statistics; (2) the wide (highly anharmonic) amplitudes of the He-dopant and He-He motions; and (3) both the weakly attractive (long-range) and the strongly repulsive (short-range) interaction between the helium atoms. Special emphasis is placed on the dependence of the cluster properties on the number of helium atoms, and on the comparison between the two helium isotopes. In particular, we analyze the analogies between quantum rings comprising electrons and (3)He atoms. The synthetic vibro-rotational Raman spectra of Cl2(X) immersed in ((3,4)He)N clusters (N≤ 4) are discussed as a function of the cluster size and the nuclear statistics. It is shown that the Coriolis couplings play a key role in modifying the spectral dopant profile in (3)He. Finally, we point out possible directions for future research using the quantum-chemistry ansatz.

Vibrational dynamics of the H5(+) and its isotopologues from multiconfiguration time-dependent Hartree calculations.

Full-dimensional multiconfiguration time-dependent Hartree (MCTDH) computations are reported for the vibrational states of the H(5)(+) and its H(4)D(+), H(3)D(2)(+), H(2)D(3)(+), HD(4)(+), D(5)(+) isotopologues employing two recent analytical potential energy surfaces of Xie et al. [J. Chem. Phys. 122, 224307 (2005)] and Aguado et al. [J. Chem. Phys. 133, 024306 (2010)]. The potential energy operators are constructed using the n-mode representation adapted to a four-combined mode cluster expansion, including up to seven-dimensional grids, chosen adequately to take advantage in representing the MCTDH wavefunction. An error analysis is performed to quantify the convergence of the potential expansion to reproduce the reference surfaces at the energies of interest. An extensive analysis of the vibrational ground state properties of these isotopes and comparison with the reference diffusion Monte Carlo results by Acioli et al. [J. Chem. Phys. 128, 104318 (2008)] are presented. It is found that these systems are highly delocalized, interconverting between equivalent minima through rotation and internal proton transfer motions even at their vibrational ground state. Isotopic substitution affects the zero-point energy and structure, showing preference in the arrangements of the H and D within the mixed clusters, and the most stable conformers of each isotopomer are the ones with the H in the central position. Vibrational excited states are also computed and by comparing the energies and structures predicted from the two surfaces, the effect of the potential topology on them is discussed.

Theoretical investigation of the He4Br2 conformers.

Full dimensional quantum dynamics calculations of the three lowest isomers of the He(4)Br(2) van der Waals molecule in its ground electronic state are reported. The calculations are performed using the multiconfiguration time-dependent Hartree (MCTDH) method and a realistic potential form that includes the sum of three body ab initio coupled-cluster single double triple [CCSD(T)] He-Br(2) interactions plus the He-He and Br-Br interactions. This potential exhibits several multiple minima, with the three lowest ones lying very close in energy, just within 2 cm(-1). Such small differences are also found in the calculated binding energies of the three most stable conformers, indicating the floppiness of the system and, thus, the need of accurate potential forms and quantum full dynamics methods to treat this kind of complexes. The 12 dimensional results reported in this work present benchmark data and, thus, can serve to evaluate approximate methods aiming to describe higher order rare gas-dihalogen (N > 4) complexes. A comparison with previous studies using different potential forms and approaches to the energetics for the He(4)Br(2) cluster is also presented.

Quantum-dynamics study of the H5+ cluster: full dimensional benchmark results on its vibrational states.

A full-dimensional quantum dynamics study is carried out for the highly fluxional H(5)(+) cation on a recent reference potential energy surface by using the multi configuration time-dependent Hartree method. With five equivalent light atoms and shallow barriers between various low-lying stationary points on the surface, the spectroscopic characterization of H(5)(+) represents a huge challenge for accurate quantum dynamics simulations. The present calculation is the first such a study on this cation, which together with its isotope analogies are of primary importance in the interstellar chemistry. The vibrational ground state properties and several vibrationally excited states corresponding to low vibrational frequency motions, not yet directly observable by the experiment, are presented and analyzed.

Full-dimensional multi configuration time dependent Hartree calculations of the ground and vibrationally excited states of He2,3Br2 clusters.

Quantum dynamics calculations are reported for the tetra-, and penta-atomic van der Waals He(N)Br(2) complexes using the multiconfiguration time-dependent Hartree (MCTDH) method. The computations are carried out in satellite coordinates, and the kinetic energy operator in this set of coordinates is given. A scheme for the representation of the potential energy surface based on the sum of the three-body HeBr(2) interactions at CSSD(T) level plus the He-He interaction is employed. The potential surfaces show multiple close lying minima, and a quantum description of such highly floppy multiminima systems is presented. Benchmark, full-dimensional converged results on ground vibrational/zero-point energies are reported and compared with recent experimental data available for all these complexes, as well as with previous variational quantum calculations for the smaller HeBr(2) and He(2)Br(2) complexes on the same surface. Some low-lying vibrationally excited eigenstates are also computed by block improved relaxation calculations. The binding energies and the corresponding vibrationally averaged structures are determined for different conformers of these complexes. Their relative stability is discussed, and contributes to evaluate the importance of the multiple-minima topology of the underlying potential surface.

Ab initio characterization of the Ne-I2 van der Waals complex: intermolecular potentials and vibrational bound states.

A theoretical study of the potential energy surface and bound states is performed for the ground state of the NeI(2) van der Waals (vdW) complex. The three-dimensional interaction energies are obtained from ab initio coupled-cluster, coupled-cluster single double (triple)/complete basis set, calculations using large basis sets, of quadruple- through quintuple-zeta quality, in conjunction with relativistic effective core potentials for the heavy iodine atoms. For the analytical representation of the surface two different schemes, based on fitting and interpolation surface generation techniques, are employed. The surface shows a double-minimum topology for linear and T-shaped configurations. Full variational quantum mechanical calculations are carried out using the model surfaces, and the vibrationally averaged structures and energetics for the NeI(2) isomers are determined. The accuracy of the potential energy surfaces is validated by a comparison between the present results and the corresponding experimental data available. In lieu of more experimental measurements, we also report our results/predictions on higher bound vibrational vdW levels, and the influence of the employed surface on them is discussed.

Quantum features of a barely bound molecular dopant: Cs2(3Σu) in bosonic helium droplets of variable size.

We present in this work the study of small (4)He(N)-Cs(2)((3)Σ(u)) aggregates (2 ≤ N ≤ 30) through combined variational, diffusion Monte Carlo (DMC), and path integral Monte Carlo (PIMC) calculations. The full surface is modeled as an addition of He-Cs(2) interactions and He-He potentials. Given the negligible strength and large range of the He-Cs(2) interaction as compared with the one for He-He, a propensity of the helium atoms to pack themselves together, leaving outside the molecular dopant is to be expected. DMC calculations determine the onset of helium gathering at N = 3. To analyze energetic and structural properties as a function of N, PIMC calculations with no bosonic exchange, i.e., Boltzmann statistics, at low temperatures are carried out. At T = 0.1 K, although acceptable one-particle He-Cs(2) distributions are obtained, two-particle He-He distributions are not well described, indicating that the proper symmetry should be taken into account. PIMC distributions at T = 1 K already compare well with DMC ones and show minor exchange effects, although binding energies are still far from having converged in terms of the number of quantum beads. As N increases, the He-He PIMC pair correlation function shows a clear tendency to coincide with the experimental boson-liquid helium one at that temperature. It supports the picture of a helium droplet which carries the molecular impurity on its surface, as found earlier for other triplet dimers.

Internal proton transfer and H2 rotations in the H5(+) cluster: a marked influence on its thermal equilibrium state.

Classical and path integral Monte Carlo (CMC, PIMC) "on the fly" calculations are carried out to investigate anharmonic quantum effects on the thermal equilibrium structure of the H5(+) cluster. The idea to follow in our computations is based on using a combination of the above-mentioned nuclear classical and quantum statistical methods, and first-principles density functional (DFT) electronic structure calculations. The interaction energies are computed within the DFT framework using the B3(H) hybrid functional, specially designed for hydrogen-only systems. The global minimum of the potential is predicted to be a nonplanar configuration of C(2v) symmetry, while the next three low-lying stationary points on the surface correspond to extremely low-energy barriers for the internal proton transfer and to the rotation of the H2 molecules, around the C2 axis of H5(+), connecting the symmetric C(2v) minima in the planar and nonplanar orientations. On the basis of full-dimensional converged PIMC calculations, results on the quantum vibrational zero-point energy (ZPE) and state of H5(+) are reported at a low temperature of 10 K, and the influence of the above-mentioned topological features of the surface on its probability distributions is clearly demonstrated.

Theoretical characterization of intermolecular vibrational states through the multi-configuration time dependent Hartree approach: the He2,3ICl clusters.

Benchmark, full-dimensional calculations on the ground and excited vibrational states for the tetra-, and penta-atomic weakly bound He(2,3)ICl complexes are reported. The representation of the potential energy surfaces includes three-body HeICl potentials parameterized to coupled-cluster singles, doubles, and perturbative triples ab initio data. These terms are important in accurately describing the interactions of such highly floppy systems. The corresponding 6D/9D computations are performed with the multi-configuration time dependent Hartree method, using natural potential fits, and a mode combination scheme to optimize the computational effort in the improved relaxation calculations. For these complexes several low-lying vibrational states are computed, and their binding energies and radial/angular probability density distributions are obtained. We found various isomers which are assigned to different structural models related with combinations of the triatomic isomers, like linear, T-shaped, and antilinear ones. Comparison of these results with recent experimental data is presented, and the quantitative deviations found with respect to the experiment are discussed.

Toward a realistic density functional theory potential energy surface for the H5+ cluster.

The potential energy surface of H(5)(+) is characterized using density functional theory. The hypersurface is evaluated at selected configurations employing different functionals, and compared with results obtained from ab initio CCSD(T) calculations. The lowest ten stationary points (minima and saddle-points) on the surface are located, and the features of the short-, intermediate-, and long-range intermolecular interactions are also investigated. A detailed analysis of the surface's topology, and comparisons with extensive CCSD(T) results, as well as a recent ab initio analytical surface, shows that density functional theory calculations using the B3(H) functional represent very well all aspects studied on the H(5)(+) potential. These include the tiny energy difference between the minimum at 1-C(2v) configuration and the 2-D(2d) one corresponding to the transition state for the proton transfer between the two equivalent C(2v) minima, and also the correct asymptotic behavior of the long-range interactions. The calculated binding energy and dissociation enthalpies compare very well with previous benchmark coupled-cluster ab initio data, and with experimental data available. Based on these results the use of such approach to perform first-principles molecular dynamics simulations could provide reliable information regarding the dynamics of protonated hydrogen clusters.

Structuring a quantum solvent around a weakly bound dopant: the He-Cs2(3Sigma(u)) complex.

The structure and energetics of (3,4)HeCs(2)((3)Sigma(u)) molecules are analyzed from first principles. Fixing the cesium dimer at its equilibrium distance, the electronic structure was determined through ab initio methods at the CCSD(T) level of theory using a large basis set to compute the interaction energies. At the T-shaped geometry, there is a shallow well with a depth of approximately 2 cm(-1) placed at R approximately 6.75 A, R being the distance from the center of mass of Cs(2) to He. That depth gradually decreases to approximately 0.75 cm(-1), while R increases to about 11.5 A at linear arrangements. A simple model of adding atom-atom Lennard-Jones potentials with well-depth and equilibrium distance parameters depending on the angular orientation was found to accurately reproduce the ab initio points. Using this analytical form, variational calculations at zero total angular momentum are performed, predicting a single bound level at approximately -0.106 (approximately -0.042) cm(-1) for the boson (fermion) species. Further calculations using Quantum Monte Carlo methods are carried out and found to be in good agreement with the variational ones. On the basis of the present results, such analytical expression could in turn be used to describe the structure and binding of larger complexes and therefore opens the possibility to further studies involving such aggregates.

Intermolecular ab initio potential and spectroscopy of the ground state of HeI2 complex revisited.

The structure, energetics, and spectroscopy of ground-state HeI(2) molecule are analyzed from first principles. Ab initio methodology at CCSD(T) level of theory was employed, and large basis sets were used to compute the interaction energies. Scalar relativistic effects were accounted for by relativistic effective core potentials for the iodine atoms. Recent experimental investigations of the HeI(2) rovibronic spectra have estimated the ground-state binding energies of 16.6 +/- 0.6 and 16.3 +/- 0.6 cm(-1) for the T-shaped and linear isomers, respectively. Given the extremely small difference between the two conformers, special attention was paid in the choice of basis sets used and the extrapolation schemes employed, as well as the fitting process for its analytical representation. The complete analytical form is provided, and variational fully quantum mechanical calculations were carried out by using the new parametrized surface, to evaluate vibrationally averaged structures and binding energies for the different conformers. The results obtained are in good accord with recent data available from experimental investigations of the He-I(2) rovibronic spectra.

Ab initio vibrational predissociation dynamics of He-I2(B) complex.

Three-dimensional quantum mechanical calculations on the vibrational predissociation dynamics of HeI2 B state complex are performed using a potential energy surface accurately fitted to unrestricted open-shell coupled cluster ab initio data, further enabling extrapolation for large I2 bond lengths. A Lanczos iterative method with an optimized complex absorbing potential is used to determine energies and lifetimes of the vibrationally predissociating He,I2(B,v') complex for v'