A mixed quantum-classical molecular dynamics study of the hydroxyl stretch in methanol/carbon tetrachloride mixtures: equilibrium hydrogen-bond structure and dynamics at the ground state and the infrared absorption spectrum

Vibrational Spectroscopic Map, Vibrational Spectroscopy, and Intermolecular Interaction

Vibrational spectroscopy is an essential tool in chemical analyses, biological assays, and studies of functional materials. Over the past decade, various coherent nonlinear vibrational spectroscopic techniques have been developed and enabled researchers to study time-correlations of the fluctuating frequencies that are directly related to solute-solvent dynamics, dynamical changes in molecular conformations and local electrostatic environments, chemical and biochemical reactions, protein structural dynamics and functions, characteristic processes of functional materials, and so on. In order to gain incisive and quantitative information on the local electrostatic environment, molecular conformation, protein structure and interprotein contacts, ligand binding kinetics, and electric and optical properties of functional materials, a variety of vibrational probes have been developed and site-specifically incorporated into molecular, biological, and material systems for time-resolved vibrational spectroscopic investigation. However, still, an all-encompassing theory that describes the vibrational solvatochromism, electrochromism, and dynamic fluctuation of vibrational frequencies has not been completely established mainly due to the intrinsic complexity of intermolecular interactions in condensed phases. In particular, the amount of data obtained from the linear and nonlinear vibrational spectroscopic experiments has been rapidly increasing, but the lack of a quantitative method to interpret these measurements has been one major obstacle in broadening the applications of these methods. Among various theoretical models, one of the most successful approaches is a semiempirical model generally referred to as the vibrational spectroscopic map that is based on a rigorous theory of intermolecular interactions. Recently, genetic algorithm, neural network, and machine learning approaches have been applied to the development of vibrational solvatochromism theory. In this review, we provide comprehensive descriptions of the theoretical foundation and various examples showing its extraordinary successes in the interpretations of experimental observations. In addition, a brief introduction to a newly created repository Web site (http://frequencymap.org) for vibrational spectroscopic maps is presented. We anticipate that a combination of the vibrational frequency map approach and state-of-the-art multidimensional vibrational spectroscopy will be one of the most fruitful ways to study the structure and dynamics of chemical, biological, and functional molecular systems in the future.

Quantitative Analysis of Infrared Spectra of Binary Alcohol + Cyclohexane Solutions with Quantum Chemical Calculations

Hydrogen bonding has profound effects on the behavior of molecules. Fourier-transform infrared spectroscopy is the technique most commonly used to qualitatively identify hydrogen-bonding moieties present in a chemical sample. However, quantitative analysis of infrared (IR) spectra is nontrivial for the hydroxyl stretching region where hydrogen bonding is most prominently expressed in organic alcohols and water. Specifically, the breadth and extreme overlap of the O-H stretching bands, and the order of magnitude variability of their IR attenuation coefficients complicates the analysis. In the present work, sequential molecular dynamics simulations and quantum mechanical calculations are used to develop a function to relate the integrated IR attenuation coefficient to the vibrational frequencies of hydroxyl bands across the O-H stretching region. This relationship is then used as a guide to develop an attenuation coefficient scaling function to quantitatively determine concentrations of alcohols in a hydrocarbon solution from experimental IR spectra by integration across the entire hydroxyl frequency range.

Nonadiabatic Dynamics via the Symmetrical Quasi-Classical Method in the Presence of Anharmonicity

The symmetrical quasi-classical (SQC) method recently proposed by Miller and Cotton allows one to simulate nonadiabatic dynamics based on an algorithm with classical-like scaling with respect to system size. This is made possible by casting the electronic degrees of freedom in terms of mapping variables that can be propagated in a classical-like manner. While SQC was shown to be rather accurate when applied to benchmark models with harmonic electronic potential energy surfaces, it was also found to become inaccurate and to suffer numerical instabilities when applied to anharmonic systems. In this paper, we propose an extended SQC (E-SQC) methodology for overcoming those discrepancies by describing the anharmonic nuclear modes, which are coupled to the electronic degrees of freedom, in terms of classical-like mapping variables. The accuracy of E-SQC relative to standard SQC is demonstrated on benchmark models with quartic and Morse potential energy surfaces.

Hydrogen bond dynamics in bulk alcohols

Hydrogen-bonded liquids play a significant role in numerous chemical and biological phenomena. In the past decade, impressive developments in multidimensional vibrational spectroscopy and combined molecular dynamics-quantum mechanical simulation have established many intriguing features of hydrogen bond dynamics in one of the fundamental solvents in nature, water. The next class of a hydrogen-bonded liquid--alcohols--has attracted much less attention. This is surprising given such important differences between water and alcohols as the imbalance between the number of hydrogen bonds, each molecule can accept (two) and donate (one) and the very presence of the hydrophobic group in alcohols. Here, we use polarization-resolved pump-probe and 2D infrared spectroscopy supported by extensive theoretical modeling to investigate hydrogen bond dynamics in methanol, ethanol, and isopropanol employing the OH stretching mode as a reporter. The sub-ps dynamics in alcohols are similar to those in water as they are determined by similar librational and hydrogen-bond stretch motions. However, lower density of hydrogen bond acceptors and donors in alcohols leads to the appearance of slow diffusion-controlled hydrogen bond exchange dynamics, which are essentially absent in water. We anticipate that the findings herein would have a potential impact on fundamental chemistry and biology as many processes in nature involve the interplay of hydrophobic and hydrophilic groups.

Histidine Orientation Modulates the Structure and Dynamics of a de Novo Metalloenzyme Active Site

The ultrafast dynamics of a de novo metalloenzyme active site is monitored using two-dimensional infrared spectroscopy. The homotrimer of parallel, coiled coil α-helices contains a His3-Cu(I) metal site where CO is bound and serves as a vibrational probe of the hydrophobic interior of the self-assembled complex. The ultrafast spectral dynamics of Cu-CO reveals unprecedented ultrafast (2 ps) nonequilibrium structural rearrangements launched by vibrational excitation of CO. This initial rapid phase is followed by much slower ∼40 ps vibrational relaxation typical of metal-CO vibrations in natural proteins. To identify the hidden coupled coordinate, small molecule analogues and the full peptide were studied by QM and QM/MM calculations, respectively. The calculations show that variation of the histidines' dihedral angles in coordinating Cu controls the coupling between the CO stretch and the Cu-C-O bending coordinates. Analysis of different optimized structures with significantly different electrostatic field magnitudes at the CO ligand site indicates that the origin of the stretch-bend coupling is not directly due to through-space electrostatics. Instead, the large, ∼3.6 D dipole moments of the histidine side chains effectively transduce the electrostatic environment to the local metal coordination orientation. The sensitivity of the first coordination sphere to the protein electrostatics and its role in altering the potential energy surface of the bound ligands suggests that long-range electrostatics can be leveraged to fine-tune function through enzyme design.

Mixed quantum-classical molecular dynamics study of the hydroxyl stretch in methanol/carbon-tetrachloride mixtures II: excited state hydrogen bonding structure and dynamics, infrared emission spectrum, and excited state lifetime

We present a mixed quantum-classical molecular dynamics study of the hydrogen-bonding structure and dynamics of a vibrationally excited hydroxyl stretch in methanol/carbon-tetrachloride mixtures. The adiabatic Hamiltonian of the quantum-mechanical hydroxyl is diagonalized on-the-fly to obtain the ground and first-excited adiabatic energy levels and wave functions which depend parametrically on the instantaneous configuration of the classical degrees of freedom. The dynamics of the classical degrees of freedom are determined by Hellmann-Feynman forces obtained by taking the expectation value of the force with respect to the ground or excited vibrational wave functions. Polarizable force fields are used which were previously shown to reproduce the experimental infrared absorption spectrum rather well, for different isotopomers and over a wide composition range [Kwac, K.; Geva, E. J. Phys. Chem. B 2011, 115, 9184]. We show that the agreement of the absorption spectra with experiment can be further improved by accounting for the dependence of the dipole moment derivatives on the configuration of the classical degrees of freedom. We find that the propensity of a methanol molecule to form hydrogen bonds increases upon photoexcitation of its hydroxyl stretch, thereby leading to a sizable red-shift of the corresponding emission spectrum relative to the absorption spectrum. Treating the relaxation from the first excited to the ground state as a nonadiabatic process, and calculating its rate within the framework of Fermi's golden rule and the harmonic-Schofield quantum correction factor, we were able to predict a lifetime which is of the same order of magnitude as the experimental value. The experimental dependence of the lifetime on the transition frequency is also reproduced. Nonlinear mapping relations between the hydroxyl transition frequency and bond length in the excited state and the electric field along the hydroxyl bond axis are established. These mapping relations make it possible to reduce the computational cost of the mixed quantum-classical treatment to that of a fully classical treatment.

Time-averaging approximation in the interaction picture: absorption line shapes for coupled chromophores with application to liquid water

The time-averaging approximation (TAA), originally developed to calculate vibrational line shapes for coupled chromophores using mixed quantum/classical methods, is reformulated. In the original version of the theory, time averaging was performed for the full one-exciton Hamiltonian, while herein the time averaging is performed on the coupling (off-diagonal) Hamiltonian in the interaction picture. As a result, the influence of the dynamic fluctuations of the transition energies is more accurately described. We compare numerical results of the two versions of the TAA with numerically exact results for the vibrational absorption line shape of the OH stretching modes in neat water. It is shown that the TAA in the interaction picture yields theoretical line shapes that are in better agreement with exact results.