New General Tools for Constrained Geometry Optimizations

From skin sensitizers to wastewater: the unknown photo-deactivation process of low-lying excited states of isothiazolinones. A non-adiabatic dynamics investigation

Isothiazolinones represent a class of heterocyclic compounds widely used in various applications, including as biocides in cosmetics, detergents, and paints, as well as in industrial wastewater treatment. Indeed, the presence of isothiazolinones in the environment and their associated potential health hazards have raised significant concerns. In this study, a non-adiabatic dynamics investigation was conducted using state-of-the-art methodologies to explore the photochemistry of isothiazolinones. A simplified model, isothiazol-3(2H)-one (ISO), was employed to represent this compound class. The study validated the model and demonstrated that ISO can return to its ground state through the cleavage of the S-N or S-C bonds, with no significant energy barrier observed. Non-adiabatic dynamics simulations provided insights into the time scales and detailed processes of isothiazolinone photodissociation. The preferred route for deactivation was found to be the cleavage of the S-N bond. This research enhances our understanding of the photodeactivation processes of isothiazolinones and their potential environmental impact.

Quantum-Chemistry Study of the Photophysical Properties of 4-Thiouracil and Comparisons with 2-Thiouracil

DNA in living beings is constantly damaged by exogenous and endogenous agents. However, in some cases, DNA photodamage can have interesting applications, as it happens in photodynamic therapy. In this work, the current knowledge on the photophysics of 4-thiouracil has been extended by further quantum-chemistry studies to improve the agreement between theory and experiments, to better understand the differences with 2-thiouracil, and, last but not least, to verify its usefulness as a photosensitizer for photodynamic therapy. This study has been carried out by determining the most favorable deactivation paths of UV-vis photoexcited 4-thiouracil by means of the photochemical reaction path approach and an efficient combination of the complete-active-space second-order perturbation theory//complete-active-space self-consistent field (CASPT2//CASSCF), (CASPT2//CASPT2), time-dependent density functional theory (TDDFT), and spin-flip TDDFT (SF-TDDFT) methodologies. By comparing the data computed herein for both 4-thiouracil and 2-thiouracil, a rationale is provided on the relatively higher yields of intersystem crossing, triplet lifetime and singlet oxygen production of 4-thiouracil, and the relatively higher yield of phosphorescence of 2-thiouracil.

Minimum-Energy Conical Intersections by Compressed Multistate Pair-Density Functional Theory

Compressed multistate pair-density functional theory (CMS-PDFT) is a multistate version of multiconfiguration pair-density functional theory that can capture the correct topology of coupled potential energy surfaces (PESs) around conical intersections. In this work, we develop interstate coupling vectors (ISCs) for CMS-PDFT in the OpenMolcas and PySCF/mrh electronic structure packages. Yet, the main focus of this work is using ISCs to calculate minimum-energy conical intersections (MECIs) by CMS-PDFT. This is performed using the projected constrained optimization method in OpenMolcas, which uses ISCs to restrain the iterations to the conical intersection seam. We optimize the S1/S0 MECIs for ethylene, butadiene, and benzene and show that CMS-PDFT gives smooth PESs in the vicinities of the MECIs. Furthermore, the CMS-PDFT MECIs are in good agreement with the MECI calculated by the more expensive XMS-CASPT2 method.

The copious photochemistry of 2,6-diaminopurine: Luminescence, triplet population, and ground state recovery

9H- and 7H-2,6-Diaminopurine (26DAP) photoinduced events in vacuum were studied at the MS-CASPT2/cc-pVDZ level of theory. The S1 1 (ππ* La ) state is initially populated evolving barrierless towards its minimum energy structure, from where two photochemical events can take place in both tautomers. The first is the return of the electronic population to the ground state via the C6 conical intersection (CI-C6). The second involves an internal conversion to the ground through the C2 conical intersection (CI-C2). According to our geodesic interpolated paths connecting the critical structures, the second route is less favorable in both tautomers, due to the presence of high energy barriers. Our calculations suggest a competition between fluorescence and ultrafast relaxation to the electronic ground state via internal conversion process. Based on our calculated potential energy surfaces and experimental excited state lifetimes from the literature, we can infer that the 7H- must have a greater fluorescence yield than the 9H-tautomer. We also explored the triplet state population mechanisms on the 7H-26DAP to understand their long-lived components observed experimentally.

Exploration of biochemical reactivity with a QM/MM growing string method

In this work, we have implemented the single-ended growing string method using a hybrid internal/Cartesian coordinate scheme within our in-house QM/MM package, QoMMMa, representing the first implementation of the growing string method in the QM/MM framework. The goal of the implementation was to facilitate generation of QM/MM reaction pathways with minimal user input, and also to improve the quality of the pathways generated as compared to the widely used adiabatic mapping approach. We have validated the algorithm against a reaction which has been studied extensively in previous computational investigations - the Claisen rearrangement catalysed by chorismate mutase. The nature of the transition state and the height of the barrier was predicted well using our algorithm, where more than 88% of the pathways generated were deemed to be of production quality. Directly compared to using adiabatic mapping, we found that while our QM/MM single-ended growing string method is slightly less efficient, it readily produces reaction pathways with fewer discontinuites and thus minimises the need for involved remapping of unsatisfactory energy profiles.

Photophysics of tz Adenine and tz Guanine fluorescent nucleobases embedded into DNA and RNA

UV-VIS photoinduced events of tz A and tz G embedded into DNA and RNA are described by combining the Extended Multi-State Second-Order Perturbation Theory (XMS-CASPT2) and electrostatic embedding molecular mechanics methods (QM/MM). Our results point out that the S1 1 (ππ* La ) state is the bright state in both environments. After the photoexcitation to the S1 1 (ππ* La ) state, the electronic population evolves barrierless towards its minimum, from where the excess of energy can be dissipated by fluorescence. As the minimum energy crossing point structure between the ground and first bright states lies in a high-energy region, the direct internal conversion to the ground state is an unviable mechanism. Other spectroscopic properties (for instance, absorption and Stokes shifts) and comparisons with photochemical properties of canonical nucleobases are also provided.

A-Band Absorption Spectrum of the ClSO Radical: Electronic Structure of the Sulfinyl Group

Sulfur oxide species (RSOx) play a critical role in many fields, ranging from biology to atmospheric chemistry. Chlorine-containing sulfur oxides may play a key role in sulfate aerosol formation in Venus' cloud layer by catalyzing the oxidation of SO to SO2 via sulfinyl radicals (RSO). We present results from the gas-phase UV-vis transient absorption spectroscopy study of the simplest sulfinyl radical, ClSO, generated from the pulsed-laser photolysis of thionyl chloride at 248 nm (at 40 Torr of N2 and 292 K). A weak absorption spectrum from 350 to 480 nm with a peak at 385 nm was observed, with partially resolved vibronic bands (spacing = 226 cm-1), and a peak cross section σ(385 nm) = (7.6 ± 1.9) × 10-20 cm2. From ab initio calculations at the EOMEE-CCSD/ano-pVQZ level, we assigned this band to 12A' ← X2A″ and 22A' ← X2A″ transitions. The spectrum was modeled as a sum of a bound-to-free transition to the 12A' state and a bound-to-bound transition to the 22A' state with similar oscillator strengths; the prediction agreed well with the observed spectrum. We attributed the vibronic structure to a progression in the bending vibration of the 22A' state. Further calculations at the XDW-CASPT2 level predicted a conical intersection between the excited 12A' and 22A' potential energy surfaces near the Franck-Condon region. The geometry of the minimum-energy conical intersection was similar to that of the ground-state geometry. The lack of structure at shorter wavelengths could be evidence of a short excited-state lifetime arising from strong vibronic coupling. From simplified molecular orbital analysis, we attributed the ClSO spectrum to transitions involving the out-of-plane π/π* orbitals along the S-O bond and the in-plane orbital possessing a σ/σ* character along the S-Cl bond. We hypothesize that these orbitals are common to other sulfinyl radicals, RSO, which would share a combination of a strong and a weak transition in the UV (near 300 nm) and visible (400-600 nm) regions.

Smooth Things Come in Threes: A Diabatic Surrogate Model for Conical Intersection Optimization

The optimization of conical intersection structures is complicated by the nondifferentiability of the adiabatic potential energy surfaces. In this work, we build a pseudodiabatic surrogate model, based on Gaussian process regression, formed by three smooth and differentiable surfaces that can adequately reproduce the adiabatic surfaces. Using this model with the restricted variance optimization method results in a notable decrease of the overall computational effort required to obtain minimum energy crossing points.

Resolving Ultrafast Photoinitiated Dynamics of the Hachimoji 5-Aza-7-Deazaguanine Nucleobase: Impact of Synthetically Expanding the Genetic Alphabet

The guanine derivative, 5-aza-7-deazaguanine (^5N7C G) has recently been proposed as one of four unnatural bases, termed Hachimoji (8-letter) to expand the genetic code. We apply steady-state and time-resolved spectroscopy to investigate its electronic relaxation mechanism and probe the effect of atom substitution on the relaxation mechanism in polar protic and polar aprotic solvents. Mapping of the excited state potential energy surfaces is performed, from which the critical points are optimized by using the state-of-art Extended Multi-State Complete Active Space Second-Order Perturbation Theory. It is demonstrated that excitation to the lowest energy ¹ ππ* state of ^5N7C G results in complex dynamics leading to ca. 10 to 30-fold slower relaxation (depending on solvent) compared to guanine. A significant conformational change occurs at the S₁ minimum, resulting in a 10-fold greater fluorescence quantum yield compared to guanine. The fluorescence quantum yield and S₁ decay lifetime increase going from water to acetonitrile to propanol. The solvent-dependent results are supported by the quantum chemical calculations showing an increase in the energy barrier between the S₁ minimum and the S₁ /S₀ conical intersection going from water to propanol. The longer lifetimes might make ^5N7C G more photochemical active to adjacent nucleobases than guanine or other nucleobases within DNA.

Normal-Mode Vibrational Analysis of Weakly Bound Oligomers at Constrained Stationary Points of Arbitrary Order

Following the full realization of the importance of noncovalent interactions, finding and characterizing stationary points (SP), of various order, for weakly bound oligomers have become important tasks for computational chemists. An efficient algorithm and an associated computer code, called oligoCGO, are described, facilitating constrained geometry optimization of oligomers of arbitrary structure and complexity and normal-mode vibrational analysis at nonstationary geometries. To minimize the adverse effects of nonzero forces on harmonic vibrational analyses at constrained stationary points (cSP), two residual gradient correction (RGC) schemes are proposed. RGC₁, for which a rigorous justification is given, is based on ignoring the remaining forces in internal-coordinate space. RGC₂ modifies the geometry of the cSP in a single Newton step and recalculates the Cartesian Hessian at this updated geometry. As demonstrated by 10 examples related to the water-water, water-methane, and methane-methane dimers as well as the methane trimer, without RGC the harmonic analysis of cSPs may result in even qualitatively incorrect results when compared to reference values obtained at the nearby unconstrained SPs (uSP). Both RGC protocols work exceedingly well, and the corrected harmonic wavenumbers of the cSPs are very close to their uSP counterparts.

Role of conical intersection seam topography in the chemiexcitation of 1,2-dioxetanes

Chemiexcitation, the generation of electronic excited states by a thermal reaction initiated on the ground state, is an essential step in chemiluminescence, and it is mediated by the presence of a conical intersection that allows a nonadiabatic transition from ground state to excited state. Conical intersections classified as sloped favor chemiexcitation over ground state relaxation. The chemiexcitation yield of 1,2-dioxetanes is known to increase upon methylation. In this work we explore to which extent this trend can be attributed to changes in the conical intersection topography or accessibility. Since conical intersections are not isolated points, but continuous seams, we locate regions of the conical intersection seams that are close to the configuration space traversed by the molecules as they react on the ground state. We find that conical intersections are energetically and geometrically accessible from the reaction trajectory, and that topographies favorable to chemiexcitation are found in all three molecules studied. Nevertheless, the results suggest that dynamic effects are more important for explaining the different yields than the static features of the potential energy surfaces.

Mechanistic photophysics and photochemistry of unnatural bases and sunscreen molecules: insights from electronic structure calculations

Photophysics and photochemistry are basic subjects in the study of light-matter interactions and are ubiquitous in diverse fields such as biology, energy, materials, and environment. A full understanding of mechanistic photophysics and photochemistry underpins many recent advances and applications. This contribution first provides a short discussion on the theoretical calculation methods we have used in relevant studies, then we introduce our latest progress on the mechanistic photophysics and photochemistry of two classes of molecular systems, namely unnatural bases and sunscreens. For unnatural bases, we disclose the intrinsic driving forces for the ultrafast population to reactive triplet states, impacts of the position and degree of chalcogen substitutions, and the effects of complex environments. For sunscreen molecules, we reveal the photoprotection mechanisms that dissipate excess photon energy to the surroundings by ultrafast internal conversion to the ground state. Finally, relevant theoretical challenges and outlooks are discussed.

Photophysical Deactivation Mechanisms of the Pyrimidine Analogue 1-Cyclohexyluracil

The photophysical relaxation mechanisms of 1-cyclohexyluracil, in vacuum and water, were investigated by employing the Multi-State CASPT2 (MS-CASPT2, Multi-State Complete Active-Space Second-Order Perturbation Theory) quantum chemical method and Dunning's cc-pVDZ basis sets. In both environments, our results suggest that the primary photophysical event is the population of the S11(ππ*) bright state. Afterwards, two likely deactivation pathways can take place, which is sustained by linear interpolation in internal coordinates defined via Z-Matrix scans connecting the most important characteristic points. The first one (Route 1) is the same relaxation mechanism observed for uracil, its canonical analogue, i.e., internal conversion to the ground state through an ethylenic-like conical intersection. The other route (Route 2) is the direct population transfer from the S11(ππ*) bright state to the T23(nπ*) triplet state via an intersystem crossing process involving the (S11(ππ*)/T23(nπ*))STCP singlet-triplet crossing point. As the spin-orbit coupling is not too large in either environment, we propose that most of the electronic population initially on the S11(ππ*) state returns to the ground following the same ultrafast deactivation mechanism observed in uracil (Route 1), while a smaller percentage goes to the triplet manifold. The presence of a minimum on the S11(ππ*) potential energy hypersurface in water can help to understand why experimentally it is noticed suppression of the triplet states population in polar protic solvent.

Electronic Structure Methods for the Description of Nonadiabatic Effects and Conical Intersections

Nonadiabatic effects are ubiquitous in photophysics and photochemistry, and therefore, many theoretical developments have been made to properly describe them. Conical intersections are central in nonadiabatic processes, as they promote efficient and ultrafast nonadiabatic transitions between electronic states. A proper theoretical description requires developments in electronic structure and specifically in methods that describe conical intersections between states and nonadiabatic coupling terms. This review focuses on the electronic structure aspects of nonadiabatic processes. We discuss the requirements of electronic structure methods to describe conical intersections and nonadiabatic couplings, how the most common excited state methods perform in describing these effects, and what the recent developments are in expanding the methodology and implementing nonadiabatic couplings.

Role of Triplet States in the Photodynamics of Aniline

The dynamics of excited heteroaromatic molecules is a key to understanding the photoprotective properties of many biologically relevant chromophores that dissipate their excitation energy nonreactively and thereby prevent the detrimental effects of ultraviolet radiation. Despite their structural variability, most substituted aromatic compounds share a common feature of a repulsive ¹πσ* potential energy surface. This surface can lead to photoproducts, and it can also facilitate the population transfer back to the ground electronic state by means of a ¹πσ*/S₀ conical intersection. Here, we explore a hidden relaxation route involving the triplet electronic state of aniline, which has recently been discovered by means of time-selected photofragment translational spectroscopy [J. Chem. Phys. 2019, 151, 141101]. By using the recently available analytical gradients for multiconfiguration pair-density functional theory, it is now possible to locate the minimum-energy crossing points between states of different spin and therefore compute the intersystem crossing rates with a multireference method, rather than with the less reliable single-reference methods. Using such calculations, we demonstrate that the population loss of aniline in the T₁(³ππ*) state is dominated by C₆H₅NH₂ → C₆H₅NH· + H· dissociation, and we explain the long nonradiative lifetimes of the T₁(³ππ*) state at the excitation wavelengths of 294-264 nm.

On the population of triplet states of 2-seleno-thymine

The population and depopulation mechanisms leading to the lowest-lying triplet states of 2-Se-Thymine were studied at the MS-CASPT2/cc-pVDZ level of theory. Several critical points on different potential energy hypersurfaces were optimized, including minima, conical intersections, and singlet-triplet crossings. The accessibility of all relevant regions on the potential energy hypersurfaces was investigated by means of minimum energy paths and linear interpolation in internal coordinates techniques. Our analysis indicates that, after the population of the bright S2 state in the Franck-Condon region, the first photochemical event is a barrierless evolution towards one of its two minima. After that, three viable photophysical deactivation paths can take place. In one of them, the population in the S2 state is transferred to the T2 state via intersystem crossing and subsequently to the T1 state by internal conversion. Alternatively, the S1 state could be accessed by internal conversion through two distinct conical intersections with S2 state followed by singlet-triplet crossing with the T2 state. The absence of a second minimum on the T1 state and a small energy barrier on pathway along the potential energy surface towards the ground state from the lowest triplet state are attributed as potential reasons to explain why the lifetime of the triplet state of 2-Se-Thymine might be reduced in comparison with its thio-analogue.

Restricted-Variance Constrained, Reaction Path, and Transition State Molecular Optimizations Using Gradient-Enhanced Kriging

Gaussian process regression has recently been explored as an alternative to standard surrogate models in molecular equilibrium geometry optimization. In particular, the gradient-enhanced Kriging approach in association with internal coordinates, restricted-variance optimization, and an efficient and fast estimate of hyperparameters has demonstrated performance on par or better than standard methods. In this report, we extend the approach to constrained optimizations and transition states and benchmark it for a set of reactions. We compare the performance of the newly developed method with the standard techniques in the location of transition states and in constrained optimizations, both isolated and in the context of reaction path computation. The results show that the method outperforms the current standard in efficiency as well as in robustness.

Solvation Structures and Deactivation Pathways of Luminescent Isothiazole-Derived Nucleobases: tz, tz, and tz

The photophysical relaxation pathways of ^tzA, ^tzG, and ^tzI luminescent nucleobases were investigated with the MS-CASPT2 quantum-chemical method and double-ζ basis sets (cc-pVDZ) in gas and condensed phases (1,4-dioxane and water) with the sequential Monte Carlo/CASPT2 and free energy gradient (FEG) methods. Solvation shell structures, in the ground and excited states, were examined with the pairwise radial distribution function (G(r)) and solute-solvent hydrogen-bond networks. Site-specific hydrogen bonding analysis evidenced relevant changes between both electronic states. The three luminescent nucleobases share a common photophysical pattern, summarized as the lowest-lying ¹(ππ*) bright state that is populated directly after the absorption of radiation and evolves barrierless to the minimum energy structure, from where the excess of energy is released by fluorescence. From the ¹(ππ*)_min region, the conical intersection with the ground state ((ππ*/GS)_CI) is not accessible due to the presence of high energetic barriers. By combining the present results with those reported earlier by us for the pyrimidine fluorescent nucleobases, we present a comprehensive description of the photophysical properties of this important class of new fluorescent nucleosides.

Modern quantum chemistry with [Open]Molcas

MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree-Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.

Multireference Electron Correlation Methods: Journeys along Potential Energy Surfaces

Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way and, therefore, can yield accurate and predictive results even when single-reference methods or multiconfigurational self-consistent field theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces. This includes the optimization of molecular geometries, such as equilibrium geometries and conical intersections and on-the-fly photodynamics simulations, both of which depend heavily on the ability of the method to properly explore the potential energy surface. Because such applications require nuclear gradients and derivative couplings, the availability of analytical nuclear gradients greatly enhances the scope of quantum chemical methods. This review focuses on the developments and advances made in the past two decades. A detailed account of the analytical nuclear gradient and derivative coupling theories is presented. Emphasis is given to the software infrastructure that allows one to make use of these methods. Notable applications of multireference electron correlation methods to chemistry, including geometry optimizations and on-the-fly dynamics, are summarized at the end followed by a discussion of future prospects.

Light-Induced Spin Crossover in an Intermediate-Spin Penta-Coordinated Iron(III) Complex

(PMe₃)₂FeCl₃ is an Fe(III) complex that exists in the intermediate-spin ground state in a distorted trigonal bipyramidal geometry. An electronic state with high-spin configuration lies in close vicinity to the ground state, making it a potential spin crossover candidate. A mechanistic account of the spin crossover from the lowest quartet state (Q₀) to the lowest sextet state (S₁) of this complex is provided by exploring both thermal and light-induced pathways. The presence of a large barrier between the two spin states suggests a possible thermal spin crossover at a rather high temperature. The light-induced spin crossover is investigated by employing complete active space self-consistent field calculations together with dynamic correlation and spin-orbit coupling for the lowest seven quartet and lowest five sextet states. The system in the Q₀ state upon light absorption is excited to the optically bright Q₄ LMCT state. By following minimum energy pathways along the electronic states, two light-induced pathways for spin crossover are identified. From the Q₄ state, the system can photo-regenerate the ground intermediate-spin state (Q₀) through an internal conversion of Q₄/Q₃ followed by Q₃/S₁ and S₁/Q₀ intersystem crossings. In an alternate route, through Q₄/S₂ intersystem crossing followed by S₂/S₁ internal conversion, the system can complete the spin crossover from the Q₀ to S₁ state.

Photochemical Relaxation Pathways of 9 H-8-Azaguanine and 8 H-8-Azaguanine

The photochemical reaction path approach, the MS-CASPT2 quantum-chemical method, and double-ζ basis sets (cc-pVDZ) were used to study 9 H-8-azaguanine and 8 H-8-azaguanine relaxation pathways. Several potential energy hypersurfaces were characterized by means of equilibrium geometries, surface crossings (conical intersections and singlet-triplet intersystem crossings), minimum energy paths, and linear interpolation in internal coordinates. The 9 H-8-azaguanine main photochemical event begins with the direct population of the ¹(ππ* L_a) state, which evolves toward a conical intersection with the ground state after surmounting a small energy barrier, explaining why it is nonfluorescent. For 8 H-8-azaguanine, two relaxation mechanisms are possible, depending on the excitation energy. If the S₁ ¹(ππ*) state is initially populated (lower-energy region), the system evolves barrierless to the S₁ ¹(ππ*)_min region, from where the excess energy is released. If the ¹(ππ* L_a) state is populated (higher-energy radiation), the system will encounter conical intersections with the S₂ ¹(n_Oπ*) and S₁ ¹(ππ*) states before evolving to the ¹(ππ* L_a)_min region, from where a conical intersection with the ground state is accessible, favoring radiationless deactivation to the ground state. However, because a fraction of the population can be transferred from ¹(ππ* L_a) to the S₁ ¹(ππ*) state, emission from the S₁ ¹(ππ*)_min region is also expected, although weaker than it would be if the S₁ ¹(ππ*) state were populated directly. Irrespective of the excitation energy, the emissive state is the same and a single fluorescence band is observed, with the strongest emission occurring upon excitation in the lower-energy region, as observed experimentally. Therefore, our computational study corroborates experimental results, attributing the emission of the neutral form of 8-azaguanine in solution to the presence of the minor 8 H-8-azaguanine tautomer, while the 9 H-8-azaguanine major tautomer is nonfluorescent.

Discovery of conical intersection mediated photochemistry with growing string methods

Conical intersections (CIs) are important features of photochemistry that determine yields and selectivity. Traditional CI optimizers require significant human effort and chemical intuition, which typically restricts searching to only a small region of the CI space. Herein, a systematic approach utilizing the growing string method is introduced to locate multiple CIs. Unintuitive MECI are found using driving coordinates that can be generated using a combinatorial search, and subsequent optimization allows reaction pathways, transition states, products, and seam-space pathways to be located. These capabilities are demonstrated by application to two prototypical photoisomerization reactions and the dimerization of butadiene. In total, many reaction pathways were uncovered, including the elusive stilbene hula-twist mechanism, and a previously unidentified product in butadiene dimerization. Overall, these results suggest that growing string methods provide a predictive strategy for exploring photochemistry.

Synthesis and Computational Study of a Pyridylcarbene Fe(II) Complex: Unexpected Effects of fac/ mer Isomerism in Metal-to-Ligand Triplet Potential Energy Surfaces

The synthesis and the steady-state absorption spectrum of a new pyridine-imidazolylidene Fe(II) complex (Fe-NHC) are presented. A detailed mechanism of the triplet metal-to-ligand charge-transfer states decay is provided on the basis of minimum energy path (MEP) calculations used to connect the lowest-lying singlet, triplet, and quintet state minima. The competition between the different decay pathways involved in the photoresponse is assessed by analyzing the shapes of the obtained potential energy surfaces. A qualitative difference between facial ( fac) and meridional ( mer) isomers' potential energy surface (PES) topologies is evidenced for the first time in iron-based complexes. Indeed, the mer complex shows a steeper triplet path toward the corresponding ³MC minimum, which lies at a lower energy as compared to the fac isomer, thus pointing to a faster triplet decay of the former. Furthermore, while a major role of the metal-centered quintet state population from the triplet ³MC region is excluded, we identify the enlargement of iron-nitrogen bonds as the main normal modes driving the excited-state decay.

Combining the Advantages of Alkene and Azo E-Z Photoisomerizations: Mechanistic Insights into Ketoimine Photoswitches

We carried out CASPT2//(TD)DFT and CASPT2//CASSCF studies on the working mechanism of imine switches, including a camphorquinone-derived ketoimine (shortened as k-Imine) switch designed by Lehn as well as a model camphorquinone alkene-imine (a-Imine) proposed in this study. Under the experimental conditions (light irradiation with 455 and 365 nm for E and Z, respectively), k-Imine is excited to the S₁:(n_N,π*) state and then decays toward a perpendicular intermediate following the C═N bond rotation coordinate. During the bond rotation, a mild energy barrier caused by the strong interaction of S₁:(n_N,π*) and S₂:(n_O,π*) states will more or less slow down the rotation speed of k-Imine. The large difference in irradiation light wavelength supports k-Imine as a two-way photoswitch. The photoisomerization of a-Imine obeys a similar but fully barrierless pattern while requiring a higher excitation energy to reach the (n_N,π*) state. The good directionality of thermal isomerization toward E(a-Imine), plus the barrierless photoisomerization, allows for the design of a thermal and photo-operated switch. For both imines, a minimal-energy crossing point (MECI) located at the perpendicular region, with low relative energy and close to the rotary path, ensures the directionality of C═N bond rotation and confirms imines as optimal candidates for photoswitches and motors.

Assessment of the Potential Energy Hypersurfaces in Thymine within Multiconfigurational Theory: CASSCF vs. CASPT2

The present study provides new insights into the topography of the potential energy hypersurfaces (PEHs) of the thymine nucleobase in order to rationalize its main ultrafast photochemical decay paths by employing two methodologies based on the complete active space self-consistent field (CASSCF) and the complete active space second-order perturbation theory (CASPT2) methods: (i) CASSCF optimized structures and energies corrected with the CASPT2 method at the CASSCF geometries and (ii) CASPT2 optimized geometries and energies. A direct comparison between these strategies is drawn, yielding qualitatively similar results within a static framework. A number of analyses are performed to assess the accuracy of these different computational strategies under study based on a variety of numerical thresholds and optimization methods. Several basis sets and active spaces have also been calibrated to understand to what extent they can influence the resulting geometries and subsequent interpretation of the photochemical decay channels. The study shows small discrepancies between CASSCF and CASPT2 PEHs, displaying a shallow planar or twisted ¹(ππ*) minimum, respectively, and thus featuring a qualitatively similar scenario for supporting the ultrafast bi-exponential deactivation registered in thymine upon UV-light exposure. A deeper knowledge of the PEHs at different levels of theory provides useful insight into its correct characterization and subsequent interpretation of the experimental observations. The discrepancies displayed by the different methods studied here are then discussed and framed within their potential consequences in on-the-fly non-adiabatic molecular dynamics simulations, where qualitatively diverse outcomes are expected.

Exploring the isomerization paths of push–pull hexatrienes

The first photo-excited singlet state of all-trans donor (amine), acceptor (cyano)-substituted hexatriene forms donor side-twisted minima and its passage towards the cis–trans–trans isomer is characterized by an S₀/S₁ conical intersection.

Mechanism of excited state deactivation of indan-1-ylidene and fluoren-9-ylidene malononitriles

A joint experimental and computational study on the non-radiative double bond isomerisation decay channel of indan-1-ylidene malononitrile and fluoren-9-ylidene malononitrile is presented in this work.

Molcas 8: New capabilities for multiconfigurational quantum chemical calculations across the periodic table

In this report, we summarize and describe the recent unique updates and additions to the Molcas quantum chemistry program suite as contained in release version 8. These updates include natural and spin orbitals for studies of magnetic properties, local and linear scaling methods for the Douglas-Kroll-Hess transformation, the generalized active space concept in MCSCF methods, a combination of multiconfigurational wave functions with density functional theory in the MC-PDFT method, additional methods for computation of magnetic properties, methods for diabatization, analytical gradients of state average complete active space SCF in association with density fitting, methods for constrained fragment optimization, large-scale parallel multireference configuration interaction including analytic gradients via the interface to the Columbus package, and approximations of the CASPT2 method to be used for computations of large systems. In addition, the report includes the description of a computational machinery for nonlinear optical spectroscopy through an interface to the QM/MM package Cobramm. Further, a module to run molecular dynamics simulations is added, two surface hopping algorithms are included to enable nonadiabatic calculations, and the DQ method for diabatization is added. Finally, we report on the subject of improvements with respects to alternative file options and parallelization.

Towards biochemically relevant QM computations on nucleic acids: controlled electronic structure geometry optimization of nucleic acid structural motifs using penalty restraint functions

Recent developments in dispersion-corrected density functional theory methods allow for the first time the description of large fragments of nucleic acids (hundreds of atoms) with an accuracy clearly surpassing the accuracy of common biomolecular force fields. Such calculations can significantly improve the description of the potential energy surface of nucleic acid molecules, which may be useful for studies of molecular interactions and conformational preferences of nucleic acids, as well as verification and parameterization of other methods. The first of such studies, however, demonstrated that successful applications of accurate QM calculations to larger nucleic acid building blocks are hampered by difficulties in obtaining geometries that are biochemically relevant and are not biased by non-native structural features. We present an approach that can greatly facilitate large-scale QM studies on nucleic acids, namely electronic structure geometry optimization of nucleic acid fragments utilizing a penalty function to restrain key internal coordinates with a specific focus on the torsional backbone angles. This work explores the viability of these restraint optimizations for DFT-D3, PM6-D3H and HF-3c optimizations on a set of examples (a UpA dinucleotide, a DNA G-quadruplex and a B-DNA fragment). Evaluation of different penalty function strengths reveals only a minor system-dependency and reasonable restraint values range from 0.01 to 0.05 Eh rad(-2) for the backbone torsions. Restraints are crucial to perform the QM calculations on biochemically relevant conformations in implicit solvation and gas phase geometry optimizations. The reasons for using restrained instead of constrained or unconstrained optimizations are explained and an open-source external optimizer is provided.

Theoretical study of the hydroxyl radical addition to uracil and photochemistry of the formed U6OH• adduct

Hydroxyl radical ((•)OH) is produced in biological systems by external or endogenous agents. It can damage DNA/RNA by attacking pyrimidine nucleobases through the addition to the C5═C6 double bond. The adduct resulting from the attachment at the C5 position prevails in the experimental measurements, although the reasons for this preference remain unclear. The first aim of this work is therefore to shed light on the comprehension of this important process. Thus, the thermal (•)OH addition to the C5═C6 double bond of uracil has been studied theoretically by using DFT, MP2, and the multiconfigurational CASPT2//CASSCF methodologies. The in-vacuo results obtained with the latter protocol plus the analysis of solvent effects support the experimental observation. A significant lower barrier height is predicted for the C5 pathway with respect to that of the C6 route. In contrast to the C5 adduct, the C6 adduct is able to absorb visible light. Hence, the second aim of the work is to study the photochemistry of this species using the CASPT2//CASSCF methodology within the framework of the photochemical reaction path approach (PRPA). The nonradiative decay to the ground state of this compound has been characterized. A photoreactive character is predicted for the C6 adduct in the excited states according to the presence of excited-state minima along the main decay channel. Finally, a new mechanism of photodissociation has been explored, which implies the photoinduced regeneration of the canonical nucleobase by irradiating with visible light, being therefore relevant in RNA protection against damage by reactive oxygen species.

Revisiting the Nonadiabatic Process in 1,2-Dioxetane

Determining the ground and excited-state decomposition mechanisms of 1,2-dioxetane is essential to understand the chemiluminescence and bioluminescence phenomena. Several experimental and theoretical studies has been performed in the past without reaching a converged description. The reason is in part associated with the complex nonadiabatic process taking place along the reaction. The present study is an extension of a previous work (De Vico, L.; Liu, Y.-J.; Krogh, J. W.; Lindh, R. J. Phys. Chem. A 2007, 111, 8013-8019) in which a two-step mechanism was established for the chemiluminescence involving asynchronous O-O' and C-C' bond dissociations. New high-level multistate multi configurational reference second-order perturbation theory calculations and ab initio molecular dynamics simulations at constant temperature are performed in the present study, which provide further details on the mechanisms and allow to rationalize further experimental observations. In particular, the new results explain the high ratio of triplet to singlet dissociation products.