Bifurcation of Excited-State Population Leads to Anti-Kasha Luminescence in a Disulfide-Decorated Organometallic Rhenium Photosensitizer.

We report a rhenium diimine photosensitizer equipped with a peripheral disulfide unit on one of the bipyridine ligands, [Re(CO)3(bpy)(S-Sbpy4,4)]+ (1+, bpy = 2,2'-bipyridine, S-Sbpy4,4 = [1,2]dithiino[3,4-c:6,5-c']dipyridine), showing anti-Kasha luminescence. Steady-state and ultrafast time-resolved spectroscopies complemented by nonadiabatic dynamics simulations are used to disclose its excited-state dynamics. The calculations show that after intersystem crossing the complex evolves to two different triplet minima: a (S-Sbpy4,4)-ligand-centered excited state (3LC) lying at lower energy and a metal-to-(bpy)-ligand charge transfer (3MLCT) state at higher energy, with relative yields of 90% and 10%, respectively. The 3LC state involves local excitation of the disulfide group into the antibonding σ* orbital, leading to significant elongation of the S-S bond. Intriguingly, it is the higher-lying 3MLCT state, which is assigned to display luminescence with a lifetime of 270 ns: a signature of anti-Kasha behavior. This assignment is consistent with an energy barrier ≥ 0.6 eV or negligible electronic coupling, preventing reaction toward the 3LC state after the population is trapped in the 3MLCT state. This study represents a striking example on how elusive excited-state dynamics of transition-metal photosensitizers can be deciphered by synergistic experiments and state-of-the-art calculations. Disulfide functionalization lays the foundation of a new design strategy toward harnessing excess energy in a system for possible bimolecular electron or energy transfer reactivity.

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)2]3+ (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)2]3+.

Efficient Reverse Intersystem Crossing in Carbene-Copper-Amide TADF Emitters via an Intermediate Triplet State.

The mechanism behind reverse intersystem crossing (rISC) in metal-based TADF emitters is still under debate. Thermal rISC necessitates small singlet/triplet energy gaps as realized in donor-acceptor systems with charge-transfer excited states. However, their associated spin-orbit couplings are too small to account for effective rISC. Here, we report the first nonadiabatic dynamics simulation of the rISC process in a carbene-copper(I)-carbazolyl TADF emitter. Efficient rISC on a picosecond time scale is demonstrated for an initial triplet minimum geometry that exhibits a perpendicular orientation of the ligands. The dynamics involves an intermediate higher-lying triplet state of metal-to-ligand charge transfer character (3 MLCT), which enables large spin-orbit couplings with the lowest singlet charge transfer state. The mechanism is completed in the S1 state, where the complex can return to a co-planar coordination geometry that presents high fluorescence efficiency.

peri-Acenoacene Ribbons with Zigzag BN-Doped Peripheries.

Here, we report the synthesis of BN-doped graphenoid nanoribbons, in which peripheral carbon atoms at the zigzag edges have been selectively replaced by boron and nitrogen atoms as BN and NBN motifs. This includes high-yielding ring closure key steps that, through N-directed borylation reaction using solely BBr3, allow the planarization of meta-oligoarylenyl precursors, through the formation of B-N and B-C bonds, to give ter-, quater-, quinque-, and sexi-arylenyl nanoribbons. X-ray single-crystal diffraction studies confirmed the formation of the BN and NBN motifs and the zigzag-edged topology of the regularly doped ribbons. Steady-state absorption and emission investigations at room temperature showed a systematic bathochromic shift of the UV-vis absorption and emission envelopes upon elongation of the oligoarylenyl backbone, with the nanoribbon emission featuring a TADF component. All derivatives displayed phosphorescence at 77 K. Electrochemical studies showed that the π-extension of the peri-acenoacene framework provokes a lowering of the first oxidative event (from 0.83 to 0.40 V), making these nanoribbons optimal candidates to engineer p-type organic semiconductors.

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 SO2, 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.

Intersystem Crossing and Triplet Dynamics in an Iron(II) N-Heterocyclic Carbene Photosensitizer.

The electronic excited states of the iron(II) complex [FeII(tpy)(pyz-NHC)]2+ [tpy = 2,2':6',2″-terpyridine; pyz-NHC = 1,1'-bis(2,6-diisopropylphenyl)pyrazinyldiimidazolium-2,2'-diylidene] and their relaxation pathways have been theoretically investigated. To this purpose, trajectory surface-hopping simulations within a linear vibronic coupling model including a 244-dimensional potential energy surface (PES) with 20 singlet and 20 triplet coupled states have been used. The simulations show that, after excitation to the lowest-energy absorption band of predominant metal-to-ligand charge-transfer character involving the tpy ligand, almost 80% of the population undergoes intersystem crossing to the triplet manifold in about 50 fs, while the remaining 20% decays through internal conversion to the electronic ground state in about 300 fs. The population transferred to the triplet states is found to deactivate into two different regions of the PESs, one where the static dipole moment is small and shows increased metal-centered character and another with a large static dipole moment, where the electron density is transferred from the tpy to pyz-NHC ligand. Coherent oscillations of 400 fs are observed between these two sets of triplet populations, until the mixture equilibrates to a ratio of 60:40. Finally, the importance of selecting suitable normal modes is highlighted-a choice that can be far from straightforward in transition-metal complexes with hundreds of degrees of freedom.

Finite-temperature Wigner phase-space sampling and temperature effects on the excited-state dynamics of 2-nitronaphthalene.

The concept of finite temperature Wigner phase-space sampling allowing the population of vibrationally excited states is introduced and employed to study temperature effects on the absorption spectrum of 2-nitronaphtalene (2NN) and its relaxation dynamics. It is found that, despite the fact that the general deactivation mechanism of 2NN after light irradiation does not change with increasing temperature, i.e., after excitation to the singlet manifold, 2NN deactivates via internal conversion in less than 100 fs and then undergoes intersystem crossing to the triplet manifold before decaying to the lowest triplet state in less than 150 fs, the intersystem crossing rate increases at higher temperatures. This increase is attributed to the appearance of a new deactivation pathway, which is not operative at lower temperatures. The present example illustrates that appropriate initial conditions beyond the idealized temperature of 0 K are indispensable to obtain reliable excited state mechanisms, which can be then related to experimental conditions.

Mechanism of Ultrafast Intersystem Crossing in 2-Nitronaphthalene.

Nitronaphthalene derivatives efficiently populate their electronically excited triplet states upon photoexcitation through ultrafast intersystem crossing (ISC). Despite having been studied extensively by time-resolved spectroscopy, the reasons behind their ultrafast ISC remain unknown. Herein, we present the first ab initio nonadiabatic molecular dynamics study of a nitronaphthalene derivative, 2-nitronaphthalene, including singlet and triplet states. We find that there are two distinct ISC reaction pathways involving different electronic states at distinct nuclear configurations. The high ISC efficiency is explained by the very small electronic and nuclear alterations that the chromophore needs to undergo during the singlet-triplet transition in the dominating ISC pathway after initial dynamics in the singlet manifold. The insights gained in this work are expected to shed new light on the photochemistry of other nitro polycyclic aromatic hydrocarbons that exhibit ultrafast intersystem crossing.

New compact density matrix averaged ANO basis sets for relativistic calculations.

When including relativistic effects in quantum chemical calculations, basis sets optimized for relativistic Hamiltonians such as the atomic natural orbital-relativistic core-correlated (ANO-RCC) basis set have to be used to avoid large errors that appear upon contraction of the basis set. While the large size of the ANO-RCC basis set in terms of primitive basis functions allows for highly accurate calculations, it also hinders its applicability to large sized systems due to the computational costs. To tackle this problem, a new compact relativistic ANO basis set, the ANO-eXtra Small (XS) basis set, is introduced for elements H-Ca. The number of primitive basis functions in ANO-XS is about half that of the ANO-RCC basis set. This greatly reduces the computational costs in the integral calculations especially when used in combination with Cholesky decomposition. At the same time, the ANO-XS basis set is able to predict molecular properties such as bond lengths and excitation energies with reasonable errors compared to the larger ANO-RCC basis set. The main intention for the ANO-XS basis set is to be used in conjunction with the ANO-RCC basis set for large systems that can be divided with regions demanding different qualities of basis sets. This is exemplified in CASPT2 calculations for an Ir(C3H4N)3 complex, where substituting the larger ANO-RCC for the compact ANO-XS basis set at the ligand atoms yields only minor differences for a large number of excited states compared to calculations employing the ANO-RCC basis set on all atoms. Thus, accurate calculations including relativistic effects for large systems become more affordable with the new ANO-XS basis set.

Communication: electron transfer mediated decay enabled by spin-orbit interaction in small krypton/xenon clusters.

In this work we study the influence of relativistic effects, in particular spin-orbit coupling, on electronic decay processes in KrXe2 clusters of various geometries. For the first time it is shown that inclusion of spin-orbit coupling has decisive influence on the accessibility of a specific decay pathway in these clusters. The radiationless relaxation process is initiated by a Kr 4s ionization followed by an electron transfer from xenon to krypton and a final second ionization of the system. We demonstrate the existence of competing electronic decay pathways depending in a subtle way on the geometry and level of theory. For our calculations a fully relativistic framework was employed where omission of spin-orbit coupling leads to closing of two decay pathways. These findings stress the relevance of an adequate relativistic description for clusters with heavy elements and their fragmentation dynamics.

Singlet and Triplet Pathways Determine the Thermal Z/E Isomerization of an Arylazopyrazole-Based Photoswitch.

Understanding the thermal isomerization mechanism of azobenzene derivatives is essential to designing photoswitches with tunable half-lives. Herein, we employ quantum chemical calculations, nonadiabatic transition state theory, and photosensitized experiments to unravel the thermal Z/E isomerization of a heteroaromatic azoswitch, the phenylazo-1,3,5-trimethylpyrazole. In contrast to the parent azobenzene, we predict two pathways to be operative at room temperature. One is a conventional ground-state reaction occurring via inversion of the aryl group, and the other is a nonadiabatic process involving intersystem crossing to the lowest-lying triplet state and back to the ground state, accompanied by a torsional motion around the azo bond. Our results illustrate that the fastest reaction rate is not controlled by the mechanism involving the lowest activation energy, but the size of the spin-orbit couplings at the crossing between the singlet and the triplet potential energy surfaces is also determinant. It is therefore mandatory to consider all of the multiple reaction pathways in azoswitches in order to predict experimental half-lives.

Can range-separated functionals be optimally tuned to predict spectra and excited state dynamics in photoactive iron complexes?

Density functional theory is an efficient computational tool to investigate photophysical and photochemical processes in transition metal complexes, giving invaluable assistance in interpreting spectroscopic and catalytic experiments. Optimally tuned range-separated functionals are particularly promising, as they were created to address some of the fundamental deficiencies present in approximate exchange-correlation functionals. In this paper, we scrutinize the selection of optimally tuned parameters and its influence on the excited state dynamics, using the example of the iron complex [Fe(cpmp)2]2+ with push-pull ligands. Various tuning strategies are contemplated based on pure self-consistent DFT protocols, as well as on the comparison with experimental spectra and multireference CASPT2 results. The two most promising sets of optimal parameters are then employed to carry out nonadiabatic surface-hopping dynamics simulations. Intriguingly, we find that the two sets lead to very different relaxation pathways and timescales. While the set of optimal parameters from one of the self-consistent DFT protocols predicts the formation of long-lived metal-to-ligand charge transfer triplet states, the set in better agreement with CASPT2 calculations leads to deactivation in the manifold of metal-centered states, in better agreement with the experimental reference data. These results showcase the complexity of iron-complex excited state landscapes and the difficulty of obtaining an unambiguous parametrization of long-range corrected functionals without experimental input.

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 and long-time excited state kinetics of an NIR-emissive vanadium(iii) complex II. Elucidating triplet-to-singlet excited-state dynamics.

We report the non-adiabatic dynamics of VIIICl3(ddpd), a complex based on the Earth-abundant first-row transition metal vanadium with a d2 electronic configuration which is able to emit phosphorescence in solution in the near-infrared spectral region. Trajectory surface-hopping dynamics based on linear vibronic coupling potentials obtained with CASSCF provide molecular-level insights into the intersystem crossing from triplet to singlet metal-centered states. While the majority of the singlet population undergoes back-intersystem crossing to the triplet manifold, 1-2% remains stable during the 10 ps simulation time, enabling the phosphorescence described in Dorn et al. Chem. Sci., 2021, DOI: 10.1039/D1SC02137K. Competing with intersystem crossing, two different relaxation channels via internal conversion through the triplet manifold occur. The nuclear motion that drives the dynamics through the different electronic states corresponds mainly to the increase of all metal-ligand bond distances as well as the decrease of the angles of trans-coordinated ligand atoms. Both motions lead to a decrease in the ligand-field splitting, which stabilizes the interconfigurational excited states populated during the dynamics. Analysis of the electronic character of the states reveals that increasing and stabilizing the singlet population, which in turn can result in enhanced phosphorescence, could be accomplished by further increasing the ligand-field strength.

The ANO-R Basis Set.

In this work, the new ANO-R basis set for all elements of the first six periods is introduced. The ANO-R basis set is an all-electron basis set that was constructed including scalar-relativistic effects of the exact-two component (X2C) Hamiltonian and modeling the atomic nucleus by a Gaussian charge distribution, which makes the basis set suitable for calculations of both light and heavy elements. For high accuracy, it takes advantage of the general contraction scheme and was developed at the CASSCF/CASPT2 level of theory. The distinguishing feature of the ANO-R basis set is its compactness in terms of both primitive and contracted basis functions, thus containing no superfluous functions for a given quality. An optimum number of primitive basis functions was selected based on studying the convergence toward the complete basis set limit for each element individually. The primitive basis sets were then contracted using the density-averaged atomic-natural-orbital (ANO) scheme, and suitable contraction levels were determined solely based on the natural orbital occupation numbers that describe the contribution of each natural orbital to the one-particle density matrix. Rather than following the common "split-valence n-tuple zeta plus polarization functions" structure, the resulting basis sets ANO-R0 to ANO-R3 possess a unique composition for each element, ensuring that no unnecessary functions are included while the basis sets are still balanced across the first six periods (H-Rn).

Nonadiabatic Dynamics Simulation Predict Intersystem Crossing in Nitroaromatic Molecules on a Picosecond Time Scale.

Previous time-resolved spectroscopic experiments and static quantum-chemical calculations attributed nitronaphthalene derivatives one of the fastest time scales for intersystem crossing within organic molecules, reaching the 100 fs mark. Nonadiabatic dynamics simulations on three nitronaphthalene derivatives challenge this view, showing that the experimentally observed ∼100 fs process corresponds to internal conversion in the singlet manifolds. Intersystem crossing, instead, takes place on a longer time scale of ∼1 ps. The dynamics simulations further reveal that the spin transitions occur via two distinct pathways with different contribution for the three systems, which are determined by electronic factors and the torsion of the nitro group. This study, therefore, indicates that the existence of sub-picosecond intersystem crossing in other nitroaromatic molecules should be questioned.

Vibrational Sampling and Solvent Effects on the Electronic Structure of the Absorption Spectrum of 2-Nitronaphthalene.

The influence of vibrational motion on electronic excited state properties is investigated for the organic chromophore 2-nitronaphtalene in methanol. Specifically, the performance of two vibrational sampling techniques - Wigner sampling and sampling from an ab initio molecular dynamics trajectory- is assessed, in combination with implicit and explicit solvent models. The effects of the different sampling/solvent combinations on the energy and electronic character of the absorption bands are analyzed in terms of charge transfer and exciton size, computed from the electronic transition density. The absorption spectra obtained using sampling techniques and its underlying properties are compared to those of the electronic excited states calculated at the Franck-Condon equilibrium geometry. It is found that the absorption bands of the vibrational ensembles are red-shifted compared to the Franck-Condon bright states, and this red-shift scales with the displacement from the equilibrium geometry. Such displacements are found larger and better described when using ensembles from the harmonic Wigner distribution than snapshots from the molecular dynamics trajectory. Particularly relevant is the torsional motion of the nitro group that quenches the charge transfer character of some of the absorption bands. This motion, however, is better described in the molecular dynamics trajectory. Thus, none of the vibrational sampling approaches can satisfactorily capture all important aspects of the nuclear motion. The inclusion of solvent also red-shifts the absorption bands with respect to the gas phase. This red-shift scales with the charge-transfer character of the bands and is found larger for the implicit than for the explicit solvent model. The advantages and drawbacks of the different sampling and solvent models are discussed to guide future research on the calculation of UV-vis spectra of nitroaromatic compounds.

The IPEA dilemma in CASPT2.

Multi-configurational second order perturbation theory (CASPT2) has become a very popular method for describing excited-state properties since its development in 1990. To account for systematic errors found in the calculation of dissociation energies, an empirical correction applied to the zeroth-order Hamiltonian, called the IPEA shift, was introduced in 2004. The errors were attributed to an unbalanced description of open-shell versus closed-shell electronic states and is believed to also lead to an underestimation of excitation energies. Here we show that the use of the IPEA shift is not justified and the IPEA should not be used to calculate excited states, at least for organic chromophores. This conclusion is the result of three extensive analyses. Firstly, we survey the literature for excitation energies of organic molecules that have been calculated with the unmodified CASPT2 method. We find that the excitation energies of 356 reference values are negligibly underestimated by 0.02 eV. This value is an order of magnitude smaller than the expected error based on the calculation of dissociation energies. Secondly, we perform benchmark full configuration interaction calculations on 137 states of 13 di- and triatomic molecules and compare the results with CASPT2. Also in this case, the excited states are underestimated by only 0.05 eV. Finally, we perform CASPT2 calculations with different IPEA shift values on 309 excited states of 28 organic small and medium-sized organic chromophores. We demonstrate that the size of the IPEA correction scales with the amount of dynamical correlation energy (and thus with the size of the system), and gets immoderate already for the molecules considered here, leading to an overestimation of the excitation energies. It is also found that the IPEA correction strongly depends on the size of the basis set. The dependency on both the size of the system and of the basis set, contradicts the idea of a universal IPEA shift which is able to compensate for systematic CASPT2 errors in the calculation of excited states.

Quenching of Charge Transfer in Nitrobenzene Induced by Vibrational Motion.

Although nitrobenzene is the smallest nitro-aromatic molecule, the nature of its electronic structure is still unclear. Most notably, the lowest-energy absorption band was assessed in numerous studies providing conflicting results regarding its charge-transfer character. In this study, we employ a combination of molecular dynamics and quantum chemical methods to disentangle the nature of the lowest-energy absorption band of nitrobenzene. Surprisingly, the charge-transfer transition from the benzene moiety to the nitro group is found to be quenched by a flow of charge into the opposite direction induced by vibrational motion. Beyond clarifying the charge-transfer character of nitrobenzene, we show that the widely used approach of analyzing the ground-state minimum-energy geometry provides completely wrong conclusions.