A study of molecular vibrational relaxation mechanism in condensed phase based upon mixed quantum-classical molecular dynamics. I. A test of IBC model for the relaxation of a nonpolar solute in nonpolar solvent at high density

In order to investigate vibrational relaxation mechanism in condensed phase, a series of mixed quantum-classical molecular dynamics calculations have been executed for nonpolar solute in nonpolar solvent and polar solute in polar solvent. In the first paper (Paper I), relaxation mechanism of I2 in Ar, where Lennard-Jones force is predominant in the interaction, is investigated as a function of density and temperature, focusing our attention on the isolated binary collision (IBC) model. The model was originally established for the relaxation in gas phase. A key question, here, is "can we apply the IBC model to the relaxation in the high-density fluid?" Analyzing the trajectory of solvent molecule as well as its interaction with the solute, we found that collisions between them may be defined clearly even in the high-density fluid. Change of the survival probability of the vibrationally first excited state on collision was traced. The change caused by collisions with a particular solvent molecule was also traced together with the interaction between them. Each collision makes a contribution to the relaxation by a stepwise change in the probability. The analysis clearly shows that the relaxation is caused by collisions even in the high-density fluid. The difference between stepwise relaxation and the continuous one found for the total relaxation in the low-density fluid and in the high-density one, respectively, was clarified to come from just the difference in frequency of the collision. The stronger the intensity of the collision is, the greater the relaxation caused by the collision is. Further, the shorter the collision time is, the greater the resultant relaxation is. The discussion is followed by the succeeding paper (Paper II), where we report that molecular mechanism of the relaxation of a polar molecule in supercritical water is significantly different from that assumed in the IBC model despite that the density dependence of the relaxation rate showed a linear correlation with the local density of water around the solute, the linear correlation being apparently in good accordance with the IBC model. The puzzle will be solved in Paper II.     

Theoretical studies of solute vibrational energy relaxation at liquid interfaces

Recent advances in the theoretical understanding of solute vibrational energy relaxation at liquid interfaces and surfaces are described. Non-equilibrium molecular dynamics simulations of the relaxation of an initially excited solute molecule are combined with equilibrium force autocorrelation calculations to gain insight into the factors that influence the vibrational relaxation rate. Diatomic and triatomic nonpolar, polar, and ionic solute molecules adsorbed at the liquid/vapor interface of several liquids as well as at the water/CCl(4) liquid/liquid interface are considered. In general, the vibrational relaxation rate is significantly slower (a factor of 3 to 4) at the liquid/vapor and liquid/liquid interface than in the bulk due to the reduced density, which gives rise to a reduced contribution of the repulsive solvent-solute forces on the vibrational mode. The surface effects on the ionic solutes are much smaller (50% or less slower relaxation relative to the bulk). This is due to the fact that ionic solutes at the interface are able to keep part of their solvation shell to a degree that depends on their size. Thus, a significant portion of the repulsive forces is maintained. A high degree of correlation is found between the peak height of the solvent-solute radial distribution function and the vibrational relaxation rate. The relaxation rate at the liquid/liquid interface strongly depends on the location of the solute across the interface and correlates with the change in the density and polarity profile of the interface.     

Vibrational relaxation of azide ions in liquid-to-supercritical water

The dynamics of vibrational energy relaxation (VER) of the aqueous azide anion was studied over a wide temperature (300 K ≤ T ≤ 663 K) and density (0.6 g cm(-3) ≤ ρ ≤ 1.0 g cm(-3)) range thereby covering the liquid and the supercritical phase of the water solvent. Femtosecond mid-infrared spectroscopy on the ν(3) band associated with the asymmetric stretching vibration of the azide anion was used to monitor the relaxation dynamics in a time-resolved fashion. The variation of the vibrational relaxation rate constant with temperature and density was found to be rather small. Surprisingly, the simple isolated binary collision model is able to fully reproduce the experimentally observed temperature and density dependence of the relaxation rate provided a local density correction around the vibrationally excited solute based on classical molecular dynamics simulations is used. The simulations further suggest that head-on collisions of the solvent with the terminal nitrogen atoms rather than side-on collisions with the central nitrogen atom of the azide govern the vibrational energy relaxation of this system. Finally, the importance of hydrogen bonding for the VER dynamics in this system is briefly discussed.     

Vibrational relaxation of OH and CH fundamentals of polar and nonpolar molecules in the condensed phase

Studies of vibrational energy flow in various polar and nonpolar molecules that follows the ultrafast excitation of the CH and OH stretch fundamentals, modeled using semiclassical methods, are reviewed. Relaxation rates are calculated using Landau-Teller theory and a time-dependent method, both of which consider a quantum mechanical solute molecule coupled to a classical bath of solvent molecules. A wide range of decay rates are observed, ranging from 1 ps for neat methanol to 50 ps for neat bromoform. In order to understand the flow rates, it is argued that an understanding of the subtle mixing between the solute eigenstates is needed and that solute anharmonicities are critical to facilitating condensed phase vibrational relaxation. The solvent-assisted shifts of the solute vibrational energy levels are seen to play a critical role of enhancing or decreasing lifetimes.     

Solvation in supercritical water

Solvation in supercritical water under equilibrium and nonequilibrium conditions is studied via molecular dynamics simulations. The influence of solute charge distributions and solvent density on the solvation structures and dynamics is examined with a diatomic probe solute molecule. It is found that the solvation structure varies dramatically with the solute dipole moment, especially in low-density water, in accord with many previous studies on ion solvation. This electrostrictive effect has important consequences for solvation dynamics. In the case of a nonequilibrium solvent relaxation, if there are sufficiently many water molecules close to the solute at the outset of the relaxation, the solvent response measured as a dynamic Stokes shift is almost completely governed by inertial rotations of these water molecules. By contrast, in the opposite case of a low local solvent density near the solute, not only rotations but also translations of water molecules play an important role in solvent relaxation dynamics. The applicability of a linear response is found to be significantly restricted at low water densities.     

On the mechanism of hydrophobic association of nanoscopic solutes

The hydration behavior of two planar nanoscopic hydrophobic solutes in liquid water at normal temperature and pressure is investigated by calculating the potential of mean force between them at constant pressure as a function of the solute-solvent interaction potential. The importance of the effect of weak attractive interactions between the solute atoms and the solvent on the hydration behavior is clearly demonstrated. We focus on the underlying mechanism behind the contrasting results obtained in various recent experimental and computational studies on water near hydrophobic solutes. The length scale where crossover from a solvent separated state to the contact pair state occurs is shown to depend on the solute sizes as well as on details of the solute-solvent interaction. We find the mechanism for attractive mean forces between the plates is very different depending on the nature of the solute-solvent interaction which has implications for the mechanism of the hydrophobic effect for biomolecules.     

Solute rotational dynamics at the water liquid/vapor interface

The rotational dynamics of a number of diatomic molecules adsorbed at different locations at the interface between water and its own vapors are studied using classical molecular dynamics computer simulations. Both equilibrium orientational and energy correlations and nonequilibrium orientational and energy relaxation correlations are calculated. By varying the dipole moment of the molecule and its location, and by comparing the results with those in bulk water, the effects of dielectric and mechanical frictions on reorientation dynamics and on rotational energy relaxation can be studied. It is shown that for nonpolar and weekly polar solutes, the equilibrium orientational relaxation is much slower in the bulk than at the interface. As the solute becomes more polar, the rotation slows down and the surface and bulk dynamics become similar. The energy relaxation (both equilibrium and nonequilibrium) has the opposite trend with the solute dipole (larger dipoles relax faster), but here again the bulk and surface results converge as the solute dipole is increased. It is shown that these behaviors correlate with the peak value of the solvent-solute radial distribution function, which demonstrates the importance of the first hydration shell structure in determining the rotational dynamics and dependence of these dynamics on the solute dipole and location.     

Vibrational energy relaxation of a diatomic molecule in a room-temperature ionic liquid

Vibrational energy relaxation (VER) dynamics of a diatomic solute in ionic liquid 1-ethyl-3-methylimidazolium hexafluorophosphate (EMI(+)PF(6) (-)) are studied via equilibrium and nonequilibrium molecular dynamics simulations. The time scale for VER is found to decrease markedly with the increasing solute dipole moment, consonant with many previous studies in polar solvents. A detailed analysis of nonequilibrium results shows that for a dipolar solute, dissipation of an excess solute vibrational energy occurs almost exclusively via the Lennard-Jones interactions between the solute and solvent, while an oscillatory energy exchange between the two is mainly controlled by their electrostatic interactions. Regardless of the anharmonicity of the solute vibrational potential, VER becomes accelerated as the initial vibrational energy increases. This is attributed primarily to the enhancement in variations of the solvent force on the solute bond, induced by large-amplitude solute vibrations. One interesting finding is that if a time variable scaled with the initial excitation energy is employed, dissipation dynamics of the excess vibrational energy of the dipolar solute tend to show a universal behavior irrespective of its initial vibrational state. Comparison with water and acetonitrile shows that overall characteristics of VER in EMI(+)PF(6) (-) are similar to those in acetonitrile, while relaxation in water is much faster than the two. It is also found that the Landau-Teller theory predictions for VER time scale obtained via equilibrium simulations of the solvent force autocorrelation function are in reasonable agreement with the nonequilibrium results.     

Molecular dynamics study on the solvent dependent heme cooling following ligand photolysis in carbonmonoxy myoglobin

The time scale and mechanism of vibrational energy relaxation of the heme moiety in myoglobin was studied using molecular dynamics simulation. Five different solvent models, including normal water, heavy water, normal glycerol, deuterated glycerol and a nonpolar solvent, and two forms of the heme, one native and one lacking acidic side chains, were studied. Structural alteration of the protein was observed in native myoglobin glycerol solution and native myoglobin water solution. The single-exponential decay of the excess kinetic energy of the heme following ligand photolysis was observed in all systems studied. The relaxation rate depends on the solvent used. However, this dependence cannot be explained using bulk transport properties of the solvent including macroscopic thermal diffusion. The rate and mechanism of heme cooling depends upon the detailed microscopic interaction between the heme and solvent. Three intermolecular energy transfer mechanisms were considered: (i) energy transfer mediated by hydrogen bonds, (ii) direct vibration-vibration energy transfer via resonant interaction, and (iii) energy transfer via vibration-translation or vibration-rotation interaction, or in other words, thermal collision. The hydrogen bond interaction and vibration-vibration interaction between the heme and solvent molecules dominates the energy transfer in native myoglobin aqueous solution and native myoglobin glycerol solutions. For modified myoglobin, the vibration-vibration interaction is also effective in glycerol solution, different from aqueous solution. Thermal collisions form the dominant energy transfer pathway for modified myoglobin in water solution, and for both native myoglobin and modified myoglobin in a nonpolar environment. For native myoglobin in a nonpolar solvent solution, hydrogen bonds between heme isopropionate side chains and nearby protein residues, absent in the modified myoglobin nonpolar solvent solution, are key interactions influencing the relaxation pathways.     

Why density augmentation occurs in dilute supercritical solutions

The conventional explanation of the density augmentation in a supercritical solvent, observed spectroscopically when a small amount of a solute was added, involved the clustering of the solvent about individual solute molecules. Here it is suggested that the augmentation is not caused by the solute, but rather it is due to the preexisting near critical fluctuations in the pure solvent and the preference of the solute for the high density regions of the supercritical solvent. It is also shown that the local composition of the solute molecules about a solute molecules is enhanced compared to its bulk composition.     

Vibrational relaxation studies in methyl acetate: role of microenvironment and hydrodynamic forces

The Raman spectra were recorded for the C=O stretching vibration of methyl acetate as solutions in various polar and non-polar solvents. The isotropic component was obtained and the vibrational relaxation rates were calculated. On the basis of dependence of isotropic Raman bandwidth on the hydrodynamic properties of the solute-solvent systems, the information was obtained on the microenvironment prevailing in the neighborhood of the C=O stretching mode. Significant correlation is observed between vibrational relaxation rate and solvent parameters namely viscosity, density, refractive index and molecular radius. Microviscosity, involving the size of solute and solvent molecules, is found to be crucial in determining the bandwidth, hence the relaxation rate. The microenvironment appears to play an important role in the vibrational relaxation process.     

Novel quantum mechanical/molecular mechanical method combined with the theory of energy representation: free energy calculation for the Beckmann rearrangement promoted by proton transfers in the supercritical water

The Beckmann rearrangement of acetone oxime promoted by proton transfers in the supercritical water has been investigated by means of the hybrid quantum mechanical/molecular mechanical approach combined with the theory of energy representation (QM/MM-ER) recently developed. The transition state (TS) structures have been explored by ab initio calculations for the reaction of hydrated acetone oxime on the assumption that the reaction is catalyzed by proton transfers along the hydrogen bonds connecting the solute and the solvent water molecules. Up to two water molecules have been considered as reactants that take part in the proton transfers. As a result of the density functional theory calculations with B3LYP functional and aug-cc-pVDZ basis set, it has been found that participation of two water molecules in the reaction reduces the activation free energy by -12.3 kcal/mol. Furthermore, the QM/MM-ER simulations have revealed that the TS is more stabilized than the reactant state in the supercritical water by 2.7 kcal/mol when two water molecules are involved in the reaction. Solvation free energies of the reactant and the TS have been decomposed into terms due to the electronic polarization of the solute, electron density fluctuation, and others to elucidate the origin of the stabilization of the TS as compared with the reactant. It has been revealed that the promotion of the chemical reaction due to the hydration mainly originates from the interaction between the nonpolarized solute and the solvent water molecules at the supercritical state.     

High-density hydration layer of lysozymes: molecular dynamics decomposition of solution scattering data

A characterization of the physical properties of protein hydration water is critical for understanding protein structure and function. Recent small-angle X-ray and neutron scattering data indicate that the density of water on the surface of lysozyme is significantly higher than in bulk water. Here, we provide an interpretation of the scattering results using a molecular dynamics simulation, which allows us to make quantitative predictions about density variations in the first hydration shell. The perturbation relative to bulk water involves statistically significant changes in the average water structure in the first hydration layer. The water density in the first hydration shell is increased by 5% with respect to the bulk. In regions of higher water density, the water dipoles align more parallel to each other and the number of hydrogen bonds per water molecule is higher. Increased water density is found for water molecules interacting with hydrogen and carbon atoms in the backbone or with nonpolar or negatively charged side-chain groups.     

Hydrophobic hydration in an orientational lattice model

To shed light on the microscopic mechanism of hydrophobic hydration, we study a simplified lattice model for water solutions in which the orientational nature of hydrogen bonding as well as the degeneracy related to proton distribution are taken into account. Miscibility properties of the model are looked at for both polar (hydrogen bonding) and nonpolar (non-hydrogen bonding) solutes. A quasichemical solution for the pure system is reviewed and extended to include the different kinds of solute. A Monte Carlo study of our model yields a novel feature for the local structure of the hydration layer: energy correlation relaxation times for solvation water are larger than the corresponding relaxation times for bulk water. This result suggests the presence of ordering of water particles in the first hydration shell. A nonassociating model solvent, represented by a lattice gas, presents opposite behavior, indicating that this effect is a result of the directionality of the interaction. In presence of polar solutes, we find an ordered mixed pseudophase at low temperatures, indicating the possibility of closed loops of immiscibility.     

The role of local structure on formic acid decomposition in supercritical water: A hybrid quantum mechanics/Monte Carlo study

We explored water-assisted decompositions of formic acid in supercritical water in terms of local structure near reactant. A hybrid quantum mechanics/molecular mechanics (QM/MM) simulation used in this paper includes QM part as first solvation shell members around the reactant. A present QM/MM approach can simulate supercritical water solution with a reasonable computational load while keeping the simulation preciseness because a density functional theory of B3LYP/6-31+G(d) level was iterated at every 1000 Monte Carlo solute moves. The formic acid converts mainly decarboxylation by water-assisted mechanism, and the coordinated water molecules play an important role for understanding supercritical water density dependence of the reaction. We analyzed a contour map based on the solute–solvent interaction energy along with the reaction pathway. Coordinated water molecule restricted the dehydration pathway by means of hydrogen bond with formic acid, however, the coordinated water promotes the decarboxylation pathway by means of stabilization of the transition state structure with one catalytic water molecule. The contour map of the pair interaction energy along the reaction path elucidates the role of local structure on reactions in supercritical water.     

Vibrational relaxation in hydrogen bonded systems: I. Theory

A theory is proposed which expresses the influence of intermolecular interaction in H-bonded solutions on the solute IR band shape in terms of parameters characterising the modulation and strength of this interaction. It is found that appreciable vibrational relaxation can occur by interaction between oscillators with different frequencies, provided that the energy gap is closed by the hydrogen bond spectral density. Isotopic substitution of the solvent provides a means for experimental verification, which is the subject of a subsequent paper.     

Vibrational energy relaxation of the bend fundamental of dilute water in liquid chloroform and d-chloroform

The population lifetimes of the bend fundamental of dilute water in liquid chloroform (8.5 ps) and d-chloroform (28.5 ps) display an interesting solvent isotope effect. As the lowest excited vibrational state of the molecule, the water bend fundamental relaxes directly to the ground state with about 1600 cm-1 of energy released to the other degrees of freedom. The strong solvent isotope effect along with the large energy gap indicates the participation of solvent vibrational modes in this vibrational energy relaxation process. We calculate the vibrational energy relaxation rates of the water bend in chloroform and d-chloroform using the Landau-Teller formula with a new potential model developed and parametrized self-consistently to describe the chloroform-water interaction. The computed values are in reasonable agreement with the experimental results, and the trend for the isotope effect is correct. It is found that energy transfer to the solvent vibrations does indeed play an important role. Nevertheless, no single dominant solvent accepting mode can be identified; the relaxation appears to involve both the bend and the C-Cl stretches, and frequency changes of all of these modes upon deuteration contribute to the observed solvent isotope effect.     

An application of the novel quantum mechanical/molecular mechanical method combined with the theory of energy representation: An ionic dissociation of a water molecule in the supercritical water

A novel quantum chemical approach recently developed has been applied to an ionic dissociation of a water molecule (2H(2)O-->H(3)O(+)+OH(-)) in ambient and supercritical water. The method is based on the quantum mechanical/molecular mechanical (QM/MM) simulations combined with the theory of energy representation (QM/MM-ER), where the energy distribution function of MM solvent molecules around a QM solute serves as a fundamental variable to determine the hydration free energy of the solute according to the rigorous framework of the theory of energy representation. The density dependence of the dissociation free energy in the supercritical water has been investigated for the density range from 0.1 to 0.6 g/cm(3) with the temperature fixed at a constant. It has been found that the product ionic species significantly stabilizes in the high density region as compared with the low density. Consequently, the dissociation free energy decreases monotonically as the density increases. The decomposition of the hydration free energy has revealed that the entropic term (-TDeltaS) strongly depends on the density of the solution and dominates the behavior of the dissociation free energy with respect to the variation of the density. The increase in the entropic term in the low density region can be attributed to the decrease in the translational degrees of freedom brought about by the aggregation of solvent water molecules around the ionic solute.