1
|
Begušić T, Blake GA. Two-dimensional infrared-Raman spectroscopy as a probe of water's tetrahedrality. Nat Commun 2023; 14:1950. [PMID: 37029146 PMCID: PMC10082090 DOI: 10.1038/s41467-023-37667-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2022] [Accepted: 03/22/2023] [Indexed: 04/09/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Tomislav Begušić
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, 91125, USA.
| | - Geoffrey A Blake
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, 91125, USA.
- Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA, 91125, USA.
| |
Collapse
|
2
|
Takahashi H, Tanimura Y. Discretized hierarchical equations of motion in mixed Liouville-Wigner space for two-dimensional vibrational spectroscopies of liquid water. J Chem Phys 2023; 158:044115. [PMID: 36725520 DOI: 10.1063/5.0135725] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Hideaki Takahashi
- Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan
| | - Yoshitaka Tanimura
- Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan
| |
Collapse
|
3
|
Sertcan B, Mousavi SJ, Iannuzzi M, Hamm P. Low-frequency anharmonic couplings in crystalline bromoform: Theory. J Chem Phys 2023; 158:014203. [PMID: 36610974 DOI: 10.1063/5.0134278] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Theoretical calculations of the low-frequency anharmonic couplings of the β-phase of crystalline bromoform are presented based on density functional theory quantum chemistry calculations. The electrical and mechanical anharmonicities between intra- and intermolecular modes are calculated, revealing that the electrical anharmonicity dominates the cross-peak intensities in the 2D Raman-THz response and crystalline, as well as liquid, bromoform. Furthermore, the experimentally observed difference in relative cross-peak intensities between the two intramolecular modes of bromoform and the intermolecular modes can be explained by the C3v-symmetry of bromoform in combination with orientational averaging. The good agreement with the experimental results provides further evidence for our interpretation that the 2D Raman-THz response of bromoform is, indeed, related to the anharmonic coupling between the intra- and intermolecular modes.
Collapse
Affiliation(s)
- Beliz Sertcan
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Seyyed Jabbar Mousavi
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Marcella Iannuzzi
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Peter Hamm
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| |
Collapse
|
4
|
Begušić T, Tao X, Blake GA, Miller TF. Equilibrium-nonequilibrium ring-polymer molecular dynamics for nonlinear spectroscopy. J Chem Phys 2022; 156:131102. [PMID: 35395895 DOI: 10.1063/5.0087156] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Tomislav Begušić
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Xuecheng Tao
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Geoffrey A Blake
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Thomas F Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
5
|
Jo JY, Tanimura Y. Full molecular dynamics simulations of molecular liquids for single-beam spectrally controlled two-dimensional Raman spectroscopy. J Chem Phys 2021; 154:124115. [PMID: 33810650 DOI: 10.1063/5.0044661] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Single-beam spectrally controlled (SBSC) two-dimensional (2D) Raman spectroscopy is a unique 2D vibrational measurement technique utilizing trains of short pulses that are generated from a single broadband pulse by pulse shaping. This approach overcomes the difficulty of 2D Raman spectroscopy in dealing with small-signal extraction and avoids complicated low-order cascading effects, thus providing a new possibility for measuring the intramolecular and intermolecular modes of molecular liquids using fifth-order 2D Raman spectroscopy. Recently, for quantitatively investigating the mode-mode coupling mechanism, Hurwitz et al. [Opt. Express 28, 3803 (2020)] have developed a new pulse design for this measurement to separate the contributions of the fifth- and third-order polarizations, which are often overlapped in the original single-beam measurements. Here, we describe a method for simulating these original measurements and the new 2D Raman measurements on the basis of a second-order response function approach. We carry out full molecular dynamics simulations for carbon tetrachloride and liquid water using an equilibrium-nonequilibrium hybrid algorithm, with the aim of explaining the key features of the SBSC 2D Raman spectroscopic method from a theoretical point of view. The predicted signal profiles and intensities provide valuable information that can be applied to 2D spectroscopy experiments, allowing them to be carried out more efficiently.
Collapse
Affiliation(s)
- Ju-Yeon Jo
- Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan
| | - Yoshitaka Tanimura
- Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan
| |
Collapse
|
6
|
Shalit A, Mousavi SJ, Hamm P. 2D Raman–THz Spectroscopy of Binary CHBr3–MeOH Solvent Mixture. J Phys Chem B 2021; 125:581-586. [DOI: 10.1021/acs.jpcb.0c08962] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- A. Shalit
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - S. J. Mousavi
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - P. Hamm
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| |
Collapse
|
7
|
Sidler D, Hamm P. A Feynman diagram description of the 2D-Raman-THz response of amorphous ice. J Chem Phys 2020; 153:044502. [PMID: 32752676 DOI: 10.1063/5.0018485] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The 2D-Raman-THz response in all possible time-orderings (Raman-THz-THz, THz-Raman-THz, and THz-THz-Raman) of amorphous water ice is calculated in two ways: from atomistic molecular dynamics simulations and with the help of a Feynman diagram model, the latter of which power-expands the potential energy surface and the dipole and polarizability surfaces up to leading order. Comparing both results allows one to dissect the 2D-Raman-THz response into contributions from mechanical anharmonicity, as well as electrical dipole and polarizability anharmonicities. Mechanical anharmonicity dominates the 2D-Raman-THz response of the hydrogen-bond stretching and hydrogen-bond bending bands of water, and dipole anharmonicity dominates that of the librational band, while the contribution of polarizability anharmonicity is comparably weak. A distinct echo of the hydrogen-bond stretching band is observed for the THz-Raman-THz pulse sequence, again dominated by mechanical anharmonicity. A peculiar mechanism is discussed, which is based on the coupling between the many normal modes within the hydrogen-bond stretching band and which will inevitably generate such an echo for an amorphous structure.
Collapse
Affiliation(s)
- David Sidler
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Peter Hamm
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| |
Collapse
|
8
|
Sun X. Hybrid equilibrium-nonequilibrium molecular dynamics approach for two-dimensional solute-pump/solvent-probe spectroscopy. J Chem Phys 2019; 151:194507. [DOI: 10.1063/1.5130926] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Affiliation(s)
- Xiang Sun
- Division of Arts and Sciences, NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China; Department of Chemistry, New York University, New York, New York 10003, USA; NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China; and State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200062, China
| |
Collapse
|
9
|
Magdău IB, Mead GJ, Blake GA, Miller TF. Interpretation of the THz-THz-Raman Spectrum of Bromoform. J Phys Chem A 2019; 123:7278-7287. [DOI: 10.1021/acs.jpca.9b05165] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Ioan B. Magdău
- Division of Chemistry & Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Griffin J. Mead
- Division of Chemistry & Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Geoffrey A. Blake
- Division of Chemistry & Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
- Division of Geological & Planetary Sciences, California Institute of Technology, Pasadena, California 91125, United States
| | - Thomas F. Miller
- Division of Chemistry & Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
10
|
Jansen TLC, Saito S, Jeon J, Cho M. Theory of coherent two-dimensional vibrational spectroscopy. J Chem Phys 2019; 150:100901. [DOI: 10.1063/1.5083966] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Thomas la Cour Jansen
- University of Groningen, Zernike Institute for Advanced Materials, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Shinji Saito
- Institute for Molecular Science, Myodaiji, Okazaki, Aichi 444-8585, Japan and The Graduate University for Advanced Studies, Myodaiji, Okazaki, Aichi 444-8585, Japan
| | - Jonggu Jeon
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science (IBS), Seoul 02841, South Korea
| | - Minhaeng Cho
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science (IBS), Seoul 02841, South Korea
- Department of Chemistry, Korea University, Seoul 02841, South Korea
| |
Collapse
|
11
|
Sidler D, Hamm P. Feynman diagram description of 2D-Raman-THz spectroscopy applied to water. J Chem Phys 2019; 150:044202. [DOI: 10.1063/1.5079497] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
Affiliation(s)
- David Sidler
- Department of Chemistry, University of Zurich, Zurich, Switzerland
| | - Peter Hamm
- Department of Chemistry, University of Zurich, Zurich, Switzerland
| |
Collapse
|
12
|
Zhang Q, Pan Z, Zhang L, Zhang R, Chen Z, Jin T, Wu T, Chen X, Zhuang W. Ion effect on the dynamics of water hydrogen bonding network: A theoretical and computational spectroscopy point of view. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2018. [DOI: 10.1002/wcms.1373] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Qiang Zhang
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of MatterChinese Academy of SciencesFuzhouChina
- Department of ChemistryBohai UniversityJinzhouChina
| | - Zhijun Pan
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of MatterChinese Academy of SciencesFuzhouChina
| | - Lu Zhang
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of MatterChinese Academy of SciencesFuzhouChina
| | - Ruiting Zhang
- School of Physics and Optoelectronic EngineeringXidian UniversityXi'anChina
| | - Zhening Chen
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of MatterChinese Academy of SciencesFuzhouChina
| | - Tan Jin
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of MatterChinese Academy of SciencesFuzhouChina
| | - Tianmin Wu
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of MatterChinese Academy of SciencesFuzhouChina
| | - Xian Chen
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of MatterChinese Academy of SciencesFuzhouChina
| | - Wei Zhuang
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of MatterChinese Academy of SciencesFuzhouChina
| |
Collapse
|
13
|
Ramesh P, Loring RF. Thermal Population Fluctuations in Two-Dimensional Infrared Spectroscopy Captured with Semiclassical Mechanics. J Phys Chem B 2018; 122:3647-3654. [DOI: 10.1021/acs.jpcb.7b12122] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Prashanth Ramesh
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York 14853, United States
| | - Roger F. Loring
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York 14853, United States
| |
Collapse
|
14
|
|
15
|
Full molecular dynamics simulations of liquid water and carbon tetrachloride for two-dimensional Raman spectroscopy in the frequency domain. Chem Phys 2016. [DOI: 10.1016/j.chemphys.2016.07.002] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
16
|
Terahertz echoes reveal the inhomogeneity of aqueous salt solutions. Nat Chem 2016; 9:273-278. [PMID: 28221356 DOI: 10.1038/nchem.2642] [Citation(s) in RCA: 81] [Impact Index Per Article: 10.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2016] [Accepted: 09/09/2016] [Indexed: 12/15/2022]
Abstract
The structural and dynamical properties of water are known to be affected by ion solvation. However, a consistent molecular picture that describes how and to what extent ions perturb the water structure is still missing. Here we apply 2D Raman-terahertz spectroscopy to investigate the impact of monatomic cations on the relaxation dynamics of the hydrogen-bond network in aqueous salt solutions. The inherent ability of multidimensional spectroscopy to deconvolute heterogeneous relaxation dynamics is used to reveal the correlation between the inhomogeneity of the collective intermolecular hydrogen-bond modes and the viscosity of a salt solution. Specifically, we demonstrate that the relaxation time along the echo direction t1 = t2 correlates with the capability of a given cation to 'structure' water. Moreover, we provide evidence that the echo originates from the water-water modes, and not the water-cation modes, which implies that cations can structure the hydrogen-bond network to a certain extent.
Collapse
|
17
|
Ito H, Hasegawa T, Tanimura Y. Effects of Intermolecular Charge Transfer in Liquid Water on Raman Spectra. J Phys Chem Lett 2016; 7:4147-4151. [PMID: 27689824 DOI: 10.1021/acs.jpclett.6b01766] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
The low-frequency vibrational spectrum of liquid water is composed of contributions from the intermolecular librational and translation modes. The existence of these two modes introduces difficulty into the simulation of experimentally obtained Raman spectra. We constructed a polarizability function for a water model that includes intramolecular charge flow (CF) effects, intermolecular charge transfer (CT) effects, and intermolecular dipole-induced-dipole (DID) effects. We computed the one-dimensional (1D) Raman and terahertz (THz) spectra with all of these effects included (CFCT-DID) and compared with experimental spectra. We find that the CFCT-DID function provides a better description of the experimental results, because the CT effects reduce the polarizability only for translational motion to which parallelly polarized (VV) and perpendicularly polarized (VH) Raman spectra are sensitive. In our calculations of two-dimensional (2D) Raman and THz-Raman spectra, we observe the enhancement of echo signals in both cases. The details of the CFCT-DID function, along with its source code, are provided in the Supporting Information.
Collapse
Affiliation(s)
- Hironobu Ito
- Department of Chemistry, Graduate School of Science, Kyoto University , Sakyoku, Kyoto 606-8502, Japan
| | - Taisuke Hasegawa
- Department of Chemistry, Graduate School of Science, Kyoto University , Sakyoku, Kyoto 606-8502, Japan
| | - Yoshitaka Tanimura
- Department of Chemistry, Graduate School of Science, Kyoto University , Sakyoku, Kyoto 606-8502, Japan
| |
Collapse
|
18
|
Perakis F, De Marco L, Shalit A, Tang F, Kann ZR, Kühne TD, Torre R, Bonn M, Nagata Y. Vibrational Spectroscopy and Dynamics of Water. Chem Rev 2016; 116:7590-607. [DOI: 10.1021/acs.chemrev.5b00640] [Citation(s) in RCA: 225] [Impact Index Per Article: 28.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Affiliation(s)
- Fivos Perakis
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Luigi De Marco
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
- Department of Chemistry and James Franck Institute, University of Chicago, 929 East 57th Street, Chicago, Illinois 60637, United States
| | - Andrey Shalit
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Fujie Tang
- International Center for Quantum Materials, Peking University, 5 Yiheyuan Road, Haidian, Beijing 100871, China
| | - Zachary R. Kann
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706, United States,
| | - Thomas D. Kühne
- Department of Chemistry, University of Paderborn, Warburger Strasse 100, D-33098 Paderborn, Germany
| | - Renato Torre
- European Lab for Nonlinear Spectroscopy and Dipartimento di Fisica e Astronomia, Università di Firenze, Via Nello Carrara 1, Sesto Fiorentino (Firenze) I-50019, Italy
| | - Mischa Bonn
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Yuki Nagata
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| |
Collapse
|
19
|
Ito H, Tanimura Y. Simulating two-dimensional infrared-Raman and Raman spectroscopies for intermolecular and intramolecular modes of liquid water. J Chem Phys 2016; 144:074201. [DOI: 10.1063/1.4941842] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
Affiliation(s)
- Hironobu Ito
- Department of Chemistry, Graduate School of Science, Kyoto University, Sakyoku, Kyoto 606-8502, Japan
| | - Yoshitaka Tanimura
- Department of Chemistry, Graduate School of Science, Kyoto University, Sakyoku, Kyoto 606-8502, Japan
| |
Collapse
|
20
|
Ito H, Jo JY, Tanimura Y. Notes on simulating two-dimensional Raman and terahertz-Raman signals with a full molecular dynamics simulation approach. STRUCTURAL DYNAMICS (MELVILLE, N.Y.) 2015; 2:054102. [PMID: 26798823 PMCID: PMC4711663 DOI: 10.1063/1.4932597] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/29/2015] [Accepted: 09/25/2015] [Indexed: 06/05/2023]
Abstract
Recent developments in two-dimensional (2D) THz-Raman and 2D Raman spectroscopies have created the possibility for quantitatively investigating the role of many dynamic and structural aspects of the molecular system. We explain the significant points for properly simulating 2D vibrational spectroscopic studies of intermolecular modes using the full molecular dynamics approach, in particular, regarding the system size, the treatment of the thermostat, and inclusion of an Ewald summation for the induced polarizability. Moreover, using the simulation results for water employing various polarization functions, we elucidate the roles of permanent and induced optical properties in determining the 2D profiles of the signal.
Collapse
Affiliation(s)
- Hironobu Ito
- Department of Chemistry, Graduate School of Science, Kyoto University , Sakyoku, Kyoto 606-8502, Japan
| | - Ju-Yeon Jo
- Department of Chemistry, Graduate School of Science, Kyoto University , Sakyoku, Kyoto 606-8502, Japan
| | - Yoshitaka Tanimura
- Department of Chemistry, Graduate School of Science, Kyoto University , Sakyoku, Kyoto 606-8502, Japan
| |
Collapse
|
21
|
Moberg DR, Alemi M, Loring RF. Thermal weights for semiclassical vibrational response functions. J Chem Phys 2015; 143:084101. [PMID: 26328812 DOI: 10.1063/1.4929377] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Semiclassical approximations to response functions can allow the calculation of linear and nonlinear spectroscopic observables from classical dynamics. Evaluating a canonical response function requires the related tasks of determining thermal weights for initial states and computing the dynamics of these states. A class of approximations for vibrational response functions employs classical trajectories at quantized values of action variables and represents the effects of the radiation-matter interaction by discontinuous transitions. Here, we evaluate choices for a thermal weight function which are consistent with this dynamical approximation. Weight functions associated with different semiclassical approximations are compared, and two forms are constructed which yield the correct linear response function for a harmonic potential at any temperature and are also correct for anharmonic potentials in the classical mechanical limit of high temperature. Approximations to the vibrational linear response function with quantized classical trajectories and proposed thermal weight functions are assessed for ensembles of one-dimensional anharmonic oscillators. This approach is shown to perform well for an anharmonic potential that is not locally harmonic over a temperature range encompassing the quantum limit of a two-level system and the limit of classical dynamics.
Collapse
Affiliation(s)
- Daniel R Moberg
- Baker Laboratory, Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Mallory Alemi
- Baker Laboratory, Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Roger F Loring
- Baker Laboratory, Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
22
|
Ikeda T, Ito H, Tanimura Y. Analysis of 2D THz-Raman spectroscopy using a non-Markovian Brownian oscillator model with nonlinear system-bath interactions. J Chem Phys 2015; 142:212421. [DOI: 10.1063/1.4917033] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Tatsushi Ikeda
- Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan
| | - Hironobu Ito
- Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan
| | - Yoshitaka Tanimura
- Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan
| |
Collapse
|
23
|
Pan Z, Wu T, Jin T, Liu Y, Nagata Y, Zhang R, Zhuang W. Low frequency 2D Raman-THz spectroscopy of ionic solution: A simulation study. J Chem Phys 2015; 142:212419. [DOI: 10.1063/1.4917260] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Zhijun Pan
- Department of Chemical Physics, University of Science and Technology of China, Hefei 230026, China
| | - Tianmin Wu
- Department of Chemical Physics, University of Science and Technology of China, Hefei 230026, China
| | - Tan Jin
- State Key Lab of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, China
| | - Yong Liu
- State Key Lab of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, China
| | - Yuki Nagata
- Department for Molecular Spectroscopy, Max-Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Ruiting Zhang
- State Key Lab of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, China
| | - Wei Zhuang
- State Key Lab of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, China
| |
Collapse
|
24
|
Hamm P. 2D-Raman-THz spectroscopy: a sensitive test of polarizable water models. J Chem Phys 2015; 141:184201. [PMID: 25399140 DOI: 10.1063/1.4901216] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
In a recent paper, the experimental 2D-Raman-THz response of liquid water at ambient conditions has been presented [J. Savolainen, S. Ahmed, and P. Hamm, Proc. Natl. Acad. Sci. U. S. A. 110, 20402 (2013)]. Here, all-atom molecular dynamics simulations are performed with the goal to reproduce the experimental results. To that end, the molecular response functions are calculated in a first step, and are then convoluted with the laser pulses in order to enable a direct comparison with the experimental results. The molecular dynamics simulation are performed with several different water models: TIP4P/2005, SWM4-NDP, and TL4P. As polarizability is essential to describe the 2D-Raman-THz response, the TIP4P/2005 water molecules are amended with either an isotropic or a anisotropic polarizability a posteriori after the molecular dynamics simulation. In contrast, SWM4-NDP and TL4P are intrinsically polarizable, and hence the 2D-Raman-THz response can be calculated in a self-consistent way, using the same force field as during the molecular dynamics simulation. It is found that the 2D-Raman-THz response depends extremely sensitively on details of the water model, and in particular on details of the description of polarizability. Despite the limited time resolution of the experiment, it could easily distinguish between various water models. Albeit not perfect, the overall best agreement with the experimental data is obtained for the TL4P water model.
Collapse
Affiliation(s)
- Peter Hamm
- Department of Chemistry, University of Zurich, Winterthurerstr. 190, CH-8057 Zürich, Switzerland
| |
Collapse
|