Elucidating Heavy Atom Tunneling Kinetics in the Cope Rearrangement of Semibullvalene.

In this work, we characterize the temperature dependence of kinetic properties in heavy atom tunneling reactions by means of molecular dynamics simulations, including nuclear quantum effects (NQEs) via Path Integral theory. To this end, we consider the prototypical Cope rearrangement of semibullvalene. The reaction was studied in the 25-300K temperature range observing that the inclusion of NQEs modifies the temperature behavior of both free energy barriers and dynamical recrossing factors with respect to classical dynamics. Notably, while in classical simulations the activation free energy shows a very weak temperature dependence, it becomes strongly dependent on temperature when NQEs are included. This temperature behavior shows a transition from a regime where the quantum effects are limited and can mainly be traced back to zero point energy, to a low temperature regime where tunneling plays a dominant role. In this regime, the free energy curve tunnels below the potential energy barrier along the reaction coordinate,  allowing much faster reaction rates. Finally, the temperature dependence of the rate constants obtained from molecular dynamics simulations was compared with available experimental data and with semi-classical transition state theory calculations, showing comparable behaviors and similar transition temperatures from  thermal to (deep) tunneling regime.

Quantum versus classical unimolecular fragmentation rate constants and activation energies at finite temperature from direct dynamics simulations.

In the present work, we investigate how nuclear quantum effects modify the temperature dependent rate constants and, consequently, the activation energies in unimolecular reactions. In the reactions under study, nuclear quantum effects mainly stem from the presence of a large zero point energy. Thus, we investigate the behavior of methods compatible with direct dynamics simulations, the quantum thermal bath (QTB) and ring polymer molecular dynamics (RPMD). To this end, we first compare them with quantum reaction theory for a model Morse potential before extending this comparison to molecular models. Our results show that, in particular in the temperature range comparable with or lower than the zero point energy of the system, the RPMD method is able to correctly capture nuclear quantum effects on rate constants and activation energies. On the other hand, although the QTB provides a good description of equilibrium properties including zero-point energy effects, it largely overestimates the rate constants. The origin of the different behaviours is in the different distance distributions provided by the two methods and in particular how they differently describe the tails of such distributions. The comparison with transition state theory shows that RPMD can be used to study fragmentation of complex systems for which it may be difficult to determine the multiple reaction pathways and associated transition states.

Does High Pressure Induce Structural Reorganization in Linear Alcohols? A Computational Answer.

We present an exhaustive computational study on the effect of high pressure on normal alcohols with alkyl chains with lengths of three-to-eight carbon atoms. 1-Propanol, 1-butanol, 1-pentanol, 1-hexanol, 1-heptanol, and 1-octanol were studied by using classical molecular dynamics simulations and applying pressures in the range of 1 to 104  bar. The results of our calculations show that high-pressure values affect the structure significantly. In particular, we have observed a marked difference in behavior for alcohols with chain lengths below six and those with more than six or seven carbon atoms, with hexanol and heptanol being boundary cases. We have named the model with the most shrunk alkyl chains as the Asclepius form inspired by the Rod of Asclepius, the universally known symbol of medicine, in which a snake is coiled around a rod.

Environmental and Nuclear Quantum Effects on Double Proton Transfer in the Guanine-Cytosine Base Pair.

In the present letter, we investigate the double proton transfer (DPT) tautomerization process in guanine-cytosine (GC) DNA base pairs. In particular, we study the influence of the biological environment on the mechanism, the kinetics and thermodynamics of such DPT. To this end, we present a molecular dynamics (MD) study in the tight-binding density functional theory framework, and compare the reactivity of the isolated GC dimer with that of the same dimer embedded in a small DNA structure. The impact of nuclear quantum effects (NQEs) is also evaluated using Path Integral based MD. Results show that in the isolated dimer, the DPT occurs via a concerted mechanism, while in the model biological environment, it turns into a stepwise process going through an intermediate structure. One of the water molecules in the vicinity of the proton transfer sites plays an important role as it changes H-bond pattern during the DPT reaction. The inclusion of NQEs has the effect of speeding up the tautomeric-to-canonical reaction, reflecting the destabilization of both the tautomeric and intermediate forms.