Nonadiabatic Simulation of Exciton Dynamics in Organic Semiconductors Using Neural Network-Based Frenkel Hamiltonian and Gradients

In this study, we present a multiscale method to simulate the propagation of Frenkel singlet excitons in organic semiconductors (OSCs). The approach uses neural network models to train a Frenkel-type Hamiltonian and its gradient, obtained by the long-range correction version of density functional tight-binding with self-consistent charges. Our models accurately predict site energies, excitonic couplings, and corresponding gradients, essential for the nonadiabatic molecular dynamics simulations. Combined with the fewest switches surface hopping algorithm, the method was applied to four representative OSCs: anthracene, pentacene, perylenediimide, and diindenoperylene. The simulated exciton diffusion constants align well with experimental and reported theoretical values and offer valuable insights into exciton dynamics in OSCs.

On the description of conical intersections between excited electronic states with LR-TDDFT and ADC(2)

Conical intersections constitute the conceptual bedrock of our working understanding of ultrafast, nonadiabatic processes within photochemistry (and photophysics). Accurate calculation of potential energy surfaces within the vicinity of conical intersections, however, still poses a serious challenge to many popular electronic structure methods. Multiple works have reported on the deficiency of methods like linear-response time-dependent density functional theory within the adiabatic approximation (AA LR-TDDFT) or algebraic diagrammatic construction to second-order [ADC(2)]-approaches often used in excited-state molecular dynamics simulations-to describe conical intersections between the ground and excited electronic states. In the present study, we focus our attention on conical intersections between excited electronic states and probe the ability of AA LR-TDDFT and ADC(2) to describe their topology and topography, using protonated formaldimine and pyrazine as two exemplar molecules. We also take the opportunity to revisit the performance of these methods in describing conical intersections involving the ground electronic state in protonated formaldimine-highlighting in particular how the intersection ring exhibited by AA LR-TDDFT can be perceived either as a (near-to-linear) seam of intersection or two interpenetrating cones, depending on the magnitude of molecular distortions within the branching space.

X-ray transient absorption reveals the 1Au (nπ*) state of pyrazine in electronic relaxation

Electronic relaxation in organic chromophores often proceeds via states not directly accessible by photoexcitation. We report on the photoinduced dynamics of pyrazine that involves such states, excited by a 267 nm laser and probed with X-ray transient absorption spectroscopy in a table-top setup. In addition to the previously characterized ¹B_2u (ππ*) (S₂) and ¹B_3u (nπ*) (S₁) states, the participation of the optically dark ¹A_u (nπ*) state is assigned by a combination of experimental X-ray core-to-valence spectroscopy, electronic structure calculations, nonadiabatic dynamics simulations, and X-ray spectral computations. Despite ¹A_u (nπ*) and ¹B_3u (nπ*) states having similar energies at relaxed geometry, their X-ray absorption spectra differ largely in transition energy and oscillator strength. The ¹A_u (nπ*) state is populated in 200 ± 50 femtoseconds after electronic excitation and plays a key role in the relaxation of pyrazine to the ground state.

Identification of an ultrafast internal conversion pathway of pyrazine by time-resolved vacuum ultraviolet photoelectron spectrum simulations

The internal conversion from the optically bright S₂ (¹B_2u, ππ*) state to the dark S₁ (¹B_3u, nπ*) state in pyrazine is a standard benchmark for experimental and theoretical studies on ultrafast radiationless decay. Since 2008, a few theoretical groups have suggested significant contributions of other dark states S₃ (¹A_u, nπ*) and S₄ (¹B_2g, nπ*) to the decay of S₂. We have previously reported the results of nuclear wave packet simulations [Kanno et al., Phys. Chem. Chem. Phys. 17, 2012 (2015)] and photoelectron spectrum calculations [Mignolet et al., Chem. Phys. 515, 704 (2018)] that support the conventional two-state picture. In this article, the two different approaches, i.e., wave packet simulation and photoelectron spectrum calculation, are combined: We computed the time-resolved vacuum ultraviolet photoelectron spectrum and photoelectron angular distribution for the ionization of the wave packet transferred from S₂ to S₁. The present results reproduce almost all the characteristic features of the corresponding experimental time-resolved spectrum [Horio et al., J. Chem. Phys. 145, 044306 (2016)], such as a rapid change from a three-band to two-band structure. This further supports the existence and character of the widely accepted pathway (S₂ → S₁) of ultrafast internal conversion in pyrazine.

Coupling electrons and vibrations in molecular quantum chemistry

We derive an electron-vibration model Hamiltonian in a quantum chemical framework and explore the extent to which such a Hamiltonian can capture key effects of nonadiabatic dynamics. The model Hamiltonian is a simple two-body operator, and we make preliminary steps at applying standard quantum chemical methods to evaluate its properties, including mean-field theory, linear response, and a primitive correlated model. The Hamiltonian can be compared to standard vibronic Hamiltonians, but it is constructed without reference to potential energy surfaces through direct differentiation of the one- and two-electron integrals at a single reference geometry. The nature of the model Hamiltonian in the harmonic and linear-coupling regime is investigated for pyrazine, where a simple time-dependent calculation including electron-vibration correlation is demonstrated to exhibit the well-studied population transfer between the S₂ and S₁ excited states.

Modeling irreversible molecular internal conversion using the time-dependent variational approach with sD2 ansatz

A model of irreversible molecular internal conversion dynamics due to molecular thermal energy dissipation to the bath is presented.

Time-resolved near-edge X-ray absorption fine structure of pyrazine from electronic structure and nuclear wave packet dynamics simulations

As a demonstration of the analysis of the electronic structure and the nuclear dynamics from time-resolved near-edge X-ray absorption fine structure (TR-NEXAFS), we present the TR-NEXAFS spectra of pyrazine following the excitation to the ¹B_2u(ππ^*) state. The spectra are calculated combining the frozen-core/core-valence separated equation-of-motion coupled cluster singles and doubles approach for the spectral signatures and the multiconfiguration time-dependent Hartree method for the wave packet propagation. The population decay from the ¹B_2u(ππ^*) state to the ¹B_3u(nπ^*) and ¹A_u(nπ^*) states, followed by oscillatory flow of population between the ¹B_3u(nπ^*) and ¹A_u(nπ^*) states, is interpreted by means of visualization of the potential energy curves and the reduced nuclear densities. By examining the electronic structure of the three valence-excited states and the final core-excited states, we observe that the population dynamics is explicitly reflected in the TR-NEXAFS spectra, especially when the heteroatoms are selected as the X-ray absorption sites. This work illustrates the feasibility of extracting fine details of molecular photophysical processes from TR-NEXAFS spectra by using currently available theoretical methods.

Mapping of Wave Packet Dynamics at Conical Intersections by Time- and Frequency-Resolved Fluorescence Spectroscopy: A Computational Study

Monitoring of wave packet dynamics at conical intersections by time- and frequency-resolved fluorescence spectroscopy has been investigated theoretically for a three-state two-mode model of a conical intersection coupled to a dissipative environment. The ideal and the actually measurable time- and frequency-gated fluorescence spectra are accurately and efficiently simulated by combining the hierarchy equations-of-motion method for dissipative quantum dynamics with the methodology of the equation-of-motion phase-matching approach for the calculation of spectroscopic signals. It is shown that time- and frequency-resolved fluorescence spectra reveal essential aspects of the wave packet dynamics at conical intersections and the effects of environment-induced dissipation. The results of the present work indicate that fluorescence up-conversion spectroscopy with femtosecond time resolution is an efficient tool for the characterization of ultrafast dynamics at conical intersections.

Imaging large amplitude out-of-plane motion in photoexcited pentafluorobenzene using time-resolved photoelectron spectroscopy: a computational study

Excited-state dynamics of pentafluorobenzene is studied in detail for a quartic vibronic coupling model including the six b1 vibrational modes of the molecule and the two lowest excited electronic states. The study analyzes the influence of the large-amplitude out-of-plane vibrational motion on the electronic dynamics and extends to the simulation of the emerging time-resolved photoelectron spectra. The mapping of coherent non-separable electron-nuclear motion into oscillatory photoelectron signals is discussed.

Probing the Interstate Coupling near a Conical Intersection by Optical Spectroscopy

Conical intersections are points where adiabatic potential energy surfaces cross. The interstate coupling between the potential energy surfaces plays a crucial role in many processes associated with conical intersections. Still no method exists to measure this coupling driving the chemical reactions between the potential energy surfaces involved. In this Letter, using a generic model for photoisomerization, we propose a novel experimental approach to estimate the coupling that mixes the electronic states near a conical intersection. The approach is based on analyzing the vibrational wavepacket of the reactant in the adiabatic ground and excited electronic states. The nuclear wavepacket dynamics are extracted from linear absorption and two-dimensional electronic spectroscopy. Comparing the frequencies of the coupling mode in the adiabatic ground and excited states from models with and without coupling between the potential energy surfaces suggests an experimental tool to determine the interstate coupling.