Unsupervised Machine Learning in the Analysis of Nonadiabatic Molecular Dynamics Simulation

The all-atomic full-dimensional-level simulations of nonadiabatic molecular dynamics (NAMD) in large realistic systems has received high research interest in recent years. However, such NAMD simulations normally generate an enormous amount of time-dependent high-dimensional data, leading to a significant challenge in result analyses. Based on unsupervised machine learning (ML) methods, considerable efforts were devoted to developing novel and easy-to-use analysis tools for the identification of photoinduced reaction channels and the comprehensive understanding of complicated molecular motions in NAMD simulations. Here, we tried to survey recent advances in this field, particularly to focus on how to use unsupervised ML methods to analyze the trajectory-based NAMD simulation results. Our purpose is to offer a comprehensive discussion on several essential components of this analysis protocol, including the selection of ML methods, the construction of molecular descriptors, the establishment of analytical frameworks, their advantages and limitations, and persistent challenges.

Deciphering Photoinduced Catalytic Reaction Mechanisms in Natural and Artificial Photosynthetic Systems on Multiple Temporal and Spatial Scales Using X-ray Probes

Utilization of renewable energies for catalytically generating value-added chemicals is highly desirable in this era of rising energy demands and climate change impacts. Artificial photosynthetic systems or photocatalysts utilize light to convert abundant CO2, H2O, and O2 to fuels, such as carbohydrates and hydrogen, thus converting light energy to storable chemical resources. The emergence of intense X-ray pulses from synchrotrons, ultrafast X-ray pulses from X-ray free electron lasers, and table-top laser-driven sources over the past decades opens new frontiers in deciphering photoinduced catalytic reaction mechanisms on the multiple temporal and spatial scales. Operando X-ray spectroscopic methods offer a new set of electronic transitions in probing the oxidation states, coordinating geometry, and spin states of the metal catalytic center and photosensitizers with unprecedented energy and time resolution. Operando X-ray scattering methods enable previously elusive reaction steps to be characterized on different length scales and time scales. The methodological progress and their application examples collected in this review will offer a glimpse into the accomplishments and current state in deciphering reaction mechanisms for both natural and synthetic systems. Looking forward, there are still many challenges and opportunities at the frontier of catalytic research that will require further advancement of the characterization techniques.

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.

Vibrational Coherences of the Photoinduced Mixed-Valent Creutz-Taube Ion Revealed by Excited State Dynamics

A recent study of photoinduced mixed-valency in the one-electron reduced form (μ-pz)[RuII(NH3)5]24+ of the Creutz-Taube ion used transient absorption spectroscopy with vis-NIR broadband detection to uncover a mixed-valent excited state with a typical intervalence charge transfer band and a nanosecond lifetime [Pieslinger et al. Angew. Chem., Int. Ed. 2022, 61, e202211747]. Herein, we use excited state dynamics simulations with implicit solvation to elucidate the electronic and vibrational evolution in the first 10 ps after the optical excitation. A manifold of excited states with weak interaction between the metal centers is populated already at time zero due to the breakdown of the Condon approximation and dominates the population of electronic states at short time scales (<0.5 ps). A long-lived vibrational wave packet mostly confined to oscillations of the metal center-bridge distances is observed. The oscillations are traced to the electronic structure properties of states with weak metal-metal coupling. The long-lived mixed-valent excited state of the Creutz-Taube ion analogue is formed vibrationally cold and has a more compact geometry. While experimentally, intersystem crossing and vibrational relaxation were deduced to be completed within 1 ps, our analysis indicates that both processes might persist at longer times.

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.

An on-the-fly deep neural network for simulating time-resolved spectroscopy: predicting the ultrafast ring opening dynamics of 1,2-dithiane

Revolutionary developments in ultrafast light source technology are enabling experimental spectroscopists to probe the structural dynamics of molecules and materials on the femtosecond timescale. The capacity to investigate ultrafast processes afforded by these resources accordingly inspires theoreticians to carry out high-level simulations which facilitate the interpretation of the underlying dynamics probed during these ultrafast experiments. In this Article, we implement a deep neural network (DNN) to convert excited-state molecular dynamics simulations into time-resolved spectroscopic signals. Our DNN is trained on-the-fly from first-principles theoretical data obtained from a set of time-evolving molecular dynamics. The train-test process iterates for each time-step of the dynamics data until the network can predict spectra with sufficient accuracy to replace the computationally intensive quantum chemistry calculations required to produce them, at which point it simulates the time-resolved spectra for longer timescales. The potential of this approach is demonstrated by probing dynamics of the ring opening of 1,2-dithiane using sulphur K-edge X-ray absorption spectroscopy. The benefits of this strategy will be more markedly apparent for simulations of larger systems which will exhibit a more notable computational burden, making this approach applicable to the study of a diverse range of complex chemical dynamics.

Atomic-scale observation of solvent reorganization influencing photoinduced structural dynamics in a copper complex photosensitizer

Photochemical reactions in solution are governed by a complex interplay between transient intramolecular electronic and nuclear structural changes and accompanying solvent rearrangements. State-of-the-art time-resolved X-ray solution scattering has emerged in the last decade as a powerful technique to observe solute and solvent motions in real time. However, disentangling solute and solvent dynamics and how they mutually influence each other remains challenging. Here, we simultaneously measure femtosecond X-ray emission and scattering to track both the intramolecular and solvation structural dynamics following photoexcitation of a solvated copper photosensitizer. Quantitative analysis assisted by molecular dynamics simulations reveals a two-step ligand flattening strongly coupled to the solvent reorganization, which conventional optical methods could not discern. First, a ballistic flattening triggers coherent motions of surrounding acetonitrile molecules. In turn, the approach of acetonitrile molecules to the copper atom mediates the decay of intramolecular coherent vibrations and induces a further ligand flattening. These direct structural insights reveal that photoinduced solute and solvent motions can be intimately intertwined, explaining how the key initial steps of light harvesting are affected by the solvent on the atomic time and length scale. Ultimately, this work takes a step forward in understanding the microscopic mechanisms of the bidirectional influence between transient solvent reorganization and photoinduced solute structural dynamics.

Mind the GAP: quantifying the breakdown of the linear vibronic coupling Hamiltonian

Excited state dynamics play a critical role across a broad range of scientific fields. Importantly, the highly non-equilibrium nature of the states generated by photoexcitation means that excited state simulations should usually include an accurate description of the coupled electronic-nuclear motion, which often requires solving the time-dependent Schrödinger equation (TDSE). One of the biggest challenges for these simulations is the requirement to calculate the PES over which the nuclei evolve. An effective approach for addressing this challenge is to use the approximate linear vibronic coupling (LVC) Hamiltonian, which enables a model potential to be parameterised using relatively few quantum chemistry calculations. However, this approach is only valid provided there are no large amplitude motions in the excited state dynamics. In this paper we introduce and deploy a metric, the global anharmonicity parameter (GAP), which can be used to assess the accuracy of an LVC potential. Following its derivation, we illustrate its utility by applying it to three molecules exhibiting different rigidity in their excited states.

Trajectory Propagation of Symmetrical Quasi-classical Dynamics with Meyer-Miller Mapping Hamiltonian Using Machine Learning

The long short-term memory recurrent neural network (LSTM-RNN) approach is applied to realize the trajectory-based nonadiabatic dynamics within the framework of the symmetrical quasi-classical dynamics method based on the Meyer-Miller mapping Hamiltonian (MM-SQC). After construction, the LSTM-RNN model allows us to propagate the entire trajectory evolutions of all involved degrees of freedoms (DOFs) from initial conditions. The proposed idea is proven to be reliable and accurate in the simulations of the dynamics of several site-exciton electron-phonon coupling models and three Tully's scattering models. It indicates that the LSTM-RNN model perfectly captures the dynamical information on the trajectory evolution in the MM-SQC dynamics. Our work proposes a novel machine learning approach in the simulation of trajectory-based nonadiabatic dynamic of complex systems with a large number of DOFs.

The principal component analysis of the ring deformation in the nonadiabatic surface hopping dynamics

The analysis of the leading active molecular motions in the on-the-fly trajectory surface hopping simulation provides the essential information to understand the geometric evolution in nonadiabatic dynamics. When the ring deformation is involved, the identification of the key active coordinates becomes challenging. A "hierarchical" protocol based on the dimensionality reduction and clustering approaches is proposed for the automatic analysis of the ring deformation in the nonadiabatic molecular dynamics. The representative system keto isocytosine is taken as the prototype to illustrate this protocol. The results indicate that the current hierarchical analysis protocol is a powerful way to clearly clarify both the major and minor active molecular motions of the ring distortion in nonadiabatic dynamics.

Towards Optimized Photoluminescent Copper(I) Phenanthroline-Functionalized Complexes: Control of the Photophysics by Symmetry-Breaking and Spin–Orbit Coupling

The electronic and structural alterations induced by the functionalization of the 1,10-phenanthroline (phen) ligand in [Cu(I) (phen-R₂)₂]⁺ complexes (R=H, CH₃, tertio-butyl, alkyl-linkers) and their consequences on the luminescence properties and thermally activated delay fluorescence (TADF) activity are investigated using the density functional theory (DFT) and its time-dependent (TD) extension. It is shown that highly symmetric molecules with several potentially emissive nearly-degenerate conformers are not promising because of low S₁/S₀ oscillator strengths together with limited or no S₁/T₁ spin–orbit coupling (SOC). Furthermore, steric hindrance, which prevents the flattening of the complex upon irradiation, is a factor of instability. Alternatively, linking the phenanthroline ligands offers the possibility to block the flattening while maintaining remarkable photophysical properties. We propose here two promising complexes, with appropriate symmetry and enough rigidity to warrant stability in standard solvents. This original study paves the way for the supramolecular design of new emissive devices.

Trend in light-induced excited-state spin trapping in Fe(II)-based spin crossover systems

A computational study of the light-induced excited spin-state trapping (LIESST) in a number of Fe(II) spin crossover complexes, coordinated by monodentate, bidentate and multidentate ligands is carried out, with the goal to uncover the trend in the low temperature relaxation rate. A nine order of magnitude change in low temperature relaxation rate is observed among the complexes. The trend is rationalized in terms of the change in metal-ligand covalency, numerically estimated by the crystal orbital Hamiltonian population, thus influencing the back donation or delocalization of the electrons from the low-lying Fe(II)-centered molecular orbital to the empty low-lying ligand-centered π* antibonding molecular orbitals.

Toward Simulation of Fe(II) Low-Spin → High-Spin Photoswitching by Synergistic Spin-Vibronic Dynamics

A new theoretical approach is presented and applied for the simulation of Fe(II) low-spin (LS, singlet, t_2g⁶e_g⁰) → high-spin (HS, quintet, t_2g⁴e_g²) photoswitching dynamics of the octahedral model complex [Fe(NCH)₆]²⁺. The utilized synergistic methodology heavily exploits the strengths of complementary electronic structure and spin-vibronic dynamics methods. Specifically, we perform 3D quantum dynamics (QD) and full-dimensional trajectory surface hopping (TSH, in conjunction with a linear vibronic coupling model), with the modes for QD selected by TSH. We follow a hybrid approach which is based on the application of time-dependent density functional theory (TD-DFT) excited-state potential energy surfaces (PESs) and multiconfigurational second-order perturbation theory (CASPT2) spin–orbit couplings (SOCs). Our method delivers accurate singlet–triplet–quintet intersystem crossing (ISC) dynamics, as assessed by comparison to our recent high-level ab initio simulations and related time-resolved experimental data. Furthermore, we investigate the capability of our simulations to identify the location of ISCs. Finally, we assess the approximation of constant SOCs (calculated at the Franck–Condon geometry), whose validity has central importance for the combination of TD-DFT PESs and CASPT2 SOCs. This efficient methodology will have a key role in simulating LS → HS dynamics for more complicated cases, involving higher density of states and varying electronic character, as well as the analysis of ultrafast experiments.

Excited-state structure of copper phenanthroline-based photosensitizers

Cu diimine complexes present a noble metal free alternative to classical Ru, Re, Ir and Pt based photosensitizers in solution photochemistry, photoelectrochemical or dye-sensitized solar cells. Optimization of these dyes requires understanding of factors governing the key photochemical properties: excited state lifetime and emission quantum yield. The involvement of exciplex formation in the deactivation of the photoexcited state is a key question. We investigate the excited-state structure of [Cu(dmp)₂]⁺ and [Cu(dsbtmp)₂]⁺ (dmp = 2,9-dimethyl-1,10-phenanthroline, dsbtmp = 2,9-di-sec-butyl-3,4,7,8-tetramethyl-1,10-phenanthroline) using pump-probe X-ray absorption spectroscopy (XAS) and DFT. Features of XAS that distinguish flattened tetrahedral site and 5-coordinated geometry with an additional solvent near Cu(II) center are identified. Pump-probe XAS demonstrates that for both complexes the excited state is 4-coordinated. For [Cu(dmp)₂]⁺ the exciplex is 0.24 eV higher in energy than the flattened triplet state, therefore it can be involved in deactivation pathways as a non-observable state that forms slower than it decays. For [Cu(dsbtmp)₂]⁺ the excited-state structure is characterized by Cu-N distances of 1.98 and 2.07 Å and minor distortions, leading to a 3 orders of magnitude longer excited-state lifetime.

Surface Hopping Dynamics on Vibronic Coupling Models

The simulation of photoinduced non-adiabatic dynamics is of great relevance in many scientific disciplines, ranging from physics and materials science to chemistry and biology. Upon light irradiation, different relaxation processes take place in which electronic and nuclear motion are intimately coupled. These are best described by the time-dependent molecular Schrödinger equation, but its solution poses fundamental practical challenges to contemporary theoretical chemistry. Two widely used and complementary approaches to this problem are multiconfigurational time-dependent Hartree (MCTDH) and trajectory surface hopping (SH). MCTDH is an accurate fully quantum-mechanical technique but often is feasible only in reduced dimensionality, in combination with approximate vibronic coupling (VC) Hamiltonians, or both (i.e., reduced-dimensional VC potentials). In contrast, SH is a quantum–classical technique that neglects most nuclear quantum effects but allows nuclear dynamics in full dimensionality by calculating potential energy surfaces on the fly. If nuclear quantum effects do not play a central role and a linear VC (LVC) Hamiltonian is appropriate—e.g., for stiff molecules that generally keep their conformation in the excited state—then it seems advantageous to combine the efficient LVC and SH techniques. In this Account, we describe how surface hopping based on an LVC Hamiltonian (SH/LVC)—as recently implemented in the SHARC surface hopping package—can provide an economical and automated approach to simulate non-adiabatic dynamics. First, we illustrate the potential of SH/LVC in a number of showcases, including intersystem crossing in SO₂, intra-Rydberg dynamics in acetone, and several photophysical studies on large transition-metal complexes, which would be much more demanding or impossible to perform with other methods. While all of the applications provide very useful insights into light-induced phenomena, they also hint at difficulties faced by the SH/LVC methodology that need to be addressed in the future. Second, we contend that the SH/LVC approach can be useful to benchmark SH itself. By the use of the same (LVC) potentials as MCTDH calculations have employed for decades and by relying on the efficiency of SH/LVC, it is possible to directly compare multiple SH test calculations with a MCTDH reference and ponder the accuracy of various correction algorithms behind the SH methodology, such as decoherence corrections or momentum rescaling schemes. Third, we demonstrate how the efficiency of SH/LVC can also be exploited to identify essential nuclear and electronic degrees of freedom to be employed in more accurate MCTDH calculations. Lastly, we show that SH/LVC is able to advance the development of SH protocols that can describe nuclear dynamics including explicit laser fields—a very challenging endeavor for trajectory-based schemes. To end, this Account compiles the typical costs of contemporary SH simulations, evidencing the great advantages of using parametrized potentials. The LVC model is a sleeping beauty that, kissed by SH, is fueling the field of excited-state molecular dynamics. We hope that this Account will stimulate future research in this direction, leveraging the advantages of the SH/VC schemes to larger extents and extending their applicability to uncharted territories.

Photoinduced Low-Spin → High-Spin Mechanism of an Octahedral Fe(II) Complex Revealed by Synergistic Spin-Vibronic Dynamics

The Fe(II) low-spin (LS; ¹A_1g, t_2g⁶e_g⁰) → high-spin (HS; ⁵T_2g, t_2g⁴e_g²) light-induced excited spin state trapping (LIESST) mechanism solely involving metal-centered states is revealed by synergistic spin-vibronic dynamics simulations. For the octahedral [Fe(NCH)₆]²⁺ complex, we identify an initial ∼100 fs ¹T_1g → ³T_2g intersystem crossing, controlled by vibronic coupling to antisymmetric Fe-N stretching motion. Subsequently, population branching into ³T_1g, ⁵T_2g (HS), and ¹A_1g (LS) is observed on a subpicosecond time scale, with the dynamics dominated by coherent Fe-N breathing wavepackets. These findings are consistent with ultrafast experiments, methodologically establish a new state of the art, and will give a strong impetus for further intriguing dynamical studies on LS → HS photoswitching.

The Quest to Simulate Excited-State Dynamics of Transition Metal Complexes

This Perspective describes current computational efforts in the field of simulating photodynamics of transition metal complexes. We present the typical workflows and feature the strengths and limitations of the different contemporary approaches. From electronic structure methods suitable to describe transition metal complexes to approaches able to simulate their nuclear dynamics under the effect of light, we give particular attention to build a bridge between theory and experiment by critically discussing the different models commonly adopted in the interpretation of spectroscopic experiments and the simulation of particular observables. Thereby, we review all the studies of excited-state dynamics on transition metal complexes, both in gas phase and in solution from reduced to full dimensionality.

Capturing photochemical and photophysical transformations in iron complexes with ultrafast X-ray spectroscopy and scattering

Light-driven chemical transformations provide a compelling approach to understanding chemical reactivity with the potential to use this understanding to advance solar energy and catalysis applications. Capturing the non-equilibrium trajectories of electronic excited states with precision, particularly for transition metal complexes, would provide a foundation for advancing both of these objectives. Of particular importance for 3d metal compounds is characterizing the population dynamics of charge-transfer (CT) and metal-centered (MC) electronic excited states and understanding how the inner coordination sphere structural dynamics mediate the interaction between these states. Recent advances in ultrafast X-ray laser science has enabled the electronic excited state dynamics in 3d metal complexes to be followed with unprecedented detail. This review will focus on simultaneous X-ray emission spectroscopy (XES) and X-ray solution scattering (XSS) studies of iron coordination and organometallic complexes. These simultaneous XES-XSS studies have provided detailed insight into the mechanism of light-induced spin crossover in iron coordination compounds, the interaction of CT and MC excited states in iron carbene photosensitizers, and the mechanism of Fe-S bond dissociation in cytochrome c.

Analysis of bath motion in MM-SQC dynamics via dimensionality reduction approach: Principal component analysis

The system-plus-bath model is an important tool to understand the nonadiabatic dynamics of large molecular systems. Understanding the collective motion of a large number of bath modes is essential for revealing their key roles in the overall dynamics. Here, we applied principal component analysis (PCA) to investigate the bath motion in the basis of a large dataset generated from the symmetrical quasi-classical dynamics method based on the Meyer-Miller mapping Hamiltonian nonadiabatic dynamics for the excited-state energy transfer in the Frenkel-exciton model. The PCA method clearly elucidated that two types of bath modes, which either display strong vibronic coupling or have frequencies close to that of the electronic transition, are important to the nonadiabatic dynamics. These observations were fully consistent with the physical insights. The conclusions were based on the PCA of the trajectory data and did not involve significant pre-defined physical knowledge. The results show that the PCA approach, which is one of the simplest unsupervised machine learning dimensionality reduction methods, is a powerful one for analyzing complicated nonadiabatic dynamics in the condensed phase with many degrees of freedom.

Ultrafast processes: coordination chemistry and quantum theory

Coordination compounds, characterized by fascinating and tunable electronic properties, are capable of binding easily to proteins, polymers, wires and DNA. Upon irradiation, these molecular systems develop functions finding applications in solar cells, photocatalysis, luminescent and conformational probes, electron transfer triggers and diagnostic or therapeutic tools. The control of these functions is activated by the light wavelength, the metal/ligand cooperation and the environment within the first picoseconds (ps). After a brief summary of the theoretical background, this perspective reviews case studies, from 1st row to 3rd row transition metal complexes, that illustrate how spin-orbit, vibronic coupling and quantum effects drive the photophysics of this class of molecules at the early stage of the photoinduced elementary processes within the fs-ps time scale range.

Photoinduced Dynamics with Constrained Vibrational Motion: FrozeNM Algorithm

Ab initio molecular dynamics (AIMD) simulation, analyzed in terms of vibrational normal modes, is a widely used technique that facilitates understanding of complex structural motions and coupling between electronic and nuclear degrees of freedom. Usually, only a subset of vibrations is directly involved in the process of interest. The impact of these vibrations can be evaluated by performing AIMD simulations by selectively freezing certain motions. Herein, we present frozen normal mode (FrozeNM), a new algorithm to apply normal-mode constraints in AIMD simulations, as implemented in the nonadiabatic excited state molecular dynamics code. We further illustrate its capacity by analyzing the impact of normal-mode constraints on the photoinduced energy transfer between polyphenylene ethynylene dendrimer building blocks. Our results show that the electronic relaxation can be significantly slowed down by freezing a well-selected small subset of active normal modes characterized by their contributions in the direction of energy transfer. The application of these constraints reduces the nonadiabatic coupling between electronic excited states during the entire dynamical simulations. Furthermore, we validate reduced dimensionality models by freezing all the vibrations, except a few active modes. Altogether, we consider FrozeNM as a useful tool that can be broadly used to underpin the role of vibrational motion in a studied process and to formulate reduced models that describe essential physical phenomena.

On the interplay of solvent and conformational effects in simulated excited-state dynamics of a copper phenanthroline photosensitizer

QM/MM direct dynamics simulations in acetonitrile reveal the interplay between solvent and conformational effects in the photoinduced ultrafast flattening of a copper photosensitizer.

Vibronic coherence evolution in multidimensional ultrafast photochemical processes

The complex choreography of electronic, vibrational, and vibronic couplings used by photoexcited molecules to transfer energy efficiently is remarkable, but an unambiguous description of the temporally evolving vibronic states governing these processes has proven experimentally elusive. We use multidimensional electronic-vibrational spectroscopy to identify specific time-dependent excited state vibronic couplings involving multiple electronic states, high-frequency vibrations, and low-frequency vibrations which participate in ultrafast intersystem crossing and subsequent relaxation of a photoexcited transition metal complex. We discover an excited state vibronic mechanism driving long-lived charge separation consisting of an initial electronically-localized vibrational wavepacket which triggers delocalization onto two charge transfer states after propagating for ~600 femtoseconds. Electronic delocalization consequently occurs through nonadiabatic internal conversion driven by a 50 cm^-1 coupling resulting in vibronic coherence transfer lasting for ~1 picosecond. This study showcases the power of multidimensional electronic-vibrational spectroscopy to elucidate complex, non-equilibrium energy and charge transfer mechanisms involving multiple molecular coordinates.

Unconventional two-step spin relaxation dynamics of [Re(CO)3(im)(phen)]+ in aqueous solution

Changes of molecular spin are ubiquitous in chemistry and biology. Among spin flip processes, one of the fastest is intersystem crossing (ISC) in transition metal complexes. Here, we investigate the spin relaxation dynamics and emission spectrum of [Re(CO)3(im)(phen)]+ (im = imidazole, phen = phenanthroline) using extensive full-dimensional excited-state dynamics simulations in explicit aqueous solution. Contrary to what has been observed in other transition metal complexes, the transition from the singlet to triplet states occurs via a two-step process, with clearly separable electronic and nuclear-driven components with two different time scales. The initially excited electronic wave function is a "molecular spin-orbit wave packet" that evolves almost instantaneously, with an 8 fs time constant, into an approximate 25 : 75 singlet-to-triplet equilibrium. Surprisingly, this ISC process is an order of magnitude faster than it was previously documented for this and other rhenium(i) carbonyl diimine complexes from emission spectra. Simulations including explicit laser field interactions evidence that few-cycle UV laser pulses are required to follow the creation and evolution of such molecular spin-orbit wave packets. The analysis of the dynamics also reveals a retarded ISC component, with a time constant of 420 fs, which can be explained invoking intramolecular vibrational energy redistribution. The emission spectrum is shown to be characterized by ISC convoluted with internal conversion and vibrational relaxation. These results provide fundamental understanding of ultrafast intersystem crossing in transition metal complexes.

Tracking multiple components of a nuclear wavepacket in photoexcited Cu(I)-phenanthroline complex using ultrafast X-ray spectroscopy

Disentangling the strong interplay between electronic and nuclear degrees of freedom is essential to achieve a full understanding of excited state processes during ultrafast nonadiabatic chemical reactions. However, the complexity of multi-dimensional potential energy surfaces means that this remains challenging. The energy flow during vibrational and electronic relaxation processes can be explored with structural sensitivity by probing a nuclear wavepacket using femtosecond time-resolved X-ray Absorption Near Edge Structure (TR-XANES). However, it remains unknown to what level of detail vibrational motions are observable in this X-ray technique. Herein we track the wavepacket dynamics of a prototypical [Cu(2,9-dimethyl-1,10-phenanthroline)₂]⁺ complex using TR-XANES. We demonstrate that sensitivity to individual wavepacket components can be modulated by the probe energy and that the bond length change associated with molecular breathing mode can be tracked with a sub-Angstrom resolution beyond optical-domain observables. Importantly, our results reveal how state-of-the-art TR-XANES provides deeper insights of ultrafast nonadiabatic chemical reactions.

Performance of trajectory surface hopping method in the treatment of ultrafast intersystem crossing dynamics

We carried out extensive studies to examine the performance of the fewest-switches surface hopping method in the description of the ultrafast intersystem crossing dynamic of various singlet-triplet (S-T) models by comparison with the results of the exact full quantum dynamics. Different implementation details and some derivative approaches were examined. As expected, it is better to perform the trajectory surface hopping calculations in the spin-adiabatic representation or by the local diabatization approach, instead of in the spin-diabatic representation. The surface hopping method provides reasonable results for the short-time dynamics in the S-T model with weak spin-orbital coupling (diabatic coupling), although it does not perform well in the models with strong spin-orbital coupling (diabatic coupling). When the system accesses the S-T potential energy crossing with rather high kinetic energy, the trajectory surface hopping method tends to produce a good description of the nonadiabatic intersystem crossing dynamics. The impact of the decoherence correction on the performance of the trajectory surface hopping is system dependent. It improves the result accuracy in many cases, while its influence may also be minor for other cases.

Computer-Aided Design of Luminescent Linear N-Heterocyclic Carbene Copper(I) Pyridine Complexes

Multireference configuration interaction methods including spin-orbit interactions have been employed to investigate the photophysical properties of various linear N-heterocyclic carbene (NHC) copper(I) pyridine complexes with the aim of designing performant thermally activated delayed fluorescence (TADF) emitters for use in organic-light-emitting diodes. Our theoretical results indicate that this structural motif is very favorable for generating excited triplet states with high quantum yield. The first excited singlet (S_MLCT) and corresponding triplet state (T_MLCT) are characterized by d_σ → π_Py metal-to-ligand charge-transfer (MLCT) excitations. Efficient intersystem crossing (ISC) and reverse ISC (rISC) between these states is mediated by a near-degenerate second triplet state (T_MLCT/LC) with large d_π → π_Py contributions. Spin-vibronic coupling is strong and is expected to play a major role in the (r)ISC processes. The calculations reveal, however, that the luminescence is effectively quenched by locally excited triplet states if the NHC ligand carries two diisopropylphenyl (DIPP) substituents. When DIPP is replaced with 1-adamantyl residues, this quenching process is suppressed and TADF in the UV spectral regime is predicted to proceed at a rate of about 1/μs. The introduction of +I substituents on the carbene and -M substituents on the pyridine allows tuning of the emission wavelength from the UV to the blue-green or green spectral region.

Surface hopping dynamics including intersystem crossing using the algebraic diagrammatic construction method

We report an implementation for employing the algebraic diagrammatic construction to second order [ADC(2)] ab initio electronic structure level of theory in nonadiabatic dynamics simulations in the framework of the SHARC (surface hopping including arbitrary couplings) dynamics method. The implementation is intended to enable computationally efficient, reliable, and easy-to-use nonadiabatic dynamics simulations of intersystem crossing in organic molecules. The methodology is evaluated for the 2-thiouracil molecule. It is shown that ADC(2) yields reliable excited-state energies, wave functions, and spin-orbit coupling terms for this molecule. Dynamics simulations are compared to previously reported results using high-level multi-state complete active space perturbation theory, showing favorable agreement.

Excited State Properties of Heteroleptic Cu(I) 4H-Imidazolate Complexes

The excited state properties of three heteroleptic copper(I) xantphos 4H-imidazolate complexes are investigated by means of femtosecond and nanosecond time-resolved transient absorption spectroscopy in dichloromethane solution. The subpicosecond spectral changes observed after excitation into the MLCT absorption band are interpreted as intersystem crossing from the singlet to the triplet manifold. This interpretation is corroborated by DFT and TD-DFT results, indicating a comparable molecular geometry in the ground state (and hence the nonrelaxed singlet state) and the excited triplet state. Population of the triplet state is followed by planarization of the N-aryl rings of the 4H-imidazolate ligand on a 10 ps time scale. The planarization strongly depends on the substitution pattern of the N-aryls and correlates with the reduced moment of inertia for the planarization motion. The triplet state subsequently decays to the ground state in about 100 ns. These results demonstrate that the excited state processes of copper(I) complexes depend on the specific ligand(s) and their substitution pattern. Thus, the work presented points to a possibility to design copper(I) complexes with specific photophysical properties.

Trajectory Surface-Hopping Dynamics Including Intersystem Crossing in [Ru(bpy)3]2

Surface-hopping dynamics coupled to linear response TDDFT and explicit nonadiabatic and spin-orbit couplings have been used to model the ultrafast intersystem crossing (ISC) dynamics in [Ru(bpy)₃]²⁺. Simulations using an ensemble of trajectories starting from the singlet metal-to-ligand charge transfer (¹MLCT) band show that the manifold of ³MLCT triplet states is first populated from high-lying singlet states within 26 ± 3 fs. ISC competes with an intricate internal conversion relaxation process within the singlet manifold to the lowest singlet state. Normal-mode analysis and principal component analysis, combined with further dynamical simulations where the nuclei are frozen, unequivocally demonstrate that it is not only the high density of states and the large spin-orbit couplings of the system that promote ISC. Instead, geometrical relaxation involving the nitrogen atoms is required to allow for state mixing and efficient triplet population transfer.

Accelerating direct quantum dynamics using graphical processing units

The direct dynamics variational multi-configurational Gaussian (DD-vMCG) method is combined with electronic structure calculations accelerated by Graphical Processing Units (GPUs). This is used to identify GPU acceleration will have a significant effect for both ground and excited state simulations.