Constructing grids for molecular quantum dynamics using an autoencoder

A challenge for molecular quantum dynamics (QD) calculations is the curse of dimensionality with respect to the nuclear degrees of freedom. A common approach that works especially well for fast reactive processes is to reduce the dimensionality of the system to a few most relevant coordinates. Identifying these can become a very difficult task, since they often are highly unintuitive. We present a machine learning approach that utilizes an autoencoder that is trained to find a low-dimensional representation of a set of molecular configurations. These configurations are generated by trajectory calculations performed on the reactive molecular systems of interest. The resulting low-dimensional representation can be used to generate a potential energy surface grid in the desired subspace. Using the G-matrix formalism to calculate the kinetic energy operator, QD calculations can be carried out on this grid. In addition to step-by-step instructions for the grid construction, we present the application to a test system.Comment: 24 pages, 6 figures, articl

Quantum Dynamics in an Explicit Solvent Environment: A Photochemical Bond Cleavage Treated with a Combined QD/MD Approach

In quantum chemistry methods to describe environmental effects on different levels of complexity are available in the common program packages. Electrostatic effects of a solvent for example can be included in an implicit or explicit way. For chemical reactions with large structural changes additional mechanical effects come into play. Their treatment within a quantum dynamical context is a major challenge, especially when excited states are involved. Recently, we introduced a method that realizes an implicit description. Here, we present an approach combining quantum dynamics and molecular dynamics. It explicitly incorporates the solvent environment, whereby the electrostatic as well as the dynamic effects are captured. This new method is demonstrated for the ultrafast photoinduced bond cleavage of diphenylmethylphosphonium ions (Ph₂CH–PPh₃⁺), a common precursor to generate reactive carbocations in solution

Quantum Dynamics of a Photochemical Bond Cleavage Influenced by the Solvent Environment: A Dynamic Continuum Approach

In every day chemistry, solvents are used to influence the outcome of chemical synthesis. Electrostatic effects stabilize polar configurations during the reaction and in addition dynamic solvent effects can emerge. How the dynamic effects intervene on the ultrafast time scale is in the focus of this theoretical study. We selected the photoinduced bond cleavage of Ph₂CH–PPh₃⁺ for which the electrostatic interactions are negligible. Elaborate ultrafast pump–probe studies already exist and serve as a reference. We compared quantum dynamical simulations with and without environment and noticed the necessity to model the influence of the solvent cage on the reactive motions of the solute. The frictional force induced by the dynamic viscosity of the solvent is implemented in the quantum mechanical formalism with a newly developed approach called the dynamic continuum ansatz. Only when the environment is included are the experimentally observed products reproduced on the subpicosecond time scale

Ultrafast Reactive Quantum Dynamics Coupled to Classical Solvent Dynamics Using an Ehrenfest Approach

The inclusion of solvent effects in the theoretical analysis of molecular processes becomes increasingly important. Currently, it is not feasible to directly include the solvent on the quantum level. We use an Ehrenfest approach to study the coupled time evolution of quantum dynamically treated solutes and classical solvents system. The classical dynamics of the solvent is coupled to the wavepacket dynamics of the solute and rotational and translational degrees of freedom of the solute are included classically. This allows quantum dynamics simulations for ultrafast processes that are decided by environment interactions without explicit separation of time scales. We show the application to the dissociation of ICN in liquid Ar as a proof of principal system and to the more applied example of uracil in water