Can Fulvene S1 Decay Be Controlled? A CASSCF Study with MMVB Dynamics

Ab Initio Multiple Spawning Nonadiabatic Dynamics with Different CASPT2 Flavors: A Fully Open-Source PySpawn/OpenMolcas Interface

We introduce an interface between PySpawn, a simulation package to run ab initio multiple spawning (AIMS) nonadiabatic dynamics, and OpenMolcas, a software package to perform multiconfigurational perturbations theory (CASPT2) electronic structure calculations. Our interface allows us to exploit all the functionalities of the two codes: the modular and efficient Python implementation of the AIMS algorithm and the extensive analysis tools offered by PySpawn, with the cutting-edge implementation of CASPT2 equations in OpenMolcas, including the recently introduced analytical gradients and different flavors. Both are fully open-source and free of charge, making the following implementation unique in the current plethora of software for nonadiabatic dynamics. This represents an important step toward a wider application of AIMS-based nonadiabatic dynamics combined with high-accuracy excited-state calculations. The importance and the need for such an implementation are demonstrated by application to the ultrafast relaxation of fulvene from S1 to S0, which is drastically affected by the potential energy surface on which the nuclear wavepacket is propagated. Additionally, the decay is influenced by the CASPT2 flavor adopted, posing interesting questions in the choice of one over the other and opening the door to deeper studies on the effect of CASPT2 formulations in nonadiabatic dynamics.

Optimizing Conical Intersections without Explicit Use of Non-Adiabatic Couplings

We present two alternative methods for optimizing minimum energy conical intersection (MECI) molecular geometries without knowledge of the derivative coupling (DC). These methods are based on the utilization of Lagrange multipliers: (i) one method uses an approximate calculation of the DC, while the other (ii) do not require the DC. Both methods use the fact that information on the DC is contained in the Hessian of the squared energy difference. Tests done on a set of small molecular systems, in comparison with other methods, show the ability of the proposed methods to optimize MECIs. Finally, we apply the methods to the furimamide molecule, to optimize and characterize its S1/S2 MECI, and to optimizing the S0/S1 MECI of the silver trimer.

Recommendations for Velocity Adjustment in Surface Hopping

This study investigates velocity adjustment directions after hopping in surface hopping dynamics. Using fulvene and a protonated Schiff base (PSB4) as case studies, we investigate the population decay and reaction yields of different sets of dynamics with the velocity adjusted in either the nonadiabatic coupling, gradient difference, or momentum directions. For the latter, in addition to the conventional algorithm, we investigated the performance of a reduced kinetic energy reservoir approach recently proposed. Our evaluation also considered velocity adjustment in the directions of approximate nonadiabatic coupling vectors. While results for fulvene are susceptible to the adjustment approach, PSB4 is not. We correlated this dependence to the topography near the conical intersections. When nonadiabatic coupling vectors are unavailable, the gradient difference direction is the best adjustment option. If the gradient difference is also unavailable, a semiempirical vector direction or the momentum direction with a reduced kinetic energy reservoir becomes an excellent option to prevent an artificial excess of back hoppings. The precise velocity adjustment direction is less crucial for describing the nonadiabatic dynamics than the kinetic energy reservoir's size.

Benchmarking non-adiabatic quantum dynamics using the molecular Tully models

On-the-fly non-adiabatic dynamics methods are becoming more important as tools to characterise the time evolution of a system after absorbing light. These methods, which calculate quantities such as state energies, gradients and interstate couplings at every time step, circumvent the requirement for pre-computed potential energy surfaces. There are a number of different algorithms used, the most common being Tully Surface Hopping (TSH), but all are approximate solutions to the time-dependent Schrödinger equation and benchmarking is required to understand their accuracy and performance. For this, a common set of systems and observables are required to compare them. In this work, we validate the on-the-fly direct dynamics variational multi-configuration Gaussian (DD-vMCG) method using three molecular systems recently suggested by Ibele and Curchod as molecular versions of the Tully model systems used to test one-dimensional non-adiabatic behaviour [Ibele et al., Phys. Chem. Chem. Phys. 2020, 22, 15183-15196]. Parametrised linear vibronic potential energy surfaces for each of the systems were also tested and compared to on-the-fly results. The molecules, which we term the Ibele-Curchod models, are ethene, DMABN and fulvene and the authors used them to test and compare several versions of the Ab Initio Multiple Spawning (AIMS) method alongside TSH. The three systems present different deactivation pathways after excitation to their ππ* bright states. When comparing DD-vMCG to AIMS and TSH, we obtain crucial differences in some cases, for which an explanation is provided by the classical nature and the chosen initial conditions of the TSH simulations.

Sampling effects in quantum mechanical/molecular mechanics trajectory surface hopping non-adiabatic dynamics

The impact of different initial conditions in non-adiabatic trajectory surface hopping dynamics within a hybrid quantum mechanical/molecular mechanics scheme is investigated. The influence of a quantum sampling, based on a Wigner distribution, a fully thermal sampling, based on classical molecular dynamics, and a quantum sampled system, but thermally equilibrated with the environment, is investigated on the relaxation dynamics of solvated fulvene after light irradiation. We find that the decay from the first singlet excited state to the ground state shows high dependency on the initial condition and simulation parameters. The three sampling methods lead to different distributions of initial geometries and momenta, which then affect the fate of the excited state dynamics. We evaluated both the effect of sampling geometries and momenta, analysing how the ultrafast decay of fulvene changes accordingly. The results are expected to be of interest to decide how to initialize non-adiabatic dynamics in the presence of the environment. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.

Fewest switches surface hopping with Baeck-An couplings

In the Baeck-An (BA) approximation, first-order nonadiabatic coupling vectors are given in terms of adiabatic energy gaps and the second derivative of the gaps with respect to the coupling coordinate. In this paper, a time-dependent (TD) BA approximation is derived, where the couplings are computed from the energy gaps and their second time-derivatives. TD-BA couplings can be directly used in fewest switches surface hopping, enabling nonadiabatic dynamics with any electronic structure methods able to provide excitation energies and energy gradients. Test results of surface hopping with TD-BA couplings for ethylene and fulvene show that the TD-BA approximation delivers a qualitatively correct picture of the dynamics and a semiquantitative agreement with reference data computed with exact couplings. Nevertheless, TD-BA does not perform well in situations conjugating strong couplings and small velocities. Considered the uncertainties in the method, TD-BA couplings could be a competitive approach for inexpensive, exploratory dynamics with a small trajectories ensemble. We also assessed the potential use of TD-BA couplings for surface hopping dynamics with time-dependent density functional theory (TDDFT), but the results are not encouraging due to singlet instabilities near the crossing seam with the ground state.

Fewest switches surface hopping with Baeck-An couplings

In the Baeck-An (BA) approximation, first-order nonadiabatic coupling vectors are given in terms of adiabatic energy gaps and the second derivative of the gaps with respect to the coupling coordinate. In this paper, a time-dependent (TD) BA approximation is derived, where the couplings are computed from the energy gaps and their second time-derivatives. TD-BA couplings can be directly used in fewest switches surface hopping, enabling nonadiabatic dynamics with any electronic structure methods able to provide excitation energies and energy gradients. Test results of surface hopping with TD-BA couplings for ethylene and fulvene show that the TD-BA approximation delivers a qualitatively correct picture of the dynamics and a semiquantitative agreement with reference data computed with exact couplings. Nevertheless, TD-BA does not perform well in situations conjugating strong couplings and small velocities. Considered the uncertainties in the method, TD-BA couplings could be a competitive approach for inexpensive, exploratory dynamics with a small trajectories ensemble. We also assessed the potential use of TD-BA couplings for surface hopping dynamics with time-dependent density functional theory (TDDFT), but the results are not encouraging due to singlet instabilities near the crossing seam with the ground state.

AIMSWISS-Ab initio multiple spawning with informed stochastic selections

Ab initio multiple spawning (AIMS) offers a reliable strategy to describe the excited-state dynamics and nonadiabatic processes of molecular systems. AIMS represents nuclear wavefunctions as linear combinations of traveling, coupled Gaussians called trajectory basis functions (TBFs) and uses a spawning algorithm to increase as needed the size of this basis set during nonadiabatic transitions. While the success of AIMS resides in this spawning algorithm, the dramatic increase in TBFs generated by multiple crossings between electronic states can rapidly lead to intractable dynamics. In this Communication, we introduce a new flavor of AIMS, coined ab initio multiple spawning with informed stochastic selections (AIMSWISS), which proposes a parameter-free strategy to beat the growing number of TBFs in an AIMS dynamics while preserving its accurate description of nonadiabatic transitions. The performance of AIMSWISS is validated against the photodynamics of ethylene, cyclopropanone, and fulvene. This technique, built upon the recently developed stochastic-selection AIMS, is intended to serve as a computationally affordable starting point for multiple spawning simulations.

Strategies for Design of Potential Singlet Fission Chromophores Utilizing a Combination of Ground-State and Excited-State Aromaticity Rules

Singlet exciton fission photovoltaic technology requires chromophores with their lowest excited states arranged so that 2E(T1) < E(S1) and E(S1) < E(T2). Herein, qualitative theory and quantum chemical calculations are used to develop explicit strategies on how to use Baird's 4n rule on excited-state aromaticity, combined with Hückel's 4n + 2 rule for ground-state aromaticity, to tailor new potential chromophores for singlet fission. We first analyze the E(T1), E(S1), and E(T2) of benzene and cyclobutadiene (CBD) as excited-state antiaromatic and aromatic archetypes, respectively, and reveal that CBD fulfills the criteria on the state ordering for a singlet fission chromophore. We then look at fulvenes, a class of compounds that can be tuned by choice of substituents from Baird-antiaromatic to Baird-aromatic in T1 and S1 and from Hückel-aromatic to Hückel-antiaromatic in S0. The T1 and S1 states of most substituted fulvenes (159 of 225) are described by singly excited HOMO → LUMO configurations, providing a rational for the simultaneous tuning of E(T1) and E(S1) along an approximate (anti)aromaticity coordinate. Key to the tunability is the exchange integral (KH,L), which ideally is constant throughout the compound class, providing a constant ΔE(S1 - T1). This leads us to a geometric model for the identification of singlet fission chromophores, and we explore what factors limit the model. Candidates with calculated E(T1) values of ∼1 eV or higher are identified among benzannelated 4nπ-electron compound classes and siloles. In brief, it is clarified how the joint utilization of Baird's 4n and Hückel's 4n + 2 rules, together with substituent effects (electronic and steric) and benzannelation, can be used to tailor new chromophores with potential use in singlet fission photovoltaics.

Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials

Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.

A molecular perspective on Tully models for nonadiabatic dynamics

We present a series of standardized molecular tests for nonadiabatic dynamics, reminiscent of the one-dimensional Tully models proposed in 1990.

An overview of the cycloaddition chemistry of fulvenes and emerging applications

The unusual electronic properties and unique reactivity of fulvenes have interested researchers for over a century. The propensity to form dipolar structures at relatively low temperatures and to participate as various components in cycloaddition reactions, often highly selectively, makes them ideal for the synthesis of complex polycyclic carbon scaffolds. As a result, fulvene cycloaddition chemistry has been employed extensively for the synthesis of natural products. More recently, fulvene cycloaddition chemistry has also found application to other areas including materials chemistry and dynamic combinatorial chemistry. This highlight article discusses the unusual properties of fulvenes and their varied cycloaddition chemistry, focussing on applications in organic and natural synthesis, dynamic combinatorial chemistry and materials chemistry, including dynamers, hydrogels and charge transfer complexes. Tables providing comprehensive directories of fulvene cycloaddition chemistry are provided, including fulvene intramolecular and intermolecular cycloadditions complete with reactant partners and their resulting cyclic adducts, which provide a useful reference source for synthetic chemists working with fulvenes and complex polycyclic scaffolds.

Proton-Gated Photoisomerization of Amino-Substituted Dibenzofulvene Rotors

Derivatives of 9-(2,2,2-triphenylethylidene)-fluorene (TEF) undergo E/Z photoisomerization and are candidates for light-powered molecular actuators and switches. The 2-substituted derivatives bearing an amino group (ATEF) or a dimethylamino group (DTEF) are weakly photoactive in the absence of acids, but protonation increases photoisomerization quantum yields by factors of 30-60. Such compounds may be useful for incorporation into pH-switchable photoactive devices.

Synthesis and Photoisomerization of Substituted Dibenzofulvene Molecular Rotors

The synthesis, spectral and structural characterization, and photoisomerization of a family of 2-substituted dibenzofulvene molecular actuators based on (2,2,2-triphenylethylidene)fluorene (TEF) are reported. The 2-substituted species investigated are nitro (NTEF), cyano (CTEF), and iodo (ITEF). X-ray structures of these three compounds and three intermediates were determined to assign alkene configuration and investigate the effects of the 2-substituents on steric gearing. The addition-elimination reaction of Z-9 with trityl anion to form Z-10 proceeded with complete retention of configuration. Rates of photoisomerization were measured at irradiation wavelengths between 266-355 nm in acetonitrile/dioxane solutions at room temperature. Photoisomerization quantum yields (φ) were calculated by means of a mathematical model that accounts for a certain degree of photodecomposition in the cases of CTEF and ITEF. Quantum yields vary significantly with substituent, having maximum values of φ=0.26 for NTEF, 0.39 for CTEF, and 0.50 for ITEF. NTEF is photochemically robust and has a large quantum yield for photoisomerization in the near-UV, making it a particularly promising drive rotor moiety for light-powered molecular devices.

Hybrid QTAIM and electrostatic potential-based quantum topology phase diagrams for water clusters

The topological diversity of sets of isomers of water clusters (W = H₂O)_n, 7 ≤ n ≤ 10, is analyzed employing the scalar fields of total electronic charge density ρ(r) and the molecular electrostatic potential (MESP).

Non-adiabatic dynamics close to conical intersections and the surface hopping perspective

Conical intersections play a major role in the current understanding of electronic de-excitation in polyatomic molecules, and thus in the description of photochemistry and photophysics of molecular systems. This article reviews aspects of the basic theory underlying the description of non-adiabatic transitions at conical intersections, with particular emphasis on the important case when the dynamics of the nuclei are treated classically. Within this classical nuclear motion framework, the main aspects of the surface hopping methodology in the conical intersection context are presented. The emerging picture from this treatment is that of electronic transitions around conical intersections dominated by the interplay of the nuclear velocity and the derivative non-adiabatic coupling vector field.

Seeking small molecules for singlet fission: a heteroatom substitution strategy

We design theoretically small molecule candidates for singlet fission chromophores, aiming to achieve a balance between sufficient diradical character and kinetic persistence. We develop a perturbation strategy based on the captodative effect to introduce diradical character into small π-systems. Specifically, this can be accomplished by replacing pairs of not necessarily adjacent C atoms with isoelectronic and isosteric pairs of B and N atoms. Three rules of thumb emerge from our studies to aid further design: (i) Lewis structures provide insight into likely diradical character; (ii) formal radical centers of the diradical must be well-separated; (iii) stabilization of radical centers by a donor (N) and an acceptor (B) is essential. Following the rules, we propose candidate molecules. Employing reliable multireference calculations for excited states, we identify three likely candidate molecules for SF chromophores. These include a benzene, a napthalene, and an azulene, where four C atoms are replaced by a pair of B and a pair of N atoms.

Photochemistry and photophysics at extended seams of conical intersection

The role of extended seams of conical intersection in excited-state mechanisms is reviewed. Seams are crossings of the potential energy surface in many dimensions where the decay from the excited to the ground state can occur, and the extended seam is composed of different segments lying along a reaction coordinate. Every segment is associated with a different primary photoproduct, which gives rise to competing pathways. This idea is first illustrated for fulvene and ethylene, and then it is used to explain more complex cases such as the dependence of the isomerisation of retinal chromophore isomers on the protein environment, the dependence of the efficiency of the azobenzene photochemical switch on the wavelength of irradiation and the direction of the isomerisation, and the coexistence of different mechanisms in the photo-induced Wolff rearrangement of diazonaphthoquinone. The role of extended seams in the photophysics of the DNA nucleobases and the relationship between two-state seams and three-state crossings is also discussed. As an outlook, the design of optical control strategies based on the passage of the excited molecule through the seam is considered, and it is shown how the excited-state lifetime of fulvene can be modulated by shaping the energy of the seam.

The response of the electronic structure to electronic excitation and double bond torsion in fulvene: a combined QTAIM, stress tensor and MO perspective

New insights into the double bond isomerization of fulvene in the ground and excited electronic states are provided by newly developed QTAIM and stress tensor tools. The S0 and S1 states follow the 'biradical' torsion model, but the double bond is stiffer in the S0 state; by contrast, the S2 state follows the 'zwitterionic' torsion. Differences are explained in terms of the ellipticity and bond critical point (BCP) stiffness for both QTAIM and the stress tensor. Overall, the wave-function based analysis is found to be in agreement with the work of Bonačić-Koutecký and Michl that the bond-twisted species can have biradical or zwitterionic character, depending on the state. Using QTAIM and the stress tensor a new understanding of bond torsion is revealed; the electronic charge density around the twisted bond is found not to rotate in concert with the nuclei of the rotated -CH2 methylene group. The ability to visualize how the bond stiffness varies between individual electronic states and how this correlates with the QTAIM and stress tensor bond stiffness is highlighted. In addition, the most and least preferred morphologies of bond-path torsion are visualized. Briefly we discuss the prospects for using this new QTAIM and stress tensor analysis for excited state chemistry.

Synthesis of functionalized fulvenes: [3 + 2] annulation of ethyl α-chlorocyclopropaneformates with 1,3-dicarbonyl compounds

An efficient method is reported to produce acidic multi-substituted fulvenes at high yields (up to 94%) through the base-promoted [3 + 2] annulation reaction of ethyl α-chlorocyclopropaneformates with 1,3-dicarbonyl compounds.

Photophysics of fulvene under the non-resonant stark effect. Shaping the conical intersection seam

We introduce a mechanistic strategy to control the excited state lifetime of fulvene based on shaping the topography of an extended seam of intersection with the non-resonant dynamic Stark effect. Fulvene has a very short excited state lifetime due to an energetically accessible seam of intersection which lies along the methylene torsion coordinate, and the initial decay occurs at the seam segment around the planar conical intersection structure. We have followed a three-step approach to simulate the control. First, we have calculated the effect of a non-resonant electric field on the potential energy surface at the ab initio level, including the field in a self-consistent way. The relative energy of the planar segment of the seam is increased by the non-resonant field. In the second step we simulate the control carrying out MCTDH quantum dynamics propagations under a static non-resonant field to derive the main control mechanisms. At moderately intense fields (epsilon < or = 0.03 a.u.) the decay is faster as compared to the field free case because the vibrational overlap between the excited and ground state vibrational functions is increased. However, at more intense fields (epsilon = 0.04 a.u.) the planar conical intersection is energetically inaccessible and the decay occurs at a slower time scale, at the segment of the seam with more twisted geometries. In the third step, the control over the dynamics is exerted with a non-resonant dynamic field. The acceleration of the decay due to the improved vibrational overlap does not occur, but the decay can be made slower with a dynamic field of 0.08 a.u. The results show the viability of our approach to control the photophysics shaping the topology of the conical intersection seam, and they prove that the extended nature of the seam is crucial for simulating and understanding the control.

Exploring the Photophysical Properties of Molecular Systems Using Excited State Accelerated ab Initio Molecular Dynamics

In the present work, we employ excited state accelerated ab initio molecular dynamics (A-AIMD) to efficiently study the excited state energy landscape and photophysical topology of a variety of molecular systems. In particular, we focus on two important challenges for the modeling of excited electronic states: (i) the identification and characterization of conical intersections and crossing seams, in order to predict different and often competing radiationless decay mechanisms, and (ii) the description of the solvent effect on the absorption and emission spectra of chemical species in solution. In particular, using as examples the Schiff bases formaldimine and salicylidenaniline, we show that A-AIMD can be readily employed to explore the conformational space around crossing seams in molecular systems with very different photochemistry. Using acetone in water as an example, we demonstrate that the enhanced configurational space sampling may be used to accurately and efficiently describe both the prominent features and line-shapes of absorption and emission spectra.

Quantum dynamics study of fulvene double bond photoisomerization: the role of intramolecular vibrational energy redistribution and excitation energy

The double bond photoisomerization of fulvene has been studied with quantum dynamics calculations using the multi-configuration time-dependent Hartree method. Fulvene is a test case to develop optical control strategies based on the knowledge of the excited state decay mechanism. The decay takes place on a time scale of several hundred femtoseconds, and the potential energy surface is centered around a conical intersection seam between the ground and excited state. The competition between unreactive decay and photoisomerization depends on the region of the seam accessed during the decay. The dynamics are carried out on a four-dimensional model surface, parametrized from complete active space self-consistent field calculations, that captures the main features of the seam (energy and locus of the seam and associated branching space vectors). Wave packet propagations initiated by single laser pulses of 5-25 fs duration and 1.85-4 eV excitation energy show the principal characteristics of the first 150 fs of the photodynamics. Initially, the excitation energy is transferred to a bond stretching mode that leads the wave packet to the seam, inducing the regeneration of the reactant. The photoisomerization starts after the vibrational energy has flowed from the bond stretching to the torsional mode. In our propagations, intramolecular energy redistribution (IVR) is accelerated for higher excess energies along the bond stretch mode. Thus, the competition between unreactive decay and isomerization depends on the rate of IVR between the bond stretch and torsion coordinates, which in turn depends on the excitation energy. These results set the ground for the development of future optical control strategies.