Surface hopping simulation of the vibrational relaxation of I2 in liquid xenon using the collective probabilities algorithm

Vibrational energy relaxation of a diatomic molecule in a superfluid helium nanodroplet: influence of the nanodroplet size, interaction energy and energy gap

The influence of the nanodroplet size, molecule-helium interaction potential energy and ν = 1 - ν = 0 vibrational energy gap on the vibrational energy relaxation (VER) of a diatomic molecule (X₂) in a superfluid helium nanodroplet [HeND or (⁴He)_N; finite quantum solvent at T = 0.37 K] has been studied using a hybrid quantum approach recently proposed by us and taking as a reference the VER results on the I₂@(⁴He)₁₀₀ doped nanodroplet (Vilà et al., Phys. Chem. Chem. Phys., 2018, 20, 118, which corresponds to the first theoretical study on the VER of molecules embedded in a HeND). This has allowed us to obtain a deeper insight into the vibrational relaxation dynamics. The nanodroplet size has a very small effect on the VER, as this process mainly depends on the interaction between the molecule and the nanodroplet first solvation shell. Regarding the interaction potential energy and the energy gap, both factors play an important and comparable role in the VER time properties (global relaxation time, lifetime and transition time). As the former becomes stronger the relaxation time properties decrease in a significant way (their inverse follows a linear dependence with respect to the ν = 1 - ν = 0 coupling term) and they also decrease in a significant manner when the energy gap diminishes (linear dependence on the ν = 1 - ν = 0 energy difference). We expect that this study will motivate further work on the vibrational relaxation process in HeNDs.

Vibrational energy relaxation dynamics of diatomic molecules inside superfluid helium nanodroplets. The case of the I2 molecule

The vibrational relaxation (VER) of a X₂ molecule in a ⁴He superfluid nanodroplet (HeND; 0.37 K) was studied adapting a quantum approach recently proposed by us. In the first theoretical study on the VER of molecules inside HeND the I₂ molecule was examined [cascade mechanism (ν → ν − 1; ν − 1 → ν − 2; …) and time scale of ns].

Hybrid quantum/classical simulations of the vibrational relaxation of the amide I mode of N-methylacetamide in D2O solution

Hybrid quantum/classical molecular dynamics (MD) is applied to simulate the vibrational relaxation (VR) of the amide I mode of deuterated N-methylacetamide (NMAD) in aqueous (D(2)O) solution. A novel version of the vibrational molecular dynamics with quantum transitions (MDQT) treatment is developed in which the amide I mode is treated quantum mechanically while the remaining degrees of freedom are treated classically. The instantaneous normal modes of the initially excited NMAD molecule (INM(0)) are used as internal coordinates since they provide a proper initial partition of the system in quantum and classical subsystems. The evolution in time of the energy stored in each individual normal mode is subsequently quantified using the hybrid quantum-classical instantaneous normal modes (INM(t)). The identities of both the INM(0)s and the INM(t)s are tracked using the equilibrium normal modes (ENMs) as templates. The results extracted from the hybrid MDQT simulations show that the quantum treatment of the amide I mode accelerates the whole VR process versus pure classical simulations and gives better agreement with experiments. The relaxation of the amide I mode is found to be essentially an intramolecular vibrational redistribution (IVR) process with little contribution from the solvent, in agreement with previous theoretical and experimental studies. Two well-defined relaxation mechanisms are identified. The faster one accounts for ≈40% of the total vibrational energy that flows through the NMAD molecule and involves the participation of the lowest frequency vibrations as short-life intermediate modes. The second and slower mechanism accounts for the remaining ≈60% of the energy released and is associated to the energy flow through specific mid-range and high-frequency modes.

Mixed Quantum-Classical Molecular Dynamics Analysis of the Molecular-Level Mechanisms of Vibrational Frequency Shifts

A detailed analysis of the origins of vibrational frequency shifts of diatomic molecules (I2 and ICl) in a rare gas (Xe) liquid is presented. Specifically, vibrationally adiabatic mixed quantum-classical molecular dynamics simulations are used to obtain the instantaneous frequency shifts and correlate the shifts to solvent configurations. With this approach, important mechanistic questions are addressed, including the following: How many solvent atoms determine the frequency shift? What solvent atom configurations lead to blue shifts, and which lead to red shifts? What is the effect of solute asymmetry? The mechanistic analysis can be generally applied and should be useful in understanding what information is provided by infrared and Raman spectra about the environment of the probed vibrational mode.

The Ehrenfest method with quantum corrections to simulate the relaxation of molecules in solution: Equilibrium and dynamics

The use of the Ehrenfest method to simulate the relaxation of molecules in solution is explored. Using the cyanide ion dissolved in water as a test model, the independent trajectory (IT) and the bundle of trajectories (BT) approximations are shown to provide very different results for the time evolution of the vibrational populations of the solute. None of these approximations reproduce the Boltzmann equilibrium vibrational populations accurately. A modification of the Ehrenfest method based on the use of quantum correction factors is thus proposed to solve this problem. The simulations carried out using the modified Ehrenfest method provide IT and BT relaxation times which are closer to each other and which agree quite well with previous hybrid perturbative results.

Fewest Switches Adiabatic Surface Hopping As Applied to Vibrational Energy Relaxation

In this contribution quantum/classical surface hopping methodology is applied to vibrational energy relaxation of a quantum oscillator in a classical heat bath. The model of a linearly damped (harmonic) oscillator is chosen which can be mapped onto the Brownian motion (Caldeira-Leggett) Hamiltonian. In the simulations Tully's fewest switches surface hopping scheme is adopted with inclusion of dephasing in the adiabatic basis using a simple decoherence algorithm. The results are compared to the predictions of a Redfield-type quantum master equation modeling using the classical heat bath force correlation function as input. Thereby a link is established between both types of quantum/classical approaches. Viewed from the latter perspective, surface hopping with dephasing may be interpreted as "on-the-fly" stochastic realization of a quantum/classical Pauli master equation.