Constructing Grids for Molecular Quantum Dynamics Using an Autoencoder

The principal component analysis of the ring deformation in the nonadiabatic surface hopping dynamics

The analysis of the leading active molecular motions in the on-the-fly trajectory surface hopping simulation provides the essential information to understand the geometric evolution in nonadiabatic dynamics. When the ring deformation is involved, the identification of the key active coordinates becomes challenging. A "hierarchical" protocol based on the dimensionality reduction and clustering approaches is proposed for the automatic analysis of the ring deformation in the nonadiabatic molecular dynamics. The representative system keto isocytosine is taken as the prototype to illustrate this protocol. The results indicate that the current hierarchical analysis protocol is a powerful way to clearly clarify both the major and minor active molecular motions of the ring distortion in nonadiabatic dynamics.

Excited-State Distortions Promote the Photochemical 4π-Electrocyclizations of Fluorobenzenes via Machine Learning Accelerated Photodynamics Simulations

Benzene fluorination increases chemoselectivities for Dewar-benzenes via 4π-disrotatory electrocyclization. However, the origin of the chemo- and regioselectivities of fluorobenzenes remains unexplained because of the experimental limitations in resolving the excited-state structures on ultrafast timescales. The computational cost of multiconfigurational nonadiabatic molecular dynamics simulations is also currently cost-prohibitive. We now provide high-fidelity structural information and reaction outcome predictions with machine-learning-accelerated photodynamics simulations of a series of fluorobenzenes, C₆ F_6-n H_n , n=0-3, to study their S₁ →S₀ decay in 4 ns. We trained neural networks with XMS-CASPT2(6,7)/aug-cc-pVDZ calculations, which reproduced the S₁ absorption features with mean absolute errors of 0.04 eV (<2 nm). The predicted nonradiative decay constants for C₆ F₄ H₂ , C₆ F₆ , C₆ F₃ H₃ , and C₆ F₅ H are 116, 60, 28, and 12 ps, respectively, in broad qualitative agreement with the experiments. Our calculations show that a pseudo Jahn-Teller distortion of fluorinated benzenes leads to an S₁ local-minimum region that extends the excited-state lifetimes of fluorobenzenes. The pseudo Jahn-Teller distortions reduce when fluorination decreases. Our analysis of the S₁ dynamics shows that the pseudo-Jahn-Teller distortions promote an excited-state cis-trans isomerization of a π_C-C bond. We characterized the surface hopping points from our NAMD simulations and identified instantaneous nuclear momentum as a factor that promotes the electrocyclizations.

Analysis of bath motion in MM-SQC dynamics via dimensionality reduction approach: Principal component analysis

The system-plus-bath model is an important tool to understand the nonadiabatic dynamics of large molecular systems. Understanding the collective motion of a large number of bath modes is essential for revealing their key roles in the overall dynamics. Here, we applied principal component analysis (PCA) to investigate the bath motion in the basis of a large dataset generated from the symmetrical quasi-classical dynamics method based on the Meyer-Miller mapping Hamiltonian nonadiabatic dynamics for the excited-state energy transfer in the Frenkel-exciton model. The PCA method clearly elucidated that two types of bath modes, which either display strong vibronic coupling or have frequencies close to that of the electronic transition, are important to the nonadiabatic dynamics. These observations were fully consistent with the physical insights. The conclusions were based on the PCA of the trajectory data and did not involve significant pre-defined physical knowledge. The results show that the PCA approach, which is one of the simplest unsupervised machine learning dimensionality reduction methods, is a powerful one for analyzing complicated nonadiabatic dynamics in the condensed phase with many degrees of freedom.

Pathways to New Applications for Quantum Control

In 1998, the first successful quantum control experiment with application to a molecular framework was conducted with a shaped laser pulse, optimizing the branching ratio between different organometallic reaction channels. This work induced a vast activity in quantum control during the next 10 years, and different optimization aims were achieved in the gas phase, liquid phase, and even in biologically relevant molecules like light-harvesting complexes. Accompanying and preceding this development were important advances in theoretical quantum control simulations. They predicted several control scenarios and explained how and why quantum control experiments work. After many successful proofs of concept in molecular science, the big challenge is to expand its huge conceptual potential of directly being able to steer nuclear and/or electronic motion to more applied implementations. In this Account, based on several recent advances, we give a personal evaluation of where the field of molecular quantum control is at the moment and especially where we think promising applications can be in the near future. One of these paths leads to synthetic chemistry. The synthesis of novel pharmaceutical compounds or natural products often involves many synthetic steps, each one devouring resources and lowering the product yield. Shaped laser pulses can possibly act as photonic reagents and shorten the synthetic route toward the desired product. Chemical synthesis usually takes place in solution, and by including explicit solvent molecules in our quantum control simulations, we were able to identify their highly inhomogeneous influence on chemical reactions and how this affects potential quantum control. More important, we demonstrated for a synthetically relevant example that these complications can be overcome in theory, and laser pulses can be optimized to initiate the desired carbon-carbon bond formation. Putting this into context with the recently emerging concept of flow chemistry, which brings several practical advantages to the application of laser pulses, we want to encourage experimental groups to exploit this concept. Another path was opened by several additions to the commonly used laser pulse optimization algorithm (optimal control theory, OCT), several of which were developed in our group. The OCT algorithm as such is a purely mathematical optimization procedure, with no direct relation to experimental requirements. This means that usually the electric fields obtained out of OCT optimizations do not resemble laser pulses that can be achieved experimentally. However, the previously mentioned additions are aimed at closing the gap toward the experiment. In a recent quantum control study of our group, these algorithmic developments came to fruition. We were able to suggest a shaped laser pulse which can induce a long-living wave packet in the excited state of uracil. This might pave the way for novel experiments dedicated to investigating the formation of biological photodamage in DNA and RNA. The pulse we suggest is surprisingly simple because of the extended OCT algorithm and fulfills all criteria to be experimentally accessible.