F. Sebastianelli

Coupled translation-rotation eigenstates of H2 in C60 and C70 on the spectroscopically optimized interaction potential: Effects of cage anisotropy on the energy level structure and assignments

We have developed a quantitatively accurate pairwise additive five-dimensional (5D) potential energy surface (PES) for H(2) in C(60) through fitting to the recently published infrared (IR) spectroscopic measurements of this system for H(2) in the vibrationally excited nu=1 state. The PES is based on the three-site H(2)-C pair potential introduced in this work, which in addition to the usual Lennard-Jones (LJ) interaction sites on each H atom of H(2) has the third LJ interaction site located at the midpoint of the H-H bond. For the optimal values of the three adjustable parameters of the potential model, the fully coupled quantum 5D calculations on this additive PES reproduce the six translation-rotation (T-R) energy levels observed so far in the IR spectra of H(2)@C(60) to within 0.6%. This is due in large part to the greatly improved description of the angular anisotropy of the H(2)-fullerene interaction afforded by the three-site H(2)-C pair potential. The same H(2)-C pair potential spectroscopically optimized for H(2)@C(60) was also used to construct the pairwise additive 5D PES of H(2) (nu=1) in C(70). This PES, because of the lower symmetry of C(70) (D(5h)) relative to that of C(60) (I(h)), exhibits pronounced anisotropy with respect to the direction of the translational motion of H(2) away from the cage center, unlike that of H(2) in C(60). As a result, the T-R energy level structure of H(2) in C(70) from the quantum 5D calculations on the optimized PES, the quantum numbers required for its assignment, and the degeneracy patterns which arise from the T-R coupling for translationally excited H(2) are all qualitatively different from those determined previously for H(2)@C(60) [M. Xu et al., J. Chem. Phys. 128, 011101 (2008).

Quantum dynamics of H2, D2, and HD in the small dodecahedral cage of clathrate hydrate: Evaluating H2-water nanocage interaction potentials by comparison of theory with inelastic neutron scattering experiments

We have performed rigorous quantum five-dimensional (5D) calculations and analysis of the translation-rotation (T-R) energy levels of one H(2), D(2), and HD molecule inside the small dodecahedral (H(2)O)(20) cage of the structure II clathrate hydrate, which was treated as rigid. The H(2)- cage intermolecular potential energy surface (PES) used previously in the molecular dynamics simulations of the hydrogen hydrates [Alavi et al., J. Chem. Phys. 123, 024507 (2005)] was employed. This PES, denoted here as SPC/E, combines an effective, empirical water-water pair potential [Berendsen et al., J. Phys. Chem. 91, 6269 (1987)] and electrostatic interactions between the partial charges placed on H(2)O and H(2). The 5D T-R eigenstates of HD were calculated also on another 5D H(2)-cage PES denoted PA-D, used by us earlier to investigate the quantum T-R dynamics of H(2) and D(2) in the small cage [Xu et al., J. Phys. Chem. B 110, 24806 (2006)]. In the PA-D PES, the hydrogen-water pair potential is described by the ab initio 5D PES of the isolated H(2)-H(2)O dimer. The quality of the SPC/E and the PA-D H(2)-cage PESs was tested by direct comparison of the T-R excitation energies calculated on them to the results of two recent inelastic neutron scattering (INS) studies of H(2) and HD inside the small clathrate cage. The translational fundamental and overtone excitations, as well as the triplet splittings of the j=0-->j=1 rotational transitions, of H(2) and HD in the small cage calculated on the SPC/E PES agree very well with the INS results and represent a significant improvement over the results computed on the PA-D PES. Our calculations on the SPC/E PES also make predictions about several spectroscopic observables for the encapsulated H(2), D(2), and HD, which have not been measured yet.

Quantum dynamics of coupled translational and rotational motions of H2 inside C60

We report rigorous quantum calculations of the translation-rotation (T-R) eigenstates of the H2 molecule in C60. The resulting level structure can be explained in terms of a few dominant features. These include the coupling between the orbital and the rotational angular momenta of H2 to give the total angular momentum lambda, and the splitting of the sevenfold degeneracy of T-R levels with lambda=3 by the nonsphericity of C60, according to the rules of the icosahedral I h group.

Ring-breaking electron attachment to uracil: Following bond dissociations via evolving resonances

Calculations are carried out at various distinct energies to obtain both elastic cross sections and S-matrix resonance indicators (poles) from a quantum treatment of the electron scattering from gas-phase uracil. The low-energy region confirms the presence of pi(*) resonances as revealed by earlier calculations and experiments which are compared with the present findings. They turn out to be little affected by bond deformation, while the transient negative ions (TNIs) associated with sigma(*) resonances in the higher energy region ( approximately 8 eV) indeed show that ring deformations which allow vibrational redistribution of the excess electron energy into the molecular target strongly affect these shape resonances: They therefore evolve along different dissociative pathways and stabilize different fragment anions. The calculations further show that the occurrence of conical intersections between sigma(*) and pi(*)-type potential energy surfaces (real parts) is a very likely mechanism responsible for energy transfers between different TNIs. The excess electron wavefunctions for such scattering states, once mapped over the molecular space, provide nanoscopic reasons for the selective breaking of different bonds in the ring region.

Quantum dynamics of small H2 and D2 clusters in the large cage of structure II clathrate hydrate: Energetics, occupancy, and vibrationally averaged cluster structures

We report diffusion Monte Carlo (DMC) calculations of the quantum translation-rotation (T-R) dynamics of one to five para-H(2) (p-H(2)) and ortho-D(2) (o-D(2)) molecules inside the large hexakaidecahedral (5(12)6(4)) cage of the structure II clathrate hydrate, which was taken to be rigid. These calculations provide a quantitative description of the size evolution of the ground-state properties, energetics, and the vibrationally averaged geometries, of small (p-H(2))(n) and (o-D(2))(n) clusters, n=1-5, in nanoconfinement. The zero-point energy (ZPE) of the T-R motions rises steeply with the cluster size, reaching 74% of the potential well depth for the caged (p-H(2))(4). At low temperatures, the rapid increase of the cluster ZPE as a function of n is the main factor that limits the occupancy of the large cage to at most four H(2) or D(2) molecules, in agreement with experiments. Our DMC results concerning the vibrationally averaged spatial distribution of four D(2) molecules, their mean distance from the cage center, the D(2)-D(2) separation, and the specific orientation and localization of the tetrahedral (D(2))(4) cluster relative to the framework of the large cage, agree very well with the low-temperature neutron diffraction experiments involving the large cage with the quadruple D(2) occupancy.

Hydrogen Molecules inside Fullerene C70: Quantum Dynamics, Energetics, Maximum Occupancy, And Comparison with C60

Recent synthesis of the endohedral complexes of C(70) and its open-cage derivative with one and two H(2) molecules has opened the path for experimental and theoretical investigations of the unique dynamic, spectroscopic, and other properties of systems with multiple hydrogen molecules confined inside a nanoscale cavity. Here we report a rigorous theoretical study of the dynamics of the coupled translational and rotational motions of H(2) molecules in C(70) and C(60), which are highly quantum mechanical. Diffusion Monte Carlo (DMC) calculations were performed for up to three para-H(2) (p-H(2)) molecules encapsulated in C(70) and for one and two p-H(2) molecules inside C(60). These calculations provide a quantitative description of the ground-state properties, energetics, and the translation-rotation (T-R) zero-point energies (ZPEs) of the nanoconfined p-H(2) molecules and of the spatial distribution of two p-H(2) molecules in the cavity of C(70). The energy of the global minimum on the intermolecular potential energy surface (PES) is negative for one and two H(2) molecules in C(70) but has a high positive value when the third H(2) is added, implying that at most two H(2) molecules can be stabilized inside C(70). By the same criterion, in the case of C(60), only the endohedral complex with one H(2) molecule is energetically stable. Our results are consistent with the fact that recently both (H(2))(n)@C(70) (n = 1, 2) and H(2)@C(60) were prepared, but not (H(2))(3)@C(70) or (H(2))(2)@C(60). The ZPE of the coupled T-R motions, from the DMC calculations, grows rapidly with the number of caged p-H(2) molecules and is a significant fraction of the well depth of the intermolecular PES, 11% in the case of p-H(2)@C(70) and 52% for (p-H(2))(2)@C(70). Consequently, the T-R ZPE represents a major component of the energetics of the encapsulated H(2) molecules. The inclusion of the ZPE nearly doubles the energy by which (p-H(2))(3)@C(70) is destabilized and increases by 66% the energetic destabilization of (p-H(2))(2)@C(60). For these reasons, the T-R ZPE has to be calculated accurately and taken into account for reliable theoretical predictions regarding the stability of the endohedral fullerene complexes with hydrogen molecules and their maximum H(2) content.

Coupled Translation-Rotation Eigenstates of H2, HD, and D2 in the Large Cage of Structure II Clathrate Hydrate: Comparison with the Small Cage and Rotational Raman Spectroscopy †

We report fully coupled quantum five-dimensional calculations of the translation-rotation (T-R) energy levels of one H(2), HD, and D(2) molecule confined inside the large hexakaidecahedral (5(12)6(4)) cage of the structure II clathrate hydrate. Highly converged T-R eigenstates have been obtained for excitation energies beyond the j = 2 rotational levels of the guest molecules, in order to allow comparison with the recent Raman spectroscopic measurements. The translationally excited T-R states are assigned with the quantum numbers n and l of the 3D isotropic harmonic oscillator. However, the translational excitations are not harmonic, since the level energies depend not only on n but also on l. For l > 1, the T-R levels having the same n,l values are split into groups of almost degenerate levels. The splitting patterns follow the predictions of group theory for the environment of T(d) symmetry, which is created by the configuration of the oxygen atoms of the large cage. The 2j + 1 degeneracy of the j = 1 and 2 rotational levels of the encapsulated hydrogen molecule is lifted entirely by the angular anisotropy of the H(2)-cage interaction potential. The patterns and magnitudes of the j = 1, 2 rotational level splittings, and the energies of the sublevels, in the large cage are virtually identical with those calculated for the small cage. This is in agreement with, and sheds light on, the observation that the S(0)(0) (j = 0-->2) bands in the rotational Raman spectra measured for simple H(2) hydrate and the binary hydrate of H(2) with tetrahydrofuran are remarkably similar with respect to their frequencies, widths, shapes, and internal structure, when the H(2) occupancy of the large cage of simple H(2) hydrate is low.

Electron scattering cross sections from HCN over a broad energy range (0.1–10 000 eV): Influence of the permanent dipole moment on the scattering process

We report theoretical integral and differential cross sections for electron scattering from hydrogen cyanide derived from two ab initio scattering potential methods. For low energies (0.1-100 eV), we have used the symmetry adapted-single centre expansion method using a multichannel scattering formulation of the problem. For intermediate and high energies (10-10,000 eV), we have applied an optical potential method based on a screening corrected independent atom representation. Since HCN is a strong polar molecule, further dipole-induced excitations have been calculated in the framework of the first Born approximation and employing a transformation to a space-fixed reference frame of the calculated K-matrix elements. Results are compared with experimental data available in the literature and a complete set of recommended integral elastic, inelastic, and total scattering cross sections is provided from 0.1 to 10,000 eV.

Dissociative Electron Attachment to Formamide: Direct and Indirect Pathways from Resonant Intermediates.

Dissociative electron attachment (DEA) to formamide (HCONH2), the smallest molecule with a peptide bond, is investigated with electron-molecule scattering calculations. At the equilibrium geometry we identify two resonances of A′′ and A′ symmetry at 3.77 and 14.90 eV, respectively. To further assess potential bond-breaking pathways for the transient negative ions (TNIs), the behavior of the resonances upon bond stretching of the C−H and C−N bond is investigated. While along the C−H dissociation coordinate neither resonance changes significantly, we find instead that both resonances are stabilized upon stretching the peptide C−N bond, with their resonance energy and width coming down rapidly, most strongly so for the A′ resonance. The A′ resonance is thus seen to disappear when the C−N bond is stretched for more than 1 A, where it presumably smoothly connects to a bound anion state, a direct DEA pathway for the A′ TNI to yield NH2− and HCO. The A′′ resonance is found instead not to be purely dissociative a...

Quantum and classical structures for 4He clusters with the H− impurity

The structural properties and the energetics of some of the smaller HenH− clusters, with n varying from 2 to 14, are examined both with classical and quantum treatments, making also a comparison with the corresponding neutral systems Hen. The results of the calculations, the physical reliability of the employed interaction modeling, and the comparison with previous results are discussed. The emerging picture shows that for clusters of this size the dopant ion H− always locates itself outside the Hen moiety.