Search

It is well-known that photolysis of pyrimidine nucleobases, such as uracil, in an aqueous environment results in the formation of hydrate as one of the main products. Although several hypotheses regarding photohydration have been proposed in the past, e.g., the zwitterionic and “hot” ground-state mechanisms, its detailed mechanism remains elusive. Here, theoretical nonadiabatic simulations of the uracil photodynamics reveal the formation of a highly energetic but kinetically stable intermediate that features a half-chair puckered pyrimidine ring and a strongly twisted intracyclic double bond. The existence and the kinetic stability of the intermediate are confirmed by a variety of computational chemistry methods. According to the simulations, the unusual intermediate is mainly formed almost immediately (∼50–200 fs) upon photoabsorption and survives long enough to engage in a hydration reaction with a neighboring water. A plausible mechanism of uracil photohydration is proposed on the basis of the modeling of nucleophilic insertion of water into the twisted double bond of the intermediate.

Singlet fission is a highly desired process in photovoltaic devices as it can significantly enhance photoelectric conversion efficiency. Exploitation of this process in photovoltaics is hindered by the lack of appropriate chromophores. We used mixed-reference spin-flip time-dependent density functional theory (TDDFT) to investigate five dicarbonyl anthracenes and phenanthrenes, with the purpose to design singlet fission chromophores. Two molecules were found to be promising candidates. For all the dicarbonyl molecules, the oxygen lone pair orbitals were found to be involved in the excited states that are relevant to singlet fission.

Multiple ERI (Electron Repulsion Integral) tensor contractions (METC) with several matrices are ubiquitous in quantum chemistry. In response theories, the contraction operation, rather than ERI computations, can be the major bottleneck, as its computational demands are proportional to the multiplicatively combined contributions of the number of excited states and the kernel pre-factors. This paper presents several high-performance strategies for METC. Optimal approaches involve either the data layout reformations of interim density and Fock matrices, the introduction of intermediate ERI quartet buffer, and loop-reordering optimization for a higher cache hit rate. The combined strategies remarkably improve the performance of the MRSF (mixed reference spin flip)-TDDFT (time-dependent density functional theory) by nearly 300%. The results of this study are not limited to the MRSF-TDDFT method and can be applied to other METC scenarios.

The molecular orbital (MO) concept is a useful tool, which relates the molecular ground-state energy with the energies (and occupations) of the individual orbitals. However, analysis of the excited states from linear response computations is performed in terms of the initial state MOs or some other forms of orbitals, e.g., natural or natural transition orbitals. Because these orbitals lack the respective energies, they do not allow developing a consistent orbital picture of the excited states. Herein, we argue that Dyson’s orbitals enable description of the response states compatible with the concepts of molecular orbital theory. The Dyson orbitals and their energies obtained by mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) for the response ground state are remarkably similar to the canonical MOs obtained by the usual DFT calculation. For excited states, the Dyson orbitals provide a chemically sensible picture of the electronic transitions, thus bridging the chasm between orbital theory and response computations.

The photochemistry of nucleobases, important for their role as building blocks of DNA, is largely affected by the electrostatic environment in which they are soaked. For example, despite the numerous studies of thymine in solution and DNA, there is still a debate on the photochemical deactivation pathways after UV absorption. Many theoretical models are oversimplified due to the lack of computationally accurate and efficient electronic structure methodologies that capture excited state electron correlation effects when nucleobases are embedded in large electrostatic media. Here, we combine mixed-reference spin-flip time-dependent density-functional theory (MRSF-TDDFT) with electrostatic embedding QM/MM using electrostatic potential fitting (ESPF) atomic charges, as a strategy to accurately and efficiently describe the electronic structure of chromophores polarized by an electrostatic medium. In particular, we develop analytic expressions for the energy and gradient of MRSF/MM based on the ESPF coupling using atom-centered grids and total charge conservation. We apply this methodology to the study of solvation effects on thymine photochemistry in water and thymine dimers in DNA. In the former, the combination of trajectory surface hopping (TSH) non-adiabatic molecular dynamics (NAMD) with MRSF/MM remarkably revealed the accelerated deactivation decay pathways, which is consistent with the experimental decay time of ~ 400 fs. The enhanced hopping rate can be explained by the preferential stabilization of corresponding conical interactions due to their increased dipole moments. Structurally, it is a consequence of characteristic methyl puckered geometries near the conical intersection region. For the thymine dimer in B-DNA, we found new photochemical pathways through conical intersections that could explain the formation of cyclobutadiene dimers and 6-4 photoproducts.

Relativistic MRSF-TDDFT is developed considering the spin-orbit coupling (SOC) within the mean-field approximation. The resulting SOC-MRSF faithfully reproduces the experiments with very high accuracy, which is also consistent with the values by four-component (4c) relativistic CASSCF and 4c-CASPT2 in the spin-orbit-energy splitting calculations of the C, Si and Ge atoms. Even for the fifth-row element Sn, the SOC-MRSF yielded accurate splittings (~3% error). In the SOC calculations of the molecular 4-thiothymine with a third-row element, SOC-MRSF values are in excellent agreement with those of SO-GMC-QDPT2 level, regardless of geometries and exchange-correlation functionals. The same SOC-MRSF predicted the anticipated chance of S1 (n pi*) -> T1 (pi pi*) intersystem crossing, even in thymine with only second-row elements. With its accuracy and practicality, thus, SOC-MRSF is a promising electronic structure protocol in challenging situations such as nonadiabatic molecular dynamics (NAMD) incorporating both internal conversions and intersystem crossings in large systems.

Internal conversion (IC) between the two lowest singlet excited states, 11Bu+ and 21Ag-, of s-trans-butadiene and s-trans-hexatriene is investigated using a series of single- and multi- reference wave function and density functional theory (DFT) methodologies. Three independent types of the equation-of-motion coupled-cluster (EOMCC) theory capable of providing an accurate and balanced description of one- as well as two-electron transitions, abbreviated as δ-CR-EOMCC(2,3), DIP-EOMCC(4h2p){No}, and DEA-EOMCC(4p2h){Nu} or DEA-EOMCC(3p1h,4p2h){Nu}, consistently predict that the 11Bu+/21Ag- crossing in both molecules occurs along the bond length alternation coordinate. However, the analogous 11Bu+ and 21Ag- potentials obtained with some multireference approaches, such as CASSCF and MRCIS(D), as well as with the linear-response formulation of time-dependent DFT (TDDFT), do not cross. Hence, caution needs to be exercised when studying the low-lying singlet excited states of polyenes with conventional multiconfigurational methods and TDDFT. The multistate many-body perturbation theory methods, such as XMCQDPT2, do correctly reproduce the curve crossing. Among the simplest and least expensive computational methodologies, the DFT approaches that incorporate the contributions of doubly excited configurations, abbreviated as MRSF (mixed reference spin-flip) TDDFT and SSR(4,4), accurately reproduce our best EOMCC results. This is highly promising for nonadiabatic molecular dynamics simulations in larger systems.

International audience; It is demonstrated that the challenging core-hole particle (CHP) orbital relaxation for core electron spectra can be readily achieved by the mixed-reference spin-flip (MRSF)–time-dependent density functional theory (TDDFT). With the additional scalar relativistic effects on K-edge excitation energies of 24 second- and 17 third-row molecules, the particular ΔCHP–MRSF(R) exhibited near perfect predictions with RMSE ∼0.5 eV, featuring a median value of 0.3 and an interquartile range of 0.4. Overall, the CHP effect is 2–4 times stronger than relativistic ones, contributing more than 20 eV in the cases of sulfur and chlorine third-row atoms. Such high precision allows to explain the splitting and spectral shapes of O, N, and C atom K-edges in the ground state of thymine with atom as well as orbital specific accuracy. The same protocol with a double hole particle relaxation also produced remarkably accurate K-edge spectra of core to valence hole excitation energies from the first (nO8π*) and second (ππ*) excited states of thymine, confirming the assignment of 1s → n excitation for the experimentally observed 526.4 eV peak. Regarding both accuracy and practicality, therefore, MRSF–TDDFT provides a promising protocol for core electron spectra of both ground and excited electronic states alike.

Oxidative-denitrogenation (ODN) of indole (IND) and methyl-substituted INDs (methyl-INDs), representative neutral nitrogen-containing compounds (NCCs), was carried out with TiO2@C and H2O2 as heterogeneous catalyst and oxidant, respectively, under ultrasound irradiation. The oxidation of INDs progressed through radical formation, as evidenced by electron spin resonance and radical scavenger experiments. The oxidized position of INDs in the ODN process was checked via characterization of the obtained products. It was observed that the oxidation finally occurred on the carbon rather than on the nitrogen atom of INDs, unlike the oxidation of basic NCCs (e.g., oxidation on the nitrogen atom, as respective N-oxides were formed) and sulfur-containing compounds. To understand the relative reactivity and oxidation position, electron density (ED) on the nitrogen atom of the studied INDs and relative stability of representative intermediates/products were calculated. It could be confirmed that ED on the nitrogen atom of the INDs is very important in the oxidation of INDs since the ODN reactivity of INDs was enhanced with increasing ED on the nitrogen atom of the investigated INDs. Moreover, theoretical analyses of the relative stability of substrate and intermediates/products (especially for IND) can explain the route for the observed final products in ODN. In other words, oxygen on the nitrogen atom, obtained via the first step of oxidation (electrophilic addition of an active oxygen atom on nitrogen), moves to the nearby carbon atom, because of the relative stability of the intermediates and products.

Computation of the ionization energies and of the respective Dyson orbitals based on the use of the extended Koopmans theorem (EKT) is implemented in connection with an ensemble density functional theory (eDFT) method, the state-interaction state-averaged spin-restricted ensemble-referenced Kohn–Sham (SI-SA-REKS or SSR) method. The new methodology enables fast computation of the ionization energies and evaluation of the respective Dyson orbitals, the square norms of which are related with the ionization probabilities, in the ground and excited electronic states of molecules. As the application of EKT recycles the intermediate quantities from the SSR analytical energy gradient, evaluation of the ionization energies and probabilities can be carried out on-the-fly during the nonadiabatic molecular dynamics simulations. This opens up a perspective for fast theoretical simulation of the time-resolved photoelectron spectroscopy observations. In the present work, the new methodology is tested in the computation of the ionization energies and Dyson orbitals of several molecules in the ground and excited electronic states, including strongly correlated species, such as the ozone molecule, dissociating chemical bonds, and conical intersections.

Internal conversion (IC) between the two lowest singlet excited states, 1

Symmetry-breaking charge transfer (SBCT) is an important process at the early stages of the photoinduced processes in multichromophore systems such as the photosynthetic apparatus. We investigated the photoinduced SBCT dynamics of 9,9'-bianthracene (BA), a representative molecule showing SBCT, by time-resolved fluorescence (TF) with the highest time-resolution and excited-state quantum mechanics/effective fragment potential molecular dynamics (MD) simulation. TF experiments show that the SBCT kinetics matches quantitatively with the solvation function excluding the initial ultrafast component that is assigned to the inertial motion of the solvent. Therefore, it is established that the SBCT of BA is coupled solely with the rotational diffusion of solvent molecules excluding the inertial motion of solvents. MD simulations show that random rotational fluctuation of solvents mostly in the first solvation shell generates a transient electric field as high as 1.0 × 109 V m-1, which provides an asymmetric environment required for the generation of a CT state in this symmetric dimer. Once the CT state is formed, the dipole moment in the solute causes further rotation of solvent molecules leading to an augmented electric field, which in turn further stabilizes the CT state prohibiting the reverse reaction.

The one-dimensional projection (ODP) approach is extended to two-dimensional umbrella sampling (TDUS) and is applied to three different complex systems in combination with a reactive force field (ReaxFF). TDUS is capable of showing detailed features of the free-energy surface (FES) of the double-proton transfer of the acetic acid dimer. It also revealed the direct relationship between the types of hydrogen bonding and binding strengths in the case of adrenaline molecular recognition by SIVSF (Serine, Isoleucine, Valine, Cysteine, and Phenylalanine). The study of polymer aggregation using TDUS shows that aggregation is preferred with a less-polar solvent, which is also consistent with the experimental observation of a tape-casting process. Therefore, TDUS can be generally useful in FES explorations from simple chemical reactions to complex processes of molecular recognition and polymer aggregation.

Time-resolved photoelectron spectra during the photochemical ring-opening reaction of 1,3-cyclohexadiene (CHD) are modeled by an ensemble density functional theory (eDFT) method. The computational methodology employed in this work is capable of correctly describing the multi-reference effects arising in the ground and excited electronic states of molecules, which is important for the correct description of the ring-opening reaction of CHD. The geometries of molecular species along the non-adiabatic molecular dynamics (NAMD) trajectories reported in a previous study of the CHD photochemical ring-opening were used in this work to calculate the ionization energies and the respective Dyson orbitals for all possible ionization channels. The obtained theoretical time-resolved spectra display decay characteristics in a reasonable agreement with the experimental observations; i.e., the decay (and rise) of the most mechanistically significant signals occurs on the timescale of 100–150 fs. This is very different from the excited state population decay characteristics (τ _(S1) = 234 ± 8 fs) obtained in the previous NAMD study. The difference between the population decay and the decay of the photoelectron signal intensity is traced back to the geometric transformation that the molecule undergoes during the photoreaction. This demonstrates the importance of including the geometric information in interpretation of the experimental observations.

With the help of quantum mechanical methods, the formation of H2SO4 by the oxidation of H2SO3 with H2O2 was studied theoretically. Both stepwise and concerted mechanisms were calculated. It was found that the direct oxidation of H2SO3 by H2O2 alone requires prohibitive activation energies of >38.6 kcal/mol. However, the addition of one water molecule exhibits a strong catalytic effect that dramatically reduces the overall reaction barrier to 6.2 (2.3 with PCM) kcal/mol. The deprotonated HSO3- species also reduces the overall reaction barrier to 5.6 (-5.8 with PCM) kcal/mol. Both of these proceed via concerted pathways. On the other hand, the stepwise mechanisms generally produce intermediates with a hydroperoxy group (-O-O-H), which is a result of a nucleophilic attack by the oxygens of H2O2. While studying the catalytic effect of water, a previously unknown hydroperoxy intermediate (HOO)S(OH)3, where sulfur is coordinated with three OH groups, was found. This work also reveals a rearrangement step of another hydroperoxy intermediate (HOO)SO2- to HSO4- that was found in earlier experimental studies. For all of the mechanisms calculated, the final H2SO4 is formed with a significant exothermicity of >60 kcal/mol. In general, even without sunlight, it was found that the formation of sulfuric acid by hydrogen peroxide can occur in a heterogeneous moisturized environment.

The mixed-reference spin-flip (MRSF) time-dependent density functional theory (TDDFT) method eliminates the notorious spin contamination of SF-TDDFT, thus enabling identification of states of proper spin-symmetry for automatic geometry optimization and molecular dynamics simulations. Here, we analyze and optimize the MRSF-TDDFT in the calculations of the vertical excitation energies (VEEs) and the singlet-triplet (ST) gaps. The dependence of the obtained VEEs and ST gaps on the intrinsic parameters of the MRSF-TDDFT method is investigated, and prescriptions for the proper use of the method are formulated. For VEEs, MRSF-TDDFT displays similar or better accuracy than SF-TDDFT (ca. 0.5 eV), while considerably outperforming the LR-TDDFT for the ST gaps. As a result, a new functional of STG1X (dubbed here), especially for ST gaps is suggested on the basis of splitting between the components of the atomic multiplets.

The mixed-reference spin-flip time-dependent density functional theory (MRSF-TD-DFT) method eliminates the erroneous spin contamination of the SF-TD-DFT methodology, while retaining the conceptual and practical simplicity of the latter. The availability of the analytic gradient of the energy of the MRSF-TD-DFT response states enables automatic geometry optimization of the targeted states. Here, we apply the new method to optimize the geometry of several S1/S0 conical intersections occurring in typical organic molecules. We demonstrate that MRSF-TD-DFT is capable of producing the correct double-cone topology of the intersections and describing the geometry of the lowest-energy conical intersections and their relative energies with accuracy matching that of the best multireference wavefunction ab initio methods. In this regard, MRSF-TD-DFT differs from many popular single-reference methods, such as, e.g., the linear response TD-DFT method, which fail to produce the correct topology of the intersections. As the new methodology completely eliminates the ambiguity with the identification of the response states as proper singlets or triplets, which is plaguing the SF-TD-DFT calculations, it can be used for automatic geometry optimization and molecular dynamic simulations not requiring constant human intervention.

Sudden polarization (SP) is one of the manifestations of electron transfer in the electronically excited states of molecules. Proposed initially to explain the unusual reactivity of photoexcited olefins, SP often occurs in the excited states of molecules possessing strongly correlated diradical ground state. Theoretical description of SP involves mixing between the singly excited and the doubly excited zwitterionic states, which makes it inaccessible with the use of the popular linear-response time-dependent density functional theory methods. In this work, an extended variant of the state-interaction state-averaged spin-restricted ensemble-referenced Kohn–Sham (SI-SA-REKS, or SSR) method is applied to study SP in a number of organic diradical systems. To this end, the analytical derivative formalism is derived and implemented for the SSR(3,2) method (see the main text for explanation of the acronym), which enables the automatic geometry optimization and obtains the relaxed density matrices as well as the electron binding energies and respective Dyson’s orbitals. Application of the new method to SP in the lowest singlet excited state of ethylene agrees with the results obtained previously with the use of multireference methods of wavefunction theory. A number of interesting manifestations of SP are observed, such as the charge transfer in photoexcited tetramethyleneethene (TME) diradical mediated by the vibrational motion and conductivity switching in the excited state of a donor–acceptor dyad placed in an external electric field.

We present a fast and accurate numerical algorithm for computing the first-order nonadiabatic coupling matrix element (NACME). The algorithm employs the truncated Leibniz formula (TLF) approximation within the finite-difference method, which makes it easily applicable in connection with any wave function-based methodology. In this work, we used the algorithm in connection with the recently developed mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT, MRSF for brevity). The accuracy is assessed for NACME between the singlet electronic states of a dissociating hydrogen molecule. It is demonstrated that an intermediate approximation, TLF(1), affords a negligible numeric error on the order of ∼10⁻¹⁰ a.u. while enabling a fast computation of NACME. As the MRSF method yields the correct description of the dissociation curves of H₂ for all the electronic states involved, the numeric TLF(1)/MRSF NACME values are in excellent agreement with the reference analytical values obtained by the full configuration interaction. For polyatomic molecules, the MRSF NAC vectors agree very closely with the MRCISD NAC vectors. Hence, the proposed protocol is a promising tool for the evaluation of NACMEs.

Non-radiative relaxation of the photoexcited thymine in the gas phase shows an unusually long excited-state lifetime, and, over the years, a number of models, i.e., S1-trapping, S2-trapping, and S1&S2-trapping, have been put forward to explain its mechanism. Here, we investigate this mechanism using non-adiabatic molecular dynamics (NAMD) simulations in connection with the recently developed mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) method. We show that the previously predicted S2-trapping model was due to an artifact caused by an insufficient account of the dynamic electron correlation. The current work supports the S1-trapping mechanism with two lifetimes, τ1 = 30 ± 1 fs and τ2 = 6.1 ± 0.035 ps, quantitatively consistent with the recent time-resolved experiments. Upon excitation to the S2 (ππ*) state, thymine undergoes an ultrafast (ca. 30 fs) S2→S1 internal conversion and resides around the minimum on the S1 (nOπ*) surface, slowly decaying to the ground state (ca. 6.1 ps). While the S2→S1 internal conversion is mediated by fast bond length alternation distortion, the subsequent S1→S0 occurs through several conical intersections, involving a slow puckering motion.

The non-adiabatic dynamics of furan excited in the ππ* state (S2 in the Franck–Condon geometry) was studied using non-adiabatic molecular dynamics simulations in connection with an ensemble density functional method. The time-resolved photoelectron spectra were theoretically simulated in a wide range of electron binding energies that covered the valence as well as the core electrons. The dynamics of the decay (rise) of the photoelectron signal were compared with the excited-state population dynamics. It was observed that the photoelectron signal decay parameters at certain electron binding energies displayed a good correlation with the events occurring during the excited-state dynamics. Thus, the time profile of the photoelectron intensity of the K-shell electrons of oxygen (decay constant of 34 ± 3 fs) showed a reasonable correlation with the time of passage through conical intersections with the ground state (47 ± 2 fs). The ground-state recovery constant of the photoelectron signal (121 ± 30 fs) was in good agreement with the theoretically obtained excited-state lifetime (93 ± 9 fs), as well as with the experimentally estimated recovery time constant (ca. 110 fs). Hence, it is proposed to complement the traditional TRPES observations with the trXPS (or trNEXAFS) measurements to obtain more reliable estimates of the most mechanistically important events during the excited-state dynamics.

A new adaptive algorithm for penalty function optimization for minimum-energy three-states conical intersections (ME3CI) is suggested. The new algorithm differs from the original penalty function algorithm by (a) removing the redundancy in the target function, (b) using an adaptive increment for the penalty function weighting factor, and (c) using tighter convergence criteria for the energy gap. The latter was introduced to guarantee convergence to a true conical intersection rather than to a narrowly avoided crossing geometry. The new algorithm was tested in the optimization of the ME3CI geometries in butadiene and malonaldehyde, where all of the previously found true ME3CI geometries were recovered. The previously found butadiene's CI

In different branches of time-dependent density functional theory (TDDFT), the static and dynamic electron correlation enters in different ways. The standard spin-conserving linear response (LR-TDDFT) methodology includes explicitly the contributions of the singly-excited configurations; however, it relies on an implicit account of the electron correlation through an (approximate) exchange-correlation (XC) functional. In the mixed-reference spin-flip TDDFT (MRSF-TDDFT), a number of doubly-excited (DE) configurations are explicitly included in the description of their response states. Here, the importance of the explicit account of DE is investigated for the lowest four excited singlet states of all-trans-polyenes up to C₂₄H₂₆. For the optically bright 1B_u⁺ state, the DE contribution in MRSF-TDDFT approaches 10% with the increasing system size. For the optically dark 2Ag⁻ state, the DE contribution increases from ca. 13% (C₄H₆) to nearly 30% (C₂₄H₂₆). An even more considerable DE contribution (∼50%) is observed in the higher 1B_u⁻ states. As LR-TDDFT misses these contributions entirely, its ability to accurately describe the excited states is limited by the XC functional. The hybrid XC functionals with a small fraction of the exact exchange, e.g., B3LYP, may mimic certain effects of DE through the self-interaction error (SIE). However, the description of the 1B_u⁺ state by LR-TDDFT remains poor. On the other hand, MRSF-TDDFT can flexibly take an implicit (through the XC functional) and an explicit (through DE) account of the electron correlation, which enables a more balanced description of various types of the excited states regardless of their character, thus reducing the chances of failure.