Yinan Shu

Configuration interaction singles natural orbitals: An orbital basis for an efficient and size intensive multireference description of electronic excited states

Multireference quantum chemical methods, such as the complete active space self-consistent field (CASSCF) method, have long been the state of the art for computing regions of potential energy surfaces (PESs) where complex, multiconfigurational wavefunctions are required, such as near conical intersections. Herein, we present a computationally efficient alternative to the widely used CASSCF method based on a complete active space configuration interaction (CASCI) expansion built from the state-averaged natural orbitals of configuration interaction singles calculations (CISNOs). This CISNO-CASCI approach is shown to predict vertical excitation energies of molecules with closed-shell ground states similar to those predicted by state averaged (SA)-CASSCF in many cases and to provide an excellent reference for a perturbative treatment of dynamic electron correlation. Absolute energies computed at the CISNO-CASCI level are found to be variationally superior, on average, to other CASCI methods. Unlike SA-CASSCF, CISNO-CASCI provides vertical excitation energies which are both size intensive and size consistent, thus suggesting that CISNO-CASCI would be preferable to SA-CASSCF for the study of systems with multiple excitable centers. The fact that SA-CASSCF and some other CASCI methods do not provide a size intensive/consistent description of excited states is attributed to changes in the orbitals that occur upon introduction of non-interacting subsystems. Finally, CISNO-CASCI is found to provide a suitable description of the PES surrounding a biradicaloid conical intersection in ethylene.

Defect-Induced Conical Intersections Promote Nonradiative Recombination.

We apply multireference electronic structure calculations to demonstrate the presence of conical intersections between the ground and the first excited electronic states of three silicon nanocrystals containing defects characteristic of the oxidized silicon surface. These intersections are accessible upon excitation at visible wavelengths and are predicted to facilitate nonradiative recombination with a rate that increases with decreasing particle size. This work illustrates a new framework for identifying defects responsible for nonradiative recombination.

Simulated evolution of fluorophores for light emitting diodes

Organic light emitting diodes based on fluorophores with a propensity for thermally activated delayed fluorescence (TADF) are able to circumvent limitations imposed on device efficiency by spin statistics. Molecules with a propensity for TADF necessarily have two properties: a small gap between the lowest lying singlet and triplet excited states and a large transition dipole moment for fluorescence. In this work, we demonstrate the use of a genetic algorithm to search a region of chemical space for molecules with these properties. This algorithm is based on a flexible and intuitive representation of the molecule as a tree data structure, in which the nodes correspond to molecular fragments. Our implementation takes advantage of hybrid parallel graphics processing unit accelerated computer clusters to allow efficient sampling while retaining a reasonably accurate description of the electronic structure (in this case, CAM-B3LYP/6-31G(∗∗)). In total, we have identified 3792 promising candidate fluorophores from a chemical space containing 1.26 × 10(6) molecules. This required performing electronic structure calculations on only 7518 molecules, a small fraction of the full space. Several novel classes of molecules which show promise as fluorophores are presented.

Communication: Non-radiative recombination via conical intersection at a semiconductor defect

Localization of electronic excitations at molecule-sized semiconductor defects often precedes non-radiative (NR) decay, and it is known that many molecules undergo NR decay via conical intersection. Herein, we report the direct simulation of fast and efficient NR decay via a conical intersection at a known semiconductor defect. It is suggested that this silicon epoxide defect may selectively quench photoluminescence (PL) in small silicon nanocrystals (band gap > ~2.8 eV), and thus influence both the observed PL yield and PL energy of oxidized silicon nanocrystals.

Reducing the propensity for unphysical wavefunction symmetry breaking in multireference calculations of the excited states of semiconductor clusters

Unphysical spatial symmetry breaking in multiconfigurational self-consistent field calculations can lead to undesirable artifacts in the potential energy surfaces and electronic properties of molecules. Herein, we report several examples of such symmetry breaking in calculations of the excited states of small semiconductor clusters and related molecules at the state-averaged complete active space self-consistent field (SA-CASSCF) level of theory. A multireference approach is proposed to reduce its incidence: the singly excited active space complete active space configuration interaction (SEAS-CASCI) method. In SEAS-CASCI, the orbitals are determined by variationally minimizing an energy expression that does not depend on the off-diagonal Hamiltonian matrix elements which drive symmetry breaking at the SA-CASSCF level of theory. By application to several highly symmetric molecules, SEAS-CASCI is demonstrated to reduce the propensity for unphysical spatial symmetry breaking and eliminate resulting errors in the potential energy surfaces and molecular properties relative to the SA-CASSCF description. The SEAS method is also found to eliminate unphysical wavefunction distortion in asymmetric molecules. Finally, SEAS-CASCI is demonstrated to accurately describe the biradicaloid region of the potential energy surface of ethylene.

Dual-Functional Tamm-Dancoff Approximation: A Convenient Density Functional Method that Correctly Describes S1/S0 Conical Intersections

Time-dependent Kohn-Sham density functional theory has been used successfully to compute vertical excitation energies, especially for large molecular systems. However, the lack of double excitation character in the excited amplitudes produced by linear response in the adiabatic approximation holds it back from broader applications in photochemistry; for example, it shows (3N - 7)-dimensional conical intersection seams (where N is the number of atoms) between ground and excited states, although the correct dimensionality is 3N - 8. In this letter, we present a new, conceptually simple, easy-to-implement, and easy-to-use way to employ time-dependent Kohn-Sham density functional theory that has global accuracy comparable with the conventional single-functional version and that recovers the double cone topology of the potential energy surfaces at S1/S0 conical intersection seams. The new method is called the dual-functional Tamm-Dancoff approximation (DF-TDA).

Understanding Nonradiative Recombination through Defect-Induced Conical Intersections

Defects are known to introduce pathways for the nonradiative recombination of electronic excitations in semiconductors, but implicating a specific defect as a nonradiative center remains challenging for both experiment and theory. In this Perspective, we present recent progress toward this goal involving the identification and characterization of defect-induced conical intersections (DICIs), points of degeneracy between the ground and first excited electronic states of semiconductor materials that arise from the deformation of specific defects. Analysis of DICIs does not require the assumption of weak correlation between the electron and hole nor of stationary nuclei. It is demonstrated that in some cases an energetically accessible DICI is present even when no midgap state is predicted by single-particle theories (e.g., density functional theory). We review recent theoretical and computational developments that enable the location of DICIs in semiconductor nanomaterials and present insights into the photoluminescence of silicon nanocrystals gleaned from DICIs.

Complete active space configuration interaction from state-averaged configuration interaction singles natural orbitals: Analytic first derivatives and derivative coupling vectors

A new complete active space configuration interaction (CASCI) method was recently introduced that uses state-averaged natural orbitals from the configuration interaction singles method (configuration interaction singles natural orbital CASCI, CISNO-CASCI). This method has been shown to perform as well or better than state-averaged complete active space self-consistent field for a variety of systems. However, further development and testing of this method have been limited by the lack of available analytic first derivatives of the CISNO-CASCI energy as well as the derivative coupling between electronic states. In the present work, we present a Lagrangian-based formulation of these derivatives as well as a highly efficient implementation of the resulting equations accelerated with graphical processing units. We demonstrate that the CISNO-CASCI method is practical for dynamical simulations of photochemical processes in molecular systems containing hundreds of atoms.

Dynamics of recombination via conical intersection in a semiconductor nanocrystal

The ultrafast dynamics of nonradiative recombination at dangling bond defects is elucidated by nanoscale multireference ab initio molecular dynamics simulations.

Dual-Functional Tamm-Dancoff Approximation with Self-Interaction-Free Orbitals: Vertical Excitation Energies and Potential Energy Surfaces Near an Intersection Seam

Recently we have developed the dual-functional Tamm-Dancoff approximation (DF-TDA) method. DF-TDA is an alternative to linear-response time-dependent density functional theory (LR-TDDFT) with the advantage of providing a correct double-cone topology of S1/S0 conical intersections. In the DF-TDA method, we employ different functionals, which are denoted G and F, for orbital optimization and Hamiltonian construction. We use the notation DF-TDA/G:F. In the current work, we propose that G be the same as F except for having 100% Hartree-Fock exchange. We use the notation F100 to denote functional F with this modification. A motivation for this is that functionals with 100% Hartree-Fock exchange are one-electron self-interaction-free. Here we validate the use of F100/M06 to compute vertical excitation energies and the global potential energy surface of ammonia near a conical intersection to further validate the F100 method for photochemical problems.