Full-dimensional multi-state simulation of the photodissociation of thioanisole

Fitting of Coupled Potential Energy Surfaces via Discovery of Companion Matrices by Machine Intelligence

Fitting coupled potential energy surfaces is a critical step in simulating electronically nonadiabatic chemical reactions and energy transfer processes. Analytic representation of coupled potential energy surfaces enables one to perform detailed dynamics calculations. Traditionally, fitting is performed in a diabatic representation to avoid fitting the cuspidal ridges of coupled adiabatic potential energy surfaces at conical intersection seams. In this work, we provide an alternative approach by carrying out fitting in the adiabatic representation using a modified version of the Frobenius companion matrices, whose usage was first proposed by Opalka and Domcke. Their work involved minimizing the errors in fits of the characteristic polynomial coefficients (CPCs) and diagonalizing the resulting companion matrix, whose eigenvalues are adiabatic potential energies. We show, however, that this may lead to complex eigenvalues and spurious discontinuities. To alleviate this problem, we provide a new procedure for the automatic discovery of CPCs and the diagonalization of a companion matrix by using a special neural network architecture. The method effectively allows analytic representation of global coupled adiabatic potential energy surfaces and their gradients with only adiabatic energy input and without experience-based selection of a diabatization scheme. We demonstrate that the new procedure, called the companion matrix neural network (CMNN), is successful by showing applications to LiH, H3, phenol, and thiophenol.

Parametrically Managed Activation Functions for Improved Global Potential Energy Surfaces for Six Coupled 5A' States and Fourteen Coupled 3A' States of O + O2

We report new potential energy surfaces for six coupled 5A' states and 14 coupled 3A' states of O3. The new surfaces are created by parametrically managed diabatization by deep neural network (PM-DDNN). The PM-DDNN method uses calculated adiabatic potential energy surfaces to discover and fit an underlying adiabatic-equivalent set of diabatic surfaces and their couplings and obtains the fit to the adiabatic surfaces by diagonalization of the diabatic potential energy matrix (DPEM). The procedure yields the adiabatic surfaces and their gradients, as well as the DPEM and its gradient. If desired one can also compute the nonadiabatic coupling due to the transformation. The present work improves on previous work by using a new coordinate to guide the decay of the neural network contribution to the many-body fit to the whole DPEM. The main objective was to obtain smoother potentials than the previous ones with better suitability for dynamics calculations, and this was achieved. Furthermore, we obtained suitably small deviations from the input reference data. For the six coupled 5A' surfaces, the 60,366 data below 10 eV are fit with a mean unsigned error (MUE) of 49 meV, and for the 14 coupled 3A' surfaces, the 76,733 data below 10 eV are fit with an MUE of 28 meV. The data below 5 eV fit even more accurately with MUEs of 37 meV (5A') and 20 meV (3A').

ChemPotPy: A Python Library for Analytic Representations of Potential Energy Surfaces and Diabatic Potential Energy Matrices

Constructing analytic representations of global and semiglobal potential energy surfaces is difficult and can be laborious, and it is even harder when one needs coupled potential energy surfaces and their electronically nonadiabatic couplings. When accomplished, however, the resulting potential functions are a valuable resource. To facilitate the convenient use of potentials that have been developed, we provide a collection of existing surfaces in a library with consistent units and formats. A potential energy surface library of this type, namely PotLib, was built more than 20 years ago. However, that library only provided pristine Fortran subroutines for each potential energy surface, and therefore, it is not as user-friendly as would be desirable. Here, we report the creation of ChemPotPy, a CHEMical library of POTential energy surfaces in PYthon. ChemPotPy is a user-friendly library for analytic representation of single-state and multistate potential energy surfaces and couplings. A given entry in the library contains an analytic potential energy function or analytic functions for a set of coupled potential energy surfaces, and depending on the case, it may also include analytic or numerical gradients, nonadiabatic coupling vectors, and/or diabatic potential energy matrices and their gradients. Only three inputs, namely, the chemical formula of the system, the name of the potential energy surface or surface set, and the Cartesian geometry, are required. ChemPotPy uses the same units for input and output quantities of all surfaces and surface sets to facilitate general interfaces with the dynamics programs. The initial version of the library contains 338 entries, and we anticipate that more will be added in the future.

Dynamic Role of the Intramolecular Hydrogen Bonding in the S1 State Relaxation Dynamics Revealed by the Direct Measurement of the Mode-Dependent Internal Conversion Rate of 2-Chlorophenol and 2-Chlorothiophenol

The dynamic role of the intramolecular hydrogen bond in the S1 relaxation of cis-2-chlorophenol (2-CP) or cis-2-chlorothiophenol (2-CTP) has been investigated in a state-specific manner. Whereas ultrafast internal conversion is dominant for 2-CP, the H-tunneling competes with internal conversion for 2-CTP even at the S1 origin. The S0-S1 internal conversion rate of 2-CTP could be directly measured from the S1 lifetimes of 2-CTP-d1 (Cl-C6H4-SD) as the D-tunneling is kinetically blocked, allowing distinct estimations of tunneling and internal conversion rates with increasing the energy. The internal conversion rate of 2-CTP increases by two times at the out-of-plane torsional mode excitation, suggesting that the internal conversion is facilitated at the nonplanar geometry. It then sharply increases at ∼600 cm-1, indicating that the S1/S0 conical intersection is readily accessible at the extended C-Cl bond length. The strength of the intramolecular hydrogen bond should be responsible for the distinct dynamic behaviors of 2-CP and 2-CTP.

πσ*-Mediated Nonadiabatic Tunneling Dynamics of Thiophenols in S1: The Semiclassical Approaches

The S-H bond tunneling predissociation dynamics of thiophenol and its ortho-substituted derivatives (2-fluorothiophenol, 2-methoxythiophenol, and 2-chlorothiphenol) in S₁ (ππ*) where the H atom tunneling is mediated by the nearby S₂ (πσ*) state (which is repulsive along the S-H bond extension coordinate) have been investigated in a state-specific way using the picosecond time-resolved pump-probe spectroscopy for the jet-cooled molecules. The effects of the specific vibrational mode excitations and the SH/SD substitutions on the S-H(D) bond rupture tunneling dynamics have been interrogated, giving deep insights into the multidimensional aspects of the S₁/S₂ conical intersection, which also shapes the underlying adiabatic tunneling potential energy surfaces (PESs). The semiclassical tunneling rate calculations based on the Wentzel-Kramers-Brillouin (WKB) approximation or Zhu-Nakamura (ZN) theory have been carried out based on the ab initio PESs calculated in the (one, two, or three) reduced dimensions to be compared with the experiment. Though the quantitative experimental results could not be reproduced satisfactorily by the present calculations, the qualitative trends among different molecules in terms of the behavior of the tunneling rate versus the (adiabatic) barrier height or the number of PES dimensions could be rationalized. Most interestingly, the H/D kinetic isotope effect observed in the tunneling rate could be much better explained by the ZN theory compared to the WKB approximation, indicating that the nonadiabatic coupling matrix elements should be invoked for understanding the tunneling dynamics taking place in the proximity of the conical intersection.

Photoionization Spectroscopic and Theoretical Study on the Molecular Structures of cis- and trans-3-Chlorothioanisole

Resonance-enhanced two-photon ionization (R2PI) and mass-analyzed threshold ionization (MATI) spectra are measured for the cis- and trans-3-chlorothioanisole (3ClTA). The first electronic excitation energy (E 1) and the adiabatic ionization energy (IE) of the cis-rotamer are determined to be 33 959±3 and 65 326±5 cm-1, respectively, and those of the trans-rotamer are determined to be 34102±3 and 65 471±5 cm-1, respectively. Density functional theory (DFT) calculations confirm that both the cis- and trans-rotamers of 3ClTA are stable and coexist in their respective S0, S1, and D0 states. Both rotamers adopt planar structures with cis- being slightly more stable than trans- in the respective S0, S1, and D0 states. The conformation, substitution, and isotope effects on the molecular structure, active vibrations, and electronic transition and ionization energies of 3ClTA are analyzed.

Nonadiabatic Dynamics Algorithms with Only Potential Energies and Gradients: Curvature-Driven Coherent Switching with Decay of Mixing and Curvature-Driven Trajectory Surface Hopping

Direct dynamics by mixed quantum-classical nonadiabatic methods is an important tool for understanding processes involving multiple electronic states. Very often, the computational bottleneck of such direct simulation comes from electronic structure theory. For example, at every time step of a trajectory, nonadiabatic dynamics requires potential energy surfaces, their gradients, and the matrix elements coupling the surfaces. The need for the couplings can be alleviated by employing the time derivatives of the wave functions, which can be evaluated from overlaps of electronic wave functions at successive time steps. However, evaluation of overlap integrals is still expensive for large systems. In addition, for electronic structure methods for which the wave functions or the coupling matrix elements are not available, nonadiabatic dynamics algorithms become inapplicable. In this work, building on recent work by Baeck and An, we propose new nonadiabatic dynamics algorithms that only require adiabatic potential energies and their gradients. The new methods are named curvature-driven coherent switching with decay of mixing (κCSDM) and curvature-driven trajectory surface hopping (κTSH). We show how powerful these new methods are in terms of computation time and accuracy as compared to previous mixed quantum-classical nonadiabatic dynamics algorithms. The lowering of the computational cost will allow longer nonadiabatic trajectories and greater ensemble averaging to be affordable, and the ability to calculate the dynamics without electronic structure coupling matrix elements extends the dynamics capability to new classes of electronic structure methods.

Femtosecond Wavepacket Dynamics Reveals the Molecular Structures in the Excited (S1) and Cationic (D0) States

Molecular structures in the electronically excited (S₁) and cationic (D₀) states of 2-fluorothioanisole (2-FTA) have been precisely refined from the real-time dynamics of the femtosecond (fs) wavepacket prepared by the coherent excitation of the Franck-Condon active out-of-plane torsional modes in the S₁ ← S₀ transition at 285 nm. The simulation to reproduce the experiment in terms of the beating frequencies gives the nonplanar geometry of 2-FTA in S₁, where the out-of-plane dihedral angle (φ) of the S-CH₃ moiety is 51° with respect to the molecular plane. The behavior of the fs wavepacket in terms of the amplitudes and phases with the change of the probe (ionization) wavelength (λ_probe = 300-330 nm) provides the otherwise veiled structure of the cationic D₀ state. While the 2-FTA cation adopts the planar geometry (φ = 0°) at the global minimum, it is found to have a vertical minimum at φ ≈ 135° from the perspective of the D₀ ← S₁ vertical transition. Ab initio calculations support the experiment quite well although the simulation using the model potentials could improve the match with the experiment, giving the new interpretation for the previously disputed photoelectron spectroscopic results.

Conservation of Angular Momentum in Direct Nonadiabatic Dynamics

Direct nonadiabatic dynamics is used to study processes involving multiple electronic states from small molecules to materials. Compared with dynamics with fitted analytical potential energy surfaces, direct dynamics is more user-friendly in that it obtains all needed energies, gradients, and nonadiabatic couplings (NACs) by electronic structure calculations. However, the NAC that is usually used does not conserve angular momentum or the center of mass in widely used mixed quantum-classical nonadiabatic dynamics algorithms, in particular, trajectory surface hopping, semiclassical Ehrenfest, and coherent switching with decay of mixing. We show that by using a projection operator to remove the translational and rotational components of the originally computed NAC, one can restore the conservation.

Multidimensional characterization of the conical intersection seam in the normal mode space

Multidimensional conical intersection seam has been characterized by utilizing the dynamic resonances in the nonadiabatic transition probability experimentally observed in the predissociation of thioanisole isotopomers.

Quantifying rival bond fission probabilities following photoexcitation: C-S bond fission in t-butylmethylsulfide

We illustrate a new, collision-free experimental strategy that allows determination of the absolute probabilities of rival bond fission processes in a photoexcited molecule - here t-butylmethylsulfide (BSM). The method combines single photon ('universal') ionization laser probe methods, simultaneous imaging of all probed fragments (multi-mass ion imaging) and the use of an appropriate internal calibrant (here dimethylsulfide). Image analysis allows quantification of the dynamics of the rival B-SM and BS-M bond fission processes following ultraviolet (UV) excitation of BSM and shows the former to be twice as probable, despite the only modest (∼2%) differences in the respective ground state equilibrium C-S bond lengths or bond strengths. Rationalising this finding should provide a stringent test of the two close-lying, coupled excited states of 1A'' symmetry accessed by UV excitation in BSM and related thioethers, of the respective transition dipole moment surfaces, and of the geometry dependent non-adiabatic couplings that enable the rival C-S bond fissions.

Experimental observation of nonadiabatic bifurcation dynamics at resonances in the continuum

The surface crossing of bound and unbound electronic states in multidimensional space often gives rise to resonances in the continuum. This situation happens in the πσ*-mediated photodissociation reaction of 2-fluorothioanisole; optically-bright bound S1 (ππ*) vibrational states of 2-fluorothioanisole are strongly coupled to the optically-dark S2 (πσ*) state, which is repulsive along the S-CH3 elongation coordinate. It is revealed here that the reactive flux prepared at such resonances in the continuum bifurcates into two distinct reaction pathways with totally different dynamics in terms of energy disposal and nonadiabatic transition probability. This indicates that the reactive flux in the Franck-Condon region may either undergo nonadiabatic transition funneling through the conical intersection from the upper adiabat, or follow a low-lying adiabatic path, along which multiple dynamic saddle points may be located. Since 2-fluorothioanisole adopts a nonplanar geometry in the S1 minimum energy, the quasi-degenerate S1/S2 crossing seam in the nonplanar geometry, which lies well below the planar S1/S2 conical intersection, is likely responsible for the efficient vibronic coupling, especially in the low S1 internal energy region. As the excitation energy increases, bound-to-continuum coupling is facilitated with the aid of intramolecular vibrational redistribution, along many degrees of freedom spanning the large structural volume. This leads to the rapid domination of the continuum character of the reactive flux. This work reports direct and robust experimental observations of the nonadiabatic bifurcation dynamics of the reactive flux occurring at resonances in the continuum of polyatomic molecules.