Electronic and Nuclear Contributions to Vibrational Stark Shifts of Hydroxyl Stretching Frequencies of Water Clusters

In spite of the importance of vibrational Stark effect (VSE) and many attempts, origin of VSE is still unclear in molecular level. Here, we studied on origin of VSE of hydroxyl stretching vibration in small water clusters (monomer, dimer, and tetramer) assuming that VSE can be separated by nuclear and electronic contribution. We calculated total Stark tuning rate (Delta mu(tot)) and its nuclear contribution (Delta mu(geom)) using the ab initio method, then the electronic contribution (Delta mu(elec)) was simply obtained by the difference, Delta mu(tot) - Delta mu(geom). In all cases, the nuclear contribution has dominant contribution to VSE. The hydroxyl stretching mode with neighboring hydrogen acceptor showed larger Delta mu(geom) than that of dangling bonds. Furthermore, the calculated Delta mu(geom) became larger in larger cluster due to the hydrogen bond network. The comparison between Stark tuning rates including and excluding anharmonicity supports the importance of potential anharmonicity in VSE, as previously reported. Interestingly, a good linear relationship is observed between the hydroxyl stretch frequency (nu(geom)) and hydroxyl bond length and also between the Stark tuning rate (Delta mu(geom),) and the change of hydroxyl bond length. Similarly, the electronic contribution of calculated frequencies and Stark tuning rate (Delta mu(elec)) showed a good linear relationship with atomic charge derived by electronic perturbation (Delta q(elec)) and change of that (Delta(Delta q(elec))), respectively.