Analysis and visualization of energy densities. II. Insights from linear-response time-dependent density functional theory calculations

Multistate Energy Decomposition Analysis of Molecular Excited States

A multistate energy decomposition analysis (MS-EDA) method is described to dissect the energy components in molecular complexes in excited states. In MS-EDA, the total binding energy of an excimer or an exciplex is partitioned into a ground-state term, called local interaction energy, and excited-state contributions that include exciton excitation energy, superexchange stabilization, and orbital and configuration-state delocalization. An important feature of MS-EDA is that key intermediate states associated with different energy terms can be variationally optimized, providing quantitative insights into widely used physical concepts such as exciton delocalization and superexchange charge-transfer effects in excited states. By introducing structure-weighted adiabatic excitation energy as the minimum photoexcitation energy needed to produce an excited-state complex, the binding energy of an exciplex and excimer can be defined. On the basis of the nature of intermolecular forces through MS-EDA analysis, it was found that molecular complexes in the excited states can be classified into three main categories, including (1) encounter excited-state complex, (2) charge-transfer exciplex, and (3) intimate excimer or exciplex. The illustrative examples in this Perspective highlight the interplay of local excitation polarization, exciton resonance, and superexchange effects in molecular excited states. It is hoped that MS-EDA can be a useful tool for understanding photochemical and photobiological processes.

Divide-and-Conquer Linear-Scaling Quantum Chemical Computations

Fragmentation and embedding schemes are of great importance when applying quantum-chemical calculations to more complex and attractive targets. The divide-and-conquer (DC)-based quantum-chemical model is a fragmentation scheme that can be connected to embedding schemes. This feature article explains several DC-based schemes developed by the authors over the last two decades, which was inspired by the pioneering study of DC self-consistent field (SCF) method by Yang and Lee (J. Chem. Phys. 1995, 103, 5674-5678). First, the theoretical aspects of the DC-based SCF, electron correlation, excited-state, and nuclear orbital methods are described, followed by the two-component relativistic theory, quantum-mechanical molecular dynamics simulation, and the introduction of three programs, including DC-based schemes. Illustrative applications confirmed the accuracy and feasibility of the DC-based schemes.

Analyzing Excitation-Energy Transfer Based on the Time-Dependent Density Functional Theory in Real Time

Excitation-energy transfer is a key step in processes such as photosynthesis that convert light into other forms of energy. Time-dependent density functional theory (DFT) in real time is ideal for the first-principles simulation of such processes due to its computational efficiency. We here demonstrate how real-time DFT can be used for analyzing excitation-energy transfer from first-principles. We discuss several measures of energy transfer that are based solely on the time-dependent density, are well founded in the DFT framework, allow for intuitive understanding and visualization, and reproduce important limiting cases of an analytical model. We demonstrate their usefulness in calculations for model systems, both with static nuclei and in the context of DFT-based Ehrenfest dynamics.

Controlling Aggregation-Induced Two-Photon Absorption Enhancement via Intermolecular Interactions

Historically, two-photon absorption (2PA) cross sections reported in the literature have been derived from solution-phase measurements. However, such techniques fail to grasp the implications of how these cross sections can be impacted by varying degrees of aggregation or in the condensed phase as bulk solids or thin films. For a precise determination of how aggregation impacts 2PA at a molecular level, computational methods present themselves as ideal. Herein, a series of quadrupolar π-conjugated dyes were simulated by molecular dynamics (MD) in the gas phase and condensed phase. In the condensed phase, their intermolecular interactions and electronic coupling behavior were fully characterized, both quantitatively and qualitatively. Using quadratic-response time-dependent density functional theory, 2PA cross sections of structures derived from MD trajectories were calculated. Comparisons are made between gas-phase and condensed-phase results, and enhancement factors are defined to show how certain dyes may experience changes in their respective 2PA cross sections as a function of aggregation. It was found that these cross sections depend heavily on conformational locking in the condensed phase and relative stacking arrangements. J-aggregates were associated with enhanced 2PA and H-aggregates with quenched 2PA activity. However, in a highly disordered aggregate, the effects of these stacking arrangements are averaged out of the bulk result, and the effects of conformational locking dominate.

Visualization of Electron Density Changes Along Chemical Reaction Pathways

We propose a simple procedure for visualizing the electron density changes (EDC) during a chemical reaction, which is based on a mapping of rectangular grid points for a stationary structure into (distorted) positions around atoms of another stationary structure. Specifically, during a small step along the minimum energy pathway (MEP), the displacement of each grid point is obtained as a linear combination of the motion of all atoms, with the contribution from each atom scaled by the corresponding Hirshfeld weight. For several reactions (identity SN2, Claisen rearrangement, Diels-Alder reaction, [3+2] cycloaddition, and phenylethyl mercaptan attack on pericosine A), our EDC plots showed an expected reduction of electron densities around severed bonds (or those with the bond-order lowered), with the opposite observed for newly-formed or enhanced chemical bonds. The EDC plots were also shown for copper triflate catalyzed N2O fragmentation, where the N-O bond weakening initially occurred on a singlet surface, but continued on a triplet surface after reaching the minimum-energy crossing point (MECP) between the two potential energy surfaces.

Evaluation of molecular photophysical and photochemical properties using linear response time-dependent density functional theory with classical embedding: Successes and challenges

Time-dependent density functional theory (TDDFT) based approaches have been developed in recent years to model the excited-state properties and transition processes of the molecules in the gas-phase and in a condensed medium, such as in a solution and protein microenvironment or near semiconductor and metal surfaces. In the latter case, usually, classical embedding models have been adopted to account for the molecular environmental effects, leading to the multi-scale approaches of TDDFT/polarizable continuum model (PCM) and TDDFT/molecular mechanics (MM), where a molecular system of interest is designated as the quantum mechanical region and treated with TDDFT, while the environment is usually described using either a PCM or (non-polarizable or polarizable) MM force fields. In this Perspective, we briefly review these TDDFT-related multi-scale models with a specific emphasis on the implementation of analytical energy derivatives, such as the energy gradient and Hessian, the nonadiabatic coupling, the spin–orbit coupling, and the transition dipole moment as well as their nuclear derivatives for various radiative and radiativeless transition processes among electronic states. Three variations of the TDDFT method, the Tamm–Dancoff approximation to TDDFT, spin–flip DFT, and spin-adiabatic TDDFT, are discussed. Moreover, using a model system (pyridine–Ag20 complex), we emphasize that caution is needed to properly account for system–environment interactions within the TDDFT/MM models. Specifically, one should appropriately damp the electrostatic embedding potential from MM atoms and carefully tune the van der Waals interaction potential between the system and the environment. We also highlight the lack of proper treatment of charge transfer between the quantum mechanics and MM regions as well as the need for accelerated TDDFT modelings and interpretability, which calls for new method developments.

Adsorption and catalytic degradation of preservative parabens by graphene-family nanomaterials

Parabens pose increasing threats to human health due to endocrine disruption activity. Adsorption and degradation of parabens by three types of graphene-family nanomaterials (GFNs) were therefore investigated. For a given paraben, the maximum adsorption capacities (Q⁰) followed the order of reduced graphene oxide (RGO) > multilayered graphene (MG) > graphene oxide (GO); for a given GFN, Q⁰ followed the order of butylparaben (BuP) > propylparaben (PrP) > ethylparaben (EtP) > methylparaben (MeP), dominated by hydrophobic interaction. MeP removal by all the three GFNs was highly enhanced (0.55-4.37 times) with the assistance of H₂O₂ due to additional catalytic degradation process, and MG showed the highest removal enhancement. ∙OH was confirmed as the dominant radicals responsible for parabens degradation. For MG and RGO, the metal impurities (Fe, Cu, Mn, and Co) initiated Fenton-like reaction with H₂O₂ to generate ∙OH. GO contained oxygen-centered free radicals, which were responsible for ∙OH formation via transferring electron to H₂O₂. Four degradation byproducts of MeP were identified, including oxalic, propanedioic, fumaric, and 2,5-dihydroxybenzoic acids. Combined with density function theory calculations, the degradation sites and pathways were identified and confirmed. These findings provide useful information on mechanistic understanding towards the adsorption and degradation of parabens by GFNs.

The role of excited-state character, structural relaxation, and symmetry breaking in enabling delayed fluorescence activity in push-pull chromophores

Thermally activated delayed fluorescence (TADF) is a current promising route for generating highly efficient light-emitting devices. However, the design process of new chromophores is hampered by the complicated underlying photophysics. In this work, four closely related donor-π-acceptor-π-donor systems are investigated, two of which were synthesised previously, with the aim of elucidating their varying effectiveness for TADF. We outline that the frontier orbitals are insufficient for discriminating between the molecules. Subsequently, a detailed analysis of the excited states at a correlated ab initio level highlights the presence of a number of closely spaced singlet and triplet states of varying character. Results from five density functionals are compared against this reference revealing dramatic changes in, both, excited state energies and wavefunctions following variations in the amount of Hartree-Fock exchange included. Excited-state minima are optimised in solution showing the crucial role of structural variations and symmetry breaking for producing a strongly emissive S₁ state. The adiabatic singlet-triplet gaps thus obtained depend strongly on the range separation parameter used in the hybrid density functional calculations. More generally, this work highlights intricate differences present between singlet and triplet excited state wavefunctions and the challenges in describing them accurately.

Elucidating the Electronic Structure of a Delayed Fluorescence Emitter via Orbital Interactions, Excitation Energy Components, Charge-Transfer Numbers, and Vibrational Reorganization Energies

Recently, Wang and co-workers carried out frontier molecule orbital engineering in the design of m-Cz-BNCz, a thermally activated delayed fluorescence (TADF) molecule that emits pure green light at an external quantum efficiency of 27%. To further understand the underlying molecular design principles, we employed four advanced electronic structure analysis tools. First, an absolutely localized molecular orbitals (ALMO-) based analysis indicates an antibonding combination between the highest occupied molecular orbitals (HOMOs) of the donor 3,6-di-tert-butylcarbazole fragment and the acceptor BNCz fragment, which raises the HOMO energy and red-shifts the fluorescence emission wavelength. Second, excitation energy component analysis reveals that the S₁-T₁ gap is dominated by two-electron components of the excitation energies. Third, charge transfer number analysis, which is extended to use fragment-based Hirshfeld weights, indicates that the S₁ and T₁ excited states of m-Cz-BNCz (within time-dependent density functional theory) have notable charge transfer characters (27% for S₁ and 12% for T₁). This provides a balance between a small single-triplet gap and a substantial fluorescence intensity. Last, a vibrational reorganization energy analysis pinpoints the torsional motion between the BNCz and Cz moieties of m-Cz-BNCz as the source for its wider emission peak than that of p-Cz-BNCz. These four types of analyses are expected to be very valuable in the study and design of other TADF and functional dye molecules.

Analysis and visualization of energy densities. I. Insights from real-time time-dependent density functional theory simulations

In this article, we report a scheme to analyze and visualize the energy density fluctuations during the real-time time-dependent density functional theory (RT-TDDFT) simulations. Using Ag4-N2 complexes as examples, it is shown that the grid-based Kohn-Sham energy density can be computed at each time step using a procedure from Nakai and coworkers. Then the instantaneous energy of each molecular fragment (such as Ag4 and N2) can be obtained by partitioning the Kohn-Sham energy densities using Becke or fragment-based Hirshfeld (FBH) scheme. A strong orientation-dependence is observed for the energy flow between the Ag4 cluster and a nearby N2 molecule in the RT-TDDFT simulations. Future applications of such an energy density analysis in electron dynamics simulations are discussed.