Low frequency 2D Raman-THz spectroscopy of ionic solution: A simulation study

Monitoring MgCl2 hydrate formation from aqueous solutions using terahertz time-domain spectroscopy

The interaction of MgCl₂ with H₂O is heavily involved in biological and chemical processes. In this work, freezing-induced hydrate formation from MgCl₂ aqueous solution was monitored using terahertz time-domain spectroscopy. At low temperatures, two phase transitions from brine to hydrate formation could be clearly observed, and the formation of hydrate was accompanied by the emergence of new THz fingerprint peaks at 1.02, 1.56, and 1.84 THz, respectively. Integrating XRD and quantum chemical calculations, we attributed the absorption peaks to the vibrational modes of the formed MgCl₂·12H₂O. This demonstrates the potential of THz spectroscopy for application in the detection of biological processes in low-temperature environments, such as cell freezing.

Full molecular dynamics simulations of molecular liquids for single-beam spectrally controlled two-dimensional Raman spectroscopy

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.

Two-dimensional terahertz spectroscopy of condensed-phase molecular systems

Nonlinear terahertz (THz) spectroscopy relies on the interaction of matter with few-cycle THz pulses of electric field amplitudes up to megavolts/centimeter (MV/cm). In condensed-phase molecular systems, both resonant interactions with elementary excitations at low frequencies such as intra- and intermolecular vibrations and nonresonant field-driven processes are relevant. Two-dimensional THz (2D-THz) spectroscopy is a key method for following nonequilibrium processes and dynamics of excitations to decipher the underlying interactions and molecular couplings. This article addresses the state of the art in 2D-THz spectroscopy by discussing the main concepts and illustrating them with recent results. The latter include the response of vibrational excitations in molecular crystals up to the nonperturbative regime of light-matter interaction and field-driven ionization processes and electron transport in liquid water.

A Feynman diagram description of the 2D-Raman-THz response of amorphous ice

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.

Terahertz echoes reveal the inhomogeneity of aqueous salt solutions

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 t₁= t₂ 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.

Notes on simulating two-dimensional Raman and terahertz-Raman signals with a full molecular dynamics simulation approach

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.