“Watching” Polaron Pair Formation from First-Principles Electron–Nuclear Dynamics

Blueshift of the CN stretching vibration of acetonitrile in solution: computational and experimental study

Acetonitrile, a polar molecule that cannot form hydrogen bonds on its own, interacts with solvent molecules mainly through the lone pair of its nitrogen atom and the π electrons of its CN triple bond [Correction added on 17 July 2024, after first online publication: Acetole has been changed to Acetonitrile in the preceeding sentence.]. Interestingly, acetonitrile exhibits an unexpected strengthening of the triple bond's force constant in an aqueous environment, leading to an upshift (blueshift) in the corresponding stretching vibration: this effect contrasts with the usual consequence of hydrogen bonding on the vibrational frequencies of the acceptor groups, that is, frequency redshift. This investigation elucidates this phenomenon using Raman spectroscopy to examine the behavior of acetonitrile in organic solvent, water, and silver ion aqueous solutions, where an even more pronounced upshift is observed. Raman spectroscopy is particularly well suited for analyzing aqueous solutions due to the minimal scattering effect of water molecules across most of the vibrational spectrum. Computational approaches, both static and dynamical, based on Density Functional Theory and hybrid functionals, are employed here to interpret these findings, and accurately reproduce the vibrational frequencies of acetonitrile in different environments. Our calculations also allow an explanation for this unique behavior in terms of electric charge displacements. On the other hand, the study of the interaction of acetonitrile with water molecules and metal ions is relevant for the use of this molecule as a solvent in both chemical and pharmaceutical applications.

Conformational and environmental effects on the electronic and vibrational properties of dyes for solar cell devices

The main challenge for solar cell devices is harvesting photons beyond the visible by reaching the red-edge (650-780 nm). Dye-sensitized solar cell (DSSC) devices combine the optical absorption and the charge separation processes by the association of a sensitizer as a light-absorbing material (dye molecules, whose absorption can be tuned and designed) with a wide band gap nanostructured semiconductor. Conformational and environmental effects (i.e., solvent, pH) can drastically influence the photophysical properties of molecular dyes. This study proposes a combined experimental and computational approach for the comprehensive investigation of the electronic and vibrational properties of a unique class of organic dye compounds belonging to the family of red-absorbing dyes, known as squaraines. Our focus lies on elucidating the intricate interplay between the molecular structure, vibrational dynamics, and optical properties of squaraines using state-of-the-art density functional theory calculations and spectroscopic techniques. Through systematic vibrational and optical analyses, we show that (i) the main absorption peak in the visible range is influenced by the conformational and protonation equilibria, (ii) the solvent polarity tunes the position of the UV-vis absorption, and (iii) the vibrational spectroscopy techniques (infrared and Raman) can be used as informative tools to distinguish between different conformations and protonation states. This comprehensive understanding offers valuable insights into the design and optimization of squaraine-based DSSCs for enhanced solar energy conversion efficiency.

Evidence of Excited-State Vibrational Mode Governing the Photorelaxation of a Charge-Transfer Complex

Modern, nonlinear, time-resolved spectroscopic techniques have opened new doors for investigating the intriguing but complex world of photoinduced ultrafast out-of-equilibrium phenomena and charge dynamics. The interaction between light and matter introduces an additional dimension, where the complex interplay between electronic and vibrational dynamics needs the most advanced theoretical-computational protocols to be fully understood on the molecular scale. In this study, we showcase the capabilities of ab initio molecular dynamics simulation integrated with a multiresolution wavelet protocol to carefully investigate the excited-state relaxation dynamics in a noncovalent complex involving tetramethylbenzene (TMB) and tetracyanoquinodimethane (TCNQ) undergoing charge transfer (CT) upon photoexcitation. Our protocol provides an accurate description that facilitates a direct comparison between transient vibrational analysis and time-resolved spectroscopic signals. This molecular level perspective enhances our understanding of photorelaxation processes confined in the adiabatic regime and offers an improved interpretation of vibrational spectra. Furthermore, it enables the quantification of anharmonic vibrational couplings between high- and low-frequency modes, specifically the TCNQ "rocking" and "bending" modes. Additionally, it identifies the primary vibrational mode that governs the adiabaticity between the ground state and the CT state. This comprehensive understanding of photorelaxation processes holds significant importance in the rational design and precise control of more efficient photovoltaic and sensor devices.

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.

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.

Bridge-Mediated Metal-to-Metal Electron and Hole Transfer in a Supermolecular Dinuclear Complex: A Computational Study Using Quantum Electron-Nuclear Dynamics

Bimetallic electron donor-acceptor complexes can facilitate electron and energy transfer with excellent structural control through synthetic design. In this work, we investigate the photochemical dynamics in a Ru-Cu bimetallic complex after photoexcitation of the Ru-centered charge transfer state. The physical underpinnings of the metal-to-metal directional charge transfer process are unraveled via analyses of the quantum electronic dynamics and electron-nuclear trajectories. The effects of molecular vibrations in the photoexcited state on the charge transfer processes are also analyzed.

Nature of the Ultrafast Interligands Electron Transfers in Dye-Sensitized Solar Cells

Charge-transfer dynamics and interligand electron transfer (ILET) phenomena play a pivotal role in dye-sensitizers, mostly represented by the Ru-based polypyridyl complexes, for TiO₂ and ZnO-based solar cells. Starting from metal-to-ligand charge-transfer (MLCT) excited states, charge dynamics and ILET can influence the overall device efficiency. In this letter, we focus on N3^4- dye ( [Ru(dcbpy)₂(NCS)₂]^4-, dcbpy = 4,4'-dicarboxy-2,2'-bipyridine) to provide a first direct observation with high time resolution (<20 fs) of the ultrafast electron exchange between bpy-like ligands. ILET is observed in water solution after photoexcitation in the ∼400 nm MLCT band, and assessment of its ultrafast time-scale is here given through a real-time electronic dynamics simulation on the basis of state-of-the-art electronic structure methods. Indirect effects of water at finite temperature are also disentangled by investigating the system in a symmetric gas-phase structure. As main result, remarkably, the ILET mechanism appears to be based upon a purely electronic evolution among the dense, experimentally accessible, MLCT excited states manifold at ∼400 nm, which rules out nuclear-electronic couplings and proves further the importance of the dense electronic manifold in improving the efficiency of dye sensitizers in solar cell devices.

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.

Direct observation of the solvent organization and nuclear vibrations of [Ru(dcbpy)2(NCS)2]4-, [dcbpy = (4,4'-dicarboxy-2,2'-bipyridine)], via ab initio molecular dynamics

Environmental effects can drastically influence the optical properties and photoreactivity of molecules, particularly in the presence of polar and/or protic solvents. In this work we investigate a negatively charged Ru(II) complex, [Ru(dcbpy)₂(NCS)₂]^4- [dcbpy = (4,4'-dicarboxy-2,2'-bipyridine)], in water solution, since this system belongs to a broader class of transition-metal compounds undergoing upon photo-excitation rapid and complex charge transfer (CT) dynamics, which can be dictated by structural rearrangement and solvent environment. Ab initio molecular dynamics (AIMD) relying on a hybrid quantum/molecular mechanics scheme is used to probe the equilibrium microsolvation around the metal complex in terms of radial distribution functions of the main solvation sites and solvent effects on the overall equilibrium structure. Then, using our AIMD-based generalized normal mode approach, we investigate how the ligand vibrational spectroscopic features are affected by water solvation, also contributing to the interpretation of experimental Infra-Red spectra. Two solvation sites are found for the ligands: the sulfur and the oxygen sites can interact on average with ∼4 and ∼3 water molecules, respectively, where a stronger interaction of the oxygen sites is highlighted. On average an overall dynamic distortion of the C₂ symmetric gas-phase structure was found to be induced by water solvation. Vibrational analysis reproduced experimental values for ligand symmetric and asymmetric stretchings, linking the observed shifts with respect to the gas-phase to a complex solvent distribution around the system. This is the groundwork for future excited-state nuclear and electronic dynamics to monitor non-equilibrium processes of CT excitation in complex environments, such as exciton migration in photovoltaic technologies.

Water-Mediated Excited State Proton Transfer of Pyranine-Acetate in Aqueous Solution: Vibrational Fingerprints from Ab Initio Molecular Dynamics

In this work, we simulate the excited state proton transfer (ESPT) reaction involving the pyranine photoacid and an acetate molecule as proton acceptor, connected by a bridge water molecule. We employ ab initio molecular dynamics combined with an hybrid quantum/molecular mechanics (QM/MM) framework. Furthermore, a time-resolved vibrational analysis based on the wavelet-transform allows one to identify two low frequency vibrational modes that are fingerprints of the ESPT event: a ring wagging and ring breathing. Their composition suggests their key role in optimizing the structure of the proton donor–acceptor couple and promoting the ESPT event. We find that the choice of the QM/MM partition dramatically affects the photoinduced reactivity of the system. The QM subspace was gradually extended including the water molecules directly interacting with the pyranine–water–acetate system. Indeed, the ESPT reaction takes place when the hydrogen bond network around the reactive system is taken into account at full QM level.

Exploring the Franck-Condon region of a photoexcited charge transfer complex in solution to interpret femtosecond stimulated Raman spectroscopy: excited state electronic structure methods to unveil non-radiative pathways

We present electronic structure methods to unveil the non-radiative pathways of photoinduced charge transfer (CT) reactions that play a main role in photophysics and light harvesting technologies. A prototypical π-stacked molecular complex consisting of an electron donor (1-chloronaphthalene, 1ClN) and an electron acceptor (tetracyanoethylene, TCNE) was investigated in dichloromethane solution for this purpose. The characterization of TCNE:π:1ClN in both its equilibrium ground and photoinduced low-lying CT electronic states was performed by using a reliable and accurate theoretical-computational methodology exploiting ab initio molecular dynamics simulations. The structural and vibrational time evolution of key vibrational modes is found to be in excellent agreement with femtosecond stimulated Raman spectroscopy experiments [R. A. Mathies et al., J. Phys. Chem. A, 2018, 122, 14, 3594], unveiling a correlation between vibrational fingerprints and electronic properties. The evaluation of nonadiabatic coupling matrix elements along generalized normal modes has made possible the interpretation on the molecular scale of the activation of nonradiative relaxation pathways towards the ground electronic state. In particular, two low frequency vibrational modes such as the out of plane bending and dimer breathing and the TCNE central C[double bond, length as m-dash]C stretching play a prominent role in relaxation phenomena from the electronic CT state to the ground state one.

Ultrafast Nonradiative Decay of a Dipolar Plasmon-like State in Naphthalene

Motivated by the uncertainty in our understanding of ultrafast plasmon decay mechanisms, we examine the effect of nuclear vibrations on the dynamical behavior of the strong plasmon-like dipole response of naphthalene, known as the β peak. The real-time time-dependent density functional (RT-TDDFT) method coupled with Ehrenfest molecular dynamics is used to describe the interconnected nuclear and electronic motion. Several vibrational modes promote drastic plasmon decay in naphthalene. The most astonishing finding of this study is that activation of one particular vibrational mode (corresponding to the B_1u representation in D_2h point group symmetry) leads to a continuous drop of the dipole response corresponding to the β peak into a totally symmetric, dark, quadrupolar electronic state. A second B_1u mode provokes the sharp plasmon-like peak to split due to the breaking of structural symmetry. Nonadiabatic coupling between a B_2g vibrational mode and the β peak (a B_1u electronic state) gives rise to a B_3u vibronic state, which can be identified as one of the p-band peaks that reside close in energy to the β peak energy. Overall, strong nonadiabatic coupling initiates plasmon decay into nearby electronic states in acenes, most importantly into dark states. These findings expand our knowledge about possible plasmon decay processes and pave the way for achieving high optical performance in acene-based materials such as graphene.

Multiresolution continuous wavelet transform for studying coupled solute-solvent vibrations via ab initio molecular dynamics

Vibrational analysis in solution and the theoretical determination of infrared and Raman spectra are of key importance in many fields of chemical interest. Vibrational band dynamics of molecules and their sensitivity to the environment can also be captured by these spectroscopies in their time dependent version. However, it is often difficult to provide an interpretation of the experimental data at the molecular scale, such as molecular mechanisms or the processes hidden behind them. In this work, we present a theoretical-computational protocol based on ab initio molecular dynamics simulations and a combination of normal-like (generalized) mode analysis of solute-solvent clusters with a wavelet transform, for the first time. The case study is the vibrational dynamics of N-methyl-acetamide (NMA) in water solution, a well-known model of hydration of peptides and proteins. Amide modes are typical bands of peptide and protein backbone, and their couplings with the environment are very challenging in terms of the accurate prediction of solvent induced intensity and frequency shifts. The contribution of water molecules surrounding NMA to the composition of generalized and time resolved modes is introduced in our vibrational analysis, showing unequivocally its influence on the amide mode spectra. It is also shown that such mode compositions need the inclusion of the first shell solvent molecules to be accurately described. The wavelet analysis is proven to be strongly recommended to follow the time evolution of the spectra, and to capture vibrational band couplings and frequency shifts over time, preserving at the same time a well-balanced time-frequency resolution. This peculiar feature also allows one to perform a combined structural-vibrational analysis, where the different strengths of hydrogen bond interactions can quantitatively affect the amide bands over time at finite temperature. The proposed method allows for the direct connection between vibrational modes and local structural changes, providing a link from the spectroscopic observable to the structure, in this case the peptide backbone, and its hydration layouts.

Time-Resolved Vibrational Analysis of Excited State Ab Initio Molecular Dynamics to Understand Photorelaxation: The Case of the Pyranine Photoacid in Aqueous Solution

We present a novel time-resolved vibrational analysis for studying photoinduced nuclear relaxation. Generalized modes velocities are defined from ab initio molecular dynamics and wavelet transformed, providing the time localization of vibrational signals in the electronic excited state. The photoexcited pyranine in aqueous solution is presented as a case study. The transient and sequential activation of the simulated vibrational signals is in good agreement with vibrational dynamics obtained from femtosecond stimulated Raman spectroscopy data.

Simulation of self-heating process on the nanoscale: a multiscale approach for molecular models of nanocomposite materials

A theoretical-computational protocol to model the Joule heating process in nanocomposite materials is presented. The proposed modeling strategy is based on post processing of trajectories obtained from large scale molecular simulations. This protocol, based on molecular models, is the first one to be applied to organic nanocomposites based on carbon nanotubes (CNT). This strategy allows to keep a microscopic explicit picture of the systems, to directly catch the molecular structure underlying the process under study and, at the same time, to include macroscopic boundary conditions fixed in the experiments. As validation and first application of the proposed strategy, a detailed investigation on CNT based organic composites is reported. The effect of CNT morphologies, concentration and working conditions on Joule heating has been modelled and compared with available experiments. Further experiments are performed also in this work to increase the number of comparisons especially in specific voltage ranges where available references from literature were missing. Simulations are in both qualitative and quantitative agreement with several experiments and trends reported in the recent literature, as well as with experiments performed in this work. The proposed approach combined with large scale hybrid particle-field molecular simulations can give insights and opens to way to a rational design of self-heating nanocomposites.

Ab-initio molecular dynamics and hybrid explicit-implicit solvation model for aqueous and nonaqueous solvents: GFP chromophore in water and methanol solution as case study

Solute-solvent interactions are proxies for understanding how the electronic density of a chromophore interacts with the environment in a more exhaustive way. The subtle balance between polarization, electrostatic, and non-bonded interactions need to be accurately described to obtain good agreement between simulations and experiments. First principles approaches providing accurate configurational sampling through molecular dynamics may be a suitable choice to describe solvent effects on solute chemical-physical properties and spectroscopic features, such as optical absorption of dyes. In this context, accurate energy potentials, obtained by hybrid implicit/explicit solvation methods along with employing nonperiodic boundary conditions, are required to represent bulk solvent around a large solute-solvent cluster. In this work, a novel strategy to simulate methanol solutions is proposed combining ab initio molecular dynamics, a hybrid implicit/explicit flexible solvent model, nonperiodic boundary conditions, and time dependent density functional theory. As case study, the robustness of the proposed protocol has been gauged by investigating the microsolvation and electronic absorption of the anionic green fluorescent protein chromophore in methanol and aqueous solution. Satisfactory results are obtained, reproducing the microsolvation layout of the chromophore and, as a consequence, the experimental trends shown by the optical absorption in different solvents.

An electron density based analysis to establish the electronic adiabaticity of proton coupled electron transfer reactions

Electrons and protons are the main actors in play in proton coupled electron transfer (PCET) reactions, which are fundamental in many biological (i.e., photosynthesis and enzymatic reactions) and electrochemical processes. The mechanism, energetics and kinetics of PCET reactions are strongly controlled by the coupling between the transferred electrons and protons. Concerted PCET reactions are classified according to the electronical adiabaticity degree of the process. To discriminate among different mechanisms, we propose a new analysis based on the use of electron density based indexes. We choose, as test case, the 3-Methylphenoxyl/phenol system in two different conformations to show how the proposed analysis is a suitable tool to discriminate between the different degree of adiabaticity of PCET processes. The very low computational cost of this procedure is extremely promising to analyze and provide evidences of PCET mechanisms ruling the reactivity of many biological and catalytic systems.

Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials

Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.

Investigating the character of excited states in TiO2 nanoparticles from topological descriptors: implications for photocatalysis

Excited state topological descriptors based on the attachment/detachment one-electron charge density are used to investigate the centroids of photoactive TiO₂ nanoclusters and nanoparticles.

Unveiling anharmonic coupling by means of excited state ab initio dynamics: application to diarylethene photoreactivity

In this work, excited state ab initio molecular dynamics together with a time resolved vibrational analysis is employed to shed light on the vibrational photoinduced dynamics of a well-known diarylethene molecule experiencing a ring opening reaction upon electronic excitation. The photoreactivity of diarylethenes is recognized to be controlled by a non-adiabatic intersection point between the ground and the first excited state surfaces. The computation of an energy scan, along a suitable reaction coordinate, allows us to identify the region of potential energy surfaces in which the ground (S0) and the first excited (S1) state are well separated. The adiabatic sampling of that region in S1 shows that in the first 3 picoseconds, the central CC bond, which is subject to break, oscillates in an antiphase with respect to the energy gap ΔE(S1 - S0). A multiresolution analysis based on the wavelet transform was then applied to the structural parameters extracted from the excited state dynamics. The wavelet maps show characteristic oscillations of the frequencies, mainly CC stretching and CCC bending localized on the central 4-ring moiety. Moreover, we have identified the main frequency (methyl wagging motion) involved in the modulation of these oscillations. The anharmonic coupling within a group of vibrational modes was therefore highlighted, in good agreement with experimental evidence. For the first time, a quantitative analysis of time resolved signals from a wavelet transform/ab initio molecular dynamics approach was performed.

Multistep and multiscale electron transfer and localization dynamics at a model electrolyte/metal interface

The lifetime, coupling, and localization dynamics of electronic states in molecular films near metal electrodes fundamentally determine their propensity to act as precursors or reactants in chemical reactions, crucial for a detailed understanding of charge transport and degradation mechanisms in batteries. In the current study, we investigate the formation dynamics of small polarons and their role as intermediate electronic states in thin films of dimethyl sulfoxide (DMSO) on Cu(111) using time- and angle-resolved two-photon photoemission spectroscopy. Upon photoexcitation, a delocalized DMSO electronic state is initially populated two monolayers from the Cu surface, becoming a small polaron on a 200 fs time scale, consistent with localization due to vibrational dynamics of the DMSO film. The small polaron is a precursor state for an extremely long-lived and weakly coupled multilayer electronic state, with a lifetime of several seconds, thirteen orders of magnitude longer than the small polaron. Although the small polaron in DMSO has a lifetime of 140 fs, its role as a precursor state for long-lived electronic states could make it an important intermediate in multistep battery reactivity.

Ultrafast polaron-pair dynamics in a poly(3-hexylthiophene-2,5-diyl) device influenced by a static electric field: insights into electric-field-related charge loss

Photogenerated polaron-pair ultrafast dynamics in poly(3-hexylthiophene)-based devices are found to be influenced by external electric fields via delayed field-induced singlet exciton dissociation, yielding a bimolecular decay contribution.

Momentum-resolved TDDFT algorithm in atomic basis for real time tracking of electronic excitation

Ultrafast electronic dynamics in solids lies at the core of modern condensed matter and materials physics. To build up a practical ab initio method for studying solids under photoexcitation, we develop a momentum-resolved real-time time dependent density functional theory (rt-TDDFT) algorithm using numerical atomic basis, together with the implementation of both the length and vector gauge of the electromagnetic field. When applied to simulate elementary excitations in two-dimensional materials such as graphene, different excitation modes, only distinguishable in momentum space, are observed. The momentum-resolved rt-TDDFT is important and computationally efficient for the study of ultrafast dynamics in extended systems.

Game of Frontier Orbitals: A View on the Rational Design of Novel Charge-Transfer Materials

Since the first application of frontier molecular orbitals (FMOs) to rationalize stereospecificity of pericyclic reactions, FMOs have remained at the forefront of chemical theory. Yet, the practical application of FMOs in the rational design and synthesis of novel charge transfer materials remains under-appreciated. In this Perspective, we demonstrate that molecular orbital theory is a powerful and universal tool capable of rationalizing the observed redox/optoelectronic properties of various aromatic hydrocarbons in the context of their application as charge-transfer materials. Importantly, the inspection of FMOs can provide instantaneous insight into the interchromophoric electronic coupling and polaron delocalization in polychromophoric assemblies, and therefore is invaluable for the rational design and synthesis of novel materials with tailored properties.

Effect of a heavy heteroatom on triplet formation and interactions in single conjugated polymer molecules and aggregates

Triplet formation and interactions with emissive singlet excitons are investigated in poly(3-hexylselenophene) (P3HS) using single molecule spectroscopy. P3HS is a heavy atom analog of the more commonly studied poly(3-hexylthiophene) (P3HT), a benchmark polymer for solar cells. P3HS tends to aggregate strongly which necessitates dilution to ultra-low levels within a solid inert host in order to resolve photophysical responses of single chains. Fluorescence excitation intensity modulation is performed on isolated P3HS chains using a sequence of rectangular pulses of varying intensities to probe the presence of spin-forbidden triplet excitons. Triplet population dynamics originating from singlet-triplet and triplet-triplet interactions appear as quenching of the initial fluorescence intensity to steady-state levels on characteristic time scales of ∼1-10 μs. Over 80% of all molecules studied display significant fluorescence intensity modulation (quenching depths >50%) indicative of efficient intersystem crossing and large triplet occupancies. Because triplets are highly localized and singlet-triplet and triplet-triplet annihilation rate constants are comparable to those of intersystem crossing, multiple triplets are present at any given time on single P3HS chains. Triplet lifetimes were estimated to be ∼4 μs (upper limit) determined from recovery to the ground electronic singlet state in the absence of light and, surprisingly, triplets vanish at the onset of P3HS aggregation. This result was unexpected since P3HS triplet formation takes place on time scales <30 ps making this process competitive with most accessible non-radiative deactivation pathways.

Ab Initio Transient Vibrational Spectral Analysis

Pump probe spectroscopy techniques have enabled the direct observation of a variety of transient molecular species in both ground and excited electronic states. Time-resolved vibrational spectroscopy is becoming an indispensable tool for investigating photoinduced nuclear dynamics of chemical systems of all kinds. On the other hand, a complete picture of the chemical dynamics encoded in these spectra cannot be achieved without a full temporal description of the structural relaxation, including the explicit time-dependence of vibrational coordinates that are substantially displaced from equilibrium by electronic excitation. Here we present a transient vibrational analysis protocol combining ab initio direct molecular dynamics and time-integrated normal modes introduced in this work, relying on the recent development of analytic time-dependent density functional theory (TDDFT) second derivatives for excited states. Prototypical molecules will be used as test cases, showing the evolution of the vibrational signatures that follow electronic excitation. This protocol provides a direct route to assigning the vibrations implicated in the (photo)dynamics of several (photoactive) systems.