Ultrafast processes: coordination chemistry and quantum theory

Ultrafast Excited-State Nonadiabatic Dynamics in Pt(II) Donor-Bridge-Acceptor Assemblies: A Quantum Approach for Optical Control

The ultrafast nonadiabatic excited state dynamics of (PTZ-N-benzyl-acetylide) (trans-bis-trimethylphosphine) Pt(II) (acetylide-NDI-bis-methyl) 1, representative of a series of Pt(II) donor-bridge-acceptor assemblies experimentally studied by the Weinstein group, University of Sheffield, is investigated by means of wavepacket propagations based on the multiconfiguration time-dependent Hartree (MCTDH) method. On the basis of electronic structure data obtained at the time-dependent density functional theory (TD-DFT) level, the subpicosecond decay is simulated by solving an 11 electronic states multimode problem, up to 18 vibrational normal modes, including both spin-orbit coupling (SOC) and vibronic coupling. A careful analysis of the results, within the diabatic representation, provides the key features of the spin-vibronic mechanism at work in this complex, distinguishing between the spin-orbit and vibronically activated ultrafast processes within the excited states manifold. The knowledge of the key active normal modes that promote selectively the population of specific electronic excited states opens a route toward optical control by selectively exciting these modes in order to drive the associated nonadiabatic processes. Relevant simulations, over 2 ps, are proposed to assess the impact of these selective vibrational excitations on the branching ratio between the primary photoproducts, namely, bridge-acceptor charge-transfer (CT) and donor-acceptor charge-separated (CS) electronic states. Whereas the excitation of the localized acetylide bridge C≡C bond stretching does not modify drastically the population of the low-lying electronic states within the first two ps, vibrational excitation of the out-of-plane twisting motion of the N-benzyl group linked to the donor entity favors the population of the 1,3CS states at the expense of the lowest 1,3CT states. This quantum study opens the route to IR optical control experiments based on the specific alteration of vibrational normal modes that activate vibronic couplings between key electronic excited states. However, the presence of critical crossings along the PES channels associated with these normal modes and the role of concurrent SOC driven ultrafast transfers of population should not be underestimated.

Tuning ultrafast time-evolution of photo-induced charge-transfer states: A real-time electronic dynamics study in substituted indenotetracene derivatives

Photo-induced charge transfer (CT) states are pivotal in many technological and biological processes. A deeper knowledge of such states is mandatory for modeling the charge migration dynamics. Real-time time-dependent density functional theory (RT-TD-DFT) electronic dynamics simulations are employed to explicitly observe the electronic density time-evolution upon photo-excitation. Asymmetrically substituted indenotetracene molecules, given their potential application as n-type semiconductors in organic photovoltaic materials, are here investigated. Effects of substituents with different electron-donating characters are analyzed in terms of the overall electronic energy spacing and resulting ultrafast CT dynamics through linear response (LR-)TD-DFT and RT-TD-DFT based approaches. The combination of the computational techniques here employed provided direct access to the electronic density reorganization in time and to its spatial and rational representation in terms of molecular orbital occupation time evolution. Such results can be exploited to design peculiar directional charge dynamics, crucial when photoactive materials are used for light-harvesting applications.

Spin-Vibronic Dynamics in Open-Shell Systems beyond the Spin Hamiltonian Formalism

Vibronic coupling has a dramatic influence over a large number of molecular processes, ranging from photochemistry to spin relaxation and electronic transport. The simulation of vibronic coupling with multireference wave function methods has been largely applied to organic compounds, and only early efforts are available for open-shell systems such as transition metal and lanthanide complexes. In this work, we derive a numerical strategy to differentiate the molecular electronic Hamiltonian in the context of multireference ab initio methods and inclusive of spin-orbit coupling effects. We then provide a formulation of open quantum system dynamics able to predict the time evolution of the electron density matrix under the influence of a Markovian phonon bath up to fourth-order perturbation theory. We apply our method to Co(II) and Dy(III) molecular complexes exhibiting long spin relaxation times and successfully validate our strategy against the use of an effective spin Hamiltonian. Our study sheds light on the nature of vibronic coupling, the importance of electronic excited states in spin relaxation, and the need for high-level computational chemistry to quantify it.

Watching the Interplay between Photoinduced Ultrafast Charge Dynamics and Nuclear Vibrations

Here is presented the ultrafast hole-electron dynamics of photoinduced metal to ligand charge-transfer (MLCT) states in a Ru(II) complex, [Ru(dcbpy)2(NCS)2]4- (dcbpy = 4,4'-dicarboxy-2,2'-bipyridine), a photoactive molecule employed in dye sensitized solar cells. Via cutting-edge computational techniques, a tailored computational protocol is here presented and developed to provide a detailed analysis of the electronic manifold coupled with nuclear vibrations to better understand the nonradiative pathways and the resulting overall dye performances in light-harvesting processes (electron injection). Thus, the effects of different vibrational modes were investigated on both the electronic levels and charge transfer dynamics through a theoretical-computational approach. First, the linear response time-dependent density functional (LR-TDDFT) formalism was employed to characterize excitation energies and spacing among electronic levels (the electronic layouts). Then, to understand the ultrafast (femtosecond) charge dynamics on the molecular scale, we relied on the nonperturbative mean-field quantum electronic dynamics via real-time (RT-) TDDFT. Three vibrational modes were selected, representative for collective nuclear movements that can have a significant influence on the electronic structure: two involving NCS- ligands and one involving dcbpy ligands. As main results, we observed that such MLCT states, under vibrational distortions, are strongly affected and a faster interligand electron transfer mechanism is observed along with an increasing MLCT character of the adiabatic electronic states approaching closer in energy due to the vibrations. Such findings can help both in providing a molecular picture of multidimensional vibro-electronic spectroscopic techniques, used to characterize ultrafast coherent and noncoherent dynamics of complex systems, and to improve dye performances with particular attention to the study of energy or charge transport processes and vibronic couplings.

Quantum Biology and the Potential Role of Entanglement and Tunneling in Non-Targeted Effects of Ionizing Radiation: A Review and Proposed Model

It is well established that cells, tissues, and organisms exposed to low doses of ionizing radiation can induce effects in non-irradiated neighbors (non-targeted effects or NTE), but the mechanisms remain unclear. This is especially true of the initial steps leading to the release of signaling molecules contained in exosomes. Voltage-gated ion channels, photon emissions, and calcium fluxes are all involved but the precise sequence of events is not yet known. We identified what may be a quantum entanglement type of effect and this prompted us to consider whether aspects of quantum biology such as tunneling and entanglement may underlie the initial events leading to NTE. We review the field where it may be relevant to ionizing radiation processes. These include NTE, low-dose hyper-radiosensitivity, hormesis, and the adaptive response. Finally, we present a possible quantum biological-based model for NTE.

An energy decomposition and extrapolation scheme for evaluating electron transfer rate constants: a case study on electron self-exchange reactions of transition metal complexes

A simple approach to the analysis of electron transfer (ET) reactions based on energy decomposition and extrapolation schemes is proposed. The present energy decomposition and extrapolation-based electron localization (EDEEL) method represents the diabatic energies for the initial and final states using the adiabatic energies of the donor and acceptor species and their complex. A scheme for the efficient estimation of ET rate constants is also proposed. EDEEL is semi-quantitative by directly evaluating the seam-of-crossing region of two diabatic potentials. In a numerical test, EDEEL successfully provided ET rate constants for electron self-exchange reactions of thirteen transition metal complexes with reasonable accuracy. In addition, its energy decomposition and extrapolation schemes provide all the energy values required for activation-strain model (ASM) analysis. The ASM analysis using EDEEL provided rational interpretations of the variation of the ET rate constants as a function of the transition metal complexes. These results suggest that EDEEL is useful for efficiently evaluating ET rate constants and obtaining a rational understanding of their magnitudes.

The three kingdoms-Photoinduced electron transfer cascades controlled by electronic couplings

Excited states are the key species in photocatalysis, while the critical parameters that govern their applications are (i) excitation energy, (ii) accessibility, and (iii) lifetime. However, in molecular transition metal-based photosensitizers, there is a design tension between the creation of long-lived excited (triplet), e.g., metal-to-ligand charge transfer (3MLCT) states and the population of such states. Long-lived triplet states have low spin-orbit coupling (SOC) and hence their population is low. Thus, a long-lived triplet state can be populated but inefficiently. If the SOC is increased, the triplet state population efficiency is improved-coming at the cost of decreasing the lifetime. A promising strategy to isolate the triplet excited state away from the metal after intersystem crossing (ISC) involves the combination of transition metal complex and an organic donor/acceptor group. Here, we elucidate the excited state branching processes in a series of Ru(II)-terpyridyl push-pull triads by quantum chemical simulations. Scalar-relativistic time-dependent density theory simulations reveal that efficient ISC takes place along 1/3MLCT gateway states. Subsequently, competitive electron transfer (ET) pathways involving the organic chromophore, i.e., 10-methylphenothiazinyl and the terpyridyl ligands are available. The kinetics of the underlying ET processes were investigated within the semiclassical Marcus picture and along efficient internal reaction coordinates that connect the respective photoredox intermediates. The key parameter that governs the population transfer away from the metal toward the organic chromophore either by means of ligand-to-ligand (3LLCT; weakly coupled) or intra-ligand charge transfer (3ILCT; strongly coupled) states was determined to be the magnitude of the involved electronic coupling.

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.

Understanding Charge Dynamics in Dense Electronic Manifolds in Complex Environments

Photoinduced charge transfer (CT) excited states and their relaxation mechanisms can be highly interdependent on the environment effects and the consequent changes in the electronic density. Providing a molecular interpretation of the ultrafast (subpicosecond) interplay between initial photoexcited states in such dense electronic manifolds in condensed phase is crucial for improving and understanding such phenomena. Real-time time-dependent density functional theory is here the method of choice to observe the charge density, explicitly propagated in an ultrafast time domain, along with all time-dependent properties that can be easily extracted from it. A designed protocol of analysis for real-time electronic dynamics to be applied to time evolving electronic density related properties to characterize both in time and in space CT dynamics in complex systems is here introduced and validated, proposing easy to be read cross-correlation maps. As case studies to test such tools, we present the photoinduced charge-transfer electronic dynamics of 5-benzyluracil, a mimic of nucleic acid/protein interactions, and the metal-to-ligand charge-transfer electronic dynamics in water solution of [Ru(dcbpy)2(NCS)2]4-, dcbpy = (4,4'-dicarboxy-2,2'-bipyridine), or "N34-", a dye sensitizer for solar cells. Electrostatic and explicit ab initio treatment of solvent molecules have been compared in the latter case, revealing the importance of the accurate modeling of mutual solute-solvent polarization on CT kinetics. We observed that explicit quantum mechanical treatment of solvent slowed down the charge carriers mobilities with respect to the gas-phase. When all water molecules were modeled instead as simpler embedded point charges, the electronic dynamics appeared enhanced, with a reduced hole-electron distance and higher mean velocities due to the close fixed charges and an artificially increased polarization effect. Such analysis tools and the presented case studies can help to unveil the influence of the electronic manifold, as well as of the finite temperature-induced structural distortions and the environment on the ultrafast charge motions.

Nonadiabatic excited-state dynamics of ReCl(CO)3(bpy) in two different solvents

We present a study of excited-states relaxation of the complex ReCl(CO)₃(bpy) (bpy = 2,2-bipyridine) using a nonadiabatic TD-DFT dynamics on spin-mixed potential energy surfaces in explicit acetonitrile (ACN) and dimethylsulfoxide (DMSO) solutions up to 800 fs. ReCl(CO)₃(bpy) belongs to a group of important photosensitizers which show ultrafast biexponential subpicosecond fluorescence decay kinetics. The choice of solvents was motivated by the different excited-state relaxation dynamics observed in subpicosecond time-resolved IR (TRIR) experiments. Simulations of intersystem crossing (ISC) showed the development of spin-mixed states in both solvents. Transformation of time-dependent populations of spin-mixed states enabled to monitor the temporal evolution of individual singlet and triplet states, fitting of bi-exponential decay kinetics, and simulating the time-resolved fluorescence spectra that show only minor differences between the two solvents. Analysis of structural relaxation and solvent reorganization employing time-resolved proximal distribution functions pointed to the factors influencing the fluorescence decay time constants. Nonadiabatic dynamics simulations of time-evolution of electronic, molecular, and solvent structures emerge as a powerful technique to interpret time-resolved spectroscopic data and ultrafast photochemical reactivity.

Spin–vibronic coupling in the quantum dynamics of a Fe(III) trigonal-bipyramidal complex

The presence of a high density of excited electronic states in the immediate vicinity of the optically bright state of a molecule paves the way for numerous photo-relaxation channels. In transition-metal complexes, the presence of heavy atoms results in a stronger spin–orbit coupling, which enables spin forbidden spin-crossover processes to compete with the spin-allowed internal conversion processes. However, no matter how effectively the states cross around the Franck–Condon region, the degree of vibronic coupling, of both relativistic and non-relativistic nature, drives the population distribution among these states. One such case is demonstrated in this work for the intermediate-spin Fe(III) trigonal-bipyramidal complex. A quantum dynamical investigation of the photo-deactivation mechanism in the Fe(III) system is presented using the multi-configurational time-dependent Hartree approach based on the vibronic Hamiltonian whose coupling terms are derived from the state-averaged complete active space self-consistent field/complete active space with second-order perturbation theory (CASPT2) calculations and spin–orbit coupling of the scalar-relativistic CASPT2 states. The results of this study show that the presence of a strong (non-relativistic) vibronic coupling between the optically bright intermediate-spin state and other low-lying states of the same spin-multiplicity overpowers the spin–orbit coupling between the intermediate-spin and high-spin states, thereby lowering the chances of spin-crossover while exhibiting ultrafast relaxation among the intermediate-spin states. In a special case, where the population transfer pathway via the non-relativistic vibronic coupling is blocked, the probability of the spin-crossover is found to increase. This suggests that a careful modification of the complex by incorporation of heavier atoms with stronger relativistic effects can enhance the spin-crossover potential of Fe(III) intermediate-spin complexes.

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.

Are luminescent Ru2+ chelated complexes selective coordinative sensors for the detection of heavy cations?

The ability of [Ru(bpy)₂(bpym)]²⁺ (bpy = 2,2'-bipyridine; bpym = 2,2'-bipyrimidine) to probe specifically heavy cations has been investigated by means of density functional theory for transition metals, group 12 elements and Pb²⁺. On the basis of the calculated Gibbs free energies of complexation in water it is shown that all reactions are favorable with negative enthalpies except for Hg²⁺, with the transition metal cations forming stable bi-metallic complexes by coordination to the bpym ligand. Comparison between the optical and photophysical properties of the Ru²⁺ probe and those of the coordination compounds does not demonstrate a high selectivity due to very similar characteristics of the absorption and emission spectra. Whereas by complexation the lowest metal-to-ligand-charge-transfer (MLCT) shoulder of [Ru(bpy)₂(bpym)]²⁺ at 462 nm is more or less shifted to the red as a function of the cation, the second MLCT band at 415 nm, less sensitive to the complexation, gains in intensity and is slightly blue-shifted. The visible MLCT emission of [Ru(bpy)₂(bpym)]²⁺ at 706 nm is altered by complexation leading to near IR (800-900 nm) emission in most of the coordination compounds. Complexation to some transition metal cations (Fe, Co, Rh and Pd) generates low-lying metal-centered (MC) excited states that quench luminescence. In contrast to the conclusion of experimental findings by Kumar et al. (Chem. Commun. 2014, 50, 8488-8490), [Ru(bpy)₂(bpym)]²⁺ cannot be proposed as a fast and selective probe for monitoring Pd²⁺ in aqueous media. Indeed, it does not possess the optical and photophysical characteristics necessary to discriminate Pd²⁺ ions over a variety of other cations.

Electronic structure of strongly correlated systems: recent developments in multiconfiguration pair-density functional theory and multiconfiguration nonclassical-energy functional theory

Strong electron correlation plays an important role in transition-metal and heavy-metal chemistry, magnetic molecules, bond breaking, biradicals, excited states, and many functional materials, but it provides a significant challenge for modern electronic structure theory. The treatment of strongly correlated systems usually requires a multireference method to adequately describe spin densities and near-degeneracy correlation. However, quantitative computation of dynamic correlation with multireference wave functions is often difficult or impractical. Multiconfiguration pair-density functional theory (MC-PDFT) provides a way to blend multiconfiguration wave function theory and density functional theory to quantitatively treat both near-degeneracy correlation and dynamic correlation in strongly correlated systems; it is more affordable than multireference perturbation theory, multireference configuration interaction, or multireference coupled cluster theory and more accurate for many properties than Kohn–Sham density functional theory. This perspective article provides a brief introduction to strongly correlated systems and previously reviewed progress on MC-PDFT followed by a discussion of several recent developments and applications of MC-PDFT and related methods, including localized-active-space MC-PDFT, generalized active-space MC-PDFT, density-matrix-renormalization-group MC-PDFT, hybrid MC-PDFT, multistate MC-PDFT, spin–orbit coupling, analytic gradients, and dipole moments. We also review the more recently introduced multiconfiguration nonclassical-energy functional theory (MC-NEFT), which is like MC-PDFT but allows for other ingredients in the nonclassical-energy functional. We discuss two new kinds of MC-NEFT methods, namely multiconfiguration density coherence functional theory and machine-learned functionals.

This feature article overviews recent work on active spaces, matrix product reference states, treatment of quasidegeneracy, hybrid theory, density-coherence functionals, machine-learned functionals, spin–orbit coupling, gradients, and dipole moments.

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.

Ultrafast photoinduced dynamics in Prussian blue analogues

A review on ultrafast photoinduced processes in molecule-based magnets with an emphasis on Prussian blue analogues.