Description of Sudden Polarization in the Excited Electronic States with an Ensemble Density Functional Theory Method

Structure of Multi-State Correlation in Electronic Systems

Beyond the Hohenberg-Kohn density functional theory for the ground state, it has been established that the Hamiltonian matrix for a finite number (N) of lowest eigenstates is a matrix density functional. Its fundamental variable─the matrix density D(r)─can be represented by, or mapped to, a set of auxiliary, multiconfigurational wave functions expressed as a linear combination of no more than N2 determinant configurations. The latter defines a minimal active space (MAS), which naturally leads to the introduction of the correlation matrix functional, responsible for the electronic correlation effects outside the MAS. In this study, we report a set of rigorous conditions in the Hamiltonian matrix functional, derived by enforcing the symmetry of a Hilbert subspace, namely the subspace invariance property. We further establish a fundamental theorem on the correlation matrix functional. That is, given the correlation functional for a single state in the N-dimensional subspace, all elements of the correlation matrix functional for the entire subspace are uniquely determined. These findings reveal the intricate structure of electronic correlation within the Hilbert subspace of lowest eigenstates and suggest a promising direction for efficient simulation of excited states.

Unraveling the effect of aromaticity for the dynamics of excited states of single benzene fluorophores

The photophysical properties of a series of recently synthesized single benzene fluorophores were investigated using ensemble density functional theory calculations. The energetic stability of the ground and excited state species were counterposed against the aromaticity index derived from local vibrational modes. It was found that the large Stokes shift of the fluorophores (up to ca. 5800 cm  - 1 ) originates from the effect of electron donating and electron withdrawing substituents rather than π -delocalization and related (anti-)aromaticity. On the basis of nonadiabatic molecular dynamics simulations, the absence of fluorescence from one of the regioisomers was explained by the occurrence of easily accessible S  1 /S  0 conical intersections below the vertical excitation energy level. It is demonstrated in the manuscript that the analysis of local mode force constants and the related aromaticity index represent a useful tool for the characterization of π -delocalization effects in π -conjugated compounds.

Formulation of transition dipole gradients for non-adiabatic dynamics with polaritonic states

A general formulation of the strong coupling between photons confined in a cavity and molecular electronic states is developed for the state-interaction state-average spin-restricted ensemble-referenced Kohn-Sham method. The light-matter interaction is included in the Jaynes-Cummings model, which requires the derivation and implementation of the analytical derivatives of the transition dipole moments between the molecular electronic states. The developed formalism is tested in the simulations of the nonadiabatic dynamics in the polaritonic states resulting from the strong coupling between the cavity photon mode and the ground and excited states of the penta-2,4-dieniminium cation, also known as PSB3. Comparison with the field-free simulations of the excited-state decay dynamics in PSB3 reveals that the light-matter coupling can considerably alter the decay dynamics by increasing the excited state lifetime and hindering photochemically induced torsion about the C=C double bonds of PSB3. The necessity of obtaining analytical transition dipole gradients for the accurate propagation of the dynamics is underlined.

Assessment of the Electron Correlation Treatment on the Quantum-Classical Dynamics of Retinal Protonated Schiff Base Models: XMS-CASPT2, RMS-CASPT2, and REKS Methods

We compare the performance of three different multiconfigurational wave function-based electronic structure methods and two implementations of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method. The study is characterized by three features: (i) it uses a small set of quantum-classical trajectories rather than potential energy surface mapping, (ii) it focuses, exclusively, on the photoisomerization of retinal protonated Schiff base models, and (iii) it probes the effect of both methyl substitution and the increase in length of the conjugate π-system. For each tested method, the corresponding analytical gradients are used to drive the quantum-classical (Tully's FSSH method) trajectory propagation, including the recent multistate XMS-CASPT2 and RMS-CASPT2 gradients. It is shown that while CASSCF, XMS-CASPT2, and RMS-CASPT2 yield consistent photoisomerization dynamics descriptions, REKS produces, in some of these systems, qualitatively different behavior that is attributed to a flatter and topographically different excited state potential energy surface. The origin of this behavior can be traced back to the effect of the employed density functional approximation. The above studies are further expanded by benchmarking, at the CASSCF and REKS levels, the electronic structure methods using a QM/MM model of the visual pigment rhodopsin.

Dark State of the Thiele Hydrocarbon: Efficient Solvatochromic Emission from a Nonpolar Centrosymmetric Singlet Diradicaloid

In this work, a comprehensive investigation of the photoinduced processes and mechanisms linked to the luminescence of a novel nonperchlorinated Thiele hydrocarbon (TTH) is presented. Despite the comparable diradical character of TTH (y0 = 0.32-0.44) and the unsubstituted Thiele hydrocarbon (TH) (y0 = 0.30), the polyhalogenated species is inert and photostable, showing an intense deep-red/near-infrared (NIR) fluorescence (photoluminescence quantum yield (PLQY) = 0.84 in toluene) even at room temperature and in the solid state (PLQY = 0.19). TTH displays a large Stokes shift (307 nm in benzonitrile) and solvatochromic behavior, which is unusual for a centrosymmetric, nonpolar, and low-conjugated species. These outstanding emission features are interpreted through quantum-chemical calculations, indicating that its fluorescence arises from the low-lying dark doubly excited zwitterionic state, typically found at low excitation energies in diradicaloids, acquiring dipole moment and intensity by state mixing via twisting around the strongly elongated exocyclic CC bonds of the excited p-quinodimethane (pQDM) core, with a mechanism similar to sudden polarization occurring in olefins. Such a mechanism is derived from ns and fs transient absorption measurements.

DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

Ensemble Density Functional Theory of Neutral and Charged Excitations : Exact Formulations, Standard Approximations, and Open Questions

Recent progress in the field of (time-independent) ensemble density-functional theory (DFT) for excited states are reviewed. Both Gross-Oliveira-Kohn (GOK) and N-centered ensemble formalisms, which are mathematically very similar and allow for an in-principle-exact description of neutral and charged electronic excitations, respectively, are discussed. Key exact results, for example, the equivalence between the infamous derivative discontinuity problem and the description of weight dependencies in the ensemble exchange-correlation density functional, are highlighted. The variational evaluation of orbital-dependent ensemble Hartree-exchange (Hx) energies is discussed in detail. We show in passing that state-averaging individual exact Hx energies can lead to severe (although solvable) v-representability issues. Finally, we explore the possibility of using the concept of density-driven correlation, which has been introduced recently and does not exist in regular ground-state DFT, for improving state-of-the-art correlation density-functional approximations for ensembles. The present review reflects the efforts of a growing community to turn ensemble DFT into a rigorous and reliable low-cost computational method for excited states. We hope that, in the near future, this contribution will stimulate new formal and practical developments in the field.