We present the computational methodology that allows rigorous and efficient nine-dimensional (9D) quantum calculations of the intermolecular vibrational states of noncovalently bound trimers of diatomic molecules, with the monomers treated as rigid. The full 9D vibrational Hamiltonian of the trimer is partitioned into a 3D “frame” (or stretching) Hamiltonian and a 6D “bend” Hamiltonian. These two Hamiltonians are diagonalized separately, and a certain number of their lowest-energy eigenstates is included in the final 9D product contracted basis in which the full 9D intermolecular vibrational Hamiltonian is diagonalized. This methodology is applied to the 9D calculations of the intermolecular vibrational levels of (HF)3, a prototypical hydrogen-bonded trimer, on the rigid-monomer version of an ab initio calculated potential energy surface (PES). They are the first to include fully the stretch-bend coupling present in the trimer. The frequencies of all bending fundamentals considered from the present 9D calculations are about 10% lower than those from the earlier quantum 6D calculations that considered only the bending modes of the HF trimer. This means that the stretch-bend coupling is strong, and it is imperative to include it in any accurate treatment of the (HF)3 vibrations aiming to assess the accuracy of the PES employed. Moreover, the 9D results are in better agreement with the limited available spectroscopic data that those from the 6D calculations. In addition, the 9D results show sensitivity to the value of the HF bond length, equilibrium or vibrationally averaged, used in the calculations. The implication is that full-dimensional 12D quantum calculations will be required to obtain definitive vibrational excitation energies for a given PES. Our study also demonstrates that the nonadditive three-body interactions are very significant in (HF)3 and have to be included in order to obtain accurate intermolecular vibrational energy levels of the trimer.
We present the computational methodology, which for the first time allows rigorous twelve-dimensional (12D) quantum calculations of the coupled intramolecular and intermolecular vibrational states of hydrogen-bonded trimers of flexible diatomic molecules. Its starting point is the approach that we introduced recently for fully coupled 9D quantum calculations of the intermolecular vibrational states of noncovalently bound trimers comprised of diatomics treated as rigid. In this paper, it is extended to include the intramolecular stretching coordinates of the three diatomic monomers. The cornerstone of our 12D methodology is the partitioning of the full vibrational Hamiltonian of the trimer into two reduced-dimension Hamiltonians, one in 9D for the intermolecular degrees of freedom (DOFs) and another in 3D for the intramolecular vibrations of the trimer, and a remainder term. These two Hamiltonians are diagonalized separately, and a fraction of their respective 9D and 3D eigenstates is included in the 12D product contracted basis for both the intra- and intermolecular DOFs, in which the matrix of the full 12D vibrational Hamiltonian of the trimer is diagonalized. This methodology is implemented in the 12D quantum calculations of the coupled intra- and intermolecular vibrational states of the hydrogen-bonded HF trimer on an ab initio calculated potential energy surface (PES). The calculations encompass the one- and two-quanta intramolecular HF-stretch excited vibrational states of the trimer and low-energy intermolecular vibrational states in the intramolecular vibrational manifolds of interest. They reveal several interesting manifestations of significant coupling between the intra- and intermolecular vibrational modes of (HF)3. The 12D calculations also show that the frequencies of the v = 1, 2 HF stretching states of the HF trimer are strongly redshifted in comparison to those of the isolated HF monomer. Moreover, the magnitudes of these trimer redshifts are much larger than that of the redshift for the stretching fundamental of the donor-HF moiety in (HF)2, most likely due to the cooperative hydrogen bonding in (HF)3. The agreement between the 12D results and the limited spectroscopic data for the HF trimer, while satisfactory, leaves room for improvement and points to the need for a more accurate PES.
more » « less- NSF-PAR ID:
- 10510672
- Publisher / Repository:
- The American Institute of Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 158
- Issue:
- 23
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
more » « less
We present fully coupled, full-dimensional quantum calculations of the inter- and intra-molecular vibrational states of HCl trimer, a paradigmatic hydrogen-bonded molecular trimer. They are performed utilizing the recently developed methodology for the rigorous 12D quantum treatment of the vibrations of the noncovalently bound trimers of flexible diatomic molecules [Felker and Bačić, J. Chem. Phys. 158, 234109 (2023)], which was previously applied to the HF trimer by us. In this work, the many-body 12D potential energy surface (PES) of (HCl)3 [Mancini and Bowman, J. Phys. Chem. A 118, 7367 (2014)] is employed. The calculations extend to the intramolecular HCl-stretch excited vibrational states of the trimer with one- and two-quanta, together with the low-energy intermolecular vibrational states in the two excited v = 1 intramolecular vibrational manifolds. They reveal significant coupling between the intra- and inter-molecular vibrational modes. The 12D calculations also show that the frequencies of the v = 1 HCl stretching states of the HCl trimer are significantly redshifted relative to those of the isolated HCl monomer. Detailed comparison is made between the results of the 12D calculations on the two-body PES, obtained by removing the three-body term from the original 2 + 3-body PES, and those computed on the 2 + 3-body PES. It demonstrates that the three-body interactions have a strong effect on the trimer binding energy as well as on its intra- and inter-molecular vibrational energy levels. Comparison with the available spectroscopic data shows that good agreement with the experiment is achieved only if the three-body interactions are included. Some low-energy vibrational states localized in a secondary minimum of the PES are characterized as well.
-
more » « less
In this work the H2O–HCN complex is quantitatively characterized in two ways. First, we report a new rigid-monomer 5D intermolecular potential energy surface (PES) for this complex, calculated using the symmetry-adapted perturbation theory based on density functional theory method. The PES is based on 2833 ab initio points computed employing the aug-cc-pVQZ basis set, utilizing the autoPES code, which provides a site-site analytical fit with the long-range region given by perturbation theory. Next, we present the results of the quantum 5D calculations of the fully coupled intermolecular rovibrational states of the H2O–HCN complex for the total angular momentum J values of 0, 1, and 2, performed on the new PES. These calculations rely on the quantum bound-state methodology developed by us recently and applied to a variety of noncovalently bound binary molecular complexes. The vibrationally averaged ground-state geometry of H2O–HCN determined from the quantum 5D calculations agrees very well with that from the microwave spectroscopic measurements. In addition, the computed ground-state rotational transition frequencies, as well as the B and C rotational constants calculated for the ground state of the complex, are in excellent agreement with the experimental values. The assignment of the calculated intermolecular vibrational states of the H2O–HCN complex is surprisingly challenging. It turns out that only the excitations of the intermolecular stretch mode can be assigned with confidence. The coupling among the angular degrees of freedom (DOFs) of the complex is unusually strong, and as a result most of the excited intermolecular states are unassigned. On the other hand, the coupling of the radial, intermolecular stretch mode and the angular DOFs is weak, allowing straightforward assignment of the excitation of the former.
-
The methodological advances made in recent years have significantly extended the range and dimensionality of noncovalently bound, hydrogen-bonded and van der Waals, molecular complexes for which full-dimensional and fully coupled quantum calculations of their rovibrational states are feasible. They exploit the unexpected implication that the weak coupling between the inter- and intramolecular rovibrational degrees of freedom (DOFs) of the complexes has for the ease of computing the high-energy eigenstates of the latter. This is done very effectively by using contracted eigenstate bases to cover both intra- and intermolecular DOFs. As a result, it is now possible to calculate rigorously all intramolecular rovibrational fundamentals, together with the low-lying intermolecular rovibrational states, of complexes involving two small molecules beyond diatomics, binary polyatomic molecule-large (rigid) molecule complexes, and endohedral complexes of light polyatomic molecules confined inside (rigid) fullerene cages. In this Perspective these advances are reviewed in considerable depth. The progress made thanks to them is illustrated by a number of representative applications. Whenever possible, direct comparison is made with the available infrared, far-infrared, and microwave spectroscopic data.more » « less
-
NMR relaxation rates for protons in liquid water and neat acetonitrile were computed based on ab initio molecular dynamics (aiMD) with forces from Kohn–Sham (KS) theory as well as force-field (FF) based classical dynamics. Intra- and intermolecular dipole–dipole contributions were separated, and nearly quantitative agreement with experiment was obtained for water. Spin-rotation (SR) contributions to the intramolecular relaxation rate in acetonitrile were computed using nuclear SR coupling tensors obtained from KS theory. Their inclusion improved the total computed intramolecular rate to within a factor of two of experiment. Insufficient sampling of rare short-time collision events between neighboring acetonitrile molecules in the simulations is hypothesized as a major source of error in the intermolecular contributions.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0156976
