Ultrafast anisotropy dynamics of water molecules dissolved in acetonitrile

Infrared pump-probe experiments are performed on isolated H(2)O molecules diluted in acetonitrile in the spectral region of the OH stretching vibration. The large separation between water molecules excludes intermolecular interactions, while acetonitrile as a solvent provides substantial hydrogen bonding. Intramolecular coupling between symmetric and asymmetric modes results in the anisotropy decay to the frequency-dependent values of approximately 0-0.2 with a 0.2 ps time constant. The experimental data are consistent with a theoretical model that includes intramolecular coupling, anharmonicity, and environmental fluctuations. Our results demonstrate that intramolecular processes are essential for the H(2)O stretching mode relaxation and therefore can compete with the intermolecular energy transfer in bulk water.     

Ultrafast excited-state electron transfer at an organic liquid/aqueous interface

Ultrafast excited-state electron transfer has been monitored at the liquid/liquid interface for the first time. Second harmonic generation (SHG) pump/probe measurements monitored the electron transfer (ET) occurring between photoexcited coumarin 314 (C314) acceptor and dimethylaniline (DMA) donor molecules. In the treatment of this problem, translational diffusion of solute molecules can be neglected since the donor DMA is one of the liquid phases of the interface. The dynamics of excited-state C314 at early times are characterized by two components with exponential time constants of 362 +/- 60 fs and 14 +/- 2 ps. The 362 fs decay is attributed to the solvation of the excited-state C314, and the 14 ps to the ET from donor to acceptor. We are able to provide conclusive evidence that the 14 ps component is the ET step by monitoring the formation of the radical DMA cation. The formation time is 16 ps in agreement with the 14 ps decay of C314*. The recombination dynamics of DMA+ plus C314- was determined to be 163 ps from the observation of the DMA+ SHG signal.     

Energy transfer in single hydrogen-bonded water molecules.

We study the structure and dynamics of hydrogen-bonded complexes of H2O/HDO and acetone dissolved in carbon tetrachloride by probing the response of the O-H stretching vibrations with linear mid-infrared spectroscopy and femtosecond mid-infrared pump-probe spectroscopy. We find that the hydrogen bonds in these complexes break and reform with a characteristic time scale of approximately 1 ps. These hydrogen-bond dynamics are observed to play an important role in the equilibration of vibrational energy over the two O-H groups of the H2O molecule. For both H2O and HDO, the O-H stretching vibrational excitation relaxes with a time constant of 6.3+/-0.3 ps, and the molecular reorientation has a time constant of 6+/-1 ps.     

Dynamical properties of hydrogen sulphide motion in its clathrate hydrate from ab initio and classical isobaric-isothermal molecular dynamics

Equilibrium ab initio (AI) and classical constant pressure-constant temperature molecular dynamics (MD) simulations have been performed to investigate the dynamical properties in (type I) hydrogen sulphide hydrate at 150 and 300 K and 1 bar, and also lower temperatures, with particular scrutiny of guest motion. The rattling motions of the guests in the large and small cavities were around 45 and 75-80 cm(-1), respectively, from AIMD, with the corresponding classical MD modes being 10-12 cm(-1) less at 150 K and around 5 cm(-1) lower at 300 K. The rattling motion in the small cavity overlapped somewhat with the translational motion of the host lattice (with modes at circa 85 and 110 cm(-1)), due in part to a smaller cage radius and more frequent occurrences of guest-host hydrogen bonding leading to greater coupling in the motion. The experimentally determined H-S stretch and H-S-H bending frequencies, in the vicinity of 2550-2620 and 1175 cm(-1) [H.R. Richardson et al., J. Chem. Phys.1985, 83, 4387] were reproduced successfully in the AIMD simulations. Consideration of Kubic harmonics for the guest molecules from AIMD revealed that a preferred orientation of the dipole-vector (or C(2)-axis) exists at 150 K vis--vis the [100] cube axis in both the small and large cavities, but is markedly more significant for the small cavity, while there is no preferred orientation at 300 K. In comparison, classical MD did not reveal any preferred orientation at either temperature, or at 75 K (closer to the AIMD simulation at 150 K vis-à-vis that approach's estimated melting point). Probing rotational dynamics of the guests reinforced this temperature effect, revealing more rapid rotational time scales at 300 K with faster decay times of dipole-vector (C(2)) and H-H-vector (C(2)(y)) being similar for each cage, at around 0.25 and 0.2 ps, respectively, versus approximately 0.45 and 0.5 ps (large) and 0.8 ps (small) at 150 K. It was found that the origin of the observed preferred orientations, especially in the small cages, at 150 K via AIMD was attributable to optimization of the dipolar interaction between the guest and outward-pointing water dipoles in the cavity, with guests "flitting" rotationally between various such configurations, forming occasionally hydrogen bonds with the host molecules.     

Solvation dynamics in Ni+ (H2O)n clusters probed with infrared spectroscopy

Infrared photodissociation spectroscopy is reported for mass-selected Ni+ (H2O)n complexes in the O-H stretching region up to cluster sizes of n = 25. These clusters fragment by the loss of one or more intact water molecules, and their excitation spectra show distinct bands in the region of the symmetric and asymmetric stretches of water. The first evidence for hydrogen bonding, indicated by a broad band strongly red-shifted from the free OH region, appears at the cluster size of n = 4. At larger cluster sizes, additional red-shifted structure evolves over a broader wavelength range in the hydrogen-bonding region. In the free OH region, the symmetric stretch gradually diminishes in intensity, while the asymmetric stretch develops into a closely spaced doublet near 3700 cm(-1). The data indicate that essentially all of the water molecules are in a hydrogen-bonded network by the size of n = 10. However, there is no evidence for the formation of clathrate structures seen recently via IR spectroscopy of protonated water clusters.     

Comprehensive investigation of vibrational relaxation of non-hydrogen-bonded water molecules in liquids

The vibrational dynamics of isolated water molecules dissolved in the nonpolar organic liquids 1,2-dichloroethane (C(2)H(4)Cl(2)) and d-chloroform (CDCl(3)) have been studied using an IR pump-probe experiment with approximately 2 ps time resolution. Analyzing transient, time, and spectrally resolved data in both the OH bending and the OH stretching region, the anharmonic constants of the bending overtone (v=2) and the bend-stretch combination modes were obtained. Based on this knowledge, the relaxation pathways of single water molecules were disentangled comprehensively, proving that the vibrational energy of H(2)O molecules is relaxing following the scheme OH stretch-->OH bend overtone-->OH bend-->ground state. A lifetime of 4.8+/-0.4 ps is determined for the OH bending mode of H(2)O in 1,2-dichloroethane. For H(2)O in CDCl(3) a numerical analysis based on rate equations suggests a bending overtone lifetime of tau(020)=13+/-5 ps. The work also shows that full 2-dimensional (pump-probe) spectral resolution with access to all vibrational modes of a molecule is required for the comprehensive analysis of vibrational energy relaxation in liquids.     

Solvation dynamics of DCM in a DPPC vesicle entrapped in a sodium silicate derived sol-gel matrix

Solvation dynamics of 4-(dicyanomethylene)-2-methyl-6(p-dimethylaminostyryl) 4H-pyran (DCM) has been studied in a dipalmitoyl-phosphatidylcholine (DPPC) vesicle entrapped in a sodium silicate derived sol-gel glass. Solvation dynamics in DPPC in a sol-gel glass is described by two components of 350 +/- 50 ps (50%) and 2300 +/- 200 ps (50%) with a total dynamic Stokes shift of 1300 cm(-1). The fast component (350 ps) is similar to the fast component in a DPPC vesicle in bulk water (320 +/- 50 ps). This component may be ascribed to the dynamics of the water molecules inside the water pool of the vesicle. However, the slow component (2300 +/- 200 ps) is about 2.5 times slower compared to the slow component of solvation dynamics of DCM in a DPPC vesicle in bulk solvent (900 +/- 100 ps). The anisotropy decay of DCM in a DPPC vesicle both in sol-gel glass and in bulk water exhibits a very fast initial decay with a large residual anisotropy, which does not decay in approximately 10 ns. The time scale of anisotropy decay is very different from that of solvation dynamics.     

Vibrational dynamics of deoxyguanosine 5'-monophosphate following UV excitation

The relaxation dynamics of the DNA nucleotide deoxyguanosine 5'-monophosphate (dGMP) following 266 nm photoexcitation has been studied by transient IR spectroscopy with femtosecond time resolution. The induced dynamics of the amide I (carbonyl) stretch, the asymmetric guanine ring stretch and the phosphate asymmetric stretch are monitored in the region 1000-1800 cm(-1). Excitation and subsequent rapid internal conversion to a "hot" ground state is reflected by depletion of the vibrational ground states of the amide I stretch and guanine ring stretch. However, the vibrational ground state of the phosphate is left unperturbed, indicating the absence of vibrational coupling between the guanine ring system and the phosphate group. The vibrational ground state of the amide I is repopulated in 2.5 ps (±0.2 ps) while it takes 3.7 ps (±0.5 ps) to repopulate the guanine ring vibration. This article discusses two possible relaxation pathways of dGMP, as well as the implications of the weak phosphate dynamics.     

Vibrational and orientational dynamics of water in aqueous hydroxide solutions

We report the vibrational and orientational dynamics of water molecules in isotopically diluted NaOH and NaOD solutions using polarization-resolved femtosecond vibrational spectroscopy and terahertz time-domain dielectric relaxation measurements. We observe a speed-up of the vibrational relaxation of the O-D stretching vibration of HDO molecules outside the first hydration shell of OH(-) from 1.7 ± 0.2 ps for neat water to 1.0 ± 0.2 ps for a solution of 5 M NaOH in HDO:H(2)O. For the O-H vibration of HDO molecules outside the first hydration shell of OD(-), we observe a similar speed-up from 750 ± 50 fs to 600 ± 50 fs for a solution of 6 M NaOD in HDO:D(2)O. The acceleration of the decay is assigned to fluctuations in the energy levels of the HDO molecules due to charge transfer events and charge fluctuations. The reorientation dynamics of water molecules outside the first hydration shell are observed to show the same time constant of 2.5 ± 0.2 ps as in bulk liquid water, indicating that there is no long range effect of the hydroxide ion on the hydrogen-bond structure of liquid water. The terahertz dielectric relaxation experiments show that the transfer of the hydroxide ion through liquid water involves the simultaneous motion of ~7 surrounding water molecules, considerably less than previously reported for the proton.     

Vibrational lifetimes of molecular adsorbates on metal surfaces

We report density functional theory calculations of electron-hole pair induced vibrational lifetimes of diatomic molecules adsorbed on metal surfaces. For CO on Cu(100), Ni(100), Ni(111), Pt(100), and Pt(111), we find that the C-O internal stretch and the bending modes have lifetimes in the 1-6 ps range, and that the CO-surface stretch and the frustrated translational modes relax more slowly, with lifetimes >10 ps for all cases except CO on Ni(111). This strong mode selectivity confirms earlier calculations for CO on Cu(100) and demonstrates that the trends carry over to other metal substrates. In contrast, for NO adsorbed on Pt(111), whereas we still find that the bending mode has the shortest lifetime, about 1.3 ps, we predict the other three modes to have almost equal lifetimes of 8-10 ps. Similarly, for CN adsorbed on Pt(111), we calculate that the internal stretching and molecule-surface stretching modes have approximately equal lifetimes of about 15 ps. Our results are in reasonable agreement with experiment, where available. We discuss some of the underlying factors that may contribute to the observed mode selectivity with adsorbed CO and the altered selectivity with NO and CN.     

Direct dynamics study of ultrafast vibrational energy relaxation in ice Ih

Car-Parrinello molecular dynamics (CPMD) and a previously developed wave packet model are used to study ultrafast relaxation in water clusters. Water clusters of 15 water molecules are used to represent ice Ih. The relaxation is studied by exciting a symmetric or an asymmetric stretch mode of the central water molecule. The CPMD results suggest that relaxation occurs within 100 fs. This is in agreement with experimental work by Woutersen and Bakker and the earlier wave packet calculations. The CPMD results further indicate that the excitation energy is transferred both intramolecularly and intermolecularly on roughly the same time scale. The intramolecular energy transfer occurs predominantly between the symmetric and asymmetric modes while the bend mode is largely left unexcited on the short time scale studied here.     

Reduced coupling of water molecules near the surface of reverse micelles

We report on vibrational dynamics of water near the surface of AOT reverse micelles studied by narrow-band excitation, mid-IR pump-probe spectroscopy. Evidence of OH-stretch frequency splitting into the symmetric and asymmetric modes is clearly observed for the interfacial H(2)O molecules. The polarization memory of interfacial waters is preserved over an exceptionally extended >10 ps timescale which is a factor of 100 longer than in bulk water. These observations point towards negligibly small intermolecular vibrational coupling between the water molecules as well as strongly reduced water rotational mobility within the interfacial water layer.     

Strong slowing down of water reorientation in mixtures of water and tetramethylurea

We use mid-infrared pump-probe spectroscopy to study the ultrafast dynamics of HDO molecules in mixtures of tetramethylurea (TMU) and water. The composition of the studied solutions ranges from pure water to an equimolar mixture of water and TMU. We find that the vibrational relaxation of the OD-stretching vibration of HDO proceeds via an intermediate level in which the molecule is more strongly hydrogen bonded than in the ground state. As the TMU concentration is increased, the lifetime of the excited state and of the intermediate increase from 1.8 to 5.2 ps and from 0.7 to 2.2 ps, respectively. The orientational relaxation data indicate that the solutions contain two types of water molecules: bulk-like molecules that have the same reorientation time constant as in the pure liquid (taurot = 2.5 ps) and molecules that are strongly immobilized (taurot > 10 ps). The immobilized water molecules turn out to be involved in the solvation of the methyl groups of the tetramethylurea molecule. The fraction of immobilized water molecules grows with increasing TMU concentration, reaching a limiting value of 60% at very high concentrations.     

Nature of biological water: a femtosecond study

The quasi-bound biological or structured water molecules in a protein play a key role in many biological processes. The dynamics of the biological water has been studied by femtosecond spectroscopy and large-scale computer simulations. Solvation dynamics of biological water displays an almost bulk-water like ultrafast component (approximately 1 ps) and a surprising slow component at the 100-1000 ps time scale. In this article, we discuss several examples of the ultraslow component, its possible origin and implications in biology. We show that the ultrafast (approximately 1 ps) component arises from an extended hydrogen bond network while the ultraslow component originates from binding of a water molecule to a biological macromolecule.     

Anomalous dielectric relaxation of water molecules at the surface of an aqueous micelle

Dielectric relaxation of aqueous solutions of micelles, proteins, and many complex systems shows an anomalous dispersion at frequencies intermediate between those corresponding to the rotational motion of bulk water and that of the organized assembly or macromolecule. The precise origin of this anomalous dispersion is not well-understood. In this work we employ large scale atomistic molecular dynamics simulations to investigate the dielectric relaxation (DR) of water molecules in an aqueous micellar solution of cesium pentadecafluorooctanoate. The simulations clearly show the presence of a slow component in the moment-moment time correlation function [PhiMW(t)] of water molecules, with a time constant of about 40 ps, in contrast to only 9 ps for bulk water. Interestingly, the orientational time correlation function [Cmu(t)] of individual water molecules at the surface exhibits a component with a time constant of about 19 ps. We show that these two time constants can be related by the well-known micro-macrorelations of statistical mechanics. In addition, the reorientation of surface water molecules exhibits a very slow component that decays with a time constant of about 500 ps. An analysis of hydrogen bond lifetime and of the rotational relaxation in the coordinate frame fixed on the micellar body seems to suggest that the 500 ps component owes its origin to the existence of an extended hydrogen bond network of water molecules at the surface. However, this ultraslow component is not found in the total moment-moment time correlation function of water molecules in the solution. The slow DR of hydration water is found to be well correlated with the slow solvation dynamics of cesium ions at the water-micelle interface.     

Molecular dynamics simulations of MRI-relevant GdIII chelates: direct access to outer-sphere relaxivity

The structure and dynamics of the surrounding water were studied through molecular dynamics (MD) simulations for several GdIII polyaminocarboxylate and polyaminophosphonate complexes in aqueous solution. The radial distribution functions (rdf) show that a few water molecules are bonded to the ligand through hydrogen bonds to hydrophilic groups such as carboxylates and phosphonates. Residence times are of the order of 20-25 ps for the polyaminocarboxylate and 56ps for the polyaminophosphonate chelates. No preferred orientation or bonding of water molecules is observed in the hydrophobic region of the anisotropic macrocyclic complexes. Our rdf allow calculation of the outer-sphere contribution to the nuclear magnetic resonance dispersion (NMRD) profiles using Freed's finite differences method, including electronic relaxation. The results show that the commonly used analytical force-free model is only an empirical relationship. When experimental outer-sphere NMRD profiles are available ([Gd(teta)]- and [Gd(dotp)]5-(teta=N,N',N",N"'-tetracarboxymethyl-1,4,8,11- tetraazacyclotetradecane; dotp = N,N',N",N"'-tetraphosphonatomethyl-1,4,7,10-tetraazacyclododecane) the calculated curves are in good agreement. In the case of [Gd(teta)]-, the comparison with the experimental NMRD profile has led us to predict a very fast electronic relaxation, which has been confirmed by the EPR spectrum.     

Excess electron and lithium atom solvation in water clusters at finite temperature: an ab initio molecular dynamics study of the structural, spectral, and dynamical behavior of (H2O)6- and Li(H2O)6

The roles of hydrogen bonds in the solvation of an excess electron and a lithium atom in water hexamer cluster at 150 K have been studied by means of ab initio molecular dynamics simulations. It is found that the hydrogen bonded structures of (H(2)O)(6)(-) and Li(H(2)O)(6) clusters are very different from each other and they dynamically evolve from one conformer to other along their simulation trajectories. The populations of the single acceptor, double acceptor, and free type water molecules are found to be significantly high unlike that in pure water clusters. Free hydrogens of these type of water molecules primarily capture the unbound electron density in these clusters. It is found that the binding motifs of the free electron evolve with time and the vertical detachment energy of (H(2)O)(6)(-) and vertical ionization energy of Li(H(2)O)(6) also change with time. Assignments of the observed peaks in vibrational power spectra are done, and we found direct correlations between the time-averaged population of water molecules in different hydrogen bonding states and the spectral features. The dynamical aspects of these clusters have also been studied through calculations of time correlations of instantaneous stretch frequencies of OH modes which are obtained from the simulation trajectories through a time series analysis.     

Vibrational spectra of methane clathrate hydrates from molecular dynamics simulation

Molecular dynamics simulations were performed on methane clathrate hydrates at ambient conditions. Thermal expansion results over the temperature range 60-300 K show that the unit cell volume increases with temperature in agreement with experiment. Power spectra were obtained at 273 K from velocity autocorrelation functions for selected atoms, and normal modes were assigned. The spectra were further classified according to individual atom types, allowing the assignment of contributions from methane molecules located in small and large cages within the structure I unit cell. The symmetric C-H stretch of methane in the small cages occurs at a higher frequency than for methane located in the large cages, with a peak separation of 14 cm(-1). Additionally, we determined that the symmetric C-H stretch in methane gas occurs at the same frequency as methane in the large cages. Results of molecular dynamics simulations indicate the use of power spectra obtained from the velocity autocorrelation function is a reliable method to investigate the vibrational behavior of guest molecules in clathrate hydrates.     

Frequency analysis of amide-linked rotaxane mimetics

A vibrational analysis of 2-fold hydrogen bonds between an isophthalic amide donor and different acceptors is presented. These systems can be considered as mimetics for the hydrogen-binding situation of numerous supramolecular compounds such as rotaxanes, catenanes, knotanes, and anion receptors. We calculated pronounced red-shifts up to 65 cm(-1) for the stretching modes of the acceptor carbonyl as well as for the donor NH2 groups, whereas we observe a blue shift for the NH2 bending modes and an additional weak hydrogen bond between the acceptor and the middle C-H group of the donor. The red and blue shifts observed for different modes in various complexes have been correlated with the binding energy of the complexes, independently. In comparison with comparable single hydrogen bonds, we find for the 2-fold hydrogen bonds smaller red shifts for the N-H stretch modes of the donor but larger red shifts for the C=O stretch mode of the acceptor. Furthermore, our results indicate that the pronounced blue shift of the C-H stretch mode is basically caused by the fact that the acceptor is fixed directly above this group due to the 2-fold hydrogen bond.