No Arabic abstract
Excited state electron and hole transfer underpin fundamental steps in processes such as exciton dissociation at photovoltaic heterojunctions, photoinduced charge transfer at electrodes, and electron transfer in photosynthetic reaction centers. Diabatic states corresponding to charge or excitation localized species, such as locally excited and charge transfer states, provide a physically intuitive framework to simulate and understand these processes. However, obtaining accurate diabatic states and their couplings from adiabatic electronic states generally leads to inaccurate results when combined with low-tier electronic structure methods, such as time dependent density functional theory (TDDFT), and exorbitant computational cost when combined with high-level wavefunction-based methods. Here we introduce a DFT-based diabatization scheme, {Delta}-ALMO(MSDFT2), which directly constructs the diabatic states using absolutely localized molecular orbitals (ALMOs). We demonstrate that our method, which combines ALMO calculations with the {Delta}SCF technique to construct electronically excited diabatic states and obtains their couplings with charge-transfer states using our MSDFT2 scheme, gives accurate results for excited state electron and hole transfer in both charged and uncharged systems that underlie DNA repair, charge separation in donor-acceptor dyads, chromophore-to-solvent electron transfer, and singlet fission. This framework for the accurate and efficient construction of excited state diabats and evaluation of their couplings directly from DFT thus offers a route to simulate and elucidate photoinduced electron and hole transfer in large disordered systems, such as those encountered in the condensed phase.
We have investigated photoinduced intramolecular electron transfer dynamics following metal-to-ligand charge-transfer (MLCT) excitation of [Fe(CN)$_4$(2,2-bipyridine)]$^{2-}$ (1), [Fe(CN)$_4$(2,3-bis(2-pyridyl)pyrazine)]$^{2-}$ (2) and [Fe(CN)$_4$(2,2-bipyrimidine)]$^{2-}$ (3) complexes in various solvents with static and time-resolved UV-visible absorption spectroscopy and Fe 2p3d resonant inelastic X-ray scattering. We observe $^3$MLCT lifetimes from 180 fs to 67 ps over a wide range of MLCT energies in different solvents by utilizing the strong solvatochromism of the complexes. Intramolecular electron transfer lifetimes governing $^3$MLCT relaxation increase monotonically and (super)exponentially as the $^3$MLCT energy is decreased in 1 and 2 by changing the solvent. This behavior can be described with non-adiabatic classical Marcus electron transfer dynamics along the indirect $^3$MLCT->$^3$MC pathway, where the $^3$MC is the lowest energy metal-centered (MC) excited state. In contrast, the $^3$MLCT lifetime in 3 changes non-monotonically and exhibits a maximum. This qualitatively different behaviour results from direct electron transfer from the $^3$MLCT to the electronic ground state (GS). This pathway involves nuclear tunnelling for the high-frequency polypyridyl skeleton mode ($hbaromega$ = 1530 cm$^{-1}$), which is more displaced for 3 than for either 1 or 2, therefore making the direct pathway significantly more efficient in 3. To our knowledge, this is the first observation of an efficient $^3$MLCT->GS relaxation pathway in an Fe polypyridyl complex. Our study suggests that further extending the MLCT state lifetime requires (1) lowering the $^3$MLCT state energy with respect to the $^3$MC state and (2) suppressing the intramolecular distortion of the electron-accepting ligand in the $^3$MLCT excited state to suppress the rate of direct $^3$MLCT->GS electron transfer.
We investigate photoinduced proton-coupled electron transfer (PI-PCET) reaction through a recently devel- oped quasi-diabatic (QD) quantum dynamics propagation scheme. This scheme enables interfacing accurate diabatic-based quantum dynamics approaches with adiabatic electronic structure calculations for on-the-fly simulations. Here, we use the QD scheme to directly propagate PI-PCET quantum dynamics with the di- abatic Partial Linearized Density Matrix (PLDM) path-integral approach with the instantaneous adiabatic electron-proton vibronic states. Our numerical results demonstrate the importance of treating proton quan- tum mechanically in order to obtain accurate PI-PCET dynamics, as well as the role of solvent fluctuation and vibrational relaxation on proton tunneling in various reaction regimes that exhibit different kinetic iso- tope effects. This work opens the possibility to study the challenging PI-PCET reactions through accurate diabatic quantum dynamics approaches combined with efficient adiabatic electronic structure calculations.
Using a set of oscillator strengths and excited-state dipole moments of near full configuration interaction (FCI) quality determined for small compounds, we benchmark the performances of several single-reference wave function methods (CC2, CCSD, CC3, CCSDT, ADC(2), and ADC(3/2)) and time-dependent density-functional theory (TD-DFT) with various functionals (B3LYP, PBE0, M06-2X, CAM-B3LYP, and $omega$B97X-D). We consider the impact of various gauges (length, velocity, and mixed) and formalisms: equation of motion (EOM) emph{vs} linear response (LR), relaxed emph{vs} unrelaxed orbitals, etc. Beyond the expected accuracy improvements and a neat decrease of formalism sensitivy when using higher-order wave function methods, the present contribution shows that, for both ADC(2) and CC2, the choice of gauge impacts more significantly the magnitude of the oscillator strengths than the choice of formalism, and that CCSD yields a notable improvement on this transition property as compared to CC2. For the excited-state dipole moments, switching on orbital relaxation appreciably improves the accuracy of both ADC(2) and CC2, but has a rather small effect at the CCSD level. Going from ground to excited states, the typical errors on dipole moments for a given method tend to roughly triple. Interestingly, the ADC(3/2) oscillator strengths and dipoles are significantly more accurate than their ADC(2) counterparts, whereas the two models do deliver rather similar absolute errors for transition energies. Concerning TD-DFT, one finds: i) a rather negligible impact of the gauge on oscillator strengths for all tested functionals (except for M06-2X); ii) deviations of ca.~0.10 D on ground-state dipoles for all functionals; iii) the better overall performance of CAM-B3LYP for the two considered excited-state properties.
Two-photon excitation is an attractive means for controlling chemistry in both space and time. Isoenergetic one- and two-photon excitations (OPE and TPE) in non-centrosymmetric molecules are often assumed to reach the same excited state and, hence, to produce similar excited-state reactivity. We compare the solvent-to-solute excited-state proton transfer of the super photobase FR0-SB following isoenergetic OPE and TPE. We find up to 62 % increased reactivity following TPE compared to OPE. From steady-state spectroscopy, we rule out the involvement of different excited states and find that OPE and TPE spectra are identical in non-polar solvents but not in polar ones. We propose that differences in the matrix elements that contribute to the two-photon absorption cross sections lead to the observed enhanced isoenergetic reactivity, consistent with the predictions of our high-level coupled-cluster-based computational protocol. We find that polar solvent configurations favor greater dipole moment change between ground and excited states, which enters the probability for two-photon excitations as the absolute value squared. This, in turn, causes a difference in the Franck-Condon region reached via TPE compared to OPE. We conclude that a new method has been found for controlling chemical reactivity via the matrix elements that affect two-photon cross sections, which may be of great utility for spatial and temporal precision chemistry.
We introduce natural transition geminals as a means to qualitatively understand a transition where double excitations are important. The first two $A_{1}$ singlet states of the CH cation are used as an initial example. We calculate these states with configuration interaction singles (CIS) and state-averaged Monte Carlo configuration interaction (SA-MCCI). For each method we compare the important natural transition geminals with the dominant natural transition orbitals. We then compare SA-MCCI and full configuration interaction (FCI) with regards to the natural transition geminals using the beryllium atom. We compare using the natural transition geminals with analyzing the important configurations in the CI expansion to give the dominant transition for the beryllium atom and the carbon dimer. Finally we calculate the natural transition geminals for two electronic excitations of formamide.