Communication: Theoretical prediction of the importance of the 3B2 state in the dynamics of sulfur dioxide

Photodissociation dynamics of SO2 via the G̃1B1 state: The O(1D2) and O(1S0) product channels

Produced by both nature and human activities, sulfur dioxide (SO2) is an important species in the earth's atmosphere. SO2 has also been found in the atmospheres of other planets and satellites in the solar system. The photoabsorption cross sections and photodissociation of SO2 have been studied for several decades. In this paper, we reported the experimental results for photodissociation dynamics of SO2 via the G̃1B1 state. By analyzing the images from the time-sliced velocity map ion imaging method, the vibrational state population distributions and anisotropy parameters were obtained for the O(1D2) + SO(X3Σ-, a1Δ, b1Σ+) and O(1S0) + SO(X3Σ-) channels, and the branching ratios for the channels O(1D2) + SO(X3Σ-), O(1D2) + SO(a1Δ), and O(1D2) + SO(b1Σ+) were determined to be ∼0.3, ∼0.6, and ∼0.1, respectively. The SO products were dominant in electronically and rovibrationally excited states, which may have yet unrecognized roles in the upper planetary atmosphere.

Rationalizing the Fluorescence Behavior of Core-Substituted Naphthalene Diimides

We study the internal conversion (IC) and intersystem crossing (ISC) pathways of low-lying excited electronic states of three core-substituted naphthalene diimides (bNDI, yNDI, and gNDI) using wavepacket simulations within the linear vibronic coupling method. Our wavepacket simulations reproduce the experimental electronic absorption spectra very well. All molecules decay rapidly to S₂ upon populating a higher dipole-allowed singlet excited-state. The S₂ → S₁ IC dynamics and singlet-triplet energy gap, spin-orbit coupling strength trends suggest a favorable S₂ → T₄ ISC in gNDI. The efficient ultrafast T₄ formation and its decay to lower triplet states make gNDI nonfluorescent. Such triplet formation pathways are not operative in both bNDI and yNDI; hence, these molecules emit fluorescence from S₁ after a slower S₂ → S₁ IC.

Intersystem crossing and internal conversion dynamics with GAIMS-TeraChem: Excited state relaxation in 2-cyclopentenone

Excited states relaxation in complex molecules often involves two types of nonradiative transitions, internal conversion (IC) and intersystem crossing (ISC). In the situations when the timescales of IC and ISC are comparable, an interplay between these two types of transitions can lead to complex nonadiabatic dynamics on multiple electronic states of different characters and spin multiplicities. We demonstrate that the generalized ab initio multiple spawning (GAIMS) method interfaced with the fast graphics processing unit-based TeraChem electronic structure code can be used to model such nonadiabatic dynamics involving both the IC and ISC transitions in molecules of moderate size. We carried out 1500 fs GAIMS simulations leading to the creation of up to 2500 trajectory basis functions to study the excited states relaxation in 2-cyclopentenone. After a vertical excitation from the ground state to the bright S₂ state, the molecule quickly relaxes to the S₁ state via conical intersection. The following relaxation proceeds along two competing pathways: one involves IC to the ground state, and the other is dominated by ISC to the low-lying triplet states. The time constants describing the population transfer between the six lowest singlet and triplet states predicted by the GAIMS dynamics are in good agreement with the characteristic times of IC and ISC obtained from the analysis of the time-resolved photoelectron spectrum.

Modeling Spin-Crossover Dynamics

In this article, we review nonadiabatic molecular dynamics (NAMD) methods for modeling spin-crossover transitions. First, we discuss different representations of electronic states employed in the grid-based and direct NAMD simulations. The nature of interstate couplings in different representations is highlighted, with the main focus on nonadiabatic and spin-orbit couplings. Second, we describe three NAMD methods that have been used to simulate spin-crossover dynamics, including trajectory surface hopping, ab initio multiple spawning, and multiconfiguration time-dependent Hartree. Some aspects of employing different electronic structure methods to obtain information about potential energy surfaces and interstate couplings for NAMD simulations are also discussed. Third, representative applications of NAMD to spin crossovers in molecular systems of different sizes and complexities are highlighted. Finally, we pose several fundamental questions related to spin-dependent processes. These questions should be possible to address with future methodological developments in NAMD.

Validating fewest-switches surface hopping in the presence of laser fields

The capability of fewest-switches surface hopping (FSSH) to describe non-adiabatic dynamics under explicit excitation with external fields is evaluated. Different FSSH parameters are benchmarked against multi-configurational time dependent Hartree (MCTDH) reference calculations using SO₂ and 2-thiocytosine as model, yet realistic, molecular systems. Qualitatively, FSSH is able to reproduce the trends in the MCTDH dynamics with (also without) an explicit external field; however, no set of FSSH parameters is ideal. The adequate treatment of the overcoherence in FSSH is revealed as the driving factor to improve the description of the excitation process with respect to the MCTDH reference. Here, two corrections were tested: the augmented-FSSH (AFSSH) correction and the energy-based decoherence correction. A dependence on the employed basis is detected in AFSSH, performing better when spin-orbit and external laser field couplings are treated as off-diagonal elements instead of projecting them onto the diagonal of the Hamilton operator. In the presence of an electric field, the excited state dynamics was found to depend strongly on the vector used to rescale the kinetic energy along after a transition between surfaces. For SO₂, recurrence of the excited wave packet throughout the duration of the applied laser pulse is observed for laser pulses (>100 fs), resulting in additional interferences missed by FSSH and only visible in variational multi-configurational Gaussian when utilizing a large number of Gaussian basis functions. This feature vanishes when going toward larger molecules, such as 2-thiocytosine, where this effect is barely visible in a laser pulse 200 fs long.

PySurf: A Framework for Database Accelerated Direct Dynamics

The greatest restriction to the theoretical study of the dynamics of photoinduced processes is computationally expensive electronic structure calculations. Machine learning algorithms have the potential to reduce the number of these computations significantly. Here, PySurf is introduced as an innovative code framework, which is specifically designed for rapid prototyping and development tasks for data science applications in computational chemistry. It comes with powerful Plugin and Workflow engines, which allows intuitive customization for individual tasks. Data is automatically stored through the database framework, which enables additional interpolation of properties in previously evaluated regions of the conformational space. To illustrate the potential of the framework, a code for nonadiabatic surface hopping simulations based on the Landau-Zener algorithm is presented here. Deriving gradients from the interpolated potential energy surfaces allows for full-dimensional nonadiabatic surface hopping simulations using only adiabatic energies (energy only). Simulations of a pyrazine model and ab initio-based calculations of the SO2 molecule show that energy-only calculations with PySurf are able to correctly predict the nonadiabatic dynamics of these prototype systems. The results reveal the degree of sophistication, which can be achieved by the database accelerated energy-only surface hopping simulations being competitive to commonly used semiclassical approaches.

Spin-inversion mechanisms in O2 binding to a model heme compound: A perspective from nonadiabatic wave packet calculations

Spin-inversion dynamics in O₂ binding to a model heme complex, which consisted of Fe(II)-porphyrin and imidazole, were studied using nonadiabatic wave packet dynamics calculations. We considered three active nuclear degrees of freedom in the dynamics, including the motions along the Fe-O distance, Fe-O-O angle, and Fe out-of-plane distance. Spin-free potential energy surfaces for the singlet, triplet, quintet, and septet states were developed using density functional theory calculations, and spin-orbit coupling elements were obtained from CASSCF-level electronic structure calculations. The spin-inversion mainly occurred between the singlet state and one of the triplet states due to large spin-orbit couplings and the contributions of other states were extremely small. The present quantum dynamics calculations suggested that the narrow crossing region model plays a dominant role in the O₂ binding dynamics. In addition, the one-dimensional Landau-Zener model underestimated the nonadiabatic transition probability.

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.

Dissipation and dynamics in ultrafast intersystem crossings

The effects of dynamics and dissipation on ultrafast intersystem crossings are studied for a dissipative two-level system coupled to a local vibronic mode. A method of amplitude damping of the wave packet is presented that accounts better for the position of the wave packet and avoids spurious transitions between potential wells. It is demonstrated that Fermi's golden rule, the typical semiquantitative approach to extract population transfer rates from potential landscapes, only holds under limited conditions. Generally, the effects of dynamics and dissipation lead to deviations from the expected exponential population transfer, strong changes in transfer times and total population transfer, and significant recurrence or "spill back" of the wave packet.

Triplet state promoted reaction of SO2 with H2O by competition between proton coupled electron transfer (pcet) and hydrogen atom transfer (hat) processes

The excited triplet electronic state of SO₂ (a³B₁) reacts with water through a proton coupled electron transfer (pcet) mechanism rather than via a conventional hydrogen atom transfer (hat) process.

Highly efficient surface hopping dynamics using a linear vibronic coupling model

We report an implementation of the linear vibronic coupling (LVC) model within the surface hopping dynamics approach and present utilities for parameterizing this model in a blackbox fashion. This results in an extremely efficient method to obtain qualitative and even semi-quantitative information about the photodynamical behavior of a molecule, and provides a new route toward benchmarking the results of surface hopping computations. The merits and applicability of the method are demonstrated in a number of applications. First, the method is applied to the SO2 molecule showing that it is possible to compute its absorption spectrum beyond the Condon approximation, and that all the main features and timescales of previous on-the-fly dynamics simulations of intersystem crossing are reproduced while reducing the computational effort by three orders of magnitude. The dynamics results are benchmarked against exact wavepacket propagations on the same LVC potentials and against a variation of the electronic structure level. Four additional test cases are presented to exemplify the broader applicability of the model. The photodynamics of the isomeric adenine and 2-aminopurine molecules are studied and it is shown that the LVC model correctly predicts ultrafast decay in the former and an extended excited-state lifetime in the latter. Futhermore, the method correctly predicts ultrafast intersystem crossing in the modified nucleobase 2-thiocytosine and its absence in 5-azacytosine while it fails to describe the ultrafast internal conversion to the ground state in the latter.

Nonadiabatic dynamics: The SHARC approach

We review the Surface Hopping including ARbitrary Couplings (SHARC) approach for excited-state nonadiabatic dynamics simulations. As a generalization of the popular surface hopping method, SHARC allows simulating the full-dimensional dynamics of molecules including any type of coupling terms beyond nonadiabatic couplings. Examples of these arbitrary couplings include spin-orbit couplings or dipole moment-laser field couplings, such that SHARC can describe ultrafast internal conversion, intersystem crossing, and radiative processes. The key step of the SHARC approach consists of a diagonalization of the Hamiltonian including these couplings, such that the nuclear dynamics is carried out on potential energy surfaces including the effects of the couplings-this is critical in any applications considering, for example, transition metal complexes or strong laser fields. We also give an overview over the new SHARC2.0 dynamics software package, released under the GNU General Public License, which implements the SHARC approach and several analysis tools. The review closes with a brief survey of applications where SHARC was employed to study the nonadiabatic dynamics of a wide range of molecular systems. This article is categorized under: Theoretical and Physical Chemistry > Reaction Dynamics and KineticsSoftware > Simulation MethodsSoftware > Quantum Chemistry.

Time-resolved photoelectron imaging with a femtosecond vacuum-ultraviolet light source: Dynamics in the A∼/B∼- and F∼-bands of SO2

Time-resolved photoelectron imaging is demonstrated using the third harmonic of a 400-nm femtosecond laser pulse as the ionization source. The resulting 133-nm pulses are combined with 266-nm pulses to study the excited-state dynamics in the A∼/B∼- and F∼-band regions of SO₂. The photoelectron signal from the molecules excited to the A∼/B∼-band does not decay for at least several picoseconds, reflecting the population of bound states. The temporal variation of the photoelectron angular distribution (PAD) reflects the creation of a rotational wave packet in the excited state. In contrast, the photoelectron signal from molecules excited to the F∼-band decays with a time constant of 80 fs. This time constant is attributed to the motion of the excited-state wave packet out of the ionization window. The observed time-dependent PADs are consistent with the F∼ band corresponding to a Rydberg state of dominant s character. These results establish low-order harmonic generation as a promising tool for time-resolved photoelectron imaging of the excited-state dynamics of molecules, simultaneously giving access to low-lying electronic states, as well as Rydberg states, and avoiding the ionization of unexcited molecules.

Excited state dynamics in SO2. I. Bound state relaxation studied by time-resolved photoelectron-photoion coincidence spectroscopy

The excited state dynamics of isolated sulfur dioxide molecules have been investigated using the time-resolved photoelectron spectroscopy and time-resolved photoelectron-photoion coincidence techniques. Excited state wavepackets were prepared in the spectroscopically complex, electronically mixed (B̃)(1)B1/(Ã)(1)A2, Clements manifold following broadband excitation at a range of photon energies between 4.03 eV and 4.28 eV (308 nm and 290 nm, respectively). The resulting wavepacket dynamics were monitored using a multiphoton ionisation probe. The extensive literature associated with the Clements bands has been summarised and a detailed time domain description of the ultrafast relaxation pathways occurring from the optically bright (B̃)(1)B1 diabatic state is presented. Signatures of the oscillatory motion on the (B̃)(1)B1/(Ã)(1)A2 lower adiabatic surface responsible for the Clements band structure were observed. The recorded spectra also indicate that a component of the excited state wavepacket undergoes intersystem crossing from the Clements manifold to the underlying triplet states on a sub-picosecond time scale. Photoelectron signal growth time constants have been predominantly associated with intersystem crossing to the (c̃)(3)B2 state and were measured to vary between 750 and 150 fs over the implemented pump photon energy range. Additionally, pump beam intensity studies were performed. These experiments highlighted parallel relaxation processes that occurred at the one- and two-pump-photon levels of excitation on similar time scales, obscuring the Clements band dynamics when high pump beam intensities were implemented. Hence, the Clements band dynamics may be difficult to disentangle from higher order processes when ultrashort laser pulses and less-differential probe techniques are implemented.

Ultrafast intersystem crossing dynamics in uracil unravelled by ab initio molecular dynamics

Ab initio molecular dynamics simulations have been performed in order to investigate the relaxation dynamics of uracil after UV excitation in gas phase. Intersystem crossing (ISC) has been included for the first time into time-dependent simulations of uracil, allowing the system to relax in the singlet as well as in the triplet states. The results show a qualitatively different picture than similar simulations that include singlet states only. The inclusion of ISC effectively quenches the relaxation to the singlet ground state and instead privileges transitions from the low-lying nπ* state (S1) to a ππ* triplet state (T2) followed by rapid internal conversion to the lowest triplet state.