Electric-dipole-coupled H2O@C60 dimer: Translation-rotation eigenstates from twelve-dimensional quantum calculations

J Chem Phys. 2017 Feb 28;146(8):084303. doi: 10.1063/1.4976526.

Abstract

We report on variational solutions to the twelve-dimensional (12D) Schrödinger equation appertaining to the translation-rotation (TR) eigenstates of H2O@C60 dimer, associated with the quantized "rattling" motions of the two encapsulated H2O molecules. Both H2O and C60 moieties are treated as rigid and the cage-cage geometry is taken to be fixed. We consider the TR eigenstates of H2O@C60 monomers in the dimer to be coupled by the electric dipole-dipole interaction between water moieties and develop expressions for computing the matrix elements of that interaction in a dimer basis composed of products of monomer 6D TR eigenstates reported by us recently [P. M. Felker and Z. Bačić, J. Chem. Phys. 144, 201101 (2016)]. We use these expressions to compute TR Hamiltonian matrices of H2O@C60 dimer for two values of the water dipole moment and for various dimer geometries. 12D TR eigenstates of the dimer are then obtained by filter diagonalization. The results reveal two classes of eigenstates, distinguished by the leading order (first or second) at which dipole-dipole coupling contributes to them. The two types of eigenstates differ in the general magnitude of their dipole-induced energy shifts and in the dependence of those shifts on the value of the water dipole moment and on the distance between the H2O@C60 monomers. The dimer results are also found to be markedly insensitive to any change in the orientations of the C60 cages. Finally, the results lend some support for the interpretation that electric dipole-dipole coupling is at least partially responsible for the apparent reduced-symmetry environment experienced by H2O in the powder samples of H2O@C60 [K. S. K. Goh et al., Phys. Chem. Chem. Phys. 16, 21330 (2014)], but only if the water dipole is taken to have a magnitude close to that of free water. The methodology developed in the paper is transferable directly to the calculation of TR eigenstates of larger H2O@C60 assemblies, that will be required for more extensive modeling of crystalline H2O@C60.