Mean first-passage times for solvated LiCN isomerization at intermediate to high temperatures.

The behavior of a particle in a solvent has been framed using stochastic dynamics since the early theory of Kramers. A particle in a chemical reaction reacts slower in a diluted solvent because of the lack of energy transfer via collisions. The flux-over-population reaction rate constant rises with increasing density before falling again for very dense solvents. This Kramers turnover is observed in this paper at intermediate and high temperatures in the backward reaction of the LiNC ⇌ LiCN isomerization via Langevin dynamics and mean first-passage times (MFPTs). It is in good agreement with the Pollak-Grabert-Hänggi (PGH) reaction rates at lower temperatures. Furthermore, we find a square root behavior of the reaction rate at high temperatures and have made direct comparisons of the methods in the intermediate- and high-temperature regimes, all suggesting increased ranges in accuracy of both the PGH and MFPT approaches.

Solvated molecular dynamics of LiCN isomerization: All-atom argon solvent versus a generalized Langevin bath.

The reaction rate rises and falls with increasing density or friction when a molecule is activated by collisions with the solvent particles. This so-called Kramers turnover has recently been observed in the isomerization reaction of LiCN in an argon bath. In this paper, we demonstrate by direct comparison with those results that a reduced-dimensional (generalized) Langevin description gives rise to similar reaction dynamics as the corresponding (computationally expensive) full molecular dynamics calculations. We show that the density distributions within the Langevin description are in direct agreement with the full molecular dynamics results and that the turnover in the reaction rates is reproduced qualitatively and quantitatively at different temperatures.

The role of the CN vibration in the activated dynamics of LiNC<−>LiCN isomerization in an argon solvent at high temperatures.

The isomerization between CN-Li and Li-CN in an argon bath provides a paradigmatic example of a reaction in a solvent with tunable coupling. In previous work, we found that the rates exhibited a turnover with the density of the argon bath in the limit that the CN bond was held fixed [P. L. Garcia-Muller, R. Hernandez, R. M. Benito, and F. Borondo, J. Chem. Phys. 137, 204301 (2012)]. Here, we report the effect of the CN bond vibration on the dynamics and the persistence of the turnover. As hypothesized earlier, the CN bond is indeed weakly coupled with the reaction path despite the presence of the argon cage.