Highly Efficient Algorithms for CIS Type Excited State Wave Function Overlaps

The Role of Aqueous Solvation on the Intersystem Crossing of Nitrophenols

The photochemistry of nitrophenols is a source of smog as nitrous acid is formed from their photolysis. Nevertheless, computational studies of the photochemistry of these widespread toxic molecules are scarce. In this work, the initial photodeactivation of ortho-nitrophenol and para-nitrophenol is modeled, both in gas phase and in aqueous solution to simulate atmospheric and aerosol environments. A large number of excited states, six for ortho-nitrophenol and 11 for para-nitrophenol, have been included and were all populated during the decay. Moreover, periodic time-dependent density functional theory (TDDFT) is used for both the explicitly included solvent and the solute. A comparison to periodic QM/MM (TDDFT/MM), with electrostatic embedding, is made, showing notable differences between the decays of solvated nitrophenols simulated with QM/MM and full (TD)DFT. A reduced intersystem crossing in aqueous solution could be observed thanks to the surface hopping approach using explicit, periodic TDDFT solvation including spin-orbit couplings.

Excited-state antiaromaticity relief drives facile photoprotonation of carbons in aminobiphenyls

A combined computational and experimental study reveals that ortho-, meta- and para-aminobiphenyl isomers undergo distinctly different photochemical reactions involving proton transfer. Deuterium exchange experiments show that the ortho-isomer undergoes a facile photoprotonation at a carbon atom via excited-state intramolecular proton transfer (ESIPT). The meta-isomer undergoes water-assisted excited-state proton transfer (ESPT) and a photoredox reaction via proton-coupled electron transfer (PCET). The para-isomer undergoes a water-assisted ESPT reaction. All three reactions take place in the singlet excited-state, except for the photoredox process of the meta-isomer, which involves a triplet excited-state. Computations illustrate the important role of excited-state antiaromaticity relief in these photoreactions.

On the propensity of formation of cyclobutane dimers in face-to-face and face-to-back uracil stacks in solution

UV irradiation of RNA leads to the formation of intra- and inter-strand crosslinks of cyclobutane type. Despite the importance of this reaction, relatively little is known about how the mutual orientation of the two bases affects the outcome of the reaction. Here we report a comparative nonadiabatic molecular dynamics study of face-to-back (F2B) and face-to-face (F2F) stacked uracil-water clusters. The computations were performed using the second-order algebraic-diagrammatic-construction (ADC(2)) method. We found that F2B stacked uracil-water clusters either relax non-reactively to the ground state by an ethylenic twist around the CC bond or remain in the lowest nπ* state in which the two bases gradually move away from each other. This finding is consistent with the low propensity for the formation of intra-strand cyclobutane dimers between adjacent RNA bases. On the contrary, in F2F stacked uracil-water clusters, in addition to non-reactive deactivation, we found a pro-reactive deactivation pathway, which may lead to the formation of cyclobutane uracil dimers in the electronic ground state. On a qualitative level, the observed photodynamics of F2F stacked uracil-water clusters explains the greater propensity of RNA to form inter-strand cyclobutane-type crosslinks.

Transition orbital projection approach for excited state tracking

Quantitively comparing the features between different electronic excited states (ESs) is a crucial task in both potential energy surface (PES) studies and excited-state fragmentation approaches. However, it is still a challenging problem in regard to the comparison of complex and highly degenerate systems. Herein, we present a transition orbital projection (TOP) method to calculate the similarity of different ESs based on the configuration vectors of two types of transition densities. It fully considers four significant problems, including phase, hole-particle bijectivity, orbital permutation, and sign of configuration coefficients. TOP state-tracking-based excited-state optimization shows high robustness in several high-symmetric systems, which are difficult to describe with traditional state-tracking approaches. The TOP state-tracking method is expected to be widely applied to the PES of photochemical reactions, ES molecular dynamics to track the diabatic states, and fragmentation approaches for local excitation of large systems.

Spin-Orbit Couplings for Nonadiabatic Molecular Dynamics at the ΔSCF Level

A procedure for the calculation of spin-orbit coupling (SOC) at the delta self-consistent field (ΔSCF) level of theory is presented. Singlet and triplet excited electronic states obtained with the ΔSCF method are expanded into a linear combination of singly excited Slater determinants composed of ground electronic state Kohn-Sham orbitals. This alleviates the nonorthogonality between excited and ground electronic states and introduces a framework, similar to the auxiliary wave function at the time-dependent density functional theory (TD-DFT) level, for the calculation of observables. The ΔSCF observables of the formaldehyde system were compared to reference TD-DFT values. Our procedure gives all components (energies, gradients, nonadiabatic couplings, and SOC terms) at the ΔSCF level of theory for conducting efficient, full-atomistic nonadiabatic molecular dynamics with intersystem crossing, particularly in condensed phase systems.

Methods to Calculate Electronic Excited-State Dynamics for Molecules on Large Metal Clusters with Many States: Ensuring Fast Overlap Calculations and a Robust Choice of Phase

We present an efficient set of methods for propagating excited-state dynamics involving a large number of configuration interaction singles (CIS) or Tamm-Dancoff approximation (TDA) single-reference excited states. Specifically, (i) following Head-Gordon et al., we implement an exact evaluation of the overlap of singly-excited CIS/TDA electronic states at different nuclear geometries using a biorthogonal basis and (ii) we employ a unified protocol for choosing the correct phase for each adiabat at each geometry. For many-electron systems, the combination of these techniques significantly reduces the computational cost of integrating the electronic Schrodinger equation and imposes minimal overhead on top of the underlying electronic structure calculation. As a demonstration, we calculate the electronic excited-state dynamics for a hydrogen molecule scattering off a silver metal cluster, focusing on high-lying excited states, where many electrons can be excited collectively and crossings are plentiful. Interestingly, we find that the high-lying, plasmon-like collective excitation spectrum changes with nuclear dynamics, highlighting the need to simulate non-adiabatic nuclear dynamics and plasmonic excitations simultaneously. In the future, the combination of methods presented here should help theorists build a mechanistic understanding of plasmon-assisted charge transfer and excitation energy relaxation processes near a nanoparticle or metal surface.

The photodissociation of solvated cyclopropanone and its hydrate explored via non-adiabatic molecular dynamics using ΔSCF

The decay of cyclopropanone is a typical example of a photodecomposition process. Ethylene and carbon monoxide are formed following the excitation to the first singlet excited state through a symmetrical or asymmetrical pathway. The results obtained with non-adiabatic molecular dynamics (NAMD) using the delta self-consistent field (ΔSCF) method correspond well to previous experimental and multireference theoretical studies carried out in the gas phase. Moreover, this efficient methodology allows NAMD simulations of cyclopropanone in aqueous solution to be performed, which reveal analogue deactivation mechanisms, but a shorter lifetime and reduced photodissociation as compared to the gas-phase. The excited state dynamics of cyclopropanone hydrate, an enzyme inhibitor, in an aqueous environment are reported as well. Cyclopropanone hydrate strongly interacts with the surrounding solvent via the formation of hydrogen bonds. Excitation to the first singlet excited state shows an asymmetric pathway with cyclopropanone hydrate and propionic acid as the main photoproducts.

The lifetime and photodissociation of cyclopropanone are reduced in aqueous solution, while the excitation of solvated cyclopropanone hydrate yields a range of photoproducts.

A Computational Approach to Nontraditional Intrinsic Luminescence: Vibrationally Resolved Absorption and Fluorescence Spectra of DABCO

Recently, so-called "nontraditional intrinsic luminescence" has been reported in several macromolecular systems. Although DABCO (1,4-diazabicyclo[2.2.2]octane) is the first system in which the effect was observed, a thorough analysis of the optical properties of the molecule, which would reveal the origin of this mysterious effect, is still pending. We perform an advanced post-Hartree-Fock treatment of the low-lying electronic states of this molecule, which need to be described with care because of their pronounced Rydberg character. We take a deeper look into the low-lying electronic transitions of DABCO targeting the explanation of the complex vibronic structures of its absorption and fluorescence spectra. Two electronic states, the ¹E'(n₊3p_xy) and ¹A₂^″(n₊3p_z) states, contribute to the absorption spectrum in the 39000-46000 cm^-1 spectral range. We also reveal the spectroscopic signature of the ¹A₂^″(n₊3p_z) state. The analyses of the contributions of individual vibrational normal modes allowed the identification of those giving rise to the complex vibronic structures of the spectra. Fluorescence emission arises from the vibronic coupling of the one-photon forbidden transition between the ¹A₁^'(n₊3s) state and the electronic ground state. The spectrum, which can be interpreted in terms of populating a few vibrational normal modes, is shifted toward visible wavelengths mostly due to the forced interaction of the lone pair electrons of the two nitrogen atoms. Our work on DABCO may help to rationalize the luminescence of more complex systems containing tertiary amine groups.

Internal conversion of singlet and triplet states employing numerical DFT/MRCI derivative couplings: Implementation, tests, and application to xanthone

We present an efficient implementation of nonadiabatic coupling matrix elements (NACMEs) for density functional theory/multireference configuration interaction (DFT/MRCI) wave functions of singlet and triplet multiplicity and an extension of the Vibes program that allows us to determine rate constants for internal conversion (IC) in addition to intersystem crossing (ISC) nonradiative transitions. Following the suggestion of Plasser et al. [J. Chem. Theory Comput. 12, 1207 (2016)], the derivative couplings are computed as finite differences of wave function overlaps. Several measures have been taken to speed up the calculation of the NACMEs. Schur's determinant complement is employed to build up the determinant of the full matrix of spin-blocked orbital overlaps from precomputed spin factors with fixed orbital occupation. Test calculations on formaldehyde, pyrazine, and xanthone show that the mutual excitation level of the configurations at the reference and displaced geometries can be restricted to 1. In combination with a cutoff parameter of t_norm = 10^-8 for the DFT/MRCI wave function expansion, this approximation leads to substantial savings of cpu time without essential loss of precision. With regard to applications, the photoexcitation decay kinetics of xanthone in apolar media and in aqueous solution is in the focus of the present work. The results of our computational study substantiate the conjecture that S₁ T₂ reverse ISC outcompetes the T₂ ↝ T₁ IC in aqueous solution, thus explaining the occurrence of delayed fluorescence in addition to prompt fluorescence.

Simulation of UV absorption spectra and relaxation dynamics of uracil and uracil-water clusters

Despite many studies, the mechanisms of nonradiative relaxation of uracil in the gas phase and in aqueous solution are still not fully resolved. Here we combine theoretical UV absorption spectroscopy with nonadiabatic dynamics simulations to identify the photophysical mechanisms that can give rise to experimentally observed decay time constants. We first compute and theoretically assign the electronic spectra of uracil using the second-order algebraic-diagrammatic-construction (ADC(2)) method. The obtained electronic states, their energy differences and state-specific solvation effects are the prerequisites for understanding the photodynamics. We then use nonadiabatic trajectory-surface-hopping dynamics simulations to investigate the photoinduced dynamics of uracil and uracil-water clusters. In contrast to previous studies, we found that a single mechanism - the ethylenic twist around the C[double bond, length as m-dash]C bond - is responsible for the ultrafast component of the nonradiative decay, both in the gas phase and in solution. Very good agreement with the experimentally determined ultrashort decay time constants is obtained.

Structural basis for unique color tuning mechanism in heliorhodopsin

Microbial rhodopsins comprise an opsin protein with seven transmembrane helices and a retinal as the chromophore. An all-trans retinal is covalently bonded to a lysine residue through the retinal Schiff base (RSB) and stabilized by a negatively charged counterion. The distance between the RSB and counterion is closely related to the light energy absorption. However, in heliorhodopsin-48C12 (HeR-48C12), while E107 acts as the counterion, E107D mutation exhibits an identical absorption spectrum to the wild-type, suggesting that the distance does not affect its absorption spectra. Here we present the 2.6 Å resolution crystal structure of the Thermoplasmatales archaeon HeR E108D mutant, which also has an identical absorption spectrum to the wild-type. The structure revealed that D108 does not form a hydrogen bond with the RSB, and its counterion interaction becomes weaker. Alternatively, the serine cluster, S78, S112, and S238 form a distinct interaction network around the RSB. The absorption spectra of the E to D and S to A double mutants suggested that S112 influences the spectral shift by compensating for the weaker counterion interaction. Our structural and spectral studies have revealed the unique spectral shift mechanism of HeR and clarified the physicochemical properties of HeRs.

Fast Nonadiabatic Molecular Dynamics via Spin-Flip Time-Dependent Density-Functional Tight-Binding Approach: Application to Nonradiative Relaxation of Tetraphenylethylene with Locked Aromatic Rings

Nonadiabatic dynamics around conical intersections between ground and excited states are crucial to understand excited-state phenomena in complex chemical systems. With this background in mind, we present an approach combining fewest-switches trajectory surface hopping and spin-flip (SF) time-dependent (TD) density-functional tight binding (DFTB), which is a simplified version of SF-TD density functional theory (DFT) with semiempirical parametrizations, for computationally efficient nonadiabatic molecular dynamics simulations. The estimated computational time of the SF-TD-DFTB approach is several orders of magnitude lower than that of SF-TD-DFT. In addition, the proposed method reproduces the time scales and quantum yields in photoisomerization reactions of azobenzene at a level comparable with conventional ab initio approaches, demonstrating reasonable accuracy. Finally, we report a practical application of the developed technique to explore the nonradiative relaxation processes of tetraphenylethylene and its derivative with torsionally locked aromatic rings and discuss the effect of locking the rings on the excited-state lifetime.

UV absorption spectra of DNA bases in the 350-190 nm range: assignment and state specific analysis of solvation effects

The theoretical assignment of electronic spectra of polyatomic molecules is a challenging problem that requires the specification of the character of a large number of electronic states. We propose a procedure for automatically determining the character of electronic transitions and apply it to the study of UV spectra of DNA bases in the gas phase and in the aqueous environment. The procedure is based on the computation of electronic wave function overlaps and accounts for an extensive sampling of nuclear geometries. Novelties of this work are the theoretical assignment of the electronic spectra of DNA bases up to 190 nm and a state specific analysis of solvation effects. By accounting for different effects contributing to the total solvent shift we obtained a good agreement between the computed and experimental spectra. Effects of vibrational averaging, temperature and solvent-induced structural changes shift excitation energies to lower values. Solvent-solute electrostatic interactions are state specific and strongly destabilize nRyd states, and to lesser extent nπ* and πRyd states. Altogether, this results in the stabilization of ππ* states and destabilization of nπ*, πRyd and nRyd states in solution.

Toward the evaluation of intersystem crossing rates with variational relativistic methods

The change in electronic state from one spin multiplicity to another, known as intersystem crossing, occurs in molecules via the relativistic phenomenon of spin-orbit coupling. Current means of estimating intersystem crossing rates rely on the perturbative evaluation of spin-orbit coupling effects. This perturbative approach, valid in lighter atoms where spin-orbit coupling is weaker, is expected to break down for heavier elements where relativistic effects become dominant. Methods which incorporate spin-orbit effects variationally, such as the exact-two-component (X2C) method, will be necessary to treat this strong-coupling regime. We present a novel procedure which produces a diabatic basis of spin-pure electronic states coupled by spin-orbit terms, generated from fully variational relativistic calculations. This method is implemented within X2C using time-dependent density-functional theory and is compared to results from a perturbative relativistic study in the weak spin-orbit coupling regime. Additional calculations on a more strongly spin-orbit-coupled [UO₂Cl₄]^2- complex further illustrate the strengths of this method. This procedure will be valuable in the estimation of intersystem crossing rates within strongly spin-coupled species.