Neural network based quasi-diabatic Hamiltonians with symmetry adaptation and a correct description of conical intersections

Neural network potentials facilitating accurate complex scaling for molecular resonances: from a model to high dimensional realistic systems

Here we propose a neural network based complex scaling (NN-CS) method for computing the complex eigenvalues (Er-iΓ/2) of molecular resonances, in which the CS of the potential part in the non-Hermitian Hamiltonian is effectively achieved by NNs. Taking a two-dimensional (2D) diabatic model including two states coupled by the conical intersection for example, the NN-CS method is shown to reproduce the eigenvalues of the resonance states quite well. Subsequently, this NN-CS method with a 2D Hamiltonian model is utilized to compute the vibronic resonances in the 1nσ*-mediated photodissociation of thioanisole based on a new NN diabatic potential energy matrix. The calculated lifetimes of the vibronic resonances are found to be in good agreement with other theoretical results and available experimental data. Finally, the NN-CS method is applied to treat a much more challenging system, namely, the resonances in the six-dimensional (6D) photodissociation continuum of NH3, due to its high dimensionalities and all three dissociative coordinates needing to be scaled in the complex scaling of the potential part. Again, the calculated energy positions and widths of the 6D resonances by the NN-CS method agree well with other theoretical results. Our calculations show that the NN-CS method is able to accurately treat the vibronic resonances involving multiple coupled electronic states and resonances in high dimensional realistic systems.

Compensation States Approach in the Hybrid Diabatization Scheme: Extension to Multidimensional Data and Properties

The diabatization of reactive systems for more than just a couple of states is a very demanding problem and generally requires advanced diabatization techniques. Especially for dissociative processes, the drastic changes in the adiabatic wave functions often would require large diabatic state bases, which quickly become impractical. Recently, we addressed this problem by the compensation states approach developed in the context of our hybrid diabatization scheme. This scheme utilizes wave function as well as energy data in combination with a diabatic potential model. In regions where the initial diabatic state basis becomes insufficient for an appropriate representation of the adiabatic states, new model states are generated. The new model states compensate for the state space not spanned by the initial diabatic basis. Such a compensation state is obtained by projecting the initial diabatic state space out of the adiabatic wave function. This yields a very efficient basis representation of the electronic Hamiltonian. The present work presents two new aspects. First, it is shown how other operators like the spin-orbit operator in the framework of the Effective Relativistic Coupling by Asymptotic Representation (ERCAR) can be evaluated in this compact model state space without losing the correct wave function information and accuracy. Second, the extension of the approach to multidimensional potential energy surface models is presented for methyl iodide including the C-I dissociation coordinate and the angular H3C-I bending coordinates.

Accurate diabatization based on combined-hyperbolic-inverse-power-representation: 1,2 2A' states of BeH2

A diabatic potential energy matrix (DPEM) for the two lowest states of BeH2+ has been constructed using the combined-hyperbolic-inverse-power-representation (CHIPR) method. By imposing symmetry constraints on the coefficients of polynomials, the complete nuclear permutation inversion symmetry is correctly preserved in the CHIPR functional form. The symmetrized CHIPR functional form is then used in the diabatization by ansatz procedure. The ab initio energies are reproduced with satisfactory accuracy. In addition, the CHIPR-based DPEM also reproduces the local topology of a conical intersection. Future work will focus on a complete four-state diabatic representation with emphasis on the long-range interactions and spin-orbit couplings, which will enable accurate quantum scattering calculations for the Be+(2P) + H2 → BeH+(X1Σ+) + H(2S) reaction.

Quantum Dynamics of Photodissociation: Recent Advances and Challenges

Recent advances in constructing accurate potential energy surfaces and nonadiabatic couplings from high-level ab initio data have revealed detailed potential landscapes in not only the ground electronic state but also excited ones. They enabled quantitatively accurate characterization of photoexcited reactive systems using quantum mechanical methods. In this Perspective, we survey the recent progress in quantum mechanical studies of adiabatic and nonadiabatic photodissociation dynamics, focusing on initial state control and product energy disposal. These new insights helped to understand quantum effects in small prototypical systems, and the results serve as benchmarks for developing more approximate theoretical methods.

Parametrically Managed Activation Function for Fitting a Neural Network Potential with Physical Behavior Enforced by a Low-Dimensional Potential

Machine-learned representations of potential energy surfaces generated in the output layer of a feedforward neural network are becoming increasingly popular. One difficulty with neural network output is that it is often unreliable in regions where training data is missing or sparse. Human-designed potentials often build in proper extrapolation behavior by choice of functional form. Because machine learning is very efficient, it is desirable to learn how to add human intelligence to machine-learned potentials in a convenient way. One example is the well-understood feature of interaction potentials that they vanish when subsystems are too far separated to interact. In this article, we present a way to add a new kind of activation function to a neural network to enforce low-dimensional constraints. In particular, the activation function depends parametrically on all of the input variables. We illustrate the use of this step by showing how it can force an interaction potential to go to zero at large subsystem separations without either inputting a specific functional form for the potential or adding data to the training set in the asymptotic region of geometries where the subsystems are separated. In the process of illustrating this, we present an improved set of potential energy surfaces for the 14 lowest ³A' states of O₃. The method is more general than this example, and it may be used to add other low-dimensional knowledge or lower-level knowledge to machine-learned potentials. In addition to the O₃ example, we present a greater-generality method called parametrically managed diabatization by deep neural network (PM-DDNN) that is an improvement on our previously presented permutationally restrained diabatization by deep neural network (PR-DDNN).

Constructing Diabatic Potential Energy Matrices with Neural Networks Based on Adiabatic Energies and Physical Considerations: Toward Quantum Dynamic Accuracy

A permutation invariant polynomial-neural network (PIP-NN) approach for constructing the global diabatic potential energy matrices (PEMs) of the coupled states of molecules is proposed. Specifically, the diabatization scheme is based merely on the adiabatic energy data of the system, which is ideally a most convenient way due to not requiring additional ab initio calculations for the data of the derivative coupling or any other physical properties of the molecule. Considering the permutation and coupling characteristics of the system, particularly in the presence of conical intersections, some vital treatments for the off-diagonal terms in diabatic PEM are essentially needed. Taking the photodissociation of H₂O(X~/B~)/NH₃(X~/A~) and nonadiabatic reaction Na(3p) + H₂ → NaH(Σ⁺) + H for example, this PIP-NN method is shown to build up the global diabatic PEMs effectively and accurately. The root-mean-square errors of the adiabatic potential energies in the fitting for three different systems are all small (<10 meV). Further quantum dynamic calculations show that the absorption spectra and product branching ratios in both H₂O(X~/B~) and NH₃(X~/A~) nonadiabatic photodissociation are well reproduced on the new diabatic PEMs, and the nonadiabatic reaction probability of Na(3p) + H₂ → NaH(Σ⁺) + H obtained on the new diabatic PEMs of the 1²A₁ and 1²B₂ states is in reasonably good agreement with previous theoretical result as well, validating this new PIP-NN method.

Accurate Full-Dimensional Global Diabatic Potential Energy Matrix for the Two Lowest-Lying Electronic States of the H + O2 ↔ HO + O Reaction

A new and more accurate diabatic potential energy matrix (DPEM) is developed for the two lowest-lying electronic states of HO₂, covering both the strong interaction region and reaction asymptotes. The ab initio calculations were performed at the Davidson corrected multireference configuration interaction level with the augmented correlation-consistent polarized valence quintuple-zeta basis set (MRCI+Q/AV5Z). The accuracy of the electronic structure calculations is validated by excellent agreement with the experimental HO₂ equilibrium geometry, fundamental vibrational frequencies, and H + O₂ ↔ OH + O reaction energy. Through the combination of an electronic angular momentum-method and a configuration interaction vector-based method, the mixing angle between the first two ²A″ states of HO₂ was successfully determined. Elements of the 2×2 DPEM were fit to neural networks with a proper account of the complete nuclear permutation inversion symmetry of HO₂. The DPEM correctly predicted the properties of conical intersection seams at linear and T-shape geometries, thus providing a reliable platform for studying both the spectroscopy of HO₂ and the nonadiabatic dynamics for the H + O₂ ↔ OH + O reaction.

Permutation invariant polynomial neural network based diabatic ansatz for the (E + A) × (e + a) Jahn-Teller and Pseudo-Jahn-Teller systems

In this work, the permutation invariant polynomial neural network (PIP-NN) approach is employed to construct a quasi-diabatic Hamiltonian for system with non-Abelian symmetries. It provides a flexible and compact NN-based diabatic ansatz from the related approach of Williams, Eisfeld, and co-workers. The example of H₃ ⁺ is studied, which is an (E + A) × (e + a) Jahn-Teller and Pseudo-Jahn-Teller system. The PIP-NN diabatic ansatz is based on the symmetric polynomial expansion of Viel and Eisfeld, the coefficients of which are expressed with neural network functions that take permutation-invariant polynomials as input. This PIP-NN-based diabatic ansatz not only preserves the correct symmetry but also provides functional flexibility to accurately reproduce ab initio electronic structure data, thus resulting in excellent fits. The adiabatic energies, energy gradients, and derivative couplings are well reproduced. A good description of the local topology of the conical intersection seam is also achieved. Therefore, this diabatic ansatz completes the PIP-NN based representation of DPEM with correct symmetries and will enable us to diabatize even more complicated systems with complex symmetries.

Diabatic States of Molecules

Quantitative simulations of electronically nonadiabatic molecular processes require both accurate dynamics algorithms and accurate electronic structure information. Direct semiclassical nonadiabatic dynamics is expensive due to the high cost of electronic structure calculations, and hence it is limited to small systems, limited ensemble averaging, ultrafast processes, and/or electronic structure methods that are only semiquantitatively accurate. The cost of dynamics calculations can be made manageable if analytic fits are made to the electronic structure data, and such fits are most conveniently carried out in a diabatic representation because the surfaces are smooth and the couplings between states are smooth scalar functions. Diabatic representations, unlike the adiabatic ones produced by most electronic structure methods, are not unique, and finding suitable diabatic representations often involves time-consuming nonsystematic diabatization steps. The biggest drawback of using diabatic bases is that it can require large amounts of effort to perform a globally consistent diabatization, and one of our goals has been to develop methods to do this efficiently and automatically. In this Feature Article, we introduce the mathematical framework of diabatic representations, and we discuss diabatization methods, including adiabatic-to-diabatic transformations and recent progress toward the goal of automatization.

High-fidelity first principles nonadiabaticity: diabatization, analytic representation of global diabatic potential energy matrices, and quantum dynamics

Nonadiabatic dynamics, which goes beyond the Born-Oppenheimer approximation, has increasingly been shown to play an important role in chemical processes, particularly those involving electronically excited states. Understanding multistate dynamics requires rigorous quantum characterization of both electronic and nuclear motion. However, such first principles treatments of multi-dimensional systems have so far been rather limited due to the lack of accurate coupled potential energy surfaces and difficulties associated with quantum dynamics. In this Perspective, we review recent advances in developing high-fidelity analytical diabatic potential energy matrices for quantum dynamical investigations of polyatomic uni- and bi-molecular nonadiabatic processes, by machine learning of high-level ab initio data. Special attention is paid to methods of diabatization, high fidelity construction of multi-state coupled potential energy surfaces and property surfaces, as well as quantum mechanical characterization of nonadiabatic nuclear dynamics. To illustrate the tremendous progress made by these new developments, several examples are discussed, in which direct comparison with quantum state resolved measurements led to either confirmation of the observation or sometimes reinterpretation of the experimental data. The insights gained in these prototypical systems greatly advance our understanding of nonadiabatic dynamics in chemical systems.

Machine Learning Approach to Calculate Electronic Couplings between Quasi-diabatic Molecular Orbitals: The Case of DNA

Diabatization of one-electron states in flexible molecular aggregates is a great challenge due to the presence of surface crossings between molecular orbital (MO) levels and the complex interaction between MOs of neighboring molecules. In this work, we present an efficient machine learning approach to calculate electronic couplings between quasi-diabatic MOs without the need of nonadiabatic coupling calculations. Using MOs of rigid molecules as references, the MOs that can be directly regarded to be quasi-diabatic in molecular dynamics are selected out, state tracked, and phase corrected. On the basis of this information, artificial neural networks are trained to characterize the structure-dependent onsite energies of quasi-diabatic MOs and the intermolecular electronic couplings. A representative sequence of DNA is systematically studied as an illustration. Smooth time evolution of electronic couplings in all base pairs is obtained with quasi-diabatic MOs. In particular, our method can calculate electronic couplings between different quasi-diabatic MOs independently, and thus, this possesses unique advantages in many applications.

Ab initio symmetric quasi-classical approach to investigate molecular Tully models

We perform on-the-fly non-adiabatic molecular dynamics simulations using the symmetrical quasi-classical (SQC) approach with the recently suggested molecular Tully models: ethylene and fulvene. We attempt to provide benchmarks of the SQC methods using both the square and triangle windowing schemes as well as the recently proposed electronic zero-point-energy correction scheme (the so-called γ correction). We use the quasi-diabatic propagation scheme to directly interface the diabatic SQC methods with adiabatic electronic structure calculations. Our results showcase the drastic improvement of the accuracy by using the trajectory-adjusted γ-corrections, which outperform the widely used trajectory surface hopping method with decoherence corrections. These calculations provide useful and non-trivial tests to systematically investigate the numerical performance of various diabatic quantum dynamics approaches, going beyond simple diabatic model systems that have been used as the major workhorse in the quantum dynamics field. At the same time, these available benchmark studies will also likely foster the development of new quantum dynamics approaches based on these techniques.

Machine Learning for Electronically Excited States of Molecules

Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.

Machine Learning for Electronically Excited States of Molecules

Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.

Molecular excited states through a machine learning lens

Theoretical simulations of electronic excitations and associated processes in molecules are indispensable for fundamental research and technological innovations. However, such simulations are notoriously challenging to perform with quantum mechanical methods. Advances in machine learning open many new avenues for assisting molecular excited-state simulations. In this Review, we track such progress, assess the current state of the art and highlight the critical issues to solve in the future. We overview a broad range of machine learning applications in excited-state research, which include the prediction of molecular properties, improvements of quantum mechanical methods for the calculations of excited-state properties and the search for new materials. Machine learning approaches can help us understand hidden factors that influence photo-processes, leading to a better control of such processes and new rules for the design of materials for optoelectronic applications.

How important are the residual nonadiabatic couplings for an accurate simulation of nonadiabatic quantum dynamics in a quasidiabatic representation?

Diabatization of the molecular Hamiltonian is a standard approach to remove the singularities of nonadiabatic couplings at conical intersections of adiabatic potential energy surfaces. In general, it is impossible to eliminate the nonadiabatic couplings entirely-the resulting "quasidiabatic" states are still coupled by smaller but nonvanishing residual nonadiabatic couplings, which are typically neglected. Here, we propose a general method for assessing the validity of this potentially drastic approximation by comparing quantum dynamics simulated either with or without the residual couplings. To make the numerical errors negligible to the errors due to neglecting the residual couplings, we use the highly accurate and general eighth-order composition of the implicit midpoint method. The usefulness of the proposed method is demonstrated on nonadiabatic simulations in the cubic Jahn-Teller model of nitrogen trioxide and in the induced Renner-Teller model of hydrogen cyanide. We find that, depending on the system, initial state, and employed quasidiabatization scheme, neglecting the residual couplings can result in wrong dynamics. In contrast, simulations with the exact quasidiabatic Hamiltonian, which contains the residual couplings, always yield accurate results.

Enabling complete multichannel nonadiabatic dynamics: A global representation of the two-channel coupled, 1,21A and 13A states of NH3 using neural networks

Global coupled three-state two-channel potential energy and property/interaction (dipole and spin-orbit coupling) surfaces for the dissociation of NH₃(Ã) into NH + H₂ and NH₂ + H are reported. The permutational invariant polynomial-neural network approach is used to simultaneously fit and diabatize the electronic Hamiltonian by fitting the energies, energy gradients, and derivative couplings of the two coupled lowest-lying singlet states as well as fitting the energy and energy gradients of the lowest-lying triplet state. The key issue in fitting property matrix elements in the diabatic basis is that the diabatic surfaces must be smooth, that is, the diabatization must remove spikes in the original adiabatic property surfaces attributable to the switch of electronic wavefunctions at the conical intersection seam. Here, we employ the fit potential energy matrix to transform properties in the adiabatic representation to a quasi-diabatic representation and remove the discontinuity near the conical intersection seam. The property matrix elements can then be fit with smooth neural network functions. The coupled potential energy surfaces along with the dipole and spin-orbit coupling surfaces will enable more accurate and complete treatment of optical transitions, as well as nonadiabatic internal conversion and intersystem crossing.

A fundamental invariant-neural network representation of quasi-diabatic Hamiltonians for the two lowest states of H3

The FI-NN approach is capable of representing highly accurate diabatic PESs with particular and complicated symmetry problems.

Advances and New Challenges to Bimolecular Reaction Dynamics Theory

Dynamics of bimolecular reactions in the gas phase are of foundational importance in combustion, atmospheric chemistry, interstellar chemistry, and plasma chemistry. These collision-induced chemical transformations are a sensitive probe of the underlying potential energy surface(s). Despite tremendous progress in past decades, our understanding is still not complete. In this Perspective, we survey the recent advances in theoretical characterization of bimolecular reaction dynamics, stimulated by new experimental observations, and identify key new challenges.

Neural Network Potential Energy Surfaces for Small Molecules and Reactions

We review progress in neural network (NN)-based methods for the construction of interatomic potentials from discrete samples (such as ab initio energies) for applications in classical and quantum dynamics including reaction dynamics and computational spectroscopy. The main focus is on methods for building molecular potential energy surfaces (PES) in internal coordinates that explicitly include all many-body contributions, even though some of the methods we review limit the degree of coupling, due either to a desire to limit computational cost or to limited data. Explicit and direct treatment of all many-body contributions is only practical for sufficiently small molecules, which are therefore our primary focus. This includes small molecules on surfaces. We consider direct, single NN PES fitting as well as more complex methods that impose structure (such as a multibody representation) on the PES function, either through the architecture of one NN or by using multiple NNs. We show how NNs are effective in building representations with low-dimensional functions including dimensionality reduction. We consider NN-based approaches to build PESs in the sums-of-product form important for quantum dynamics, ways to treat symmetry, and issues related to sampling data distributions and the relation between PES errors and errors in observables. We highlight combinations of NNs with other ideas such as permutationally invariant polynomials or sums of environment-dependent atomic contributions, which have recently emerged as powerful tools for building highly accurate PESs for relatively large molecular and reactive systems.

Exclusive Neural Network Representation of the Quasi-Diabatic Hamiltonians Including Conical Intersections

We propose a numerically simple and straightforward, yet accurate and efficient neural networks-based fitting strategy to construct coupled potential energy surfaces (PESs) in a quasi-diabatic representation. The fundamental invariants are incorporated to account for the complete nuclear permutation inversion symmetry. Instead of derivative couplings or interstate couplings, a so-called modified derivative coupling term is fitted by neural networks, resulting in accurate description of near degeneracy points, such as the conical intersections. The adiabatic energies, energy gradients, and derivative couplings are well reproduced, and the vanishing of derivative couplings as well as the isotropic topography of adiabatic and diabatic energies in asymptotic regions are automatically satisfied. All of these features of the coupled global PESs are requisite for accurate dynamics simulations. Our approach is expected to be very useful in developing highly accurate coupled PESs in a quasi-diabatic representation in an efficient machine learning-based way.

High-Fidelity Potential Energy Surfaces for Gas-Phase and Gas-Surface Scattering Processes from Machine Learning

In this Perspective, we review recent advances in constructing high-fidelity potential energy surfaces (PESs) from discrete ab initio points, using machine learning tools. Such PESs, albeit with substantial initial investments, provide significantly higher efficiency than direct dynamics methods and/or high accuracy at a level that is not affordable by on-the-fly approaches. These PESs not only are a necessity for quantum dynamical studies because of delocalization of wave packets but also enable the study of low-probability and long-time events in (quasi-)classical treatments. Our focus here is on inelastic and reactive scattering processes, which are more challenging than bound systems because of the involvement of continua. Relevant applications and developments for dynamical processes in both the gas phase and at gas-surface interfaces are discussed.

Combining SchNet and SHARC: The SchNarc Machine Learning Approach for Excited-State Dynamics

In recent years, deep learning has become a part of our everyday life and is revolutionizing quantum chemistry as well. In this work, we show how deep learning can be used to advance the research field of photochemistry by learning all important properties-multiple energies, forces, and different couplings-for photodynamics simulations. We simplify such simulations substantially by (i) a phase-free training skipping costly preprocessing of raw quantum chemistry data; (ii) rotationally covariant nonadiabatic couplings, which can either be trained or (iii) alternatively be approximated from only ML potentials, their gradients, and Hessians; and (iv) incorporating spin-orbit couplings. As the deep-learning method, we employ SchNet with its automatically determined representation of molecular structures and extend it for multiple electronic states. In combination with the molecular dynamics program SHARC, our approach termed SchNarc is tested on two polyatomic molecules and paves the way toward efficient photodynamics simulations of complex systems.

Construction of Quasi-diabatic Hamiltonians That Accurately Represent ab Initio Determined Adiabatic Electronic States Coupled by Conical Intersections for Systems on the Order of 15 Atoms. Application to Cyclopentoxide Photoelectron Detachment in the Full 39 Degrees of Freedom

We present, for systems of moderate dimension, a fitting framework to construct quasi-diabatic Hamiltonians that accurately represent ab initio adiabatic electronic structure data including the effects of conical intersections. The framework introduced here minimizes the difference between the fit prediction and the ab initio data obtained in the adiabatic representation, which is singular at a conical intersection seam. We define a general and flexible merit function to allow arbitrary representations and propose a representation to measure the fit-ab initio difference at geometries near electronic degeneracies. A fit Hamiltonian may behave poorly in insufficiently sampled regions, in which case a machine learning theory analysis of the fit representation suggests a regularization to address the deficiency. Our fitting framework including the regularization is used to construct the full 39-dimensional coupled diabatic potential energy surfaces for cyclopentoxy relevant to cyclopentoxide photoelectron detachment.

The generality of the GUGA MRCI approach in COLUMBUS for treating complex quantum chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

Accurate Neural Network Representation of the Ab Initio Determined Spin-Orbit Interaction in the Diabatic Representation Including the Effects of Conical Intersections

A method for fitting ab initio determined spin-orbit coupling interactions, in the Breit-Pauli approximation, based on quasidiabatic representations using neural network fits is reported. The algorithm generalizes our recently reported neural network approach for representing the dipole interaction. The S₀, S₁, and T₁ states of formaldehyde are used as an example. First, the two singlet states S₀ and S₁ are diabatized with a modified Boys Localization diabatization method. Second, the spin-orbit coupling between singlet and triplet states is transformed to the diabatic representation. This removes the discontinuities in the adiabatic representation. The diabatized spin-orbit couplings are then fit with smooth neural network functions. The analytic representation of spin-orbit coupling interactions in a diabatic basis by neural networks will make accurate full-dimensional quantum dynamical treatment of both internal conversion and intersystem crossing possible, which will help us to gain better understanding of both processes.