Chirality-selective vibrational sum frequency generation (chiral SFG) spectroscopy has emerged as a powerful technique for the study of biomolecular hydration water due to its sensitivity to the induced chirality of the first hydration shell. Thus far, water O–H vibrational bands in phase-resolved heterodyne chiral SFG spectra have been fit using one Lorentzian function per vibrational band, and the resulting fit has been used to infer the underlying frequency distribution. Here, we show that this approach may not correctly reveal the structure and dynamics of hydration water. Our analysis illustrates that the chiral SFG responses of symmetric and asymmetric O–H stretch modes of water have opposite phase and equal magnitude and are separated in energy by intramolecular vibrational coupling and a heterogeneous environment. The sum of the symmetric and asymmetric responses implies that an O–H stretch in a heterodyne chiral SFG spectrum should appear as two peaks with opposite phase and equal amplitude. Using pairs of Lorentzian functions to fit water O–H stretch vibrational bands, we improve spectral fitting of previously acquired experimental spectra of model β-sheet proteins and reduce the number of free parameters. The fitting allows us to estimate the vibrational frequency distribution and thus reveals the molecular interactions of water in hydration shells of biomolecules directly from chiral SFG spectra.
more »
« less
Modeling Anharmonic Effects in the Vibrational Spectra of High-Frequency Modes
High-resolution vibrational spectra of C–H, O–H, and N–H stretches depend on both molecular conformation and environment as well as provide a window into the frequencies of many other vibrational degrees of freedom as a result of mode mixing. We review current theoretical strategies that are being deployed to both aid and guide the analysis of the data that are encoded in these spectra. The goal is to enhance the power of vibrational spectroscopy as a tool for probing conformational preferences, hydrogen bonding effects away from equilibrium, and energy flow pathways. Recent years have seen an explosion of new methods and strategies for solving the nuclear Schrödinger equation. Rather than attempt a comprehensive review, this work highlights specific molecular systems that we have chosen as representing bonding motifs that are important to chemistry and biology. We focus on the choices theoretical chemists make regarding the level of electronic structure theory, the representation of the potential energy surface, the selection of coordinates, preferences in basis sets, and methods of solution.
more »
« less
- Award ID(s):
- 1900095
- PAR ID:
- 10451954
- Date Published:
- Journal Name:
- Annual Review of Physical Chemistry
- Volume:
- 74
- Issue:
- 1
- ISSN:
- 0066-426X
- Page Range / eLocation ID:
- 219 to 244
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Flavins have emerged as central to electron bifurcation, signaling, and countless enzymatic reactions. In bifurcation, two electrons acquired as a pair are separated in coupled transfers wherein the energy of both is concentrated on one of the two. This enables organisms to drive demanding reactions based on abundant low‐grade chemical fuel. To enable incorporation of this and other flavin capabilities into designed materials and devices, it is essential to understand fundamental principles of flavin electronic structure that make flavins so reactive and tunable by interactions with protein. Emerging computational tools can now replicate spectra of flavins and are gaining capacity to explain reactivity at atomistic resolution, based on electronic structures. Such fundamental understanding can moreover be transferrable to other chemical systems. A variety of computational innovations have been critical in reproducing experimental properties of flavins including their electronic spectra, vibrational signatures, and nuclear magnetic resonance (NMR) chemical shifts. A computational toolbox for understanding flavin reactivity moreover must be able to treat all five oxidation and protonation states, in addition to excited states that participate in flavoprotein's light‐driven reactions. Therefore, we compare emerging hybrid strategies and their successes in replicating effects of hydrogen bonding, the surrounding dielectric, and local electrostatics. These contribute to the protein's ability to modulate flavin reactivity, so we conclude with a survey of methods for incorporating the effects of the protein residues explicitly, as well as local dynamics. Computation is poised to elucidate the factors that affect a bound flavin's ability to mediate stunningly diverse reactions, and make life possible. This article is categorized under:Structure and Mechanism > Computational Biochemistry and BiophysicsElectronic Structure Theory > Combined QM/MM MethodsTheoretical and Physical Chemistry > Spectroscopymore » « less
-
DL_POLY Quantum 2.1 is introduced here as a highly modular, sustainable, and scalable general-purpose molecular dynamics (MD) simulation software for large-scale long-time MD simulations of condensed phase and interfacial systems with the essential nuclear quantum effects (NQEs) included. The new release improves upon version 2.0 through the introduction of several emerging real-time path integral (PI) methods, including fast centroid molecular dynamics (f-CMD) and fast quasi-CMD (f-QCMD) methods, as well as our recently introduced hybrid CMD (h-CMD) method for the accurate and efficient simulation of vibrational infrared spectra. Several test cases, including liquid bulk water at 300 K and ice Ih at 150 K, are used to showcase the performance of different implemented PI methods in simulating the infrared spectra at both ambient conditions and low temperatures where NQEs become more apparent. Additionally, using different salt-in-water (i.e., dilute) and water-in-salt (i.e., concentrated) lithium bis(trifluoromethanesulfonyl)imide (Li-TFSI) aqueous electrolyte solutions, we demonstrate the applicability of our recently introduced h-CMD method implemented in DL_POLY Quantum 2.1 for the large scale simulation of infrared (IR) spectra of complex heterogeneous systems. We show that h-CMD can overcome the curvature problem of CMD and the artificial broadening of T-RPMD for the accurate simulation of the vibrational spectra of complex, heterogeneous systems with NQEs included.more » « less
-
Abstract The ability to examine the vibrational spectra of liquids with nanometer spatial resolution will greatly expand the potential to study liquids and liquid interfaces. In fact, the fundamental properties of water, including complexities in its phase diagram, electrochemistry, and bonding due to nanoscale confinement are current research topics. For any liquid, direct investigation of ordered liquid structures, interfacial double layers, and adsorbed species at liquid–solid interfaces are of interest. Here, a novel way of characterizing the vibrational properties of liquid water with high spatial resolution using transmission electron microscopy is reported. By encapsulating water between two sheets of boron nitride, the ability to capture vibrational spectra to quantify the structure of the liquid, its interaction with the liquid‐cell surfaces, and the ability to identify isotopes including H2O and D2O using electron energy‐loss spectroscopy is demonstrated. The electron microscope used here, equipped with a high‐energy‐resolution monochromator, is able to record vibrational spectra of liquids and molecules and is sensitive to surface and bulk morphological properties both at the nano‐ and micrometer scales. These results represent an important milestone for liquid and isotope‐labeled materials characterization with high spatial resolution, combining nanoscale imaging with vibrational spectroscopy.more » « less
-
null (Ed.)The vibrational predissociation of NeBr2 has been studied using a variety of theoretical and experimental methods, producing a large number of results. It is therefore a useful system for comparing different theoretical methods. Here, we apply the trajectory surface hopping (TSH) method that consists of propagating the dynamics of the system on a potential energy surface (PES) corresponding to quantum molecular vibrational states with possibility of hopping towards other surfaces until the van der Waals bond dissociates. This allows quantum vibrational effects to be added to a classical dynamics approach. We have also incorporated the kinetic mechanism for a better compression of the evolution of the complex. The novelty of this work is that it allows us to incorporate all the surfaces for (v=16,17,…,29) into the dynamics of the system. The calculated lifetimes are similar to those previously reported experimentally and theoretically. The rotational distribution, the rotational energy and jmax are in agreement with other works, providing new information for this complex.more » « less
An official website of the United States government

