Nonadiabatic Effects on Defect Diffusion in Silicon-Doped Nanographenes.

Single atom impurities in graphene, substitutional silicon defects in particular, have been observed to diffuse under electron beam irradiation. However, the relative importance of elastic and inelastic scattering in facilitating their mobility remains unclear. Here, we employ excited-state electronic structure calculations to explore potential inelastic effects, and find an electronically nonadiabatic excited-state silicon diffusion pathway involving "softened" Si-C bonding that presents an ∼2 eV lower diffusion barrier than the ground-state pathway. Beam-induced transition rates to this state indicate that the excited-state pathway is accessible through irradiation of the defect site. However, even in the limit of fully elastic scattering, upward nonadiabatic transitions are also possible along the diffusion coordinate, increasing the diffusion barrier and further demonstrating the potential for electronic nonadiabaticity to influence beam-induced atomic transformations in materials. We also propose some experimentally testable signatures of such excited-state pathways.

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 B1u representation in D2h 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 B1u mode provokes the sharp plasmon-like peak to split due to the breaking of structural symmetry. Nonadiabatic coupling between a B2g vibrational mode and the ß peak (a B1u electronic state) gives rise to a B3u 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.

Role of Vibrational Dynamics on Excited-State Electronic Coherence in a Binuclear Platinum Complex.

Long-lived quantum coherence between excited electronic states can enable highly efficient energy and charge transport processes in chemical systems. Recent pump-probe experiments on binuclear platinum complexes identified persistent, periodic beating of transient absorption anisotropy signals, indicating long excited-state coherence lifetimes. Our previous simulations of the electronic dynamics of these complexes indicate that coherence lifetimes are sensitive to the balance between competitive electronic couplings. The complexes with shorter Pt-Pt distances underwent no appreciable dephasing in the limit of static nuclei, motivating the inclusion of nuclear motion into our simulations. The tert-butyl-substituted complex is studied in this work using the Ehrenfest method for mixed quantum-classical dynamics to investigate the role of vibrational dynamics on a complex shown to support long coherence lifetimes. Results indicate that the inclusion of excited-state vibrations drives a rapid collapse of the two-state coherence prior to the experimentally determined intersystem crossing. This further suggests singlet excited-state coherences may not be prerequisites for long-lived triplet coherences.

Ab Initio Excited-State Transient Raman Analysis.

Time-resolved Raman spectroscopy has proven useful for studying the formation of polarons in conjugated polymers, verifying the presence of reactive intermediates in photochemical reactions, investigating nonradiative transitions in the short lifetime of the photoexcited species, and resolving electron-phonon coupling strengths and exciton dissociation in crystalline materials. In this paper, we present an excited state transient Raman analysis protocol combining ab initio direct molecular dynamics, transient excited state Hessian, and excited state nonresonant Raman activities evaluations. Prototypical molecules are used as test cases, showing the evolution of the transient Raman signatures that follow electronic excitation. This protocol provides a direct route to assigning the vibrations implicated in the (photo)dynamics of several (photoactive) systems, complementary to the transient infrared analysis.

Can Excited State Electronic Coherence Be Tuned via Molecular Structural Modification? A First-Principles Quantum Electronic Dynamics Study of Pyrazolate-Bridged Pt(II) Dimers.

Materials and molecular systems exhibiting long-lived electronic coherence can facilitate coherent transport, opening the door to efficient charge and energy transport beyond traditional methods. Recently, signatures of a possible coherent, recurrent electronic motion were identified in femtosecond pump-probe spectroscopy experiments on a binuclear platinum complex, where a persistent periodic beating in the transient absorption signal's anisotropy was observed. In this study, we investigate the excitonic dynamics that underlie the suspected electronic coherence for a series of binuclear platinum complexes exhibiting a range of interplatinum distances. Results suggest that the long-lived coherence can only result when competitive electronic couplings are in balance. At longer Pt-Pt distances, the electronic couplings between the two halves of the binuclear system weaken, and exciton localization and recombination is favored on short time scales. For short Pt-Pt distances, electronic couplings between the states in the coherent superposition are stronger than the coupling with other excitonic states, leading to long-lived coherence.

"Watching" Polaron Pair Formation from First-Principles Electron-Nuclear Dynamics.

The formation of polaron pairs is one of the important photophysical processes that take place after the excitation in semiconducting organic polymers. First-principles Ehrenfest excited-state dynamics is a unique tool to investigate ultrafast photoinduced charge carrier dynamics and related nonequilibrium processes involving correlated electron-nuclear dynamics. In this work the formation of polaron pairs and their dynamical evolution in an oligomer of seven thiophene units is investigated with a combined approach of first-principles exciton-nuclear dynamics and wavelet analysis. The real-time formation of a polaron pair can be observed in the dipole evolution during the excited-state dynamics. The possible driving force of the polaron pair formation is investigated through qualitative correlation between the structural dynamics and the dipole evolution. The time-dependent characteristics and spectroscopic consequences of the polaron pair formation are probed using the wavelet analysis.

Time-dependent non-equilibrium dielectric response in QM/continuum approaches.

The Polarizable Continuum Models (PCMs) are some of the most inexpensive yet successful methods for including the effects of solvation in quantum-mechanical calculations of molecular systems. However, when applied to the electronic excitation process, these methods are restricted to dichotomously assuming either that the solvent has completely equilibrated with the excited solute charge density (infinite-time limit), or that it retains the configuration that was in equilibrium with the solute prior to excitation (zero-time limit). This renders the traditional PCMs inappropriate for resolving time-dependent solvent effects on non-equilibrium solute electron dynamics like those implicated in the instants following photoexcitation of a solvated molecular species. To extend the existing methods to this non-equilibrium regime, we herein derive and apply a new formalism for a general time-dependent continuum embedding method designed to be propagated alongside the solute's electronic degrees of freedom in the time domain. Given the frequency-dependent dielectric constant of the solvent, an equation of motion for the dielectric polarization is derived within the PCM framework and numerically integrated simultaneously with the time-dependent Hartree fock/density functional theory equations. Results for small molecular systems show the anticipated dipole quenching and electronic state dephasing/relaxation resulting from out-of-phase charge fluctuations in the dielectric and embedded quantum system.

Ab initio two-component Ehrenfest dynamics.

We present an ab initio two-component Ehrenfest-based mixed quantum/classical molecular dynamics method to describe the effect of nuclear motion on the electron spin dynamics (and vice versa) in molecular systems. The two-component time-dependent non-collinear density functional theory is used for the propagation of spin-polarized electrons while the nuclei are treated classically. We use a three-time-step algorithm for the numerical integration of the coupled equations of motion, namely, the velocity Verlet for nuclear motion, the nuclear-position-dependent midpoint Fock update, and the modified midpoint and unitary transformation method for electronic propagation. As a test case, the method is applied to the dissociation of H2 and O2. In contrast to conventional Ehrenfest dynamics, this two-component approach provides a first principles description of the dynamics of non-collinear (e.g., spin-frustrated) magnetic materials, as well as the proper description of spin-state crossover, spin-rotation, and spin-flip dynamics by relaxing the constraint on spin configuration. This method also holds potential for applications to spin transport in molecular or even nanoscale magnetic devices.

From charge-transfer to a charge-separated state: a perspective from the real-time TDDFT excitonic dynamics.

In-chain donor/acceptor block copolymers comprised of alternating electron rich/poor moieties are emerging as promising semiconducting chromophores for use in organic photovoltaic devices. The mobilities of charge carriers in these materials are experimentally probed using gated organic field-effect transistors to quantify electron and hole mobilities, but a mechanistic understanding of the relevant charge diffusion pathways is lacking. To elucidate the mechanisms of electron and hole transport following excitation to optically accessible low-lying valence states, we utilize mean-field quantum electronic dynamics in the TDDFT formalism to explicitly track the evolution of these photo-accessible states. From the orbital pathway traversed in the dynamics, p- and n-type conductivities can be distinguished. The electronic dynamics of the studied polymers show the time-resolved transitions between the initial photoexcited state, a tightly-bound excitonic state that is dark to the ground state, and a partially charge separated state indicated by long-lived, out-of-phase charge oscillations along the polymer backbone. The frequency of these charge oscillations yields an insight into the characteristic mobilities of charge carriers in these materials. When the barycenters of the electron and hole densities are followed during the dynamics, a pseudo-classical picture for the translation of charge carrier densities along the polymer backbone emerges that clarifies a crucial aspect in the design of efficient organic photovoltaic materials.

Non-adiabatic molecular dynamics investigation of photoionization state formation and lifetime in Mn²âº-doped ZnO quantum dots.

The unique electronic structure of Mn(2+)-doped ZnO quantum dots gives rise to photoionization states that can be used to manipulate the magnetic state of the material and to generate zero-reabsorption luminescence. Fast formation and long non-radiative decay of this photoionization state is a necessary requirement for these important applications. In this work, surface hopping based non-adiabatic molecular dynamics are used to demonstrate the fast formation of a metal-to-ligand charge transfer state in a Mn(2+)-doped ZnO quantum dot. The formation occurs on an ultrafast timescale and is aided by the large density of states and significant mixing of the dopant Mn(2+) 3dt2 levels with the valence-band levels of the ZnO lattice. The non-radiative lifetime of the photoionization states is also investigated.

Imaging Strain-Localized Single-Photon Emitters in Layered GaSe below the Diffraction Limit.

Nanoscale strain control of exciton funneling is an increasingly critical tool for the scalable production of single photon emitters (SPEs) in two-dimensional materials. However, conventional far-field optical microscopies remain constrained in spatial resolution by the diffraction limit and thus can provide only a limited description of nanoscale strain localization of SPEs. Here, we quantify the effects of nanoscale heterogeneous strain on the energy and brightness of GaSe SPEs on nanopillars with correlative cathodoluminescence, photoluminescence, and atomic force microscopy, supported by density functional theory simulations. We report the strain-localized SPEs have a broad range of emission wavelengths from 620 to 900 nm. We reveal substantial strain-controlled SPE wavelength tunability over a ∼100 nm spectral range and 2 orders of magnitude enhancement in the SPE brightness at the pillar center due to Type-I exciton funneling. In addition, we show that radiative biexciton cascade processes contribute to observed CL photon superbunching. Also, the GaSe SPEs show excellent stability, where their properties remain unchanged after electron beam exposure. We anticipate that this comprehensive study on the nanoscale strain control of two-dimensional SPEs will provide key insights to guide the development of truly deterministic quantum photonics.

Extracting Inelastic Scattering Cross Sections for Finite and Aperiodic Materials from Electronic Dynamics Simulations.

Explicit time-dependent electronic structure theory methods are increasingly prevalent in the areas of condensed matter physics and quantum chemistry, with the broad-band optical absorptivity of molecular and small condensed-phase systems nowadays routinely studied with such approaches. In this paper, it is demonstrated that electronic dynamics simulations can similarly be employed to study cross sections for the scattering-induced electronic excitations probed in nonresonant inelastic X-ray scattering and momentum-resolved electron energy loss spectroscopies. A method is put forth for evaluating the electronic dynamic structure factor, which involves the application of a momentum boost-type perturbation and transformation of the resulting reciprocal space density fluctuations into the frequency domain. Good agreement is first demonstrated between the dynamic structure factor extracted from these electronic dynamics simulations and the corresponding transition matrix elements from linear response theory. The method is then applied to some extended (quasi)one-dimensional systems, for which the wave vector becomes a good quantum number in the thermodynamic limit. Finally, the dispersion of many-body excitations in a series of hydrogen-terminated graphene flakes (and twisted bilayers thereof) is investigated to highlight the utility of the presented approach for capturing morphology-dependent effects in the inelastic scattering cross sections of nanostructured and/or noncrystalline materials.

Understanding Beam-Induced Electronic Excitations in Materials.

Optical electronic absorption spectroscopy and corresponding selection rules have played a prominent role in the development of the theory of electronic structure of materials. In this work, we expand the modern toolbox of chemistry and materials science to include electron and ion microscopies and spectroscopies, allowing spatially resolved interrogation of materials' atomic and electronic structures by beams of charged particles. Specifically, we formulate and discuss the selection rules for electronic excitation due to the interaction between materials and beams of charged particles. We show that transition probabilities for point charge-induced electronic excitations depend strongly on the position of the externally charged particles and can significantly deviate from those derived from the electric dipole (long-wavelength) approximation. We present and implement expressions within the linear response time-dependent density functional theory (TD-DFT) framework for rates of transition between the ground and excited states induced by an external point charge. The point charge-induced transition rates for particular electronic excitations from linear response TD-DFT were validated through comparison to excited state populations from real time TD-DFT simulations following an impulsive point charge perturbation, then evaluated on a three-dimensional grid to map their spatial dependence for a small polybenzoid. This method, when combined with information about excited state energy gradients, represents a first step toward an ab initio framework for probing the structural response of materials under irradiation by charged particles resulting from inelastic scattering. Engineering electron beam interaction with matter to manipulate single atoms and localized electronic states offers a revolutionary new regime beyond laser excitation.

The Role of Excited-State Proton Relays in the Photochemical Dynamics of Water Nanodroplets.

In this work, we applied nonadiabatic excited-state molecular dynamics in tandem with ab initio electronic structure theory to illustrate a complete mechanistic landscape underpinning the ultraviolet absorption-initiated photochemical dynamics in water nanodroplets. The goal is to understand the nonequilibrium excited-state molecular dynamics initiated by the relaxation of a solvated photoelectron and consequential photochemical processes. The lowest-lying excited state shows the proton dissociation for a single water molecule forming intermediate hydronium complexes through a proton relay. At approximately 100 fs, the proton relay process gives rise to the relaxation of the excited state accompanied by a rapid increase in the nonadiabatic coupling strength with the ground state, and the nanodroplet nonradiatively decays. The nonadiabatic transition to the ground state produces excited vibrational states that facilitate the recombination of the dissociated proton and hydroxyl group, eventually leading to the desorption of water molecules from the nanodroplet. Additionally, lifetimes of transient photochemical events are also resolved for the relaxation of a solvated electron, excited-state proton relay, and nonradiative transition.

Coupling Real-Time Time-Dependent Density Functional Theory with Polarizable Force Field.

Real-time time-dependent density functional theory (RT-TDDFT) is a powerful tool for obtaining spectroscopic observables and understanding complex, time-dependent properties. Currently, performing RT-TDDFT calculations on large, fully quantum mechanical systems is not computationally feasible. Previously, polarizable mixed quantum mechanical and molecular mechanical (QM/MMPol) models have been successful in providing accurate, yet efficient, approximations to a fully quantum mechanical system. Here we develop a coupling scheme between induced dipole based QM/MMPol and RT-TDDFT. Our approach is validated by comparing calculated spectra with both real-time and linear-response TDDFT calculations. The model developed within provides an accurate method for performing RT-TDDFT calculations on extended systems while accounting for mutual polarization between the quantum mechanical and molecular mechanical regions.

Can Quantized Vibrational Effects Be Obtained from Ehrenfest Mixed Quantum-Classical Dynamics?

We explore the question of whether mean-field or "Ehrenfest" mixed quantum-classical dynamics is capable of capturing the quantized vibrational features in photoabsorption spectra that result from infrared and Raman-active vibrational transitions. We show that vibrational and electronic absorption spectra can indeed be obtained together within a single Ehrenfest simulation. Furthermore, the electronic transitions show new sidebands that are absent in electronic dynamics simulations with fixed nuclei. Inspection of the electronic sidebands reveals that the spacing corresponds to vibrational frequencies of totally symmetric vibrational modes of the ground electronic state. A simple derivation of the time-evolving dipole in the presence of external fields and vibrational motion shows the origin of these features, demonstrating that mixed quantum-classical Ehrenfest dynamics is capable of producing infrared, Raman, and electronic absorption spectra from a single simulation.

Direct ab Initio (Meta-)Surface-Hopping Dynamics.

Tractable methods for studying the molecular dynamics of chemical processes driven by electronic nonadiabaticity are highly sought after to provide insight into, for example, photochemical reaction mechanisms, molecular collisions, and thermalized electronic band structures. Starting from the time-dependent Schrödinger equation for a many-body system, a direct ab initio trajectory surface-hopping (TSH) method relying on an analytical treatment of nonadiabatic couplings between electronic states is developed in this work. An approach that combines time-dependent perturbation theory and explicit time evolution via TSH to expedite calculation of nonadiabatic transition rates, namely, meta-surface-hopping dynamics, is presented, and an extrapolatory approach using time-dependent perturbation theory for recovering unbiased transition rates is assessed. The meta-surface-hopping method is applied to the problem of estimating nonradiative relaxation rates of a photoexcited iminium ion, CH2NH2⁺, and evidence for internal consistency of the combined dynamics/perturbation theory approach is presented.

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.