Molecular dynamics integration and molecular vibrational theory. III. The infrared spectrum of water

The new symplectic molecular dynamics (MD) integrators presented in the first paper of this series were applied to perform MD simulations of water. The physical properties of a system of flexible TIP3P water molecules computed by the new integrators, such as diffusion coefficients, orientation correlation times, and infrared (IR) spectra, are in good agreement with results obtained by the standard method. The comparison between the new integrators' and the standard method's integration time step sizes indicates that the resulting algorithm allows a 3.0 fs long integration time step as opposed to the standard leap-frog Verlet method, a sixfold simulation speed-up. The accuracy of the method was confirmed, in particular, by computing the IR spectrum of water in which no blueshifting of the stretching normal mode frequencies is observed as occurs with the standard method.     

Theoretical calculations of infrared absorption, vibrational circular dichroism, and two-dimensional vibrational spectra of acetylproline in liquids water and chloroform

Infrared absorption, vibrational circular dichroism, and two-dimensional infrared pump-probe and photon echo spectra of acetylproline solutions are theoretically calculated and directly compared with experiments. In order to quantitatively determine interpeptide interaction-induced amide I mode frequency shifts, high-level quantum chemistry calculations were performed. The solvatochromic amide I mode frequency shift and fluctuation were taken into account by carrying out molecular dynamics simulations of acetylproline dissolved in liquids water and chloroform and by using the extrapolation method developed recently. We first studied correlation time scales of the two amide I vibrational frequency fluctuations, cross correlation between the two fluctuating local mode frequencies, ensemble averaged conformations of the acetylproline molecule in liquids water and chloroform. The corresponding conformations of the acetylproline in liquids water and chloroform are close to the ideal 3(10) helix and the C(7) structure, respectively. A few methods proposed to determine the angle between the two transition dipoles associated with the amide I vibrations were tested and their limitations are discussed.     

Potential of hydrogen bond in water. Comparison of the theory with vibrational spectra and results of molecular dynamics simulations

Potential of hydrogen bond is the function which relates its energy to geometrical parameters of hydrogen bridge: its length R(O…O) and angles between direction O…O and OH group [φ (H-O…O)] and/or lone pair of proton accepting oxygen atom [χ(-O…O)]. Previously we have suggested an approach to design such potentials based on the approximate numerical solution of a reverse problem of the spectrum band shape in the frames of the fluctuation theory of hydrogen bonding. In the given work this method is applied to construction of the two-parameter potentials that depend on parameters {R(O…O), φ (H-O…O} or {φ (H-O…O), χ (-O…O)}. Using them, the spectra of OH vibrations of HOD molecules in a liquid phase are calculated theoretically in good agreement with experiment in the temperature range up to 200 °C. Distributions of angles P(φ, T), P(χ, T), and energies P(E) are calculated also. The same distributions and spectra are independently calculated on the basis of the geometrical parameters of the hydrogen bridges obtained from molecular dynamics models of water. The shapes of the spectral contours and their temperature evolution calculated for computer models turned out to be similar to experimental ones only for the potential that includes the length of H-bond R(O…O).     

Simulation of vibrational spectra of polymers by molecular dynamics calculations

Molecular dynamics is used to study the vibrational spectra (infrared and Raman) of polyethylene and polytetrafluoroethylene. The calculated frequencies are in good agreement with results of normal mode analyses. The calculated intensities lead to the conclusion that, while satisfactory relative intensities can be expected from simple models, quantitative determination of the intensities necessitates the use of complete sets of electro-optical parameters.     

Experimental, theoretical calculations of the vibrational spectra and conformational analysis of 2,4-di-tert-butylphenol

In the present work, we reported a combined experimental and theoretical study on conformational stability, molecular structure and vibrational spectra of 2,4-di-tert-butylphenol (2,4-DTBP). The FT-IR (400-4000cm(-1)) and FT-Raman spectra (50-3500cm(-1)) of 2,4-DTBP were recorded. The molecular geometry, harmonic vibrational frequencies and bonding features of 2,4-DTBP in the ground-state have been calculated by using the density functional BLYP/B3LYP methods. The energy calculated by time-dependent density functional theory (TD-DFT) result complements with the experimental findings. The calculated highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) energies show that charge transfer occurs within the molecule. Finally the calculation results were compared with measured infrared and Raman spectra of the title compound which showed good agreement with observed spectra.     

Vibrational spectra,conformational equilibria,ab initio calculations and vibrational assignments of ethylmethylfluorosilane

The infrared spectra of ethylmethylfluorosilane (CH₃SiHFCH₂CH₃) have been recorded as a vapour, liquid and solid at 78 K in the 4000–50 cm⁻¹ range and isolated in an argon matrix at ca. 5 K. Infrared spectra of two different solid phases were obtained after annealing to temperatures of 120 and 130 K, and recooling to 78 K. Although the IR spectra were quite similar in the MIR region, certain differences were noted in the FIR region below 400 cm⁻¹. The most stable conformer MeMe was present after annealing to 130 K, but three bands belonging to MeH were detected after annealing to 120 K. Various infrared bands changed intensity when the argon matrix was annealed to temperatures between 20 and 35 K, and some of these were related to changes in the conformational abundance.Raman spectra of the liquid were recorded at room temperature and at various temperatures between 295 and 153 K. Spectra of an amorphous and annealed solid were recorded at 78 K. In the variable temperature Raman spectra, various bands changed in intensity and were interpreted in terms of conformational equilibria between the three possible conformers. Complete assignments were made for all the bands of the most stable conformer MeMe. From various bands assigned to the three conformers, the conformational enthalpy difference ΔH from MeMe to the intermediate energy conformer MeH was found to be 0.5 kJ mol⁻¹ and to the highest conformer MeF was 0.7 kJ mol⁻¹. At ambient temperature this leads to 39% MeMe, 32% MeH and 29% of the MeF conformer in the liquid.Ab initio calculations in the RHF, MP2, DFT approximations and very accurate G2 calculations were carried out. With one exception, the MeMe conformer had the lowest enthalpy in all these calculations, the MeH had the intermediate and the MeF the highest enthalpy, and the calculations were in good agreement with the measurements.     

Vibrational spectra, structure, and theoretical calculations of 2-fluoro- and 3-fluoropyridine

The infrared and Raman spectra of liquid and vapor-phase 2-fluoropyridine and 3-fluoropyridine have been recorded and assigned. Ab initio and DFT calculations were carried out to compute the molecular structures and to verify the vibrational assignments. The observed and calculated spectra agree extremely well. The ring bond distances of the fluoropyridines are very similar to those of pyridine except for a shortening of the C-N(F) bond in 2-fluoropyridine. The C-F bond stretching frequencies are similar to that in fluorobenzene reflecting the influence of the ring π bonding.     

Molecular dynamics simulations for CO2 spectra. III. Permanent and collision-induced tensors contributions to light absorption and scattering

Classical molecular dynamics simulations have been performed for gaseous CO(2) starting from an accurate anisotropic intermolecular potential. Through calculations of the evolutions of the positions and orientations of a large number of molecules, the time evolutions of the permanent and collision-induced electric dipole vector and polarizability tensor are obtained. These are computed from knowledge of static molecular parameters taking only the leading induction terms into account. The Laplace transforms of the auto-correlation functions of these tensors then directly yield the light absorption and scattering spectra. These predictions are, to our knowledge, the first in which the contributions of permanent and collision-induced tensors are simultaneously taken into account for gaseous CO(2), without any adjusted parameter. Comparisons between computations and measurements are made for absorption in the region of the ν(3) infrared band and for depolarized Rayleigh scattering in the roto-translational band. They demonstrate the quality of the model over spectral ranges from the band center to the far wings where the spectrum varies by several orders of magnitude. The contributions of the permanent and interaction-induced (dipole and polarizability) tensors are analyzed for the first time, through the purely permanent (allowed), purely induced, and cross permanent∕induced components of the spectra. It is shown that, while the purely induced contribution is negligible when compared to the collision-broadened allowed component, the cross term due to interferences between permanent and induced tensors significantly participates to the wings of the bands. This successfully clarifies the long lasting, confusing situation for the mechanisms governing the wings of the CO(2) spectra considered in this work.     

Vibrational and surface-enhanced vibrational spectra of 6-nitrochrysene

The infrared and Raman spectra of solids and thin solid films of 6-nitrochrysene, its electronic spectra, and resonance Raman scattering (RRS) obtained with UV-laser excitation at 325 nm are reported. The vibrational assignment is supported by ab initio computations at the B3LYP/6-311G(d, p) level of theory. The molecular organization in nanometric films evaporated onto smooth metal surfaces of silver and copper was probed using reflection-absorption infrared spectroscopy (RAIRS). The results of the surface-enhanced Raman scattering (SERS) and surface-enhanced infrared absorption (SEIRA) obtained from nanometric films evaporated onto silver island films are also discussed. It was found that the molecule efficiently interacts with silver island film surfaces, and that the interaction leads to extensive photochemical reaction at the metal surface under laser illumination.     

Molecular dynamics simulations of beta-cyclodextrin-aziadamantane complexes in water

Force-field-based atomistic simulations of host-guest supramolecular complexes between beta-cyclodextrin and several aziadamantane derivatives have been analyzed with respect to relative orientation and interaction energies, explicitly considering solvent (water) molecules. For each case, the calculations revealed two stable orientations of the guest within the host that are different in interaction energy. Fluctuation of and correlation between characteristic properties were analyzed. Among other things, it turned out that orientation angle and inclusion depth are clearly correlated. In addition, for the unsubstituted aziadamantane, the enthalpy of complex formation was calculated and compared to experimental results.     

Vibrational spectra, conformations, ab initio calculations and vibrational assignments of 3-pentyn-2-ol

The infrared spectra of 3-pentyn-2-ol, CH₃CCCH(OH)CH₃, have been recorded as a vapour and liquid at ambient temperature, as a solid at 78 K in the 4000–50 cm⁻¹ range and isolated in an argon matrix at ca. 5 K. Infrared spectra of the solid phase at 78 K were obtained before and after annealing to temperatures of 120 and 130 K. The IR spectra of the solid were quite similar to that of the liquid.

Raman spectra of the liquid were recorded at room temperature and at various temperatures between 295 and 153 K. Spectra of an amorphous and annealed solid were recorded at 78 K. In the variable temperature Raman spectra, some bands changed in relative intensity and were interpreted in terms of conformational equilibria between the three possible conformers. Complete assignments were made for all the bands of the most stable conformer in which OH is oriented anti to C₁(a_Me). From various bands assigned to a second conformer in which OH is oriented anti to H_gem(a_H), the conformational enthalpy differences was found to be between 0.4 and 0.8 kJ mol⁻¹. The highest energy conformer with OH anti to C₃(a_C) was not detected.

Quantum-chemical calculations have been carried out at the MP2 and B3LYP levels with a variety of basis sets. Except for small basis set calculations for which the a_H conformer had slightly lower energy, all the calculations revealed that a_Me was the low energy conformer. The B3LYP/cc-pVTZ calculations suggested the a_Me conformer as more stable by 0.8 and 8.3 kJ mol⁻¹ relative to a_H an a_C, respectively. Vibrational wavenumbers and infrared and Raman band intensities for two of the three conformers are reported from B3LYP/cc-pVTZ calculations.     

Multidimensional infrared spectroscopy of water. I. Vibrational dynamics in two-dimensional IR line shapes

In this and the following paper, we describe the ultrafast structural fluctuations and rearrangements of the hydrogen bonding network of water using two-dimensional (2D) infrared spectroscopy. 2D IR spectra covering all the relevant time scales of molecular dynamics of the hydrogen bonding network of water were studied for the OH stretching absorption of HOD in D2O. Time-dependent evolution of the 2D IR line shape serves as a spectroscopic observable that tracks how different hydrogen bonding environments interconvert while changes in spectral intensity result from vibrational relaxation and molecular reorientation of the OH dipole. For waiting times up to the vibrational lifetime of 700 fs, changes in the 2D line shape reflect the spectral evolution of OH oscillators induced by hydrogen bond dynamics. These dynamics, characterized through a set of 2D line shape analysis metrics, show a rapid 60 fs decay, an underdamped oscillation on a 130 fs time scale induced by hydrogen bond stretching, and a long time decay constant of 1.4 ps. 2D surfaces for waiting times larger than 700 fs are dominated by the effects of vibrational relaxation and the thermalization of this excess energy by the solvent bath. Our modeling based on fluctuations with Gaussian statistics is able to reproduce the changes in dispersed pump-probe and 2D IR spectra induced by these relaxation processes, but misses the asymmetry resulting from frequency-dependent spectral diffusion. The dynamical origin of this asymmetry is discussed in the companion paper.     

Vibrational dynamics of DNA. II. Deuterium exchange effects and simulated IR absorption spectra

In Paper I, we studied vibrational properties of normal bases, base derivatives, Watson-Crick base pairs, and multiple layer base pair stacks in the frequency range of 1400-1800 cm(-1). However, typical IR absorption spectra of single- and double-stranded DNA have been measured in D(2)O solution. Consequently, the more relevant bases and base pairs are those with deuterium atoms in replacement with labile amino hydrogen atoms. Thus, we have carried out density functional theory vibrational analyses of properly deuterated bases, base pairs, and stacked base pair systems. In the frequency range of interest, both aromatic ring deformation modes and carbonyl stretching modes appear to be strongly IR active. Basis mode frequencies and vibrational coupling constants are newly determined and used to numerically simulate IR absorption spectra. It turns out that the hydration effects on vibrational spectra are important. The numerically simulated vibrational spectra are directly compared with experiments. Also, the (18)O-isotope exchange effect on the poly(dG):poly(dC) spectrum is quantitatively described. The present calculation results will be used to further simulate two-dimensional IR photon echo spectra of DNA oligomers in the companion Paper III.     

Vibrational spectra from atomic fluctuations in dynamics simulations. II. Solvent-induced frequency fluctuations at femtosecond time resolution

The midinfrared (MIR) spectra of molecules in polar solvents exhibit inhomogeneously broadened bands whose spectral positions are shifted as compared to the gas phase. The shifts are caused by interactions with structured solvation shells and the broadenings by fluctuations of these interactions. The MIR spectra can be calculated from hybrid molecular dynamics (MD) simulations, which treat the solute molecule by density functional theory and the solvent by molecular mechanics by the so-called instantaneous normal mode analysis (INMA) or by Fourier transforming the time correlation function (FTTCF) of the molecular dipole moment. In Paper I of this work [M. Schmitz and P. Tavan, J. Chem. Phys. 121, 12233 (2004)] we explored an alternative method based on generalized virial (GV) frequencies noting, however, that GV systematically underestimates frequencies. As shown by us these artifacts are caused by solvent-induced fluctuations of the (i) equilibrium geometry, (ii) force constants, and (iii) normal mode directions as well as by (iv) diagonal and (v) off-diagonal anharmonicities. Here we now show, by analyzing the time scales of fluctuations and sample MD trajectories of formaldehyde in the gas phase and in water, that all these sources of computational artifacts can be made visible by a Fourier analysis of the normal coordinates. Correspondingly, the error sources (i) and (iii)-(v) can be removed by bandpass filtering, as long as the spectral signatures of the respective effects are well separated from the fundamental band. Furthermore, the artifacts arising from effect (ii) can be strongly diminished by a time-resolved version of the GV approach (TF-GV). The TF-GV method then yields for each mode j a trajectory of the vibrational frequency omega(j)(tmid R:tau) at a time resolution tau>tau(j), which is only limited by the corresponding oscillation time tau(j)=2pi/omega(j) and, thus, is in the femtosecond range. A correlation analysis of these trajectories clearly separates the librational motions from the conformational dynamics of the solvation shells and yields the inhomogeneously broadened MIR spectra, if the theory of motional narrowing is properly included. The MIR spectrum of formaldehyde in solution obtained by TF-GV agrees very well with the FTTCF result, if one applies the so-called "harmonic approximation" quantum correction factor and a temperature scaling to the FTTCF intensities. Also for INMA an excellent agreement is achieved if one disregards a slight INMA overestimate of linewidths.     

Vibrational dynamics of DNA. I. Vibrational basis modes and couplings

Carrying out density functional theory calculations of four DNA bases, base derivatives, Watson-Crick (WC) base pairs, and multiple-layer base pair stacks, we studied vibrational dynamics of delocalized modes with frequency ranging from 1400 to 1800 cm(-1). These modes have been found to be highly sensitive to structure fluctuation and base pair conformation of DNA. By identifying eight fundamental basis modes, it is shown that the normal modes of base pairs and multilayer base pair stacks can be described by linear combinations of these vibrational basis modes. By using the Hessian matrix reconstruction method, vibrational coupling constants between the basis modes are determined for WC base pairs and multilayer systems and are found to be most strongly affected by the hydrogen bonding interaction between bases. It is also found that the propeller twist and buckle motions do not strongly affect vibrational couplings and basis mode frequencies. Numerically simulated IR spectra of guanine-cytosine and adenine-thymine bases pairs as well as of multilayer base pair stacks are presented and described in terms of coupled basis modes. It turns out that, due to the small interlayer base-base vibrational interactions, the IR absorption spectrum of multilayer base pair system does not strongly depend on the number of base pairs.     

Vibrational spectra and quantum chemical calculations of 3,4-diaminobenzoic acid.

The Fourier Transform Raman and Fourier Transform infrared spectra of 3,4-diaminobenzoic acid (3,4-DABA) were recorded in the solid phase. Geometry optimizations were done without any constraint and harmonic-vibrational wave numbers and several thermodynamic parameters were calculated for the minimum energy conformer at ab initio and DFT levels invoking 6-311++G(d,p) basis set. The results were compared with the experimental values with the help of specific scaling procedures, the observed vibrational wavenumbers in FT-IR and FT-Raman spectra were analyzed and assigned to different normal modes of the molecule. Most of the modes have wavenumbers in the expected range, the error obtained was in general very low. The appropriate theoretical spectrograms for the IR and Raman spectra of the title molecule were also constructed.     

Molecular dynamics simulations for water and ions in protein crystals

The spatial and temporal properties of water and ions in bionanoporous materials-protein crystals-have been investigated using molecular dynamics simulations. Three protein crystals are considered systematically with different morphologies and chemical topologies: tetragonal lysozyme, orthorhombic lysozyme, and tetragonal thermolysin. It is found that the thermal fluctuations of C(alpha) atoms in the secondary structures of protein molecules are relatively weak due to hydrogen bonding. The solvent-accessible surface area per residue is nearly identical in the three protein crystals; the hydrophobic and hydrophilic residues in each crystal possess approximately the same solvent-accessible surface area. Water distributes heterogeneously and has different local structures within the biological nanopores of the three protein crystals. The mobility of water and ions in the crystals is enhanced as the porosity increases and also by the fluctuations of protein atoms particularly in the two lysozyme crystals. Anisotropic diffusion is found preferentially along the pore axis, as experimentally observed. The anisotropy of the three crystals increases in the order: tetragonal thermolysin < tetragonal lysozyme < orthorhombic lysozyme.     

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.     

Coupled calculations of vibrational frequencies and intensities. III. IR and Raman spectra of ethylene oxide and ethylene sulfide

The experimental IR and Raman spectra of ethylene oxide have been reinvestigated with particular attention to the intensities. The absolute IR intensities have been measured for the gaseous state. The spectra have been simulated by using a normal coordinate analysis coupled with a CNDO determination of the intensities. The intensity calculation using polarization functions appears to be more reliable than the standard version. Furthermore, the force field has been extended for ethylene sulfide.