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 spectroscopic techniques provide 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, 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.
more »
« less
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.
more » « less- Award ID(s):
- 1665467
- PAR ID:
- 10365555
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 156
- Issue:
- 13
- ISSN:
- 0021-9606
- Page Range / eLocation ID:
- Article No. 131102
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Abstract Quantum entanglement has emerged as a great resource for spectroscopy and its importance in two-photon spectrum and microscopy has been demonstrated. Current studies focus on the two-photon absorption, whereas the Raman spectroscopy with quantum entanglement still remains elusive, with outstanding issues of temporal and spectral resolutions. Here we study the new capabilities provided by entangled photons in coherent Raman spectroscopy. An ultrafast frequency-resolved Raman spectroscopy with entangled photons is developed for condensed-phase molecules, to probe the electronic and vibrational coherences. Using quantum correlation between the photons, the signal shows the capability of both temporal and spectral resolutions not accessible by either classical pulses or the fields without entanglement. We develop a microscopic theory for this Raman spectroscopy, revealing the electronic coherence dynamics even at timescale of 50fs. This suggests new paradigms of optical signals and spectroscopy, with potential to push detection below standard quantum limit.
-
The structure–function relationships of biomolecules have captured the interest and imagination of the scientific community and general public since the field of structural biology emerged to enable the molecular understanding of life processes. Proteins that play numerous functional roles in cellular processes have remained in the forefront of research, inspiring new characterization techniques. In this review, we present key theoretical concepts and recent experimental strategies using femtosecond stimulated Raman spectroscopy (FSRS) to map the structural dynamics of proteins, highlighting the flexible chromophores on ultrafast timescales. In particular, wavelength-tunable FSRS exploits dynamic resonance conditions to track transient-species-dependent vibrational motions, enabling rational design to alter functions. Various ways of capturing excited-state chromophore structural snapshots in the time and/or frequency domains are discussed. Continuous development of experimental methodologies, synergistic correlation with theoretical modeling, and the expansion to other nonequilibrium, photoswitchable, and controllable protein systems will greatly advance the chemical, physical, and biological sciences.more » « less
-
With the objective of understanding microscopic principles governing thermal energy flow in nanojunctions, we study phononic heat transport through metal-molecule-metal junctions using classical molecular dynamics (MD) simulations. Considering a single-molecule gold-alkanedithiol-gold junction, we first focus on aspects of method development and compare two techniques for calculating thermal conductance: (i) The Reverse Nonequilibrium MD (RNEMD) method, where heat is inputted and extracted at a constant rate from opposite metals. In this case, the thermal conductance is calculated from the nonequilibrium temperature profile that is created at the junction. (ii) The Approach-to-Equilibrium MD (AEMD) method, with the thermal conductance of the junction obtained from the equilibration dynamics of the metals. In both methods, simulations of alkane chains of a growing size display an approximate length-independence of the thermal conductance, with calculated values matching computational and experimental studies. The RNEMD and AEMD methods offer different insights, and we discuss their benefits and shortcomings. Assessing the potential application of molecular junctions as thermal diodes, alkane junctions are made spatially asymmetric by modifying their contact regions with the bulk, either by using distinct endgroups or by replacing one of the Au contacts with Ag. Anharmonicity is built into the system within the molecular force-field. We find that, while the temperature profile strongly varies (compared with the gold-alkanedithiol-gold junctions) due to these structural modifications, the thermal diode effect is inconsequential in these systems—unless one goes to very large thermal biases. This finding suggests that one should seek molecules with considerable internal anharmonic effects for developing nonlinear thermal devices.more » « less
-
more » « less
DL_POLY Quantum 2.0, a vastly expanded software based on DL_POLY Classic 1.10, is a highly parallelized computational suite written in FORTRAN77 with a modular structure for incorporating nuclear quantum effects into large-scale/long-time molecular dynamics simulations. This is achieved by presenting users with a wide selection of state-of-the-art dynamics methods that utilize the isomorphism between a classical ring polymer and Feynman’s path integral formalism of quantum mechanics. The flexible and user-friendly input/output handling system allows the control of methodology, integration schemes, and thermostatting. DL_POLY Quantum is equipped with a module specifically assigned for calculating correlation functions and printing out the values for sought-after quantities, such as dipole moments and center-of-mass velocities, with packaged tools for calculating infrared absorption spectra and diffusion coefficients.
