Redistribution of carbonyl stretch mode energy in isolated and solvated N-methylacetamide: kinetic energy spectral density analyses

Effects of interfaces on structure and dynamics of water droplets on a graphene surface: A molecular dynamics study

The structure and dynamics of water droplets on a bilayer graphene surface are investigated using molecular dynamics simulations. The effects of solid/water and air/water interfaces on the local structure of water droplets are analyzed in terms of the hydrogen bond distribution and tetrahedral order parameter. It is found that the local structure in the core region of a water droplet is similar to that in liquid water. On the other hand, the local structure of water molecules at the solid/water and air/water interfaces, referred to as the interface and surface regions, respectively, consists mainly of three-coordinated molecules that are greatly distorted from a tetrahedral structure. This study reveals that the dynamics in different regions of the water droplets affects the intermolecular vibrational density of states: It is found that in the surface and interface regions, the intensity of vibrational density of states at ∼50 cm^-1 is enhanced, whereas those at ∼200 and ∼500 cm^-1 are weakened and redshifted. These changes are attributed to the increase in the number of molecules having fewer hydrogen bonds in the interface and surface regions. Both single-molecule and collective orientation relaxations are also examined. Single-molecule orientation relaxation is found to be marginally slower than that in liquid water. On the other hand, the collective orientation relaxation of water droplets is found to be significantly faster than that of liquid water because of the destructive correlation of dipole moments in the droplets. The negative correlation between distinct dipole moments also yields a blueshifted libration peak in the absorption spectrum. It is also found that the water-graphene interaction affects the structure and dynamics of the water droplets, such as the local water structure, collective orientation relaxation, and the correlation between dipole moments. This study reveals that the water/solid and water/air interfaces strongly affect the structure and intermolecular dynamics of water droplets and suggests that the intermolecular dynamics, such as energy relaxation dynamics, in other systems with interfaces are different from those in liquid water.

Vibrational energy relaxation of interfacial OH on a water-covered α-Al2O3(0001) surface: a non-equilibrium ab initio molecular dynamics study

Vibrational relaxation of adsorbates is a sensitive tool to probe energy transfer at gas/solid and liquid/solid interfaces. The most direct way to study relaxation dynamics uses time-resolved spectroscopy. Here we report on a non-equilibrium ab initio molecular dynamics (NE-AIMD) methodology to model vibrational relaxation of OH vibrations on a hydroxylated, water-covered α-Al₂O₃(0001) surface. In our NE-AIMD approach, after exciting selected O-H bonds their coupling to surface phonons and to the water adlayer is analyzed in detail, by following both the energy flow in time, as well as the time-evolution of Vibrational Density of States (VDOS) curves. The latter are obtained from Time-dependent Correlation Functions (TCFs) and serve as prototypical, generic representatives of time-resolved vibrational spectra. As most important results, (i) we find a few-picosecond lifetime of the excited modes and (ii) identify both hydrogen-bonded aluminols and water molecules in the adsorbed water layer as main dissipative channels, while the direct coupling to Al₂O₃ surface phonons is of minor importance on the timescales of interest. Our NE-AIMD/TCF methodology is powerful for complex adsorbate systems, in principle even reacting ones, and opens a way towards time-resolved vibrational spectroscopy.

Vibrational Energy Relaxation of Deuterium Fluoride in d-Dichloromethane: Insights from Different Potentials

Vibrationally excited deuterium fluoride (DF) formed by fluorine atom reaction with a solvent was found (Science, 2015, 347, 530) to relax rapidly (less than 10 ps) in acetonitrile-d₃ (CD₃CN) and dichloromethane-d₂ (CD₂Cl₂). However, insights into how CD₂Cl₂ facilitates this energy relaxation have so far been lacking, given the weak interaction between DF and a single CD₂Cl₂. In this work, we report the results of reactive simulations with a two-state reactive empirical valence bond (EVB) potential to study the energy deposited into nascent DF after transition-state passage and of nonequilibrium molecular dynamics simulations using multiple different potential energy functions to model the relaxation dynamics. For these second simulations, we used the standard Merck molecular force field (MMFF) potential, an MMFF-based covalent-ionic empirical valence bond (EVB) potential (EVB^CI), a newly developed potential [referred to as MMFF(rDF)] which extends upon the MMFF potential by making the DF/CD₂Cl₂ interaction depend on the value of the D-F bond stretching coordinate and by taking the anisotropic charge distribution of the solvent molecules into account, the polarizable atomic multipole optimized energetics for biomolecular applications (AMOEBA) potential, and the quantum mechanics/molecular mechanics (QM/MM) potential. The relaxation is revealed to be highly sensitive to the potential used. Neither standard MMFF nor EVB^CI reproduces the experimentally observed rapid relaxation dynamics, and they also fail to provide a good description of the interaction potential between DF and CD₂Cl₂ as calculated using CCSD(T)-F12. This is attributed to the use of a point-charge model for the solute and to failing to model the anisotropic electrostatic properties of CD₂Cl₂. The MMFF(rDF), AMOEBA, and QM/MM potentials all reproduce the CCSD(T)-F12 two-body DF---CD₂Cl₂ interaction potential rather well but only with the QM/MM approach is fast vibrational relaxation obtained (lifetimes of ∼288, ∼186, and ∼8 ps, respectively), which we attribute to differences in the solute-solvent local structure. With QM/MM, a unique "many-body" interaction pattern in which DF is in close contact with two solvent Cl atoms and more than three solvent D atoms is found, but this structure is not seen with other potentials. The QM/MM dynamics also display enhanced solute-solvent interactions with vibrationally excited DF that induce a DF band redshift and hence a resonant overlap with solvent C-D modes, which facilitate the intermolecular energy transfer. Our work also suggests that potentials used to model energy relaxation need to capture the fine structure of solute-solvent interactions and not just the two-body part.

Ab initio Modeling of the Vibrational Sum-Frequency Generation Spectrum of Interfacial Water

Understanding the structural and dynamical features of interfacial water is of greatest interest in physics, chemistry, biology, and materials science. Vibrational sum-frequency generation (SFG) spectroscopy, which is sensitive to the molecular orientation and dynamics on the surfaces or at the interfaces, allows one to study a wide variety of interfacial systems. The structural and dynamical features of interfacial water at the air/water interface have been extensively investigated by SFG spectroscopy. However, the interpretations of the spectroscopic features have been under intense debate. Here, we report a simulated SFG spectrum of the air/water interface based on ab initio molecular dynamics simulations, which covers the OH stretching, bending, and libration modes of interfacial water. Quantitative agreement between our present simulations and the most recent experimental studies ensures that ab initio simulations predict unbiased structural features and electrical properties of interfacial systems. By utilizing the kinetic energy spectral density (KESD) analysis to decompose the simulated spectra, the spectroscopic features can then be assigned to specific hydrogen-bonding configurations of interfacial water molecules.

Hydrogen bonding and vibrational energy relaxation of interfacial water: A full DFT molecular dynamics simulation

The air-water interface has been a subject of extensive theoretical and experimental studies due to its ubiquity in nature and its importance as a model system for aqueous hydrophobic interfaces. We report on the structure and vibrational energy transfer dynamics of this interfacial water system studied with equilibrium and non-equilibrium molecular dynamics simulations employing a density functional theory -based description of the system and the kinetic energy spectral density analysis. The interfacial water molecules are found to make fewer and weaker hydrogen (H)-bonds on average compared to those in the bulk. We also find that (i) the H-bonded OH groups conjugate to the free OH exhibit rather low vibrational frequencies (3000-3500 cm^-1); (ii) the presence of a significant fraction (>10%) of free and randomly oriented water molecules at the interface ("labile water"), neither of whose OH groups are strong H-bond donors; (iii) the inertial rotation of free OH groups, especially from the labile water, contribute to the population decay of excited free OH groups with comparable rate and magnitude as intramolecular energy transfer between the OH groups. These results suggest that the labile water, which might not be easily detectable by the conventional vibrational sum frequency generation method, plays an important role in the surface water dynamics.

The dynamics of peptide-water interactions in dialanine: An ultrafast amide I 2D IR and computational spectroscopy study

We present a joint experimental and computational study of the dynamic interactions of dialanine (Ala-Ala) with water, comparing the results of ultrafast 2D IR and infrared transient absorption spectroscopy of its amide I vibration with spectra modeled from molecular dynamics (MD) simulations. The experimental data are analyzed to describe vibrational frequency fluctuations, vibrational energy relaxation, and chemical exchange processes. The origin of these processes in the same underlying fluctuating forces allows a common description in terms of the fluctuations and conformational dynamics of the peptide and associated solvent. By comparing computational spectroscopy from MD simulations with multiple force fields and water models, we describe how the dynamics of water hydrogen bond fluctuations and switching processes act as a source of friction that governs the dephasing and vibrational relaxation, and provide a description of coupled water and peptide motions that give rise to spectroscopic exchange processes.

Vibrational relaxation and energy migration of N-methylacetamide in water: the role of nonbonded interactions

Nonequilibrium molecular dynamics (MD) simulations together with physics-based force fields are used to follow energy flow between vibrationally excited N-methylacetamide (NMA) and water. The simulations are carried out with a previously validated force field for NMA, based on a multipolar representation of the electrostatics, and with a new fluctuating point charge model. For the water solvent, a flexible and a rigid model was employed to distinguish between the role of inter- and intramolecular degrees of freedom. On a 10 ps time scale about 90% of the available energy goes into the solvent. The remaining energy resides within internal NMA-degrees of freedom from where energy flow takes place on longer time scales. The total amount of energy transferred to the solvent on the 10 ps time scale does not depend on whether the water molecules are rigid or flexible during the simulations. Vibrational energy relaxation time scales include two regimes: one on the several 100 fs time scale and a longer one, ranging from 6 to 10 ps. This longer time scale agrees with previous simulations but overestimates the experimentally determined relaxation time by a factor of 2, which can be explained by the classical treatment of the vibrations. Including a previously determined quantum correction factor brings the long time scale into quite favorable agreement with experiment. Coupling to the bending vibration of the water molecules in H-bonding contact with the excited C═O chromophore is substantial. The equilibrium and nonequilibrium distribution of the bending angles of the water molecules in contact with the local oscillator are non-Gaussian, and one approaches the other on the subpicosecond time scale. Analysis of the water velocity distribution suggests that the C═O vibrational energy relaxes into the solvent water shells in an impulsive fashion on a picosecond time scale.

Simultaneous spectral and temporal analyses of kinetic energies in nonequilibrium systems: theory and application to vibrational relaxation of O-D stretch mode of HOD in water

A time series of kinetic energies (KE) from classical molecular dynamics (MD) simulation contains fundamental information on system dynamics. It can also be analyzed in the frequency domain through Fourier transformation (FT) of velocity correlation functions, providing energy content of different spectral regions. By limiting the FT time span, we have previously shown that spectral resolution of KE evolution is possible in the nonequilibrium situations [Jeon and Cho, J. Chem. Phys. 2011, 135, 214504]. In this paper, we refine the method by employing the concept of instantaneous power spectra, extending it to reflect an instantaneous time-correlation of velocities with those in the future as well as with those in the past, and present a new method to obtain the instantaneous spectral density of KE (iKESD). This approach enables the simultaneous spectral and temporal resolution of KE with unlimited time precision. We discuss the formal and novel properties of the new iKESD approaches and how to optimize computational methods and determine parameters for practical applications. The method is specifically applied to the nonequilibrium MD simulation of vibrational relaxation of the OD stretch mode in a hydrated HOD molecule by employing a hybrid quantum mechanical/molecular mechanical (QM/MM) potential. We directly compare the computational results with the OD band population relaxation time profiles extracted from the IR pump-probe measurements for 5% HOD in water. The calculated iKESD yields the OD bond relaxation time scale ∼30% larger than the experimental value, and this decay is largely frequency-independent if the classical anharmonicity is accounted for. From the integrated iKESD over intra- and intermolecular bands, the major energy transfer pathways were found to involve the HOD bending mode in the subps range, then the internal modes of the solvent until 5 ps after excitation, and eventually the solvent intermolecular modes. Also, strong hydrogen-bonding of HOD is found to significantly hinder the initial intramolecular energy transfer process.

Hydration Effect on Amide I Infrared Bands in Water: An Interpretation Based on an Interaction Energy Decomposition Scheme

The sensitivity of some infrared bands to the local environment can be exploited to shed light on the structure and the dynamics of biological systems. In particular, the amide I band, which is specifically related to vibrations within the peptide bonds, can give information on the ternary structure of proteins, and can be used as a probe of energy transfer. In this work, we propose a model to quantitatively interpret the frequency shift on the amide I band of a model peptide induced by the formation of hydrogen bonds in the first solvation shell. This method allows us to analyze to what extent the electrostatic interaction, electronic polarization and charge transfer affect the position of the amide I band. The impact of the anharmoniticy of the pontential energy surface on the hydration induced shift is elucidated as well.