Semiclassical multistate Liouville dynamics in the adiabatic representation

Diabatic and adiabatic representations: Electronic structure caveats

Simulations of electronically nonadiabatic processes may employ either the adiabatic or diabatic representation. Direct dynamics calculations are usually carried out in the adiabatic basis because the energy, force, and state coupling can be evaluated directly by many electronic structure methods. However, although its straightforwardness is appealing, direct dynamics is expensive when combined with quantitatively accurate electronic structure theories. This generates interest in analytically fitted surfaces to cut the expense, but the cuspidal ridges of the potentials and the singularities and vector nature of the couplings at high-dimensional, nonsymmetry-determined intersections in the adiabatic representation make accurate fitting almost impossible. This motivates using diabatic representations, where the surfaces are smooth and the couplings are also smooth and-importantly-scalar. In a recent previous work, we have developed a method called diabatization by deep neural network (DDNN) that takes advantage of the smoothness and nonuniqueness of diabatic bases to obtain them by machine learning. The diabatic potential energy matrices (DPEMs) learned by the DDNN method yield not only diabatic potential energy surfaces (PESs) and couplings in an analytic form useful for dynamics calculations, but also adiabatic surfaces and couplings in the adiabatic representation can be calculated inexpensively from the transformation. In the present work, we show how to extend the DDNN method to produce good approximations to global permutationally invariant adiabatic PESs simultaneously with DPEMs. The extended method is called permutationally restrained DDNN.

An energy-based method is proposed for the diabatization of the OH(2Pi)+F(2P)-->O(3P)+HF(1Sigma+) reaction. It is demonstrated that the diabatic representation obtained is regularized, i.e., the residual derivative couplings do not present singularities at the conical intersections appearing along the reaction path. This method only requires the knowledge of the 1,2 3A" and 1 3A' eigenvalues and does not require any adjustable parameter. Thus, many convergence problems arising in other derivative-based diabatization methods are avoided, and the description of the configuration space along the reaction path is enormously simplified. Three-dimensional coupled diabatic energy surfaces are obtained by an interpolation procedure using approximately 4000 accurate ab initio points. The angular resolved photodetachment cross sections are obtained in the diabatic and adiabatic representations using a wave packet method. An excellent agreement is obtained with recent experimental data [D. M. Neumark, Phys. Chem. Chem. Phys. 7, 433 (2005)] for high electron kinetic energies where only the triplet electronic states contribute.

The Born-Oppenheimer approximation assumes nuclear motion evolves on single, uncoupled potential energy surfaces, widely used to solve the time-independent Schrödinger equation for atomistic systems. However, for near-degenerate same-symmetry electronic states, avoided crossings in the potential energy curves occur and non-adiabatic couplings (NACs) become significant. In such cases, the adiabatic approximation is unsuitable for high-resolution spectroscopy. A unitary transformation to the diabatic representation can eliminate NACs, resulting in smooth molecular property curves that may cross. Computing this adiabatic-to-diabatic transformation (AtDT) is desirable but non-analytic for multi-state coupled systems, necessitating numerical solutions. It remains unclear if current methods yield numerically exact AtDTs ensuring rovibronic energy level equivalence between adiabatic and diabatic pictures. We demonstrate (for the first time) numerically exact equivalence of adiabatic and diabatic representations for -state diatomic molecules using abinitio data for , CH, and a model 10-state system. We show how the equivalence can be efficiently used to assess the importance of non-adiabatic effects and the impact of omitting them when computing rovibronic energies of diatomic molecules. The adiabatic and diabatic representations of the spectroscopic model, including all coupling terms, have been implemented in the diatomic code Duo.

The nonadiabatic transition probability 𝒫 between two electronic states is studied in the semiclassical approximation. Extensive use is made of the concept of global function introduced in the diabatic representation by Delos and Thorson (their t (s) function) and by ourselves in the adiabatic representation (the Massey function, T (S)). Consideration of the shapes of these global functions results in a classification of the different models which have been proposed earlier by Landau and Zener, by Rosen and Zener, by Delos and Thorson, by Nikitin, and by Bandrauk. The relationships which exist between all of these models are displayed. For all of them, the transition probability 𝒫 can be expressed as a simple analytic expression involving two parameters, ξ and ϑo. ξ is the Massey parameter. ϑo measures the initial decoupling, i.e., the component of the adiabatic states on the basis set of the diabatic states at the initial time (corresponding to the initial excitation in a photochemistry experiment). For a given value of ϑo, the different models predict relationships between 𝒫 and ξ that do not differ greatly, so that a kind of ’’universal function’’ relating 𝒫 and ξ can be defined. The previous considerations are applied to three specific examples of nonadiabatic coupling encountered in the case of H2, N2+ and CH2+. For each of them, the global functions t (s) and T (S) are calculated ab initio. The relationship between their shape and the nonadiabatic transition probability is discussed.

A semiclassical surface hopping expansion of the propagator is developed for a general representation of the “fast” variable quantum states. The representation can be the adiabatic or diabatic representation or any representation between these two. A particular representation is defined, which is optimal in the sense that it minimizes the integrated interstate coupling. The coupling is integrated over a suitable classical trajectory in this definition. Calculations for a simple one-dimensional curve crossing model problem show that the use of this optimal representation can significantly reduce the importance of multihop terms in the expansion. An approximation to this optimal representation is proposed, which is much simpler to implement numerically. Calculations for the model curve crossing problem demonstrate that this approximate optimal representation provides integrated couplings that are very close to those obtained for the optimal representation. These results suggest that this approximate optimal representation provides a computationally attractive representation for use with semiclassical surface hopping methods, when studying problems with curve crossings.

The vibrational spectroscopy of adsorbates is becoming an important investigation tool for catalysis and material science. This paper presents a semiclassical molecular dynamics method able to reproduce the vibrational energy levels of systems composed by molecules adsorbed on solid surfaces. Specifically, we extend our divide-and-conquer semiclassical method for power spectra calculations to gas-surface systems and interface it with plane-wave electronic structure codes. The Born-Oppenheimer classical dynamics underlying the semiclassical calculation is full dimensional, and our method includes not only the motion of the adsorbate but also those of the surface and the bulk. The vibrational spectroscopic peaks related to the adsorbate are accounted together with the most coupled phonon modes to obtain spectra amenable to physical interpretations. We apply the method to the adsorption of CO, NO, and H2O on the anatase-TiO2 (101) surface. We compare our semiclassical results with the single-point harmonic estimates and the classical power spectra obtained from the same trajectory employed in the semiclassical calculation. We find that CO and NO anharmonic effects of fundamental vibrations are similarly reproduced by the classical and semiclassical dynamics and that H2O adsorption is fully and properly described in its overtone and combination band relevant components only by the semiclassical approach.

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 S0, S1, and T1 states of formaldehyde are used as an example. First, the two singlet states S0 and S1 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.

An analytical polynomial formulation of the exponential of fully symmetric matrices is proposed. This polynomial formula is derived using the Cayley - Hamilton theorem and the polynomial Euclidean division. The application field of these developments is time-dependent wavepacket calculations involving multiple potential energy surface crossing. The use of this polynomial formulation in time-dependent calculations is directly connected to the choice of the split operator, proposed by Feit and Fleck (1982 J. Comput. Phys. 47 412) as a numerical propagation scheme. We show in this paper that this analytical polynomial formula is numerically more efficient than the scheme proposed by Almeida and Metiu (1988 Chem. Phys. Lett. 146 47) where the exponential of the kinetic operator is evaluated in the momentum and diabatic representation and the exponential of the potential operator is evaluated in the coordinate and adiabatic representation. The use of an optical potential ( negative imaginary potential) to absorb the wavepacket near the edges of the grid point is quite simple to implement as it appears as a multiplicative factor in front of the polynomial expression. Our scheme is exemplified on the photodissociation of a diatomic molecule involving three potential energy curves.

Exact nonadiabatic quantum evolution preserves many geometric properties of the molecular Hilbert space. In the first paper of this series ["Paper I," S. Choi and J. Vaníček, J. Chem. Phys. 150, 204112 (2019)], we presented numerical integrators of arbitrary-order of accuracy that preserve these geometric properties exactly even in the adiabatic representation, in which the molecular Hamiltonian is not separable into kinetic and potential terms. Here, we focus on the separable Hamiltonian in diabatic representation, where the split-operator algorithm provides a popular alternative because it is explicit and easy to implement, while preserving most geometric invariants. Whereas the standard version has only second-order accuracy, we implemented, in an automated fashion, its recursive symmetric compositions, using the same schemes as in Paper I, and obtained integrators of arbitrary even order that still preserve the geometric properties exactly. Because the automatically generated splitting coefficients are redundant, we reduce the computational cost by pruning these coefficients and lower memory requirements by identifying unique coefficients. The order of convergence and preservation of geometric properties are justified analytically and confirmed numerically on a one-dimensional two-surface model of NaI and a three-dimensional three-surface model of pyrazine. As for efficiency, we find that to reach a convergence error of 10-10, a 600-fold speedup in the case of NaI and a 900-fold speedup in the case of pyrazine are obtained with the higher-order compositions instead of the second-order split-operator algorithm. The pyrazine results suggest that the efficiency gain survives in higher dimensions.

The current paper reports a global investigation of all excited states below the ionic limit Fr+Na- of FrNa molecule following diabatic and adiabatic representations. The adiabatic and diabatic potential energy curves (PECs) for Σ+, Π, and Δ symmetries have been calculated for a dense grid of internuclear distances. The transition and permanent dipole moments (TDM and PDM) have been reported for both representations. Regarding the pseudopotential approach and full valence configuration interaction (FCI), the ab initio computation has been performed. Furthermore, the diabatization method was achieved by the use of the variational effective Hamiltonian theory (VEH). This latter served to assess the nonadiabatic coupling between the treated adiabatic states and to eliminate it employing a suitable unitary transformation matrix. The detailed computation of the PECs and PDM and TDM curves of the ground and excited states play a key role to coordinate experimental efforts in order to create cold and ultracold molecules in their ground stable rovibrational levels.

We present a comprehensive comparison of the following mixed quantum-classical methods for calculating electronic transition rates: (1) nonequilibrium Fermi's golden rule, (2) mixed quantum-classical Liouville method, (3) mean-field (Ehrenfest) mixed quantum-classical method, and (4) fewest switches surface-hopping method (in diabatic and adiabatic representations). The comparison is performed on the Garg-Onuchic-Ambegaokar benchmark charge-transfer model, over a broad range of temperatures and electronic coupling strengths, with different nonequilibrium initial states, in the normal and inverted regimes. Under weak to moderate electronic coupling, the nonequilibrium Fermi's golden rule rates are found to be in good agreement with the rates obtained via the mixed quantum-classical Liouville method that coincides with the fully quantum-mechanically exact results for the model system under study. Our results suggest that the nonequilibrium Fermi's golden rule can serve as an inexpensive yet accurate alternative to Ehrenfest and the fewest switches surface-hopping methods.

A theoretical description of multi-channel fragmentation of the Na 3+ cluster-ion in collision with a He target

Semiempirical potential energy surfaces for Br–HBr are constructed, based on empirical Kr–HBr, Kr–Br and Kr–Kr potentials. The electrostatic forces due to the interaction of the quadrupole of the Br atom with the multipoles of the HBr molecule are included, and make a substantial contribution to the well depths. The potentials are presented in diabatic and adiabatic representations, with and without spin–orbit coupling. When spin–orbit coupling is included, the adiabatic well depth is 342 cm−1, at a linear Br–H–Br configuration. Fully coupled bound-state calculations are carried out, and the ground state is found to be bound by 269 cm−1. If the complex can be observed spectroscopically, it will provide a sensitive probe of the potential energy surfaces in the entrance and exit valleys of the Br+HBr reaction. In addition, the complex is likely to be important as a product of photodissociation of HBr in HBr dimer.

The extended Lagrangian formulation of time-reversible Born-Oppenheimer molecular dynamics [A. M. N. Niklasson, C. J. Tymczak, and M. Challacombe, Phys. Rev. Lett. 100, 123004 (2008); Phys. Rev. Lett. 97, 123001 (2006)] enables the use of geometric integrators in the propagation of both the nuclear and the electronic degrees of freedom on the Born-Oppenheimer potential energy surface. Different symplectic integrators up to the sixth order have been adapted and optimized in the framework of ab initio self-consistent-field theory. It is shown how the accuracy can be significantly improved compared to a conventional Verlet integration at the same level of computational cost, in particular, for the case of very high accuracy requirements.