Local Substructures of Water Studied by Transient Hole-Burning Spectroscopy in the Infrared: Dynamics and Temperature Dependence

Watching Hydrogen Bonds Break: A Transient Absorption Study of Water

Ultrafast infrared transient absorption measurements of the complete hydroxyl OD stretching mode spectrum of HOD in water, from 100 fs to tens of picoseconds, observe hydrogen bond breaking and monitor the equilibration of the hydrogen bond network in water. In addition, the vibrational lifetime, the time constant for hydrogen bond breaking, and the rate of orientational relaxation are determined. The reactant and photoproduct spectra of the hydrogen bond breaking process are identified by decomposing the transient spectra into two components, the initial spectrum associated with vibrational excited states (reactants) and the long-time spectrum associated with broken hydrogen bonds (photoproducts). By properly taking into account the perturbation of the reactant spectrum decay by the growth of the photoproduct spectrum, it is found that the vibrational relaxation (1.45 ps) and orientational relaxation (1.53 ps) are wavelength independent and, therefore, independent of the degree of hydrogen bonding. Energy deposited into water by vibrational relaxation does not immediately break a hydrogen bond by predissociation nor produce a thermally equilibrated hydrogen bond distribution at an elevated temperature. Following deposition of energy by vibrational relaxation, the hydrogen bond breaking time is 800 fs, and there is a transient period of several picoseconds during which the hydrogen bond distribution is not in thermal equilibrium.

Collective Solvent Coordinates for the Infrared Spectrum of HOD in D2O Based on an ab Initio Electrostatic Map

An ab initio MP2 vibrational Hamiltonian of HOD in an external electrostatic potential parametrized by the electric field and its gradient-tensor is constructed. By combining it with the fluctuating electric field induced by the D(2)O solvent obtained from molecular dynamics simulations, we calculate the infrared absorption of the O-H stretch. The resulting solvent shift and infrared line shape for three force fields (TIP4P, SPC/E, and SW) are in good agreement with the experiment. A collective coordinate response for the solvent effect is constructed by identifying the main electrostatic field and gradient components contributing to the line shape. This allows a realistic stochastic Liouville equation simulation of the line shapes which is not restricted to Gaussian frequency fluctuations.

Spectral diffusion in a fluctuating charge model of water

We apply the combined electronic structure/molecular dynamics approach of Corcelli, Lawrence, and Skinner [J. Chem. Phys. 120, 8107 (2004)] to the fluctuating charge (SPC-FQ) model of liquid water developed by Rick, Stuart, and Berne [J. Chem. Phys. 101, 6141 (1994)]. For HOD in H(2)O the time scale for the long-time decay of the OD stretch frequency time-correlation function, which corresponds to the time scale for hydrogen-bond rearrangement in the liquid, is about 1.5 ps. This result is significantly longer than the 0.9 ps decay previously calculated for the nonpolarizable SPC/E water model. Our results for the SPC-FQ model are in better agreement with recent vibrational echo experiments.

Ultrafast vibrational spectroscopy of water and aqueous N-methylacetamide: Comparison of different electronic structure/molecular dynamics approaches

Kwac and Cho [J. Chem. Phys. 119, 2247 (2003)] have recently developed a combined electronic structure/molecular dynamics approach to vibrational spectroscopy in liquids. The method involves fitting ab initio vibrational frequencies for a solute in a cluster of solvent molecules to a linear combination of the electrostatic potentials on the solute atoms due to the charges on the solvent molecules. These authors applied their method to the N-methylacetamide-D/D(2)O system. We (S. A. Corcelli, C. P. Lawrence, and J. L. Skinner, [J. Chem. Phys. 120, 8107 (2004)]) have recently explored a closely related method, where instead of the electrostatic potential, the solute vibrational frequencies are fit to the components of the electric fields on the solute atoms due to the solvent molecules. We applied our method to the HOD/D(2)O and HOD/H(2)O systems. In order to make a direct comparison of these two approaches, in this paper we apply their method to the water system, and our method to the N-methylacetamide system. For the water system we find that the electric field method is superior to the potential approach, as judged by comparison with experiments for the absorption line shape. For the N-methylacetamide system the two methods are comparable.

The vibrational Stokes shift of water (HOD in D2O)

The vibrational Stokes shift of the OH stretching transition nu(OH) of water is the shift between the ground-state absorption and the excited-state (v=1) emission. A recent measurement on HOD in D(2)O solvent [S. Woutersen and H. J. Bakker, Phys. Rev. Lett. 83, 2077 (1999)] of a 70 cm(-1) redshift, and a subsequent calculation of a 57 cm(-1) redshift using equilibrium molecular dynamics simulations [C. P. Lawrence and J. L. Skinner, J. Chem. Phys. 117, 8847 (2002)] were in good agreement. We now report extensive measurements of the vibrational Stokes shift in HOD/D(2)O using an ultrafast IR pump, Raman probe method. The vibrational Stokes shift is seen to depend on the pump pulse frequency and on time delay; by varying these parameters it can be made to range from 112 to -32 cm(-1) (negative values indicate a blueshift in the excited state). The equilibrium vibrational Stokes shift is actually a negative rather than a positive quantity. Possible reasons for the disagreement between experiment and theory are briefly discussed.

Combined electronic structure/molecular dynamics approach for ultrafast infrared spectroscopy of dilute HOD in liquid H2O and D2O

We present a new approach that combines electronic structure methods and molecular dynamics simulations to investigate the infrared spectroscopy of condensed phase systems. This approach is applied to the OH stretch band of dilute HOD in liquid D2O and the OD stretch band of dilute HOD in liquid H2O for two commonly employed models of water, TIP4P and SPC/E. Ab initio OH and OD anharmonic transition frequencies are calculated for 100 HOD x (D2O)n and HOD x(H2O)n (n = 4-9) clusters randomly selected from liquid water simulations. A linear empirical relationship between the ab initio frequencies and the component of the electric field from the solvent along the bond of interest is developed. This relationship is used in a molecular dynamics simulation to compute frequency fluctuation time-correlation functions and infrared absorption line shapes. The normalized frequency fluctuation time-correlation functions are in good agreement with the results of previous theoretical approaches. Their long-time decay times are 0.5 ps for the TIP4P model and 0.9 ps for the SPC/E model, both of which appear to be somewhat too fast compared to recent experiments. The calculated line shapes are in good agreement with experiment, and improve upon the results of previous theoretical approaches. The methods presented are simple, and transferable to more complicated systems.

Higher-order optical correlation spectroscopy in liquids

Linear optical spectroscopies have long been used to study the behavior of liquids. Laser technology has progressed to the point that it has become possible to perform nonlinear optical experiments that probe higher-order correlation functions in liquids, opening a new window into our understanding of the microscopic details of solution-phase processes. Here we review advances that have been made in recent years in employing higher-order electronic and vibrational spectroscopies to study liquid-state dynamics and structure.