Vibrational spectroscopic determination of local solvent electric field, solute-solvent electrostatic interaction energy, and their fluctuation amplitudes

IR probes have been extensively used to monitor local electrostatic and solvation dynamics. Particularly, their vibrational frequencies are highly sensitive to local solvent electric field around an IR probe. Here, we show that the experimentally measured vibrational frequency shifts can be inversely used to determine local electric potential distribution and solute-solvent electrostatic interaction energy. In addition, the upper limits of their fluctuation amplitudes are estimated by using the vibrational bandwidths. Applying this method to fully deuterated N-methylacetamide (NMA) in D(2)O and examining the solvatochromic effects on the amide I' and II' mode frequencies, we found that the solvent electric potential difference between O(═C) and D(-N) atoms of the peptide bond is about 5.4 V, and thus, the approximate solvent electric field produced by surrounding water molecules on the NMA is 172 MV/cm on average if the molecular geometry is taken into account. The solute-solvent electrostatic interaction energy is estimated to be -137 kJ/mol, by considering electric dipole-electric field interaction. Furthermore, their root-mean-square fluctuation amplitudes are as large as 1.6 V, 52 MV/cm, and 41 kJ/mol, respectively. We found that the water electric potential on a peptide bond is spatially nonhomogeneous and that the fluctuation in the electrostatic peptide-water interaction energy is about 10 times larger than the thermal energy at room temperature. This indicates that the peptide-solvent interactions are indeed important for the activation of chemical reactions in aqueous solution.     

On the influence of the water electrostatic field on the amide group vibrational frequencies

For clusters of N-methylacetamide and water molecules the performance of the fixed-charged approximation was tested against continuum and explicit ab initio models. The dispersion of the vibrational frequencies when constant electrostatic potential was maintained at the solute atoms was compared to the distribution caused by geometry fluctuations.     

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.     

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.     

Geometry and solvent dependence of the electronic spectra of the amide group and consequences for peptide circular dichroism

The influence of geometry variations and solvent environment of N-methylacetamide on its energies and absorption intensities was systematically analyzed with the aid of the time-dependent density functional theory (TD DFT). Selective and often complicated reactions of individual electronic levels on the perturbations were found important for the resultant spectral profile. For example, the n-pi band position varied by tens of nanometers due to the C=O bond length oscillations, while it was rather unaffected by surrounding water. On the contrary, pi-pi type transition energies and intensities were broadly dispersed by the aqueous environment but exhibited a modest coordinate dependence. A simple electrostatic model used previously for absorption in the IR region (J. Chem. Phys. 2005, 122, 144501) explained these changes only partially. Additionally, electronic transfer between the solute and the solvent had to be considered for faithful modeling of the ultraviolet light absorption. The inclusion of the environment and dynamics in the modeling then provided more accurate positions, intensities, and realistic inhomogeneous widths of spectral lines. These factors were found important for absorption and circular dichroism spectra of larger peptides and proteins. This was demonstrated with a combined DFT/coupled oscillator model providing principal features observed in electronic circular dichroism spectra of standard peptide conformations.     

Heat capacity effects associated with the hydrophobic hydration and interaction of simple solutes: a detailed structural and energetical analysis based on molecular dynamics simulations

We examine the SPCE [H. J. C. Berendsen et al., J. Chem. Phys. 91, 6269 (1987)] and TIP5P [M. W. Mahoney and W. L. Jorgensen, J. Chem. Phys 112, 8910 (2000)] water models using a temperature series of molecular-dynamics simulations in order to study heat-capacity effects associated with the hydrophobic hydration and interaction of xenon particles. The temperature interval between 275 and 375 K along the 0.1-MPa isobar is studied. For all investigated models and state points we calculate the excess chemical potential for xenon employing the Widom particle insertion technique. The solvation enthalpy and excess heat capacity is obtained from the temperature dependence of the chemical potentials and, alternatively, directly by Ewald summation, as well as a reaction field based method. All three methods provide consistent results. In addition, the reaction field technique allows a separation of the solvation enthalpy into solute/solvent and solvent/solvent parts. We find that the solvent/solvent contribution to the excess heat capacity is dominating, being about one order of magnitude larger than the solute/solvent part. This observation is attributed to the enlarged heat capacity of the water molecules in the hydration shell. A detailed spatial analysis of the heat capacity of the water molecules around a pair of xenon particles at different separations reveals that even more enhanced heat capacity of the water located in the bisector plane between two adjacent xenon atoms is responsible for the maximum of the heat capacity found for the desolvation barrier distance, recently reported by Shimizu and Chan [J. Am. Chem. Soc. 123, 2083 (2001)]. The about 60% enlarged heat capacity of water in the concave part of the joint xenon-xenon hydration shell is the result of a counterplay of strengthened hydrogen bonds and an enhanced breaking of hydrogen bonds with increasing temperature. Differences between the two models with respect to the heat capacity in the xenon-xenon contact state are attributed to the different water model bulk heat capacities, and to the different spatial extension of the structure effect introduced by the hydrophobic particles. Similarities between the different states of water in the joint xenon-xenon hydration shell and the properties of stretched water are discussed.     

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.     

Molecular dynamics simulations and instantaneous normal-mode analysis of the vibrational relaxation of the C-H stretching modes of N-methylacetamide-d in liquid deuterated water

Nonequilibrium molecular dynamics (MD) simulations and instantaneous normal mode (INMs) analyses are used to study the vibrational relaxation of the C-H stretching modes (ν(s)(CH?)) of deuterated N-methylacetamide (NMAD) in aqueous (D2O) solution. The INMs are identified unequivocally in terms of the equilibrium normal modes (ENMs), or groups of them, using a restricted version of the recently proposed Min-Cost assignment method. After excitation of the parent ν(s)(CH?) modes with one vibrational quantum, the vibrational energy is shown to dissipate through both intramolecular vibrational redistribution (IVR) and intermolecular vibrational energy transfer (VET). The decay of the vibrational energy of the ν(s)(CH?) modes is well fitted to a triple exponential function, with each characterizing a well-defined stage of the entire relaxation process. The first, and major, relaxation stage corresponds to a coherent ultrashort (τ(rel) = 0.07 ps) energy transfer from the parent ν(s)(CH?) modes to the methyl bending modes δ(CH?), so that the initially excited state rapidly evolves into a mixed stretch-bend state. In the second stage, characterized by a time of 0.92 ps, the vibrational energy flows through IVR to a number of mid-range-energy vibrations of the solute. In the third stage, the vibrational energy accumulated in the excited modes dissipates into the bath through an indirect VET process mediated by lower-energy modes, on a time scale of 10.6 ps. All the specific relaxation channels participating in the whole relaxation process are properly identified. The results from the simulations are finally compared with the recent experimental measurements of the ν(s)(CH?) vibrational energy relaxation in NMAD/D?O(l) reported by Dlott et al. (J. Phys. Chem. A 2009, 113, 75.) using ultrafast infrared-Raman spectroscopy.     

The behavior of fluids near solutes: a density functional theory and computer simulation study

The density distribution of solvent near a solute particle is studied using density functional theory and Monte Carlo simulation. The fluid atoms interact with each other via a hard sphere plus Yukawa potential, and interact with the solute via a hard sphere potential. For small solute sizes, the solvent displays liquidlike ordering near the particle. When the solute become larger, a drying transition is observed at state points near the coexistence conditions of the solvent. These predictions are similar to those of a recent theory for the hydrophobic effect by Lum, Chandler, and Weeks [J. Phys. Chem. 103, 4570 (1999)], although a comparison with simulations shows that the theory of this work is quantitatively more accurate. The connection between density functional methods and the LCW approach is also established.     

Hofmeister anionic effects on hydration electric fields around water and peptide

Specific ion effects on water dynamics and local solvation structure around a peptide are important in understanding the Hofmeister series of ions and their effects on protein stability in aqueous solution. Water dynamics is essentially governed by local hydrogen-bonding interactions with surrounding water molecules producing hydration electric field on each water molecule. Here, we show that the hydration electric field on the OD bond of HOD molecule in water can be directly estimated by measuring its OD stretch infrared (IR) radiation frequency shift upon increasing ion concentration. For a variety of electrolyte solutions containing Hofmeister anions, we measured the OD stretch IR bands and estimated the hydration electric field on the OD bond to be about a hundred MV∕cm with standard deviation of tens of MV∕cm. As anion concentration increases from 1 to 6 M, the hydration electric field on the OD bond decreases by about 10%, indicating that the local H-bond network is partially broken by dissolved ions. However, the measured hydration electric fields on the OD bond and its fluctuation amplitudes for varying anions are rather independent on whether the anion is a kosmotrope or a chaotrope. To further examine the Hofmeister effects on H-bond solvation structure around a peptide bond, we examined the amide I' and II' mode frequencies of N-methylacetamide in various electrolyte D(2)O solutions. It is found that the two amide vibrational frequencies are not affected by ions, indicating that the H-bond solvation structure in the vicinity of a peptide remains the same irrespective of the concentration and character of ions. The present experimental results suggest that the Hofmeister anionic effects are not caused by direct electrostatic interactions of ions with peptide bond or water molecules in its first solvation shell. Furthermore, even though the H-bond network of water is affected by ions, thus induced change of local hydration electric field on the OD bond of HOD is not in good correlation with the well-known Hofmeister series. We anticipate that the present experimental results provide an important clue about the Hofmeister effect on protein structure and present a discussion on possible alternative mechanisms.     

A theoretical study of P4O10: vibrational analysis, infrared and Raman spectra

The normal mode frequencies and the corresponding vibrational assignments of tetraphosphorus decaoxide (P4O10) in tetrahedral (Td) symmetry are examined theoretically and experimentally. The Gaussian 98 set of quantum chemistry codes at the HF/6-311G*, MP2/6-311G*, and DFT/B3LYP/6-311G* levels of theory are used. By comparison to experimental normal mode frequencies deduced by Gilliam et al. [J. Phys. Chem. B 107 (2003) 2892], Chapman [Spectrochim. Acta A, 24 (1968) 1687], Beattie et al. [J. Chem. Soc. A (1970) 449], Konings et al. [J. Mol. Spectrosc. 152 (1992) 29] and the present work, correction factors for predominant vibrational motions are determined and compared. Normal modes were decomposed into five non-redundant motions (P-O stretch, P=O stretch, P-O-P bend, P-O-P wag, and P=O wag). Standard deviations found for the HF, MP2, and DFT corrected frequencies compared to experiment are 9, 5, and 4 cm(-1), respectively. Electron distribution for selected molecular orbitals is considered.     

Extension of the density functional derived scaled quantum mechanical force field procedure

The density functional derived scaled quantum mechanical (SQM) force field method of Rauhut and Pulay [J. Phys. Chem. 99 (1995) 3093] has been extended. The original procedure (employing B3-LYP/6-31G* computations and 11 transferable scale factors for the different kinds of internal coordinates) was capable to reproduce the vibrational fundamentals of 31 simple organic (H, C, N, O) molecules with a total mean deviation of about 13 cm(-1). The present Density Functional Theory based SQM force field method is an extension of the original one: with the help of 20 transferable scale factors can reproduce the fundamentals of 20 inorganic, organic and organosilicon molecules containing nonmetallic first and second-row atoms with a total mean deviation of 10.8 cm(-1). The transferability and reliability of the new set of scale factors are demonstrated on the examples of the a priori SQM vibrational spectra of cis and gauche-cyclopropylchlorosilane.     

Reconsideration on hydrogen bond strengthening or cleavage of photoexcited coumarin 102 in aqueous solvent: a DFT/TDDFT study

In this work, the excited-state hydrogen bonding dynamics of photoexcited coumarin 102 in aqueous solvent is reconsidered. The electronically excited states of the hydrogen bonded complexes formed by coumarin 102 (C102) chromophore and the hydrogen donating water solvent have been investigated using the time-dependent density functional theory method. Two intermolecular hydrogen bonds between C102 and water molecules are considered. The previous works (Wells et al., J Phys Chem A 2008, 112, 2511) have proposed that one intermolecular hydrogen bond would be strengthened and the other one would be cleaved upon photoexcitation to the electronically excited states. However, our theoretical calculations have demonstrated that both the two intermolecular hydrogen bonds between C102 solute and H(2)O solvent molecules are significantly strengthened in electronically excited states by comparison with those in ground state. Hence, we have confirmed again that intermolecular hydrogen bonds between C102 chromophore and aqueous solvents are strengthened not cleaved upon electronic excitation, which is in accordance with Zhao's works.     

Vibronic coupling in asymmetric bichromophores: Theory and application to diphenylmethane

The theory for modeling vibronic interactions in bichromophores was introduced in sixties by Witkowski and Moffitt [J. Chem. Phys. 33, 872 (1960)] and extended by Fulton and Gouterman [J. Chem. Phys. 35, 1059 (1961)]. The present work describes extension of this vibronic model to describe bichromophores with broken vibrational symmetry such as partly deuterated molecules. Additionally, the model is extended to include inter-chromophore vibrational modes. The model can treat multiple vibrational modes by employing Lanczos diagonalization procedure of sparse matrices. The developed vibronic model is applied to simulation of vibronic spectra of flexible bichromophore diphenylmethane and compared to high-resolution experimental spectra [J. A. Stearns, N. R. Pillsbury, K. O. Douglass, C. W. Mu?ller, T. S. Zwier, and D. F. Plusquellic, J. Chem. Phys. 129, 224305 (2008)].     

A comparative study of the properties of polar and nonpolar solvent/solute/polystyrene solutions in microwave fields via molecular dynamics

The influence of an applied microwave field on the dynamics of methylamine-dichloromethane (DCM) mixtures bound within atactic polystyrene (a-PS) over a range of polymer densities from 30 to 94 wt % polymer was examined using atomistic molecular dynamics simulations. This study is an extension of previous studies on methylamine transport in relatively polar polystyrene solutions of methanol and dimethylformamide [M. J. Purdue et al., J. Chem. Phys. 124, 204904 (2006)]. A direct comparison is made across the three types of polystyrene solutions. Consideration is given to both solvent and reagent transport within the polymer solutions under zero-field conditions and in an external electromagnetic field in the canonical ensemble (NVT) at 298.0 K. Various frequencies up to 10(4) GHz and a rms electric field intensity of 0.1 VA were applied. The simulation studies were validated by comparison of the simulated zero-field self-diffusion coefficients of DCM in a-PS with those obtained using pulsed-gradient spin-echo NMR spectrometry. Athermal effects of microwave fields on solute transport behavior within polymer solutions are discussed.     

Partial Hessian vibrational analysis: the localization of the molecular vibrational energy and entropy

The partial Hessian vibrational analysis (PHVA), in which only a subblock of the Hesssian matrix is diagonalized to yield vibrational frequencies for partially optimized systems, is extended to the calculation of vibrational enthalpy and entropy changes for chemical reactions. The utility of this method is demonstrated for various deprotonation reactions by reproducing full HVA values to within 0.1–0.4 kcal/mol, depending on the number atoms included in the PHVA. When combined with the hybrid effective fragment potential method [Gordon MS, et al. (2001) J Phys Chem A 105:293–307], the PHVA method can provide (harmonic) free-energy changes for localized chemical reactions in very large systems. Received: 21 September 2001 / Accepted: 30 October 2001 / Published online: 22 March 2002     

Vibrational energy relaxation of azulene studied by the transient grating method. I. Supercritical fluids

The vibrational energy dissipation process of the ground-state azulene in supercritical xenon, carbon dioxide, and ethane has been studied by the transient grating spectroscopy. In this method, azulene in these fluids was photoexcited by two counterpropagating subpicosecond laser pulses at 570 nm, which created a sinusoidal pattern of vibrationally hot ground-state azulene inside the fluids. The photoacoustic signal produced by the temperature rise of the solvent due to the vibrational energy relaxation of azulene was monitored by the diffraction of a probe pulse. The temperature-rise time constants of the solvents were determined at 383 and 298 K from 0.7 to 2.4 in rho(r), where rho(r) is the reduced density by the critical density of the fluids, by the fitting of the acoustic signal based on a theoretical model equation. In xenon, the temperature-rise time constant was almost similar to the vibrational energy-relaxation time constant of the photoexcited solute determined by the transient absorption measurement [D. Schwarzer, J. Troe, M. Votsmeier, and M. Zerezke, J. Chem. Phys. 105, 3121 (1996)] at the same reduced density irrespective of the solvent temperature. On the other hand, the temperature-rise time constants in ethane were larger than the vibrational energy-relaxation time constants by a factor of about 2. In carbon dioxide, the difference was small. From these results, the larger time constants of the solvent temperature rise than those of the vibrational energy relaxation in ethane and carbon dioxide were interpreted in terms of the vibrational-vibrational (V-V) energy transfer between azulene and solvent molecules and the vibrational-translational (V-T) energy transfer between solvent molecules. The contribution of the V-V energy transfer process against the V-T energy transfer process has been discussed.     

Dependence of amide vibrations on hydrogen bonding

The effect of hydrogen bonding on the amide group vibrational spectra has traditionally been rationalized by invoking a resonance model where hydrogen bonding impacts the amide functional group by stabilizing its [(-)O-C=NH (+)] structure over the [O=C-NH] structure. However, Triggs and Valentini's UV-Raman study of solvation and hydrogen bonding effects on epsilon-caprolactum, N, N-dimethylacetamide (DMA), and N-methylacetamide (NMA) ( Triggs, N. E.; Valentini, J. J. J. Phys. Chem. 1992, 96, 6922-6931) casts doubt on the validity of this model by demonstrating that, contrary to the resonance model prediction, carbonyl hydrogen bonding does not impact the AmII' frequency of DMA. In this study, we utilize density functional theory (DFT) calculations to examine the impact of hydrogen bonding on the C=O and N-H functional groups of NMA, which is typically used as a simple model of the peptide bond. Our calculations indicate that, as expected, the hydrogen bonding frequency dependence of the AmI vibration predominantly derives from the C=O group, whereas the hydrogen bonding frequency dependence of the AmII vibration primarily derives from N-H hydrogen bonding. In contrast, the hydrogen bonding dependence of the conformation-sensitive AmIII band derives equally from both C=O and N-H groups and thus, is equally responsive to hydrogen bonding at the C=O or N-H site. Our work shows that a clear understanding of the normal mode composition of the amide vibrations is crucial for an accurate interpretation of the hydrogen bonding dependence of amide vibrational frequencies.     

Study on the vibrational energy relaxation of p-nitroaniline, N,N-dimethyl-p-nitroaniline, and azulene by the transient grating method

The vibrational energy dissipation processes of the electronic ground states of p-nitroaniline and N,N-dimethyl-p-nitroaniline have been studied by transient grating spectroscopy with subpicosecond laser pulses. The rise time of the acoustic signal produced by the energy dissipation process of the hot ground state molecule was monitored. The acoustic signal was analyzed by an equation including the acoustic damping. The solvent temperature rise times in various solvents have been determined. The acoustic signals of azulene in previous papers [Y. Kimura et al., J. Chem. Phys. 123, 054512 (2005); 123, 054513 (2005)] were also reanalyzed using this equation. The temperature rise times in all cases are longer than the vibrational energy relaxation times of the solutes determined by the transient absorption measurements. The difference is discussed in terms of the energy transfer pathways from the solute to the solvent. We concluded that both the hydrogen bonding between the solute and the solvent and the lower frequency modes of the solutes play important roles in determining the energy transfer pathway from the solute to the solvent.