Bridging theory with experiment: a benchmark study of thermally averaged structural and effective spectroscopic parameters of the water molecule

J Phys Chem A. 2009 Oct 29;113(43):11665-78. doi: 10.1021/jp902690k.

Abstract

Extending our previous study on the equilibrium structures of the major isotopologues of the water molecule (Csaszar et al. J. Chem. Phys. 2005, 122, 214305), temperature-dependent averaged structural parameters (for example, r(g)- and r(a)-type distances, their related root-mean-square amplitudes, and moments corresponding to the probability distribution functions of interatomic distances), effective rotational constants, and low-order vibration-rotation interaction constants have been determined for two major symmetric isotopologues of water, H(2)(16)O and D(2)(16)O. The nuclear motion treatments employed full quantum mechanical variational procedures which utilized the accurate adiabatic semiglobal PESs of water isotopologues named CVRQD (Barletta et al. J. Chem. Phys. 2006, 125, 204307). The temperature-dependent molecular structural parameters are based on expectation value computations and Boltzmann averaging in the temperature range 0-1500 K. The precise computed average internuclear, inverse internuclear, rms amplitude, and anharmonicity parameters could support a future gas electron diffraction (GED) investigation, though water isotopologues are far from being ideal species for GED analyses. Using a clearly defined and general formalism applicable to molecules of any size, we have evaluated vibrationally averaged effective rotational constants as expectation values using inertia tensor formulas in the Eckart frame for vibrational states of H(2)(16)O and D(2)(16)O. While such variationally determined rotational constants do not correspond strictly to constants resulting from fits performed by spectroscopists, the expected good agreement is found for the A and B rotational constants for both isotopologues. Low-order vibration-rotation interaction constants, the so-called alpha- and gamma-constants, have also been determined from the computed rotational constants; the latter were derived for the first time.