Conical intersection seams in spin-orbit coupled systems with an even number of electrons: A numerical study based on neural network fit surfaces.

Abstract

In this work, we consider the existence and topography of seams of conical intersections (CIs) for two key singlet-triplet systems, including a uniformly scaled spin-orbit interaction. The basic one triplet and one singlet state system denoted as (S0,T1) and the two singlets and one triplet system denoted as (S0,S1,T1) are treated. Essential to this analysis are realistic electronic structure data taken from a recently reported neural network fit for the 1,21A and 13A states of NH3, including Hsf (spin-free) and Hso (spin-orbit) surfaces derived from high quality ab initio wavefunctions. Three types of seams for the (S0,S1,T1) system are reported, which depend on the choice of the electronic Hamiltonian, He. The nonrelativistic CI seam [He = Hsf, (S0,S1)], the energy minimized nonrelativistic singlet-triplet intersection seam [He = Hsf, (S0,T1)], and the fully relativistic seam in the spin-diabatic representation (He = Htot = Hsf + Hso) are reported as functions of R(N-H). The derivative couplings are computed using He = Htot and Hsf from the fit data. The line integral of the derivative coupling is employed to juxtapose the geometric phase in the relativistic, He = Htot, and nonrelativistic, He = Hsf, cases. It is found for the (S0,T1) system that there is no CI in the spin-adiabatic representation, while for the (S0,S1,T1) system, CI can only be formed for two pairs of spin-adiabatic electronic states. The geometric phase effect thus needs to be handled with care when it comes to spin-nonconserving dynamics simulations.