The on-the-fly surface-hopping program system Newton-X: Application to ab initio simulation of the nonadiabatic photodynamics of benchmark systems

Revisiting the intricate photodissociation mechanism of ammonia along the minor NH + H2 pathway

In this manuscript we revisit and extend the analysis of the internal conversion and intersystem crossing dynamics of ammonia originating from excitation to its first excited singlet state S1 and leading to the major NH2 + H and minor NH + H2 pathways as investigated by Wang et al. (Phys. Chem. Chem. Phys. 2022, 24, 15060). To perform extensive simulations, we use the machine-learned interatomic potentials, developed by the same authors (Y. Wang et al., J. Chem. Phys. 2021, 154, 094121) interfaced with the SHARC program package. Based on the analysis of 50 000 coupled trajectories, we fit a kinetic model for the major and minor reaction channels that also provides the minimal amount of time required for the dissociation to occur. The model predicts that the time constant associated with the rare pathway is two orders of magnitude larger than that of the frequent reaction, leading to an extrapolated quantum yield of 0.64% for the NH + H2 photoproduct. This quantum yield is in agreement with available experimental measurements carried out with an excitation pulse of 193 nm, around which we excite ammonia. A comprehensive analysis of the electronic states and structures involved in the minor channel reveals that dissociation occurs in a concerted manner via three main mechanisms. The majority of the trajectories (98%) undergo nonradiative relaxation to the electronic ground state, from where 29% directly dissociate. Additionally, we observe reverse internal conversion (58%) as well as intersystem crossing (10%), as operative pathways responsible of the rare photodissociation reaction. These findings provide valuable insights into the dynamics of ammonia photodissociation, particularly its less-studied fragmentation pathway.

Theoretical insights on the double ESIPT mechanism and fluorescence properties of H2BIo chromophore

2-{[3-(1H-benzoimidazol-2-yl)-2-hydroxy-5-methylbenzylidene] amino}-benzoic acid (H2BIo) based on proton transfer can serve as the fluorescent probe for detecting heavy metal ions. The excited-state intramolecular proton transfer (ESIPT) reaction mechanism of the H2BIo chromophore with an intramolecular asymmetric double hydrogen bond in different solvents are investigated. The reaction barrier of the ESIPT along hydrogen bond O1-H2···N3 is higher than that of ESIPT along O4-H5···N6, which indicates that the double ESIPT is a stepwise process. The time-evolving non-adiabatic excited-state dynamic simulations shows that the sequence ESIPT reactions on a time scale: the ESIPT along O1-H2···N3 is faster than the ESIPT along O4-H5···N6. The analyses of electron structure and spectra indicate that the double ESIPT couples with electron transfer, significantly enhances the fluorescence signal, thereby improving the performance of the fluorescent probe in detecting heavy metal ions.

SPADE 1.0: A Simulation Package for Non-Adiabatic Dynamics in Extended Systems

Nonadiabatic molecular dynamics (NAMD) simulations are crucial for revealing the underlying mechanisms of photochemical and photophysical processes. Typical NAMD simulation software packages rely on on-the-fly ab initio electronic structure and nonadiabatic coupling calculations, and thus become challenging when dealing with large complex systems. We here introduce a new Simulation Package for non-Adiabatic Dynamics in Extended systems (SPADE), which is designed to address the limitations of traditional surface hopping methods in dealing with these problems. By design, SPADE enables the users to define arbitrary quasi-diabatic Hamiltonians through parametrized functions and incorporates a variety of algorithms (e.g., global flux hopping probabilities, complex crossing and decoherence corrections), which can realize efficient and reliable NAMD simulations without using nonadiabatic couplings at all. All the employed methods and expressions for diabatic Hamiltonian matrix elements can be flexibly set through the input files. SPADE is mainly written in Fortran based on a modular design and has a great capacity for further implementation of new methods. SPADE can be used to simulate both model and atomistic systems as long as proper Hamiltonians are provided. As demonstrations, a series of representative models are studied to show the main features and capabilities.

Polarized molecular wires for efficient photo-generation of free electric charge carriers

In this study, five all-organic polarized molecular wires (PMWs) for photovoltaic applications are introduced and investigated by means of theoretical chemistry methods. Structurally based on polarized pyrrole and isoindole moieties, the proposed systems demonstrate efficient and ultrafast charge separation upon light absorption. The relevant PMW properties are evaluated using semi-empirical and ab initio quantum-chemical calculations, as well as nonadiabatic molecular dynamics simulations. For each system, the individual mechanism of charge carrier separation is identified and characterized, with all proposed wires exhibiting remarkable charge carrier separation efficiencies and rates. Furthermore, the designed PMW structures enable their straightforward incorporation into more extended molecular frameworks.

Static and Dynamic Studies of Excitation in a Fullerene-Anthracene Complex

We investigate the static and dynamic aspects of an anthracene-fullerene complex. Our detailed investigation includes (a) systematic scrutiny of the complex in its stable state through the examination of its chemical reactivity parameters and absorption spectra and (b) exploration of its dynamic behavior in excited states at a femtosecond time scale. To achieve this, we employ a combination of ab initio molecular dynamics (AIMD) and time-dependent density functional theory (TDDFT) to study not only the system's behavior in the excited state but also the temporal evolution of chemical reactivity parameters when it moves in a particular excited state. It shows greater reactivity in the excited state as compared to that in the ground state. Interestingly, our findings reveal that the complex can even switch between excited states during its movement in certain trajectories. Accordingly, we conduct an extensive examination of the interactions among coupling components, which coincides with the occurrence of trajectory surface hopping.

Theoretical Insights into Ultrafast Separation of Photogenerated Charges in a Push-Pull Polarized Molecular Triad

Herein, we propose a purely-organic donor-acceptor (D-A) molecular triad, with a light-absorbing polarized molecular wire (PMW) used as a central linkage, as a proof of concept for the possible future applications of the D-PMW-A arrangement in molecular photovoltaics. This work builds upon our earlier study on the PMW unit itself, which proved to be highly promising for the ultrafast photogeneration of free charge carriers. Quantum-chemical calculations performed for the D-PMW-A triad at a semi-empirical level of theory reveal a large electric dipole moment of the system, and show strong charge-transfer (CT) character of its lowest-energy excited electronic states, including theS 1 ${S_1 }$ , which favours efficient dissociation of an exciton initially formed upon the absorption of light. The confirmation for this effect was found with nonadiabatic molecular dynamics simulations, revealing an ultrafast relaxation from higher, bright excited states toS 1 ${S_1 }$ , completed on a subpicosecond timescale. The architecture of the proposed molecular triad enables its electronic coupling to the surrounding environment through chemical bonds, or noncovalent stacking interactions, which might open way for synthesis of a new class of D-PMW-A efficient molecular organic photovoltaic materials.

Coherent Excitation of the CH Stretching Vibrations in C2H4 +: The Role of the Derivative Coupling Studied by the Quantum Ehrenfest Method

We report nonadiabatic dynamics computations on C2H4 + initiated on a coherent superposition of the five lowest cationic states, employing the Quantum Ehrenfest method. In addition to the totally symmetric carbon-carbon double bond stretch and carbon-hydrogen stretches, we see that the three non-totally symmetric modes become stimulated; torsion and three different CH stretching patterns. Thus, a coherent superposition of states, of the type involved in an attochemistry experiment, leads to the stimulation of specific non-totally symmetric motions. The computations were also performed on the specific combination of the A and C states. In each case normal mode 9 (cis-asymmetric H2CCH2 stretch), out of the set of non-totally-symmetric normal modes, dominates. Thus, we can steer the nuclear motion along specific non-totally symmetric normal modes using a defined coherent superposition.

Electronic Spectroscopy of the Halogenated Criegee Intermediate, ClCHOO: Experiment and Theory

A chlorine-substituted Criegee intermediate, ClCHOO, is photolytically generated using a diiodo precursor, detected by VUV photoionization at 118 nm, and spectroscopically characterized via ultraviolet-visible (UV-vis)-induced depletion of m/z = 80 under jet cooled conditions. UV-vis excitation resonant with a π* ← π transition yields a significant ground state depletion, indicating a strong electronic transition and rapid photodissociation. The broad absorption spectrum peaks at 350 nm and is attributed to contributions from both syn (∼70%) and anti (∼30%) conformers of ClCHOO based on spectral simulations using a nuclear ensemble method. Electronic structure theory shows significant differences in the vertical excitation energies of the two conformers (330 and 371 nm, respectively) as well as their relative stabilities in the ground and excited electronic states associated with the π* ← π transition. Natural bond orbital analysis reveals significant and nonintuitive nonbonding contributions to the relative stabilities of the syn and anti conformers.

Real-Time CASSCF (Ehrenfest) Modeling of Electron Dynamics in Organic Semiconductors. Dynamics Reaction Paths Driven by Quantum Coherences. Application to a Radical Organic Semiconductor

We present a strategy for the modeling of charge carrier dynamics in organic semiconductors using conventional quantum chemistry methods, including the analytic gradient for nuclear motion. The theoretical approach uses real-time CASSCF (Ehrenfest) all-electron dynamics coupled to classical nuclear dynamics for the special case of a small number (4-8) of molecular units. The objective is to obtain mechanistic/atomistic insight at the electronic structure level, relating to spin density dynamics, to the effect of crystal structure (e.g., slippage between spin/charge carriers), and to ferromagnetic and antiferromagnetic effects. The initial conditions for our simulations use the equilibrium structures of all the molecular units. At this geometry, a localized hole on one of the units corresponds to a coherent superposition of adiabatic states. We thus generate a dynamics reaction path driven by quantum coherences. Our aim is to inform experiment and to compare with parametrized theoretical models. The methodology is demonstrated for a perfectly π-stacked ethylene model (up to 8 eclipsed molecular units) for both hole transfer and localized exciton transfer. An application for hole transfer is presented for bisdithiazolyl (S,S) and bisdiselenazolyl (Se,Se) radicals for the special case of ferromagnetic coupling. For these examples, the embedded pyridine radical model organic chromophore (up to 6 eclipsed π-stacked molecular units) has been studied on its own as well as the target bisdithiazolyl (S,S) and bisdiselenazolyl (Se,Se) systems. A significant difference between these systems and the ethylene and pyridine stacks is that the (S,S) and (Se,Se) systems exhibit molecular slippage rather than being perfectly eclipsed. This slippage may result from crystal defects or intermolecular vibrations. For the model systems, the electron dynamics is dominated by the initial and final molecular units, irrespective of the length of the chain. The intervening units act as a "superexchange bridge". Our simulations reveal that, in the presence of slippage, charge migration cannot propagate across the entire system; instead, the coherence length is limited to 3 molecular units. The results also suggest that the mechanism of charge transport is different for bisdiselenazolyl (Se,Se) (superexchange-like A -[B]→ C) and bisdithiazolyl (S,S) (direct A → C). An analysis of the spin density suggests that, in the charge carrier dynamics, the additional charge carried by the Se versus S in the "scaffold" is small. Since we use a small number of molecular units, the coupled nuclear dynamics is seen to be complementary to the electron dynamics (i.e., creating a hole causes bond length contraction while filling a hole with an electron lengthens the bond). In all the cases studied, the mechanism of charge mobility is wave-like, rather than hopping, because we use the time dependent Schrödinger equation to propagate the electronic wave function.

Non-adiabatic Couplings in Surface Hopping with Tight Binding Density Functional Theory: The Case of Molecular Motors

Nonadiabatic molecular dynamics (NAMD) has become an essential computational technique for studying the photophysical relaxation of molecular systems after light absorption. These phenomena require approximations that go beyond the Born-Oppenheimer approximation, and the accuracy of the results heavily depends on the electronic structure theory employed. Sophisticated electronic methods, however, make these techniques computationally expensive, even for medium size systems. Consequently, simulations are often performed on simplified models to interpret the experimental results. In this context, a variety of techniques have been developed to perform NAMD using approximate methods, particularly density functional tight binding (DFTB). Despite the use of these techniques on large systems, where ab initio methods are computationally prohibitive, a comprehensive validation has been lacking. In this work, we present a new implementation of trajectory surface hopping combined with DFTB, utilizing nonadiabatic coupling vectors. We selected the methaniminium cation and furan systems for validation, providing an exhaustive comparison with the higher-level electronic structure methods. As a case study, we simulated a system from the class of molecular motors, which has been extensively studied experimentally but remains challenging to simulate with ab initio methods due to its inherent complexity. Our approach effectively captures the key photophysical mechanism of dihedral rotation after the absorption of light. Additionally, we successfully reproduced the transition from the bright to dark states observed in the time-dependent fluorescence experiments, providing valuable insights into this critical part of the photophysical behavior in molecular motors.

Nonadiabatic Excited-State Molecular Dynamics with an Explicit Solvent: NEXMD-SANDER Implementation

In this article, the nonadiabatic excited-state Molecular dynamics (NEXMD) package is linked with the SANDER package, provided by AMBERTOOLS. The combination of these software packages enables the simulation of photoinduced dynamics of large multichromophoric conjugated molecules involving several coupled electronic excited states embedded in an explicit solvent by using the quantum/mechanics/molecular mechanics (QM/MM) methodology. The fewest switches surface hopping algorithm, as implemented in NEXMD, is used to account for quantum transitions among the adiabatic excited-state simulations of the photoexcitation and subsequent nonadiabatic electronic transitions, and vibrational energy relaxation of a substituted polyphenylenevinylene oligomer (PPV3-NO2) in vacuum and methanol as an explicit solvent has been used as a test case. The impact of including specific solvent molecules in the QM region is also analyzed. Our NEXMD-SANDER QM/MM implementation provides a useful computational tool to simulate qualitatively solvent-dependent effects, like electron transfer, stabilization of charge-separated excited states, and the role of solvent reorganization in the molecular optical properties, observed in solution-based spectroscopic experiments.

Effect of load-resisting force on photoisomerization mechanism of a single second generation light-driven molecular rotary motor

A pivotal aspect of molecular motors is their capability to generate load capacity from a single entity. However, few studies have directly characterized the load-resisting force of a single light-driven molecular motor. This research provides a simulation analysis of the load-resisting force for a highly efficient, second-generation molecular motor developed by Feringa et al. We investigate the M-to-P photoinduced nonadiabatic molecular dynamics of 9-(2,3-dihydro-2-methyl-1H-benz[e]inden-1-ylidene)-9H-fluorene utilizing Tully's surface hopping method at the semi-empirical OM2/MRCI level under varying load-resisting forces. The findings indicate that the quantum yield remains relatively stable under forces up to 0.003 a.u., with the photoisomerization mechanism functioning typically. Beyond this threshold, the quantum yield declines, and an alternative photoisomerization mechanism emerges, characterized by an inversion of the central double bond's twisting direction. The photoisomerization process stalls when the force attains a critical value of 0.012 a.u. Moreover, the average lifetime of the excited state oscillates around that of the unperturbed system. The quantum yield and mean lifetime of the S1 excited state in the absence of external force are recorded at 0.54 and 877.9 fs, respectively. In addition, we analyze a time-dependent fluorescence radiation spectrum, confirming the presence of a dark state and significant vibrations, as previously observed experimentally by Conyard et al.

Time-resolved photoelectron diffraction imaging of methanol photodissociation involving molecular hydrogen ejection

Imaging ultrafast atomic and molecular hydrogen motion with femtosecond time resolution is a challenge for ultrafast spectroscopy due to the low mass and small scattering cross section of the moving neutral hydrogen atoms and molecules. Here, we propose time- and momentum-resolved photoelectron diffraction (TMR-PED) as a way to overcome limitations of existing methodologies and illustrate its performance using a prototype molecular dissociation process involving the sequential ejection of a neutral hydrogen molecule and a proton from the methanol dication. By combining state-of-the-art molecular dynamics and electron-scattering methods, we show that TMR-PED allows for direct imaging of hydrogen atoms in action. More specifically, the fingerprint of hydrogen dynamics reflects the time evolution of polarization-averaged molecular-frame photoelectron angular distributions (PA-MFPADs) as would be recorded in X-ray pump/X-ray probe experiments with few-femtosecond resolution. We present the results of two precursor experiments that support the feasibility of this approach.

Theoretical Investigation on Proton Transfer Directionality and Dynamics Behavior of 3-(Benzo[d]thiazol-2-yl)-2-hydroxy-5-methoxybenzaldehyde with Two Asymmetric Proton Acceptors

A detailed theoretical investigation on the excited state intramolecular proton transfer (ESIPT) directionality and dynamics behavior of 3-(benzo[d]thiazol-2-yl)-2-hydroxy-5-methoxybenzaldehyde (BTHMB) with two unsymmetric proton acceptors (N and O2) has been performed. The hydrogen bond O1-H···N in BTHMB-a formed by the O1-H group with the N atom or O1-H···O2 in BTHMB-b formed by the O1-H group with the O2 atom is enhanced upon photoexcitation, and the strength of the O1-H···N bond is stronger, which will drive the O1-H proton to the N atom. Potential energy curves further confirm that ESIPT occurs in the N atom because of the smaller energy barrier (0.39 kcal/mol). Results of dynamics simulations manifest that no surface hopping exists between the S0 and S1 states within 300 fs, and ESIPT time constants of BTHMB-a and BTHMB-b are 48 and 151 fs, respectively. While the reverse ESIPT is observed in BTHMB-b at 294 fs, implying that the O1-H proton is transferred to the N atom instead of the O2 atom. The consistency of the calculated absorption (390 nm) and fluorescence spectra (443 and 602 nm) of BTHMB-a with the experimental values (390, 410, and 605 nm) confirms this conclusion again. The charge distribution analysis shows that the charge on the proton acceptors increases, and the O2 atom has higher electronegativity because it has more negative charges. The minimum surface electrostatic potential on the N atom in BTHMB-b correlating with the pKb value is -47.38 kcal/mol, indicating that the N atom has strong basicity. Therefore, the basicity of the N atom dominates the ESIPT process rather than the electronegativity of the O2 atom.

New theoretical insights on the nonradiative relaxation mechanism of the core structure of mycosporines: The amino-cyclohexenone central template

We present a comprehensive computational study describing the excited state dynamics and consequent photostability of amino-cyclohexenone (ACyO), the central template of mycosporine systems, widely recognized for their photoprotection of aquatic species. Photoexcitation to the first excited electronic state (S1, 1nπ*) of ACyO is considered an optically dark transition, while photoexcitation to the second excited electronic state (S21ππ*) is an optically bright 1ππ* transition and largely responsible for UV absorption properties of this molecule. We show that following initial photoexcitation to S2, ACyO relaxes via two competing deactivation mechanisms, each mediated by an S1/S0 conical intersection, which directs the excited state population to the electronic ground state (S0). Our ab initio computational results are supported with nonadiabatic dynamics simulation results, yielding an excited state lifetime of ∼280 fs for this system in vacuo. These results explain the inherent photostability of this core structure, commonplace in a wide range of microorganisms.

Transient-absorption spectroscopy of dendrimers via nonadiabatic excited-state dynamics simulations

The efficiency of light-harvesting and energy transfer in multi-chromophore ensembles underpins natural photosynthesis. Dendrimers are highly branched synthetic multi-chromophoric conjugated supra-molecules that mimic these natural processes. After photoexcitation, their repeated units participate in a number of intramolecular electronic energy relaxation and redistribution pathways that ultimately funnel to a sink. Here, a model four-branched dendrimer with a pyrene core is theoretically studied using nonadiabatic molecular dynamics simulations. We evaluate excited-state photoinduced dynamics of the dendrimer, and demonstrate on-the-fly simulations of its transient absorption pump-probe (TA-PP) spectra. We show how the evolutions of the simulated TA-PP spectra monitor in real time photoinduced energy relaxation and redistribution, and provide a detailed microscopic picture of the relevant energy-transfer pathways. To the best of our knowledge, this is the first of this kind of on-the-fly atomistic simulation of TA-PP signals reported for a large molecular system.

Ultrafast Hot Exciton Nonadiabatic Excited-State Dynamics in Green Fluorescent Protein Chromophore Analogue

The ultrafast high-energy nonadiabatic excited-state dynamics of the benzylidenedimethylimidazolinone chromophore dimer has been investigated using an electronic structure method coupled with on-the-fly quantitative wave function analysis to gain insight into the photophysics of hot excitons in biological systems. The dynamical simulation provides a rationalization of the behavior of the exciton in a dimer after the photoabsorption of light to higher-energy states. The results suggest that hot exciton localization within the manifold of excited states is caused by the hindrance of torsional rotation due to imidazolinone (I) or phenolate (P) bonds i.e., ΦI- or ΦP-dihedral rotation, in the monomeric units of a dimer. This hindrance arises due to weak π-π stacking interaction in the dimer, resulting in an energetically uphill excited-state barrier for ΦI- and ΦP-twisted rotation, impeding the isomerization process in the chromophore. Thus, this study highlights the potential impact of the weak π-π interaction in regulating the photodynamics of the green fluorescent protein chromophore derivatives.

New insights into the photocyclization reaction of a popular diarylethene switch: a nonadiabatic molecular dynamics study

Diarylethene (DAE) molecular switches have continued to attract the attention of researchers for over 20 years. Their remarkable photophysical properties endow them with countless applications in photonics and molecular technologies. However, despite extensive experimental and theoretical research, the mechanism of DAE photoswitching is not yet fully rationalized. In this work, we investigate the ring closure dynamics of a popular DAE switch, 1,2-bis(3-methyl-5-phenyl-2 thienyl)perfluorocyclopentene (PT), using nonadiabatic molecular dynamics (NAMD) simulations. Employing the fewest switches surface hopping protocol, along with the semi-empirical multireference ODM2/MRCI-SD method, we investigate possible reaction pathways for this photoprocess, as well as their timescales and resulting photoproducts. Furthermore, using a dynamic configuration-space sampling procedure, we elucidate the role of triplet states in the photocyclization of PT, supporting available experimental data for the closely related DMPT molecule, which indicate an ultrafast intersystem crossing (ISC) transition competing with the singlet-driven photoswitching reaction. Our findings not only corroborate experimental studies on DAE switches, but also provide new mechanistic insights into the potential use in the rational design of DAE switches tailored for specific technological applications.

The Role of Aqueous Solvation on the Intersystem Crossing of Nitrophenols

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

The photochemistry of Rydberg-excited cyclobutanone: Photoinduced processes and ground state dynamics

Owing to ring strain, cyclic ketones exhibit complex excited state dynamics with multiple competing photochemical channels active on the ultrafast timescale. While the excited state dynamics of cyclobutanone after π* ← n excitation into the lowest-energy excited singlet (S1) state has been extensively studied, the dynamics following 3s ← n excitation into the higher-lying singlet Rydberg (S2) state are less well understood. Herein, we employ fully quantum multiconfigurational time-dependent Hartree (MCTDH) simulations using a model Hamiltonian as well as "on-the-fly" trajectory-based surface-hopping dynamics (TSHD) simulations to study the relaxation dynamics of cyclobutanone following 3s ← n excitation and to predict the ultrafast electron diffraction scattering signature of these relaxation dynamics. Our MCTDH and TSHD simulations indicate that relaxation from the initially-populated singlet Rydberg (S2) state occurs on the timescale of a few hundreds of femtoseconds to a picosecond, consistent with the symmetry-forbidden nature of the state-to-state transition involved. There is no obvious involvement of excited triplet states within the timeframe of our simulations (<2 ps). After non-radiative relaxation to the electronic ground state (S0), vibrationally hot cyclobutanone has sufficient internal energy to form multiple fragmented products including C2H4 + CH2CO (C2; 20%) and C3H6 + CO (C3; 2.5%). We discuss the limitations of our MCTDH and TSHD simulations, how these may influence the excited state dynamics we observe, and-ultimately-the predictive power of the simulated experimental observable.

Isomerization of surface functionalized SWCNTs and the critical influence on photoluminescence: static calculations and excited-state dynamics simulations

Single-walled carbon nanotubes (SWCNTs) functionalized with sparse surface chemical groups are promising for a variety of optical applications such as quantum information and bio-imaging. However, the luminescence efficiencies and stability, two key aspects, undoubtedly govern their practical usage. Herein, we assess the surface migration of oxygen and triazine groups on as-modified SWCNT fragments by adopting transition state theory and explore the de-excitation of oxygen-functionalized SWCNT fragments by performing non-adiabatic excited-state dynamics simulations. According to the predicted moderate or even small reaction barriers, the migration of both oxygen and triazine groups is feasible from an sp3 defect configuration forming an energetically more stable sp2 configuration at moderate or even room temperatures. Such isomerization leads to drastically different light emission capabilities as indicated by the large or zero oscillator strengths. During the dynamics simulations, the lowest excited singlet (S1) state rapidly decays in energy within 20 fs and then fluctuates until the end, providing insights into the emission mechanism of SWCNTs. This study highlights the potential intrinsic limitations of surface-functionalized SWCNTs for luminescence applications.

Massively scalable workflows for quantum chemistry: BigChem and ChemCloud

Electronic structure theory, i.e., quantum chemistry, is the fundamental building block for many problems in computational chemistry. We present a new distributed computing framework (BigChem), which allows for an efficient solution of many quantum chemistry problems in parallel. BigChem is designed to be easily composable and leverages industry-standard middleware (e.g., Celery, RabbitMQ, and Redis) for distributed approaches to large scale problems. BigChem can harness any collection of worker nodes, including ones on cloud providers (such as AWS or Azure), local clusters, or supercomputer centers (and any mixture of these). BigChem builds upon MolSSI packages, such as QCEngine to standardize the operation of numerous computational chemistry programs, demonstrated here with Psi4, xtb, geomeTRIC, and TeraChem. BigChem delivers full utilization of compute resources at scale, offers a programable canvas for designing sophisticated quantum chemistry workflows, and is fault tolerant to node failures and network disruptions. We demonstrate linear scalability of BigChem running computational chemistry workloads on up to 125 GPUs. Finally, we present ChemCloud, a web API to BigChem and successor to TeraChem Cloud. ChemCloud delivers scalable and secure access to BigChem over the Internet.

Can classical mechanics sense conical intersection?

Conical intersection (CI) leads to fast electronic energy transfer. However, Hamm and Stock [Phys. Rev. Lett. 109, 173201 (2012)] showed the existence of a vibrational CI and its role in vibrational energy relaxation. In this paper, we further investigate the vibrational energy relaxation using an isolated model Hamiltonian system of four vibrational modes with two distinctively different timescales (two fast modes and two slow modes). We show that the excitation of the slow modes plays a crucial role in the energy relaxation mechanism. We also analyze the system from a mixed quantum-classical (surface hopping method) and a completely classical point of view. Notably, surface hopping and even classical simulations also capture fast energy relaxation, which is a signature of CI's existence.

Nonadiabatic dynamics with classical trajectories: The problem of an initial coherent superposition of electronic states

Advances in coherent light sources and development of pump-probe techniques in recent decades have opened the way to study electronic motion in its natural time scale. When an ultrashort laser pulse interacts with a molecular target, a coherent superposition of electronic states is created and the triggered electron dynamics is coupled to the nuclear motion. A natural and computationally efficient choice to simulate this correlated dynamics is a trajectory-based method where the quantum-mechanical electronic evolution is coupled to a classical-like nuclear dynamics. These methods must approximate the initial correlated electron-nuclear state by associating an initial electronic wavefunction to each classical trajectory in the ensemble. Different possibilities exist that reproduce the initial populations of the exact molecular wavefunction when represented in a basis. We show that different choices yield different dynamics and explore the effect of this choice in Ehrenfest, surface hopping, and exact-factorization-based coupled-trajectory schemes in a one-dimensional two-electronic-state model system that can be solved numerically exactly. This work aims to clarify the problems that standard trajectory-based techniques might have when a coherent superposition of electronic states is created to initialize the dynamics, to discuss what properties and observables are affected by different choices of electronic initial conditions and to point out the importance of quantum-momentum-induced electronic transitions in coupled-trajectory schemes.

Photo-oxidation of methanol in complexes with pyrido[2,3-b]pyrazine: a nonadiabatic molecular dynamics study

Excited-state Proton-Coupled Electron Transfer (PCET) constitutes a key step in the photo-oxidation of small, electron-rich systems possessing acidic hydrogen atoms, such as water or alcohols, which can play a vital role in green hydrogen production. In this contribution, we employ ab initio quantum-chemical methods and on-the-fly nonadiabatic molecular dynamics simulations to study the mechanism and the photodynamics of PCET in 1 : 1 pyrido[2,3-b]pyrazine complexes with methanol. We find the process to be ultrafast and efficient when the intramolecular hydrogen bond is formed with one of the β-positioned nitrogen atoms. The complex exhibiting a hydrogen bond with an isolated nitrogen site, on the contrary, shows much lower reactivity. We explain this effect with the stabilization of the reactive ππ* charge-transfer electronic state in the former case.

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.

Capturing the generation and structural transformations of molecular ions

Molecular ions are ubiquitous and play pivotal roles1-3 in many reactions, particularly in the context of atmospheric and interstellar chemistry4-6. However, their structures and conformational transitions7,8, particularly in the gas phase, are less explored than those of neutral molecules owing to experimental difficulties. A case in point is the halonium ions9-11, whose highly reactive nature and ring strain make them short-lived intermediates that are readily attacked even by weak nucleophiles and thus challenging to isolate or capture before they undergo further reaction. Here we show that mega-electronvolt ultrafast electron diffraction (MeV-UED)12-14, used in conjunction with resonance-enhanced multiphoton ionization, can monitor the formation of 1,3-dibromopropane (DBP) cations and their subsequent structural dynamics forming a halonium ion. We find that the DBP+ cation remains for a substantial duration of 3.6 ps in aptly named 'dark states' that are structurally indistinguishable from the DBP electronic ground state. The structural data, supported by surface-hopping simulations15 and ab initio calculations16, reveal that the cation subsequently decays to iso-DBP+, an unusual intermediate with a four-membered ring containing a loosely bound17,18 bromine atom, and eventually loses the bromine atom and forms a bromonium ion with a three-membered-ring structure19. We anticipate that the approach used here can also be applied to examine the structural dynamics of other molecular ions and thereby deepen our understanding of ion chemistry.

A first principles derivation of energy-conserving momentum jumps in surface hopping simulations

The fewest switches surface hopping (FSSH) method proposed by Tully in 1990 [Tully, J. Chem. Phys. 93, 1061 (1990)]-along with its many later variations-forms the basis for most practical simulations of molecular dynamics with electronic transitions in realistic systems. Despite its popularity, a rigorous formal derivation of the algorithm has yet to be achieved. In this paper, we derive the energy-conserving momentum jumps employed by FSSH from the perspective of quantum trajectory surface hopping (QTSH) [Martens, J. Phys. Chem. A 123, 1110 (2019)]. In the limit of localized nonadiabatic transitions, simple mathematical and physical arguments allow the FSSH algorithm to be derived from first principles. For general processes, the quantum forces characterizing the QTSH method provide accurate results for nonadiabatic dynamics with rigorous energy conservation, at the ensemble level, within the consistency of the underlying stochastic surface hopping without resorting to the artificial momentum rescaling of FSSH.

Twisting Aromaticity and Photoinduced Dynamics in Hexapole Helicenes

Curved aromatic molecules are attractive electronic materials, where an additional internal strain uniquely modifies their structure, aromaticity, dynamics, and optical properties. Helicenes are examples of such twisted conjugated systems. Herein, we analyze the photoinduced dynamics in different stereoisomers of a hexapole helicene by using nonadiabatic excited-state molecular dynamics simulations. We explore how changes in symmetry and structural distortion modulate the intramolecular energy redistribution. We find that distinct helical assembly leads to different rigid distorted structures that in turn impact the nonradiative energy relaxation and ultimately formation of the self-trapped exciton. Subsequently, the value of the twisting angles relative to the central triphenylene core structure controls the global molecular aromaticity and electronic localization during the internal conversion process. Our work sheds light on how the future synthesis of novel curved aromatic compounds can be directed to attain specific desired electronic properties through the modulation of their twisted aromaticity.

NEXMD v2.0 Software Package for Nonadiabatic Excited State Molecular Dynamics Simulations

We present NEXMD version 2.0, the second release of the NEXMD (Nonadiabatic EXcited-state Molecular Dynamics) software package. Across a variety of new features, NEXMD v2.0 incorporates new implementations of two hybrid quantum-classical dynamics methods, namely, Ehrenfest dynamics (EHR) and the Ab-Initio Multiple Cloning sampling technique for Multiconfigurational Ehrenfest quantum dynamics (MCE-AIMC or simply AIMC), which are alternative options to the previously implemented trajectory surface hopping (TSH) method. To illustrate these methodologies, we outline a direct comparison of these three hybrid quantum-classical dynamics methods as implemented in the same NEXMD framework, discussing their weaknesses and strengths, using the modeled photodynamics of a polyphenylene ethylene dendrimer building block as a representative example. We also describe the expanded normal-mode analysis and constraints for both the ground and excited states, newly implemented in the NEXMD v2.0 framework, which allow for a deeper analysis of the main vibrational motions involved in vibronic dynamics. Overall, NEXMD v2.0 expands the range of applications of NEXMD to a larger variety of multichromophore organic molecules and photophysical processes involving quantum coherences and persistent couplings between electronic excited states and nuclear velocity.

Modeling the Ground- and Excited-State Unimolecular Decay of the Simplest Fluorinated Criegee Intermediate, HFCOO, Formed from the Ozonolysis of Hydrofluoroolefin Refrigerants

Hydrofluoroolefins (HFO) are fourth-generation refrigerants designed to function as efficient refrigerants with no ozone depletion potential and zero global warming potential. Despite extensive studies on their chemical and physical properties, the ground- and excited-state chemistry of their atmospheric oxidation products is less well understood. This study focuses on the ground- and excited-state chemistry of the simplest fluorinated Criegee intermediate (CI), fluoroformaldehyde oxide (HFCOO), which is the simplest fluorinated CI formed from the ozonolysis of HFOs. HFCOO contains syn- and anti-conformers, which have Boltzmann populations of, respectively, 87 and 13% at 298 K. For both conformers, the calculated ground-state reaction energy profiles associated with cyclization to form fluorodioxirane is lower than the equivalent unimolecular decay path in the simplest CI, H2COO, with anti-HFCOO returning a barrier height more than half of that of H2COO. The excited-state dynamics reveal that photoexcitation to the bright S2 state of syn-HFCOO and anti-HFCOO is expected to undergo a prompt O-O fission─with the former conformer expected to dissociate with an almost unity quantum yield and to form both O (1D) + HFCO (S0) and O (3P) + HFCO (T1) products. In contrast, photoexcitation of anti-HFCOO is expected to undergo an O-O bond fission with a non-unity quantum yield. The fraction of photoexcited anti-HFCOO that dissociates is predicted to exclusively form O (1D) + HFCO (S0) products, which is in sharp contrast to H2COO. The wider implications of our results are discussed from both physical and atmospheric chemistry perspectives.

Influence of Mode-Specific Excitation on the Nonadiabatic Dynamics of Methyl Nitrate (CH3ONO2)

The impact of mode-specific vibrational excitations on initial-preparation conditions was studied by examining the excited-state population decay rates in the nonadiabatic dynamics of methyl nitrate (CH₃ONO₂). In particular, exciting a few specific modes by adding a single quantum of energy clearly decelerated the nonadiabatic dynamics population decay rates. The underlying reason for this slower population decay was explained by analyzing the profiles of the excited-state potential energy surfaces in the Franck-Condon regions and the topology of the S₁/S₀ conical intersection. This study not only provides physical insights into the key mechanisms controlling nonadiabatic dynamics but also shows the possibility of controlling nonadiabatic dynamics via mode-specific vibrational excitations.

Role of Ultrafast Internal Conversion and Intersystem Crossing in the Nonadiabatic Relaxation Dynamics of ortho-Nitrobenzaldehyde

ortho-Nitrobenzaldehyde (oNBA) is a well-known photoactivated acid and a prototypical photolabile nitro-aromatic compound. Despite extensive investigations, the ultrafast relaxation dynamics of oNBA is still not properly understood, especially concerning the role of the triplet states. In this work, we provide an in-depth picture of this dynamics by combining single- and multireference electronic structure methods with potential energy surface exploration and nonadiabatic dynamics simulations using the Surface Hopping including ARbitary Couplings (SHARC) approach. Our results reveal that the initial decay from the bright ππ* state to the S₁ minimum is barrierless. It involves three changes in electronic structure from ππ* (ring) to nπ* (nitro group), to nπ* (aldehyde group), and then to another nπ* (nitro group). The decay of the ππ* takes 60-80 fs and can be tracked with time-resolved luminescence spectroscopy, where we predict for the first time a short-lived coherence of the luminescence energy with a 25 fs period. Intersystem crossing can occur already during the S₄ → S₁ deactivation cascade but also from S₁, with a time constant of about 2.4 ps and such that first a triplet ππ* state localized on the nitro group is populated. The triplet population first evolves into an nπ* and then quickly undergoes hydrogen transfer to form a biradical intermediate, from where the ketene is eventually produced. The majority of the excited population decays from S₁ through two conical intersections of equal utilization, a previously unreported one involving a scissoring motion of the nitro group that leads back to the oNBA ground state and the one involving hydrogen transfer that leads to the ketene intermediate.

Effects of solvent polarity on the novel excited-state intramolecular thiol proton transfer and photophysical property compared with the oxygen proton transfer

Recently, the dual-fluorescent phenomena of excited state intramolecular thiol proton transfer (ESIPT) for 3-thiolflavone derivative (3NTF) were reported by Chou and coworkers for the first time [J. Am. Chem. Soc. 143 (2021) 12715-12724], which opened a new chapter in the field of ESIPT. Based on density functional theory (DFT) and time-dependent density functional theory (TDDFT), the proton transfer processes of 3NTF in toluene, dichloromethane and acetonitrile were studied. By optimizing the structure of the ground (S0) state and first excited (S1) state of 3NTF in different solvents, the hydrogen-bond parameters and proton-transfer potential energy curves were calculated. It was shown that although photo-excitation enhanced the intramolecular hydrogen bonding strength and thus promoted the occurrence of ESIPT, the solvent polarities inhibited the enhancement of the hydrogen bond of S1 state, which was not conducive to ESIPT. The electron spectra analyses were consistent with experimental data, which confirmed the rationality of molecular configurations. The time-evolved excited state dynamics simulation was performed based on the optimized structure of 3NTF, indicating that the ESIPT was an ultrafast photochemical reaction less than 180 fs. Moreover, we compared the potential energy surfaces of ESIPT, electronic structures based on natural transition orbitals (NTOs) method and electron-hole isosurfaces for the 3NTF and the traditional flavone molecule (3NHF), concluded that the unusually large Stokes shift fluorescence of 3NTF was mainly caused by the coupling of ESIPT and twisting intramolecular charge transfer (TICT), and the 3NTF isomer had the more nπ* character in the electron transition process. The nπ* ICT significantly increased with the decrease of solvent polarities, affecting the molecular photophysical properties, this made it more widely used in biomedical, photochemical, materials science and other fields.

Vibrational Funnels for Energy Transfer in Organic Chromophores

Photoinduced intramolecular energy transfers in multichromophoric molecules involve nonadiabatic vibronic channels that act as energy transfer funnels. They commonly take place through specific directions of motion dictated by the nonadiabatic coupling vectors. Vibrational funnels may support persistent coherences between electronic states and sometimes delineate the presence of minor alternative energy transfer pathways. The ultimate confirmation of their role on the interchromophoric energy transfer can be achieved by performing nonadiabatic excited-state molecular dynamics simulations by selectively freezing the nuclear motions in question. Our results point out this strategy as a useful tool to identify and evaluate the impact of these vibrational funnels on the energy transfer processes and guide the in silico design of materials with tunable properties and enhanced functionalities. Our work encourages applications of this methodology to different chemical and biochemical processes such as reactive scattering and protein conformational changes, to name a few.

Ultrafast Electronic Relaxation in Aqueous [Fe(bpy)3]2+: A Surface Hopping Study

Trajectory surface hopping simulations are performed to better understand the electronic relaxation dynamics of [Fe(bpy)₃]²⁺ in aqueous solution. Specifically, the ultrafast relaxation from the photoexcited singlet metal-to-ligand charge-transfer (MLCT) to the metastable quintet metal-centered (MC) states is simulated through the surface hopping method, where the MLCT and MC states of [Fe(bpy)₃]²⁺ in aqueous solution are computed by using a model electronic Hamiltonian developed previously. As a result, most of the trajectories are interpreted to show the sequential relaxation pathways via the triplet MC states, though some are the direct pathway from MLCT to the quintet MC states. Even though the triplet MC states are involved in the relaxation, the population transfer to the singlet MC ground state is very small, and the population of the quintet MC states reaches more than ∼96%, reasonably consistent with the unity quantum efficiency discussed experimentally.

Impact of the core on the inter-branch exciton exchange in dendrimers

Organic dendrimers with π conjugated systems are capable of capturing solar energy as a renewable source for human use. Nonetheless, further study regarding the relationship between the structure and the energy transfer mechanism in these types of molecules is still necessary. In this work, nonadiabatic excited state molecular dynamics (NEXMD) were carried out to study the intra- and inter-branch exciton migration in two tetra-branched dendrimers, C(dSSB)₄ and Ad(BuSSB)₄, which differ in their respective carbon and adamantane core. Both systems undergo a ladder decay mechanism between excited states, with back-and-forth transitions between S₁ and S₂. Despite presenting very similar absorption-emission spectra, differences in the photoinduced energy relaxation are observed. The size of the core impacts the inter-branch energy exchange and transient exciton localization/delocalization, which ultimately condition the relative energy relaxation rates, being faster in Ad(BuSSB)₄ with respect to C(dSSB)₄. Nevertheless, the photoinduced processes lead to a progressive final exciton-self-trapping in one of the branches of both dendrimers, which is a desirable feature in organic photovoltaic applications. Our results can inspire the design of more efficient dendrimers with the desired magnitude of inter-branch exciton exchange and localization/delocalization according to changes in their core.

Simultaneously improving the efficiencies of photo- and thermal isomerization of an oxindole-based light-driven molecular rotary motor by a structural redesign

We designed a novel highly efficient light-driven molecular rotary motor theoretically by using electronic structure calculations and nonadiabatic dynamics simulations, and it showed excellent performance for both photo- and thermal isomerization processes simultaneously. By the small structural modification based on 3-(2,7-dimethyl-2,3-dihydro-1H-inden-1-ylidene)-1-methylindolin-2-one (DDIYM) synthesized by Feringa et al. recently, an oxindole-based light-driven molecular rotary motor, 3-(1,5-dimethyl-4,5-dihydrocyclopenta[b]pyrrol-6(1H)-ylidene)-1-methylindolin-2-one (DDPYM), is proposed, which displays a significant electronic push-pull character and weak steric hindrance for double-bond isomerization. The newly designed motor DDPYM shows a remarkable improvement of the quantum yield for both EP → ZM and ZP → EM photoisomerization processes, compared to the original motor DDIYM. Furthermore, the rotary motion in photoisomerization processes of DDPYM behaves more like a pure axial rotational motion approximately, while that of DDIYM is an obvious precessional motion. The weakness of the steric hindrance reduces the energy barriers of the thermal helix EM → EP and ZM → ZP inversion steps, and would accelerate two ground-state isomerization steps significantly. Our results confirm the feasibility of simultaneously improving the efficiencies of photo- and thermal isomerization of oxindole-based light-driven molecular rotary motors and this design idea sheds light on the future development of more efficient molecular motors.

UV photodissociation dynamics of the acetone oxide Criegee intermediate: experiment and theory

The photodissociation dynamics of the dimethyl-substituted acetone oxide Criegee intermediate [(CH₃)₂COO] is characterized following electronic excitation to the bright ¹ππ* state, which leads to O (¹D) + acetone [(CH₃)₂CO, S₀] products. The UV action spectrum of (CH₃)₂COO recorded with O (¹D) detection under jet-cooled conditions is broad, unstructured, and essentially unchanged from the corresponding electronic absorption spectrum obtained using a UV-induced depletion method. This indicates that UV excitation of (CH₃)₂COO leads predominantly to the O (¹D) product channel. A higher energy O (³P) + (CH₃)₂CO (T₁) product channel is not observed, although it is energetically accessible. In addition, complementary MS-CASPT2 trajectory surface-hopping (TSH) simulations indicate minimal population leading to the O (³P) channel and non-unity overall probability for dissociation (within 100 fs). Velocity map imaging of the O (¹D) products is utilized to reveal the total kinetic energy release (TKER) distribution upon photodissociation of (CH₃)₂COO at various UV excitation energies. Simulation of the TKER distributions is performed using a hybrid model that combines an impulsive model with a statistical component, the latter reflecting the longer-lived (>100 fs) trajectories identified in the TSH calculations. The impulsive model accounts for vibrational activation of (CH₃)₂CO arising from geometrical changes between the Criegee intermediate and the carbonyl product, indicating the importance of CO stretch, CCO bend, and CC stretch along with activation of hindered rotation and rock of the methyl groups in the (CH₃)₂CO product. Detailed comparison is also made with the TKER distribution arising from photodissociation dynamics of CH₂OO upon UV excitation.

The vertical excitation energies and a lifetime of the two lowest singlet excited states of the conjugated polyenes from C2 to C22: Ab initio, DFT, and semiclassical MNDO-MD simulations

Electronic excited states in the series of polyene molecules were explored. Optimal ground-state geometry was used for the evaluation of vertical excitation energies. Results of a chosen set of functionals were compared to post-HF methods (EOM-CCSD, NEVPT2, CASPT2, and MRCI). In addition, the semiempirical OM2/MNDO method using MRCISD computational level was confronted with the above-mentioned techniques. Despite the fact that the first excited state has a significant double-excitation character some functionals were able to qualitatively determine the correct state order (where the lowest excited state has a A g - character). The most successful functionals in transition energies predictions were PBE, TPSS and BLYP in Tamm-Dancoff approach (TDA), which had the smallest root-mean-square deviation (RMSD) scoring towards the experimental values. Regarding RMSD scoring, the OM2/MNDO method performed fairly well, too. Besides absorption spectra, lifetimes of the first two excited states were estimated based on a stochastic approach exploring a swarm of OM2/MNDO hopping dynamics using the Tully fewest switch algorithm for each molecule. The longest lifetime of the first excited state (S₁ ) was found for decapentaene (about 5 ps). Further elongation of the conjugated chain caused a mild decrease of this value to ca 1.5 ps for docosaundecaene.

Photodynamics of the Molecular Ruby [Cr(ddpd)2]3

The introduction of strong-field ligands can enable luminescence in first-row transition-metal complexes. In this way, earth-abundant near-infrared emitters can be obtained using early 3d metals. A prime example is the molecular ruby [Cr(ddpd)₂]³⁺ (ddpd = N,N'-dimethyl-N,N'-dipyridin-2-ylpyridine-2,6-diamine) that can achieve high phosphorescence quantum yields at room temperature in aqueous solution. To understand these remarkable properties, here, we simulate its photodynamics in water using trajectory surface hopping on linear vibronic coupling potentials parametrized from multiconfigurational CASSCF/CASPT2 calculations. We find that after excitation to the second absorption band, a relaxation cascade through metal-centered states occurs. After an initial back-and-forth intersystem crossing with higher-lying doublet states, the complex relaxes through a manifold of quartet metal-centered states to the low-lying doublet metal-centered states which are responsible for the experimentally observed emission. These electronic processes are driven by an elongation of the Cr-ligand bond lengths as well as the twisting motion of the trans-coordinated pyridine units in the ddpd ligands. The low-lying doublet states are reached within 1-2 ps and are close in geometry to the doublet minima, thus explaining the high phosphorescence quantum yield of the molecular ruby [Cr(ddpd)₂]³⁺.

Investigating the factors that influence sacrificial hydrogen evolution activity for three structurally-related molecular photocatalysts: thermodynamic driving force, excited-state dynamics, and surface interaction with cocatalysts

The design of molecular organic photocatalysts for reactions such as water splitting requires consideration of factors that go beyond electronic band gap and thermodynamic driving forces. Here, we carried out a theoretical investigation of three molecular photocatalysts (1-3) that are structurally similar but that show different hydrogen evolution activities (25, 23 & 0 μmol h^-1 for 1-3, respectively). We used density functional theory (DFT) and time-dependent DFT calculations to evaluate the molecules' optoelectronic properties, such as ionization potential, electron affinity, and exciton potentials, as well as the interaction between the molecular photocatalysts and an idealized platinum cocatalyst surface. The 'static' picture thus obtained was augmented by probing the nonadiabatic dynamics of the molecules beyond the Born-Oppenheimer approximation, revealing a different picture of exciton recombination and relaxation for molecule 3. Our results suggest that slow exciton recombination, fast relaxation to the lowest-energy excited state, and a shorter charge transfer distance between the photocatalyst and the metal cocatalyst are important features that contribute to the photocatalytic hydrogen evolution activity of 1 and 2, and may partly rationalize the observed inactivity of 3, in addition to its lower light absorption profile.

Unraveling the Mechanistic Pathway for the Dual Fluorescence in Green Fluorescent Protein (GFP) Chromophore Analogue: A Detailed Theoretical Investigation

The photophysical properties of the para-sulfonamide (p-TsABDI) analogue of the green fluorescent protein (GFP) chromophore with both proton donating and accepting sites have been exploited in polar solvents to understand the origin of the unusual dual fluorescence nature of the chromophore. In the polar solvents, the compound undergoes structural rearrangement upon photoexcitation, leading to the ultrafast excited-state intermolecular proton transfer (ESIPT) phenomenon at the S₁ surface. In this work, we employed both the static electronic structure calculations and on-the-fly molecular dynamics simulation to unravel the underlying reason for this peculiar behavior of the p-TsABDI analogue in polar solvents. To represent this adequately and provide extensive information on the ESIPT mechanism mediated by the solvent molecules, we considered explicit solvent molecules using the integral equation formalism variant of polarizable continuum (IEFPCM) model. From the static calculation analysis, we can conclude that the dual emissive behavior of the compound is ascribed to the proton transfer (PT) phenomena in the excited-state. However, based on the static calculation exclusively, it is hard to ascertain the mechanistic pathway of the PT phenomena. Hence, to investigate the dynamics and reaction mechanism for the ESIPT process, we performed the on-the-fly dynamics simulation for p-TsABDI in solvent clusters. Dynamics simulation results reveal that, based on the time lag between all the proton transfer processes, the ESIPT mechanism occurs in a stepwise manner from the benzylidene moiety of the chromophore to its imidazolinone moiety. However, the nonexistence of crossings between the S₁- and S₀-states confirms the PT characteristics of the reactions.

Newton-X Platform: New Software Developments for Surface Hopping and Nuclear Ensembles

Newton-X is an open-source computational platform to perform nonadiabatic molecular dynamics based on surface hopping and spectrum simulations using the nuclear ensemble approach. Both are among the most common methodologies in computational chemistry for photophysical and photochemical investigations. This paper describes the main features of these methods and how they are implemented in Newton-X. It emphasizes the newest developments, including zero-point-energy leakage correction, dynamics on complex-valued potential energy surfaces, dynamics induced by incoherent light, dynamics based on machine-learning potentials, exciton dynamics of multiple chromophores, and supervised and unsupervised machine learning techniques. Newton-X is interfaced with several third-party quantum-chemistry programs, spanning a broad spectrum of electronic structure methods.

Exciton-vibrational dynamics induces efficient self-trapping in a substituted nanoring

Cycloparaphenylenes, being the smallest segments of carbon nanotubes, have emerged as prototypes of the simplest carbon nanohoops. Their unique structure-dynamics-optical properties relationships have motivated a wide variety of synthesis of new related nanohoop species. Studies of how chemical changes, introduced in these new materials, lead to systems with new structural, dynamics and optical properties, expand their functionalities for optoelectronics applications. Herein, we study the effect that conjugation extension of a cycloparaphenylene through the introduction of a satellite tetraphenyl substitution has on its structural and dynamical properties. Our non-adiabatic excited state molecular dynamics simulations suggest that this substitution accelerates the electronic relaxation from the high-energy band to the lowest excited state. This is partially due to efficient conjugation achieved between specific phenyl units as introduced by the tetraphenyl substitution. We observe a particular exciton redistribution during relaxation, in which the tetraphenyl substitution plays a significant role. As a result, an efficient inter-band energy transfer takes place. Besides, the observed phonon-exciton interplay induces a significant exciton self-trapping. Our results encourage and guide the future studies of new phenyl substitutions in carbon nanorings with desired optoelectronic properties.

Rapid Intersystem Crossing Induced by Ultrafast Excited-State Intramolecular Proton Transfer in 3-Mercaptopyran-4-one

Investigation into the photoinduced processes of 3-mercaptopyran-4-one is carried out using trajectory-based surface hopping simulations. Excitation into the near-degenerate higher singlet excited states reveals rapid internal conversion (IC) into S₁ on a sub-50 fs timescale. Excited-state intramolecular proton transfer (ESIPT) also takes place simultaneously with IC. We observe that following tautomerization, the molecule has multiple relaxation pathways. A channel exists for it to nonradiatively decay into the tautomer ground-state or undergo rapid intersystem crossing (ISC) into the close-lying higher triplet state, which ultimately decays into T₁. The simulations show that ISC is significantly enhanced after ESIPT, which is studied by tracking the changes in energy gaps and associated spin-orbit coupling elements.

Comparing the Excited State Dynamics of CH2OO, the Simplest Criegee Intermediate, Following Vertical versus Adiabatic Excitation

Ab initio molecular dynamics studies of CH₂OO molecules following excitation to the minimum-energy geometry of the strongly absorbing S₂ (¹ππ*) state reveal a much richer range of behaviors than just the prompt O-O bond fission, with unity quantum yield and retention of overall planarity, identified in previous vertical excitation studies from the ground (S₀) state. Trajectories propagated for 100 fs from the minimum-energy region of the S₂ state show a high surface hopping (nonadiabatic coupling) probability between the near-degenerate S₂ and S₁ (¹nπ*) states at geometries close to the S₂ minimum, which enables population transfer to the optically dark S₁ state. Greater than 80% of the excited population undergoes O-O bond fission on the S₂ or S₁ potential energy surfaces (PESs) within the analysis period, mostly from nonplanar geometries wherein the CH₂ moiety is twisted relative to the COO plane. Trajectory analysis also reveals recurrences in the O-O stretch coordinate, consistent with the resonance structure observed at the red end of the parent S₂-S₀ absorption spectrum, and a small propensity for out-of-plane motion after nonadiabatic coupling to the S₁ PES that enables access to a conical intersection between the S₁ and S₀ states and cyclization to dioxirane products.

Effect of Temperature on Photoisomerization Dynamics of a Newly Designed Two-Stroke Light-Driven Molecular Rotary Motor

The working mechanism of conventional light-driven molecular rotary motors, especially Feringa-type motors, contains two photoisomerization steps and two thermal helix inversion steps. Due to the existence of a thermal helix inversion step, both the ability to work at lower temperatures and the rotation speed are limited. In this work, a two-stroke light-driven molecular rotary motor, 2-(1,5-dimethyl-4,5-dihydrocyclopenta[b]pyrrol-6(1H)-ylidene)-1,2-dihydro-3H-pyrrol-3-one (DDPY), is proposed, which is capable of performing unidirectional and repetitive rotation by only two photoisomerization (EP→ZP and ZP→EP) steps. With trajectory surface-hopping simulation at the semi-empirical OM2/MRCI level, the EP→ZP and ZP→EP nonadiabatic dynamics of DDPY were systematically studied at different temperatures. Both EP→ZP and ZP→EP photoisomerizations are on an ultrafast timescale (ca. 200–300 fs). The decay mode of EP→ZP photoisomerization is approximately bi-exponential, while that of ZP→EP photoisomerization is found to be periodic. For EP and ZP isomers of DDPY, after the S0→S1 excitation, the dynamical processes of nonadiabatic decay are both followed by twisting about the central C=C double bond and the pyramidalization of the C atom at the stator-axle linkage. The effect of temperature on the nonadiabatic dynamics of EP→ZP and ZP→EP photoisomerizations of DDPY has been systematically investigated. The average lifetimes of the S1 excited state and quantum yields for both EP→ZP and ZP→EP photoisomerization are almost temperature-independent, while the corresponding unidirectionality of rotation is significantly increased (e.g., 74% for EP→ZP and 72% for ZP→EP at 300 K vs 100% for EP→ZP and 94% for ZP→EP at 50 K) with lowering the temperature.