The ethanol electro-oxidation catalyzed by Pd in alkaline environment involves several intermediate reaction steps promoted by the hydroxyl radical, OH*. In this work, we report on dynamical paths of the...

We report on an extended hydrodynamic modeling of the friction tensorial properties of flexible molecules including all types of natural, Z‐Matrix like, internal coordinates. We implement the new methodology by extending and updating the software DiTe [Barone et al. J. Comput. Chem. 30, 2 (2009)]. DiTe (DIffusion TEnsor) implements a hydrodynamic modeling of the generalized translational, rotational, and configurational friction and diffusion tensors of flexible molecules in which flexibility is described in terms of dihedral angles. The new tool, DiTe2, has been renewed to include also stretching and bending types of internal mobility. Furthermore, DiTe2 is able to calculate the friction and diffusion tensors along collective (or reaction) coordinates defined as linear combinations of the internal natural ones. A number of tests are reported to show the new features of DiTe2. As leitmotiv for the tests, the calmodulin protein is taken into consideration, described both at all‐atom and coarse‐grained levels. © 2018 Wiley Periodicals, Inc.

This work outlines the development and application of a multiscale computational protocol to evaluate reaction rates of elementary reactions in internal natural coordinates.

This work follows a companion article, which will be referred to as Paper I [Campeggio et al., J. Chem. Phys. 158, 244104 (2023)] in which a quantum-stochastic Liouville equation for the description of the quantum–classical dynamics of a molecule in a dissipative bath has been formulated in curvilinear internal coordinates. In such an approach, the coordinates of the system are separated into three subsets: the quantum coordinates, the classical relevant nuclear degrees of freedom, and the classical irrelevant (bath) coordinates. The equation has been derived in natural internal coordinates, which are bond lengths, bond angles, and dihedral angles. The resulting equation needs to be parameterized. In particular, one needs to compute the potential energy surfaces, the friction tensor, and the rate constants for the nonradiative jumps among the quantum states. While standard methods exist for the calculation of energy and dissipative properties, an efficient evaluation of the transition rates needs to be developed. In this paper, an approximated treatment is introduced, which leads to a simple explicit formula with a single adjustable parameter. Such an approximated expression is compared with the exact calculation of transition rates obtained via molecular dynamics simulations. To make such a comparison possible, a simple sandbox system has been used, with two quantum states and a single internal coordinate (together with its conjugate momentum). Results show that the adjustable parameter, which is an effective decoherence time, can be parameterized from the effective relaxation times of the autocorrelation functions of the conjugated momenta of the relevant nuclear coordinates.

Multiscale methods are powerful tools to describe large and complex systems. They are based on a hierarchical partitioning of the degrees of freedom (d.o.f.) of the system, allowing one to treat each set of d.o.f. in the most computationally efficient way. In the context of coupled nuclear and electronic dynamics, a multiscale approach would offer the opportunity to overcome the computational limits that, at present, do not allow one to treat a complex system (such as a biological macromolecule in explicit solvent) fully at the quantum mechanical level. Based on the pioneering work of Kapral and Ciccotti [R. Kapral and G. Ciccotti, J. Chem. Phys.110, 8919 (1999)], this work is intended to present a nonadiabatic theory that describes the evolution of electronic populations coupled with the dynamics of the nuclei of a molecule in a dissipative environment (condensed phases). The two elements of novelty that are here introduced are (i) the casting of the theory in the natural, internal coordinates, that are bond lengths, bond angles, and dihedral angles; (ii) the projection of those nuclear d.o.f. that can be considered at the level of a thermal bath, therefore leading to a quantum-stochastic Liouville equation. Using natural coordinates allows the description of structure and dynamics in the way chemists are used to describe molecular geometry and its changes. The projection of bath coordinates provides an important reduction of complexity and allows us to formulate the approach that can be used directly in the statistical thermodynamics description of chemical systems.

We estimate the kinetic constants of a series of archetypal SN2 reactions, i.e., the nucleophilic substitutions of halides in halomethane. A parameter free, multiscale approach recently developed [Campeggio et al.,...