Theory of coherent two-dimensional vibrational spectroscopy

Artificial-Intelligence-Enhanced On-the-Fly Simulation of Nonlinear Time-Resolved Spectra

Time-resolved spectroscopy is an important tool for unraveling the minute details of structural changes in molecules of biological and technological significance. The nonlinear femtosecond signals detected for such systems must be interpreted, but it is a challenging task for which theoretical simulations are often indispensable. Accurate simulations of transient absorption or two-dimensional electronic spectra are, however, computationally very expensive, prohibiting the wider adoption of existing first-principles methods. Here, we report an artificial-intelligence-enhanced protocol to drastically reduce the computational cost of simulating nonlinear time-resolved electronic spectra, which makes such simulations affordable for polyatomic molecules of increasing size. The protocol is based on the doorway-window approach for the on-the-fly surface-hopping simulations. We show its applicability for the prototypical molecule of pyrazine for which it produces spectra with high precision with respect to ab initio reference while cutting the computational cost by at least 95% compared to pure first-principles simulations.

Conformational Dynamics in Proteins: Entangled Slow Fluctuations and Nonequilibrium Reaction Events

Proteins exhibit conformational fluctuations and changes over various time scales, ranging from rapid picosecond-scale local atomic motions to slower microsecond-scale global conformational transformations. In the presence of these intricate fluctuations, chemical reactions occur and functions emerge. These conformational fluctuations of proteins are not merely stochastic random motions but possess distinct spatiotemporal characteristics. Moreover, chemical reactions do not always proceed along a single reaction coordinate in a quasi-equilibrium manner. Therefore, it is essential to understand spatiotemporal conformational fluctuations of proteins and the conformational change processes associated with reactions. In this Perspective, we shed light on the complex dynamics of proteins and their role in enzyme catalysis by presenting recent results regarding dynamic couplings and disorder in the conformational dynamics of proteins and rare but rapid enzymatic reaction events obtained from molecular dynamics simulations.

Quantum Biology and the Potential Role of Entanglement and Tunneling in Non-Targeted Effects of Ionizing Radiation: A Review and Proposed Model

It is well established that cells, tissues, and organisms exposed to low doses of ionizing radiation can induce effects in non-irradiated neighbors (non-targeted effects or NTE), but the mechanisms remain unclear. This is especially true of the initial steps leading to the release of signaling molecules contained in exosomes. Voltage-gated ion channels, photon emissions, and calcium fluxes are all involved but the precise sequence of events is not yet known. We identified what may be a quantum entanglement type of effect and this prompted us to consider whether aspects of quantum biology such as tunneling and entanglement may underlie the initial events leading to NTE. We review the field where it may be relevant to ionizing radiation processes. These include NTE, low-dose hyper-radiosensitivity, hormesis, and the adaptive response. Finally, we present a possible quantum biological-based model for NTE.

Experimental Access to Mode-Specific Coupling between Quantum Molecular Vibrations and Classical Bath Modes

The interaction of quantum-mechanical systems with a fluctuating thermal environment (bath) is fundamental to molecular mechanics and energy transport/dissipation. Its complete picture requires mode-specific measurements of this interaction and an understanding of its nature. Here, we present a combined experimental and theoretical study providing detailed insights into the coupling between a high-frequency vibrational two-level system and thermally excited terahertz modes. Experimentally, two-dimensional terahertz-infrared-visible spectroscopy reports directly on the coupling between quantum oscillators represented by CH3 streching vibrations in liquid dimethyl sulfoxide and distinct low-frequency modes. Theoretically, we present a mixed quantum-classical formalism of the sample response to enable the simultaneous quantum description of high-frequency oscillators and a classical description of the bath. We derive the strength and nature of interaction and find different coupling between CH3 stretch and low-frequency modes. This general approach enables quantitative and mode-specific analysis of coupled quantum and classical dynamics in complex chemical systems.

Molecular Dynamics of Artificially Pair-Decoupled Systems: An Accurate Tool for Investigating the Importance of Intramolecular Couplings

We propose a numerical technique to accurately simulate the vibrations of organic molecules in the gas phase, when pairs of atoms (or, in general, groups of degrees of freedom) are artificially decoupled, so that their motion is instantaneously decorrelated. The numerical technique we have developed is a symplectic integration algorithm that never requires computation of the force but requires estimates of the Hessian matrix. The theory we present to support our technique postulates a pair-decoupling Hamiltonian function, which parametrically depends on a decoupling coefficient α ∈ [0, 1]. The closer α is to 0, the more decoupled the selected atoms. We test the correctness of our numerical method on small molecular systems, and we apply it to study the vibrational spectroscopic features of salicylic acid at the Density Functional Theory ab initio level on a fitted potential. Our pair-decoupled simulations of salicylic acid show that decoupling hydrogen-bonded atoms do not significantly influence the frequencies of stretching modes, but enhance enormously the out-of-plane wagging and twisting motions of the hydroxyl and carboxyl groups to the point that the carboxyl and hydroxyl groups may overcome high potential energy barriers and change the salicylic acid conformation after a short simulation time. In addition, we found that the acidity of salicylic acid is more influenced by the dynamical couplings of the proton of the carboxylic group with the carbon ring than with the hydroxyl group.

Mixed Quantum-Classical Liouville Equation Treatment of Electronic Spectroscopy of Condensed Systems: Harmonic and Anharmonic Electron-Phonon Coupling

This Review integrates the use of electronic optical response function theory and the mixed quantum-classical (MQC) Liouville equation (MQCLE), thereby leading to electronic spectroscopy in MQC media. It further sheds light on the applicability, utility, and efficiency of the mixed quantum-classical dynamics (MQCD) formalism, which starts off with the MQCLE, in probing spectroscopy and dynamics of condensed systems, whereby quantum mechanics and classical mechanics are combined systematically. The author has been exploring and implementing MQCD to investigate electron-phonon coupling effects on electronic dephasing in harmonic and anharmonic systems by calculating linear and nonlinear optical transition analytically and numerically dipole moment time correlation functions in an MQC environment, thereby presenting an in depth spectral profile analysis and their shape and symmetry. The distinctive capability of the MQC time correlation functions is that ergodicity and stationarity properties are inherently satisfied as part of the mixed quantum-classical dynamics (MQCD) framework, unlike classical correlation functions. While some research groups have applied MQCLE to calculate vibrational spectra to study hydrogen-bonded complexes in a MQC environment and other groups calculated Optical response function to probe electron transfer dynamics using the basis mapping technique, the approach, purpose, rigor, applications, and path to the end results reported herein are different. Finally, the same framework is employed to study dissipative systems in the MQC limit, whereby the zero-phonon line adopts the correct width and eliminates its asymmetry. While the full quantum mechanical model, like the multimode Brownian oscillator (MBO) model, yields the correct width and inaccurate shape in the low-temperature limit, the MQCD formalism seems to produce an accurate zero-phonon profile. Nonlinear optical signals are also reviewed in MQC media to show the applicability and utility of this approach. The vibronic optical response functions developed here will account for geometry change, frequency change, and anharmonicity upon electronic excitation to accurately probe electronic dephasing, electron-phonon coupling, shape, and symmetry of profiles and present differences and similarities to the MBO model on pure electronic dephasing. Frequency change and anharmonicity are vitally crucial for accurately assessing electron-phonon coupling upon electronic excitation. This is an additional unique result obtained by the author to further demonstrate the applicability and utility of this approach over other approximation schemes in probing electronic dephasing including that of the MBO model.

Ultrafast Dynamics Revealed with Time-Resolved Scanning Tunneling Microscopy: A Review

A scanning tunneling microscope (STM) capable of performing pump-probe spectroscopy integrates unmatched atomic-scale resolution with high temporal resolution. In recent years, the union of electronic, terahertz, or visible/near-infrared pulses with STM has contributed to our understanding of the atomic-scale processes that happen between milliseconds and attoseconds. This time-resolved STM (TR-STM) technique is evolving into an unparalleled approach for exploring the ultrafast nuclear, electronic, or spin dynamics of molecules, low-dimensional structures, and material surfaces. Here, we review the recent advancements in TR-STM; survey its application in measuring the dynamics of three distinct systems, nucleus, electron, and spin; and report the studies on these transient processes in a series of materials. Besides the discussion on state-of-the-art techniques, we also highlight several emerging research topics about the ultrafast processes in nanoscale objects where we anticipate that the TR-STM can help broaden our knowledge.

Ab initio simulation of peak evolutions and beating maps for electronic two-dimensional signals of a polyatomic chromophore

By employing the doorway-window (DW) on-the-fly simulation protocol, we performed ab initio simulations of peak evolutions and beating maps of electronic two-dimensional (2D) spectra of a polyatomic molecule in the gas phase. As the system under study, we chose pyrazine, which is a paradigmatic example of photodynamics dominated by conical intersections (CIs). From the technical perspective, we demonstrate that the DW protocol is a numerically efficient methodology suitable for simulations of 2D spectra for a wide range of excitation/detection frequencies and population times. From the information content perspective, we show that peak evolutions and beating maps not only reveal timescales of transitions through CIs but also pinpoint the most relevant coupling and tuning modes active at these CIs.

Two-dimensional infrared-Raman spectroscopy as a probe of water's tetrahedrality

Two-dimensional spectroscopic techniques combining terahertz (THz), infrared (IR), and visible pulses offer a wealth of information about coupling among vibrational modes in molecular liquids, thus providing a promising probe of their local structure. However, the capabilities of these spectroscopies are still largely unexplored due to experimental limitations and inherently weak nonlinear signals. Here, through a combination of equilibrium-nonequilibrium molecular dynamics (MD) and a tailored spectrum decomposition scheme, we identify a relationship between the tetrahedral order of liquid water and its two-dimensional IR-IR-Raman (IIR) spectrum. The structure-spectrum relationship can explain the temperature dependence of the spectral features corresponding to the anharmonic coupling between low-frequency intermolecular and high-frequency intramolecular vibrational modes of water. In light of these results, we propose new experiments and discuss the implications for the study of tetrahedrality of liquid water.

Equivalence of quantum and classical third order response for weakly anharmonic coupled oscillators

Two-dimensional (2D) infrared (IR) spectra are commonly interpreted using a quantum diagrammatic expansion that describes the changes to the density matrix of quantum systems in response to light-matter interactions. Although classical response functions (based on Newtonian dynamics) have shown promise in computational 2D IR modeling studies, a simple diagrammatic description has so far been lacking. Recently, we introduced a diagrammatic representation for the 2D IR response functions of a single, weakly anharmonic oscillator and showed that the classical and quantum 2D IR response functions for this system are identical. Here, we extend this result to systems with an arbitrary number of bilinearly coupled, weakly anharmonic oscillators. As in the single-oscillator case, quantum and classical response functions are found to be identical in the weakly anharmonic limit or, in experimental terms, when the anharmonicity is small relative to the optical linewidth. The final form of the weakly anharmonic response function is surprisingly simple and offers potential computational advantages for application to large, multi-oscillator systems.

Direct observation of conformations of a high-mobility n-type low-bandgap copolymer in solutions and solid films

The aggregation morphologies of conjugated polymers in solutions and solid films are important for their optoelectronic applications. Due to the amorphous state of the polymers, it remains a great challenge to determine their conformations in either liquids or solids. Herein, a ps/fs synchronized 2D IR technique is applied to investigate the molecular conformations of a high-mobility n-type low-bandgap copolymer, N2200, dissolved in CHCl₃ and CCl₄, and in solid films cast from both solutions by the vibrational cross-angle method. In CCl₄, the polymer forms more aggregates and folds more and the backbone dihedral angle of C-C(NDI)/C-S(Thiophene) of its average conformation is about 10° more distorted than that in CHCl₃ and the most stable conformation for a free molecule. Anti-intuitively, the solid films cast from both solutions have the same molecular conformation, and the conformation is similar to that of the polar CHCl₃ rather than the conformation of the less polar CCl₄. The results imply that the interaction between the polymer backbones is probably stronger than its interaction with CCl₄, which can naturally guide the rearrangement of polymer chains during the evaporation of solvent molecules. This work also implies that the balance and competition between the polymer/polymer interaction and the polymer/solvent interaction seem to be the dominant factors responsible for what morphology can form in a solid film cast from solution. It is not always true that different molecular conformations must exist in solid films grown from different solutions with different polarity or different extents of aggregates with different conformations.

Simulation of two-dimensional infrared Raman spectroscopy with application to proteins

Two-dimensional infrared Raman spectroscopy is a powerful technique for studying the structure and interaction in molecular and biological systems. Here, we present a new implementation of the simulation of the two-dimensional infrared Raman signals. The implementation builds on the numerical integration of the Schrödinger equation approach. It combines the prediction of dynamics from molecular dynamics with a map-based approach for obtaining Hamiltonian trajectories and response function calculations. The new implementation is tested on the amide-I region for two proteins, where one is dominated by α-helices and the other by β-sheets. We find that the predicted spectra agree well with experimental observations. We further find that the two-dimensional infrared Raman spectra at least of the studied proteins are much less sensitive to the laser polarization used compared to conventional two-dimensional infrared experiments. The present implementation and findings pave the way for future applications for the interpretation of two-dimensional infrared Raman spectra.

Discretized hierarchical equations of motion in mixed Liouville-Wigner space for two-dimensional vibrational spectroscopies of liquid water

A model of a bulk water system describing the vibrational motion of intramolecular and intermolecular modes is constructed, enabling analysis of its linear and nonlinear vibrational spectra as well as the energy transfer processes between the vibrational modes. The model is described as a system of four interacting anharmonic oscillators nonlinearly coupled to their respective heat baths. To perform a rigorous numerical investigation of the non-Markovian and nonperturbative quantum dissipative dynamics of the model, we derive discretized hierarchical equations of motion in mixed Liouville-Wigner space, with Lagrange-Hermite mesh discretization being employed in the Liouville space of the intramolecular modes and Lagrange-Hermite mesh discretization and Hermite discretization in the Wigner space of the intermolecular modes. One-dimensional infrared and Raman spectra and two-dimensional terahertz-infrared-visible and infrared-infrared-Raman spectra are computed as demonstrations of the quantum dissipative description provided by our model.

Infrared Diffusion-Ordered Spectroscopy Reveals Molecular Size and Structure

Inspired by ideas from NMR, we have developed Infrared Diffusion-Ordered Spectroscopy (IR-DOSY), which simultaneously characterizes molecular structure and size. We rely on the fact that the diffusion coefficient of a molecule is determined by its size through the Stokes-Einstein relation, and achieve sensitivity to the diffusion coefficient by creating a concentration gradient and tracking its equilibration in an IR-frequency resolved manner. Analogous to NMR-DOSY, a two-dimensional IR-DOSY spectrum has IR frequency along one axis and diffusion coefficient (or equivalently, size) along the other, so the chemical structure and the size of a compound are characterized simultaneously. In an IR-DOSY spectrum of a mixture, molecules with different sizes are nicely separated into distinct sets of IR peaks. Extending this idea to higher dimensions, we also perform 3D-IR-DOSY, in which we combine the conformation sensitivity of femtosecond multi-dimensional IR spectroscopy with size sensitivity.

Equation-of-Motion Methods for the Calculation of Femtosecond Time-Resolved 4-Wave-Mixing and N-Wave-Mixing Signals

Femtosecond nonlinear spectroscopy is the main tool for the time-resolved detection of photophysical and photochemical processes. Since most systems of chemical interest are rather complex, theoretical support is indispensable for the extraction of the intrinsic system dynamics from the detected spectroscopic responses. There exist two alternative theoretical formalisms for the calculation of spectroscopic signals, the nonlinear response-function (NRF) approach and the spectroscopic equation-of-motion (EOM) approach. In the NRF formalism, the system-field interaction is assumed to be sufficiently weak and is treated in lowest-order perturbation theory for each laser pulse interacting with the sample. The conceptual alternative to the NRF method is the extraction of the spectroscopic signals from the solutions of quantum mechanical, semiclassical, or quasiclassical EOMs which govern the time evolution of the material system interacting with the radiation field of the laser pulses. The NRF formalism and its applications to a broad range of material systems and spectroscopic signals have been comprehensively reviewed in the literature. This article provides a detailed review of the suite of EOM methods, including applications to 4-wave-mixing and N-wave-mixing signals detected with weak or strong fields. Under certain circumstances, the spectroscopic EOM methods may be more efficient than the NRF method for the computation of various nonlinear spectroscopic signals.

Quantum-chemical calculation of two-dimensional infrared spectra using localized-mode VSCF/VCI

Computational protocols for the simulation of two-dimensional infrared (2D IR) spectroscopy usually rely on vibrational exciton models which require an empirical parameterization. Here, we present an efficient quantum-chemical protocol for predicting static 2D IR spectra that does not require any empirical parameters. For the calculation of anharmonic vibrational energy levels and transition dipole moments, we employ the localized-mode vibrational self-consistent field (L-VSCF)/vibrational configuration interaction (L-VCI) approach previously established for (linear) anharmonic theoretical vibrational spectroscopy [P. T. Panek and C. R. Jacob, ChemPhysChem 15, 3365-3377 (2014)]. We demonstrate that with an efficient expansion of the potential energy surface using anharmonic one-mode potentials and harmonic two-mode potentials, 2D IR spectra of metal carbonyl complexes and dipeptides can be predicted reliably. We further show how the close connection between L-VCI and vibrational exciton models can be exploited to extract the parameters of such models from those calculations. This provides a novel route to the fully quantum-chemical parameterization of vibrational exciton models for predicting 2D IR spectra.

Molecular photothermal effects on time-resolved IR spectroscopy

Time-resolved IR pump-probe (IR-PP) and two-dimensional IR (2D-IR) spectroscopy are valuable techniques for studying various ultrafast chemical and biological processes in solutions. The time-dependent changes of nonlinear IR signals reflecting fast molecular processes such as vibrational energy transfer and chemical exchange provide invaluable information on the rates and mechanisms of solvation dynamics and structural transitions of multi-species vibrationally interacting molecular systems. However, due to the intrinsic difficulties in distinguishing the contributions of molecule-specific processes to the time-resolved IR signals from those resulting from local heating, it becomes challenging to interpret time-resolved IR-PP and 2D-IR spectra exhibiting transient growing-in spectral components and cross-peaks unambiguously. Here, theoretical considerations of various effects of vibrational coupling, energy transfer, chemical exchange, the generation of hot ground states, molecular photothermal process, and their combinations on the lineshapes and time-dependent intensities of IR-PP spectra and 2D-IR diagonal and cross-peaks are presented. We anticipate that the present work will help researchers using IR pump-probe and 2D-IR techniques to distinguish local heating-induced photothermal signals from genuine nonlinear IR signals.

2D spectroscopies from condensed phase dynamics: Accessing third-order response properties from equilibrium multi-time correlation functions

The third-order response lies at the heart of simulating and interpreting nonlinear spectroscopies ranging from two dimensional infrared (2D-IR) to 2D electronic (2D-ES), and 2D sum frequency generation (2D-SFG). The extra time and frequency dimensions in these spectroscopies provides access to rich information on the electronic and vibrational states present, the coupling between them, and the resulting rates at which they exchange energy that are obscured in linear spectroscopy, particularly for condensed phase systems that usually contain many overlapping features. While the exact quantum expression for the third-order response is well established it is incompatible with the methods that are practical for calculating the atomistic dynamics of large condensed phase systems. These methods, which include both classical mechanics and quantum dynamics methods that retain quantum statistical properties while obeying the symmetries of classical dynamics, such as LSC-IVR, Centroid Molecular Dynamics (CMD) and Ring Polymer Molecular Dynamics (RPMD) naturally provide short-time approximations to the multi-time symmetrized Kubo transformed correlation function. Here, we show how the third-order response can be formulated in terms of equilibrium symmetrized Kubo transformed correlation functions. We demonstrate the utility and accuracy of our approach by showing how it can be used to obtain the third-order response of a series of model systems using both classical dynamics and RPMD. In particular, we show that this approach captures features such as anharmonically induced vertical splittings and peak shifts while providing a physically transparent framework for understanding multidimensional spectroscopies.

Linear and Non-Linear Population Retrieval with Femtosecond Optical Pumping of Molecular Crystals for the Generalised Uniaxial and Biaxial Systems

Femtosecond optical measurements of photoexcitable molecular crystals carry ultrafast dynamics information with structural sensitivity. The creation and detection of transient dynamics depend on the optical parameters, as well as the explicit molecular structure, crystal symmetry, crystal orientation, polarisation of the photoexciting beam, and interaction geometry. In order to retrieve the linear and non-linear population transfer in photoexcited crystals, excitation theory is combined here with the calculation of birefringence decomposition and is shown for both the generalised uniaxial and biaxial systems. A computational tool was constructed based on this treatment to allow modelling of electric field decomposition, dipole projections, and non-linear excitation population levels. This is available open source and with a GUI for ease of use. Such work has applications in two areas of ultrafast science: multidimensional optical crystallography and femtosecond time-resolved X-ray crystallography.

Equilibrium-nonequilibrium ring-polymer molecular dynamics for nonlinear spectroscopy

Two-dimensional Raman and hybrid terahertz-Raman spectroscopic techniques provide invaluable insight into molecular structures and dynamics of condensed-phase systems. However, corroborating experimental results with theory is difficult due to the high computational cost of incorporating quantum-mechanical effects in the simulations. Here, we present the equilibrium-nonequilibrium ring-polymer molecular dynamics (RPMD), a practical computational method that can account for nuclear quantum effects on the two-time response function of nonlinear optical spectroscopy. Unlike a recently developed approach based on the double Kubo transformed (DKT) correlation function, our method is exact in the classical limit, where it reduces to the established equilibrium-nonequilibrium classical molecular dynamics method. Using benchmark model calculations, we demonstrate the advantages of the equilibrium-nonequilibrium RPMD over classical and DKT-based approaches. Importantly, its derivation, which is based on the nonequilibrium RPMD, obviates the need for identifying an appropriate Kubo transformed correlation function and paves the way for applying real-time path-integral techniques to multidimensional spectroscopy.

AIM: A Mapping Program for Infrared Spectroscopy of Proteins

Here, we present a new analysis program, AIM, that allows extracting the vibrational amide-I Hamiltonian using molecular dynamics trajectories for protein infrared spectroscopy modeling. The constructed Hamiltonians can be used as input for spectral calculations allowing the calculation of infrared absorption spectra, vibrational circular dichroism, and two-dimensional infrared spectra. These spectroscopies allow the study of the structure and dynamics of proteins. We will explain the essence of how AIM works and give examples of the information and spectra that can be obtained with the program using the Trypsin Inhibitor as an example. AIM is freely available from GitHub, and the package contains a demonstration allowing easy introduction to the use of the program.

Calculating Multidimensional Optical Spectra from Classical Trajectories

Multidimensional optical spectra are measured from the response of a material system to a sequence of laser pulses and have the capacity to elucidate specific molecular interactions and dynamics whose influences are absent or obscured in a conventional linear absorption spectrum. Interpretation of complex spectra is supported by theoretical modeling of the spectroscopic observable, requiring implementation of quantum dynamics for coupled electrons and nuclei. Performing numerically correct quantum dynamics in this context may pose computational challenges, particularly in the condensed phase. Semiclassical methods based on calculating classical trajectories offer a practical alternative. Here I review the recent application of some semiclassical, trajectory-based methods to nonlinear molecular vibrational and electronic spectra.

Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Simulation of Condensed-Phase Spectroscopy with Near-Term Digital Quantum Computers

Spectroscopy is an indispensable tool for understanding the structures and dynamics of molecular systems. However, computational modeling of spectroscopy is challenging due to the exponential scaling of computational complexity with system sizes unless drastic approximations are made. Quantum computers could potentially overcome these classically intractable computational tasks, but the existing approaches using quantum computers to simulate spectroscopy can only handle isolated and static molecules. In this work, we develop a workflow that combines multi-scale modeling and a time-dependent variational quantum algorithm to compute the linear spectroscopy of systems interacting with their condensed-phase environment via the relevant time correlation function. We demonstrate the feasibility of our approach by numerically simulating the UV-vis absorption spectra of organic semiconductors. We show that our dynamical approach captures several spectral features that are otherwise overlooked by static methods. Our method can be directly used for other linear condensed-phase spectroscopy and could potentially be extended to nonlinear multi-dimensional spectroscopy.

Computational spectroscopy of complex systems

Numerous linear and non-linear spectroscopic techniques have been developed to elucidate structural and functional information of complex systems ranging from natural systems, such as proteins and light-harvesting systems, to synthetic systems, such as solar cell materials and light-emitting diodes. The obtained experimental data can be challenging to interpret due to the complexity and potential overlapping spectral signatures. Therefore, computational spectroscopy plays a crucial role in the interpretation and understanding of spectral observables of complex systems. Computational modeling of various spectroscopic techniques has seen significant developments in the past decade, when it comes to the systems that can be addressed, the size and complexity of the sample types, the accuracy of the methods, and the spectroscopic techniques that can be addressed. In this Perspective, I will review the computational spectroscopy methods that have been developed and applied for infrared and visible spectroscopies in the condensed phase. I will discuss some of the questions that this has allowed answering. Finally, I will discuss current and future challenges and how these may be addressed.

Multimode two-dimensional vibronic spectroscopy. II. Simulating and extracting vibronic coupling parameters from polarization-selective spectra

Experimental demonstrations of polarization-selection two-dimensional Vibrational-Electronic (2D VE) and 2D Electronic-Vibrational (2D EV) spectroscopies aim to map the magnitudes and spatial orientations of coupled electronic and vibrational coordinates in complex systems. The realization of that goal depends on our ability to connect spectroscopic observables with molecular structural parameters. In this paper, we use a model Hamiltonian consisting of two anharmonically coupled vibrational modes in electronic ground and excited states with linear and bilinear vibronic coupling terms to simulate polarization-selective 2D EV and 2D VE spectra. We discuss the relationships between the linear vibronic coupling and two-dimensional Huang-Rhys parameters and between the bilinear vibronic coupling term and Duschinsky mixing. We develop a description of the vibronic transition dipoles and explore how the Hamiltonian parameters and non-Condon effects impact their amplitudes and orientations. Using simulated polarization-selective 2D EV and 2D VE spectra, we show how 2D peak positions, amplitudes, and anisotropy can be used to measure parameters of the vibronic Hamiltonian and non-Condon effects. This paper, along with the first in the series, provides the reader with a detailed description of reading, simulating, and analyzing multimode, polarization-selective 2D EV and 2D VE spectra with an emphasis on extracting vibronic coupling parameters from complex spectra.

Observation of Ultrafast Coherence Transfer and Degenerate States with Polarization-Controlled Two-Dimensional Electronic Spectroscopy

Optical spectroscopy is a powerful tool to interrogate quantum states of matter. We present simulation results for the cross-polarized two-dimensional electronic spectra of the light-harvesting system LH2 of purple bacteria. We identify a spectral feature on the diagonal, which we assign to ultrafast coherence transfer between degenerate states. The implication for the interpretation of previous experiments on different systems and the potential use of this feature are discussed. In particular, we foresee that this kind of feature will be useful for identifying mixed degenerate states and for identifying the origin of symmetry breaking disorder in systems like LH2. Furthermore, this may help identify both vibrational and electronic states in biological systems such as proteins and solid-state materials such as hybrid perovskites.

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.

Two-dimensional Raman spectroscopy of Lennard-Jones liquids via ring-polymer molecular dynamics

The simulation of multidimensional vibrational spectroscopy of condensed-phase systems including nuclear quantum effects is challenging since full quantum-mechanical calculations are still intractable for large systems comprising many degrees of freedom. Here, we apply the recently developed double Kubo transform (DKT) methodology in combination with ring-polymer molecular dynamics (RPMD) for evaluating multi-time correlation functions [K. A. Jung et al., J. Chem. Phys. 148, 244105 (2018)], providing a practical method for incorporating nuclear quantum effects in nonlinear spectroscopy of condensed-phase systems. We showcase the DKT approach in the simulation of the fifth-order two-dimensional (2D) Raman spectroscopy of Lennard-Jones liquids as a prototypical example, which involves nontrivial nonlinear spectroscopic observables of systems described by anharmonic potentials. Our results show that the DKT can faithfully reproduce the 2D Raman response of liquid xenon at high temperatures, where the system behaves classically. In contrast, liquid neon at low temperatures exhibits moderate but discernible nuclear quantum effects in the 2D Raman response compared to the responses obtained with classical molecular dynamics approaches. Thus, the DKT formalism in combination with RPMD simulations enables simulations of multidimensional optical spectroscopy of condensed-phase systems that partially account for nuclear quantum effects.

Tautomeric Effect of Histidine on β-Sheet Formation of Amyloid Beta 1-40: 2D-IR Simulations

Histidine state (protonated or δ or ε tautomer) has been considered the origin of abnormal misfolding and aggregation of β-amyloid (Aβ). Our previous studies reported that the δδδ isomer of Aβ (1-40) has a greater propensity for β-sheet conformation compared to other isomers. However, direct proof of the tautomeric effect has not been reported. In this context, we calculated histidine site-specific two-dimensional infrared spectroscopy of the δδδ, εεε, and πππ (all protonated histidine) systems within the framework of classical molecular dynamics simulations aiming at connecting our previous results with the current experimental observations. Our results showed that β-sheet formation is favored for the δδδ and πππ tautomers compared with the εεε tautomer, consistent with our previous studies. This result was further supported by contact map analyses and the strength of dipole coupling between the amide-I bonds of each residue. The two-dimensional infrared diagonal trace for each tautomer included three distinctive spectrally resolvable peaks near 1680, 1686, and 1693 cm^-1, as was also observed for histidine dipeptides. However, the peak positions at His6, His13, and His14 did not show a consensus trend with the histidine or protonation state but were instead affected by the presence of surrounding hydrogen bonds. Our study provides a deeper insight into the influence of tautomerism and protonation of histidine residues in Aβ (1-40) on amyloid misfolding and provides a connection between our previous simulations and experimental observations.