Isomer dependence of HF vibrational frequency shift for ArnHF (n=4–14) van der Waals clusters: Quantum five‐dimensional bound state calculations

Automated Construction of Potential Energy Surfaces Suitable to Describe van der Waals Complexes with Highly Excited Nascent Molecules: The Rotational Spectra of Ar-CS(v) and Ar-SiS(v)

Some reactions produce extremely hot nascent products which nevertheless can form sufficiently long-lived van der Waals (vdW) complexes-with atoms or molecules from a bath gas-as to be observed via microwave spectroscopy. Theoretical calculations of such unbound resonance states can be much more challenging than ordinary bound-state calculations depending on the approach employed. One encounters not just the floppy, and perhaps multiwelled potential energy surface (PES) characteristic of vdWs complexes, but in addition, one must contend with excitation of the intramolecular modes and its corresponding influence on the PES. Straightforward computation of the (resonance) rovibrational levels of interest, involves the added complication of the unbound nature of the wave function, often treated with techniques such as introducing a complex absorbing potential. Here, we have demonstrated that a simplified approach of making a series of vibrationally effective PESs for the intermolecular coordinates-one for each reaction product vibrational quantum number of interest-can produce vdW levels for the complex with spectroscopic accuracy. This requires constructing a series of appropriately weighted lower-dimensional PESs for which we use our freely available PES-fitting code AUTOSURF. The applications of this study are the Ar-CS and Ar-SiS complexes, which are isovalent to Ar-CO and Ar-SiO, the latter of which we considered in a previously reported study. Using a series of vibrationally effective PESs, rovibrational levels and predicted microwave transition frequencies for both complexes were computed variationally. A series of shifting rotational transition frequencies were also computed as a function of the diatom vibrational quantum number. The predicted transitions were used to guide and inform an experimental effort to make complementary observations. Comparisons are given for the transitions that are within the range of the spectrometer and were successfully recorded. Calculations of the rovibrational level pattern agree to within 0.2% with experimental measurements.

Dynamics of H2 adsorbed in porous materials as revealed by computational analysis of inelastic neutron scattering spectra

This perspective article reviews the different types of quantum and classical mechanical methods that have been implemented to interpret the INS spectra for H₂ adsorbed in porous materials.

HD in C₆₀: theoretical prediction of the inelastic neutron scattering spectrum and its temperature dependence

We report rigorous quantum calculations of the inelastic neutron scattering (INS) spectra of HD@C₆₀, over a range of temperatures from 0 to 240 K and for two incident neutron wavelengths used in recent experimental investigations. The computations were performed using our newly developed methodology, which incorporates the coupled five-dimensional translation-rotation (T-R) eigenstates of the guest molecule as the initial and final states of the INS transitions, and yields highly detailed spectra. Depending on the incident neutron wavelength, the number of computed INS transitions varies from almost 500 to over 2000. The low-temperature INS spectra display the fingerprints of the coupling between the translational and rotational motions of the entrapped HD molecule, which is responsible for the characteristic splitting patterns of the T-R energy levels. INS transitions from the ground T-R state of HD to certain sublevels of excited T-R multiplets have zero intensity and are absent from the spectra. This surprising finding is explained by the new INS selection rule introduced here. The calculated spectra exhibit strong temperature dependence. As the temperature increases, numerous new peaks appear, arising from the transitions originating in excited T-R states which become populated. Our calculations show that the higher temperature features typically comprise two or more transitions close in energy and with similar intensities, interspersed with numerous other transitions whose intensities are negligible. This implies that accurately calculated energies and intensities of INS transitions which our methodology provides will be indispensable for reliable interpretation and assignment of the experimental spectra of HD@C₆₀ and related systems at higher temperatures.

Wave function delocalization and large-amplitude vibrations of helium on corrugated aromatic microsurfaces: tetracene.He and pentacene.He van der Waals complexes

We report accurate quantum three-dimensional calculations of highly excited intermolecular vibrational states of the van der Waals (vdW) complexes tetracene.He and pentacene.He in the S1 excited electronic state. The aromatic molecules were taken to be rigid and the intermolecular potential energy surfaces (IPESs) were modeled as a sum of atom-atom Lennard-Jones pair potentials. The IPESs are corrugated in the direction of the long (x) axis of the aromatic molecules, due to the presence of the symmetrically equivalent global double minimum for tetracene.He, and a triple minimum (central global minimum and two equivalent local minima) for pentacene.He, on each side of the aromatic plane. Both IPESs have two additional minor equivalent local minima further away from the center of the molecule. The vdW vibrational states analyzed in this work cover about 80% of the well depths of the IPESs. The mode coupling is generally weak for those states whose out-of-plane (z) mode is unexcited. However, the z-mode fundamental is strongly coupled to the short-axis (y) in-plane mode, so that the pure z-mode excitation could not be identified. The He atom exhibits large in-plane spatial delocalizaton already in the ground vdW vibrational state, which increases rapidly upon the excitation of the in-plane x and y modes, with little hindrance by the corrugation of the aromatic microsurfaces. For the vdW vibrational energies considered, the He atom spatial delocalization reaches Deltax and Deltay values of approximately 5 and 4 A, respectively, and is limited only by the finite size of the aromatic substrates. Side-crossing delocalization of the wave functions on both sides of the molecular plane is found at excitation energies >30 cm(-1), giving rise to the energy splittings of the pairs of states symmetric/antisymmetric with respect to the aromatic plane; the splittings show strong vdW vibrational mode specificity.

One and Two Hydrogen Molecules in the Large Cage of the Structure II Clathrate Hydrate: Quantum Translation−Rotation Dynamics Close to the Cage Wall

We have performed a rigorous theoretical study of the quantum translation-rotation (T-R) dynamics of one and two H2 and D2 molecules confined inside the large hexakaidecahedral (5(12)6(4)) cage of the sII clathrate hydrate. For a single encapsulated H2 and D2 molecule, accurate quantum five-dimensional calculations of the T-R energy levels and wave functions are performed that include explicitly, as fully coupled, all three translational and the two rotational degrees of freedom of the hydrogen molecule, while the cage is taken to be rigid. In addition, the ground-state properties, energetics, and spatial distribution of one and two p-H2 and o-D2 molecules in the large cage are calculated rigorously using the diffusion Monte Carlo method. These calculations reveal that the low-energy T-R dynamics of hydrogen molecules in the large cage are qualitatively different from that inside the small cage, studied by us recently. This is caused by the following: (i) The large cage has a cavity whose diameter is about twice that of the small cage for the hydrogen molecule. (ii) In the small cage, the potential energy surface (PES) for H2 is essentially flat in the central region, while in the large cage the PES has a prominent maximum at the cage center, whose height exceeds the T-R zero-point energy of H2/D2. As a result, the guest molecule is excluded from the central part of the large cage, its wave function localized around the off-center global minimum. Peculiar quantum dynamics of the hydrogen molecule squeezed between the central maximum and the cage wall manifests in the excited T-R states whose energies and wave functions differ greatly from those for the small cage. Moreover, they are sensitive to the variations in the hydrogen-bonding topology, which modulate the corrugation of the cage wall.

Hydrogen Molecule in the Small Dodecahedral Cage of a Clathrate Hydrate: Quantum Five-Dimensional Calculations of the Coupled Translation−Rotation Eigenstates

We report quantum five-dimensional (5D) calculations of the energy levels and wave functions of the hydrogen molecule, para-H2 and ortho-H2, confined inside the small dodecahedral (H2O)20 cage of the sII clathrate hydrate. All three translational and the two rotational degrees of freedom of H2 are included explicitly, as fully coupled, while the cage is treated as rigid. The 5D potential energy surface (PES) of the H2-cage system is pairwise additive, based on the high-quality ab initio 5D (rigid monomer) PES for the H2-H2O complex. The bound state calculations involve no dynamical approximations and provide an accurate picture of the quantum 5D translation-rotation dynamics of H2 inside the cage. The energy levels are assigned with translational (Cartesian) and rotational quantum numbers, based on calculated root-mean-square displacements and probability density plots. The translational modes exhibit negative anharmonicity. It is found that j is a good rotational quantum number, while the threefold degeneracy of the j = 1 level is lifted completely. There is considerable translation-rotation coupling, particularly for excited translational states.

ArnHF van der Waals clusters revisited: II. Energetics and HF vibrational frequency shifts from diffusion Monte Carlo calculations on additive and nonadditive potential-energy surfaces for n=1-12

The ground-state energies and HF vibrational frequency shifts of Ar(n)HF clusters have been calculated on the nonadditive potential-energy surfaces (PESs) for n=2-7 and on the pairwise-additive PESs for the clusters with n=1-12, using the diffusion Monte Carlo (DMC) method. For n>3, the calculations have been performed for the lowest-energy isomer and several higher-lying isomers which are the closest in energy. They provide information about the isomer dependence of the HF redshift, and enable direct comparison with the experimental data recently obtained in helium nanodroplets. The agreement between theory and experiment is excellent, in particular, for the nonadditive DMC redshifts. The relative, incremental redshifts are reproduced accurately even at the lower level of theory, i.e., the DMC and quantum five-dimensional (rigid Ar(n)) calculations on the pairwise-additive PESs. The nonadditive interactions make a significant contribution to the frequency shift, on the order of 10%-12%, and have to be included in the PESs in order for the theory to yield accurate magnitude of the HF redshift. The energy gaps between the DMC ground states of the cluster isomers are very different from the energy separation of their respective minima on the PES, due to the considerable variations in the intermolecular zero-point energy of different Ar(n)HF isomers.

HF in clusters of molecular hydrogen. I. Size evolution of quantum solvation by parahydrogen molecules

We present a theoretical study of the quantum solvation of the HF molecule by a small number of parahydrogen molecules, having n = 1-13 solvent particles. The minimum-energy cluster structures determined for n = 1-12 have all of the H(2) molecules in the first solvent shell. The first solvent shell closes at n = 12 and its geometry is icosahedral, with the HF molecule at the center. The quantum-mechanical ground-state properties of the clusters are calculated exactly using the diffusion Monte Carlo method. The zero-point energy of (p-H(2))(n)HF clusters is unusually large, amounting to 86% of the potential well depth for n > 7. The radial probability distribution functions (PDFs) confirm that the first solvent shell is complete for n = 12, and that the 13th p-H(2) molecule begins to fill the second solvent shell. The p-H(2) molecules execute large-amplitude motions and are highly mobile, making the solvent cage exceptionally fluxional. The anisotropy of the solvent, very pronounced for small clusters, decreases rapidly with increasing n, so that for n approximately 8-9 the solvent environment is practically isotropic. The analysis of the pair angular PDF reveals that for a given n, the parahydrogen solvent density around the HF is modulated in a pattern which clearly reflects the lowest-energy cluster configuration. The rigidity of the solvent clusters displays an interesting size dependence, increasing from n = 6 to 9, becoming floppier for n = 10, and increasing again up to n = 12, as the solvent shell is filled. The rigidity of the solvent cage appears to reach its maximum for n = 12, the point at which the first solvent shell is closed.

Stereographic projections path integral for inertia ellipsoids: Applications to Arn–HF clusters

The DeWitt formula for inertia ellipsoids mapped by stereographic projection coordinates is developed. We discover that by remapping the quaternion parameter space with stereographic projections, considerable simplification of the differential geometry for the inertia ellipsoid with spherical symmetry takes place. The metric tensor is diagonal and contains only one independent element in that case. We find no difficulties testing and implementing the DeWitt formula for the inertia ellipsoids of asymmetric tops mapped by stereographic projections. The path integral algorithm for the treatment of Rm x S2 manifolds based on a mixture of Cartesian and stereographic projection coordinates is tested for small Arn-HF clusters in the n = 2 to n = 5 range. In particular, we determine the quantum effects of the red shift and the isomerization patterns at finite temperatures. Our findings are consistent with previously reported computations and experimental data for small Arn-HF clusters.