Exploring the effects of quantum decoherence on the excited-state dynamics of molecular systems

Fewest-switches surface hopping (FSSH) is employed in order to investigate the nonadiabatic excited-state dynamics of thiophene and related compounds and hence to establish a connection between the electronic system, the critical points in configuration space and the deactivation dynamics. The potential-energy surfaces of the studied molecules were calculated with complete active space self-consistent field and time-dependent density-functional theory. They are analyzed thoroughly to locate and optimize minimum-energy conical intersections, which are essential to the dynamics of the system. The influence of decoherence on the dynamics is examined by employing different decoherence schemes. We find that irrespective of the employed decoherence algorithm, the population dynamics of thiophene give results which are sound with the expectations grounded on the analysis of the potential-energy surface. A more detailed look at single trajectories as well as on the excited-state lifetimes, however, reveals a substantial dependence on how decoherence is accounted for. In order to connect these findings, we describe how ensemble averaging cures some of the overcoherence problems of uncorrected FSSH. Eventually, we identify carbon–sulfur bond cleavage as a common feature accompanying electronic transitions between different states in the simulations of all thiophene-related compounds studied in this work, which is of interest due to their relevance in organic photovoltaics.

Control of nuclear dynamics in the benzene cation by electronic wavepacket composition

The study of coupled electron-nuclear dynamics driven by coherent superpositions of electronic states is now possible in attosecond science experiments. The objective is to understand the electronic control of chemical reactivity. In this work we report coherent 8-state non-adiabatic electron-nuclear dynamics simulations of the benzene radical cation. The computations were inspired by the extreme ultraviolet (XUV) experimental results in which all 8 electronic states were prepared with significant population. Our objective was to study the nuclear dynamics using various bespoke coherent electronic state superpositions as initial conditions in the Quantum-Ehrenfest method. The original XUV measurements were supported by Multi-configuration time-dependent Hartree (MCTDH) simulations, which suggested a model of successive passage through conical intersections. The present computations support a complementary model where non-adiabatic events are seen far from a conical intersection and are controlled by electron dynamics involving non-adjacent adiabatic states. It proves to be possible to identify two superpositions that can be linked with two possible fragmentation paths.

Strong Field Ionization-Photofragmentation on Ultrafast Evolution of Electronic States of Toluene Cations

Ultrafast time-resolved strong field ionization-photofragmentation (SFI-PF) has emerged as a useful method for investigation of dynamics of molecular cations. Here we perform a SFI-PF study on the electronic states of toluene cations. By measuring the ion yields as a function of delay time, we obtain the transients of both parent and daughter ions, which show ultrafast decays and out-of-phase oscillations. The results provide the first experimental evidence of D₁-D₀ ultrafast relaxation of toluene cations occurring in about 530 fs and indicate coincident resonance between the vibrational states in D₁ and D₀ leading to oscillations with a period of about 2.05 ps. Our study should shed some light on the ultrafast photochemistry involving complex molecular cations.

CAS without SCF-Why to use CASCI and where to get the orbitals

The complete active space self-consistent field (CASSCF) method has seen broad adoption due to its ability to describe the electronic structure of both the ground and excited states of molecules over a broader swath of the potential energy surface than is possible with the simpler Hartree-Fock approximation. However, it also has a reputation for being unwieldy, computationally costly, and un-black-box. Here, we discuss a class of alternatives, complete active space configuration interaction (CASCI) methods, paying particular attention to their application to electronic excited states. The goal of this Perspective is fourfold. First, we argue that CASCI is not merely an approximation to CASSCF, in that it can be designed to have important qualitative advantages over CASSCF. Second, we present several insights drawn from our experience experimenting with different schemes for computing orbitals to be employed in CASCI. Third, we argue that CASCI is well suited for application to nanomaterials. Finally, we reason that, with the rise in new low-scaling approaches for describing multireference systems, there is a greater need than ever to develop new methods for defining orbitals that provide an efficient and accurate description of both static correlation and electronic excitations in a limited active space.

Waveform control of molecular dynamics close to a conical intersection

Conical intersections are ubiquitous in chemical systems but, nevertheless, extraordinary points on the molecular potential energy landscape. They provide ultra-fast radiationless relaxation channels, their topography influences the product branching, and they equalize the timescales of the electron and nuclear dynamics. These properties reveal optical control possibilities in the few femtosecond regime. In this theoretical study, we aim to explore control options that rely on the carrier envelope phase of a few-cycle IR pulse. The laser interaction creates an electronic superposition just before the wave packet reaches the conical intersection. The imprinted phase information is varied by the carrier envelope phase to influence the branching ratio after the conical intersection. We test and analyze this scenario in detail for a model system and show to what extent it is possible to transfer this type of control to a realistic system like uracil.

Analytic gradients and derivative couplings for configuration interaction with all single excitations and one double excitation-En route to nonadiabatic dynamics

We present analytic gradients and derivative couplings for the simplest possible multireference configuration interaction method, CIS-1D, an electronic structure Ansatz that includes all single excitations and one lone double excitation on top of a Hartree-Fock reference state. We show that the resulting equations are numerically stable and require the evaluation of a similar number of integrals as compared to standard CIS theory; one can easily differentiate the required frontier orbitals (h and l) with minimal cost. The resulting algorithm has been implemented within the Q-Chem electronic structure package and should be immediately useful for understanding photochemistry with S₀-S₁ crossings.

Role of Multistate Intersections in Photochemistry

It has been generally accepted that the intersection of potential energy surfaces can facilitate nonadiabatic transitions and plays a crucial role in photochemistry. Although most previous studies have focused on the conical intersection of two electronic states, multistate intersections are common in polyatomic molecules, and their key roles in photochemistry have been uncovered by electronic structure calculations and nonadiabatic dynamics simulations. In this Perspective, the algorithms for searching two- or three-state intersections are first examined with an emphasis on the latest development in a general algorithm for location of multistate intersections. Then, we focus on intersystem crossing (ISC) that occurs in the region of multistate intersection, paying more attention to how the state-specific spin-orbit coupling interaction influences nonadiabatic ISC processes. Finally, the interweaving of nonadiabatic dynamics simulation and electronic structure calculation has been recognized as a correct way to ascertain the vital roles of multistate intersections in photochemical reactions.

A systematic model study quantifying how conical intersection topography modulates photochemical reactions

Despite their important role in photochemistry and expected presence in most polyatomic molecules, conical intersections have been thoroughly characterized in a comparatively small number of systems. Conical intersections can confer molecular photoreactivity or photostability, often with remarkable efficacy, due to their unique structure: at a conical intersection, the adiabatic potential energy surfaces of two or more electronic states are degenerate, enabling ultrafast decay from an excited state without radiative emission, known as nonadiabatic transfer. Furthermore, the precise conical intersection topography determines fundamental properties of photochemical processes, including excited-state decay rate, efficacy, and molecular products that are formed. However, these relationships have yet to be defined comprehensively. In this article, we use an adaptable computational model to investigate a variety of conical intersection topographies, simulate resulting nonadiabatic dynamics, and calculate key photochemical observables. We varied the vibrational mode frequencies to modify conical intersection topography systematically in four primary classes of conical intersections and quantified the resulting rate, total yield, and product yield of nonadiabatic decay. The results reveal that higher vibrational mode frequencies reduce nonadiabatic transfer, but increase the transfer rate and resulting photoproduct formation. These trends can inform progress toward experimental control of photochemical reactions or tuning of molecules' photochemical properties based on conical intersections and their topography.

Effects of conical intersections on hyperfine quenching of hydroxyl OH in collision with an ultracold Sr atom

The effect of conical intersections (CIs) on electronic relaxation, transitions from excited states to ground states, is well studied, but their influence on hyperfine quenching in a reactant molecule is not known. Here, we report on ultracold collision dynamics of the hydroxyl free-radical OH with Sr atoms leading to quenching of OH hyperfine states. Our quantum-mechanical calculations of this process reveal that quenching is efficient due to anomalous molecular dynamics in the vicinity of the conical intersection at collinear geometry. We observe wide scattering resonance features in both elastic and inelastic rate coefficients at collision energies below \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k_{\text {B}}\times 10 \, \hbox {mK}$$\end{document}kB×10mK. They are identified as either p- or d-wave shape resonances. We also describe the electronic potentials relevant for these non-reactive collisions, their diabatization procedure, as well as the non-adiabatic coupling between the diabatic potentials near the CIs.

On the Intrinsically Low Quantum Yields of Pyrimidine DNA Photodamages: Evaluating the Reactivity of the Corresponding Minimum Energy Crossing Points

The low quantum yield of photoformation of cyclobutane pyrimidine dimers and pyrimidine-pyrimidone (6-4) adducts in DNA bases is usually associated with the presence of more favorable nonreactive decay paths and with the unlikeliness of exciting the system in a favorable conformation. Here, we prove that the ability of the reactive conical intersection to bring the system either back to the absorbing conformation or to the photoproduct must be considered as a fundamental factor in the low quantum yields of the mentioned photodamage. In support of the proposed model, the one order of magnitude difference in the quantum yield of formation of the cyclobutane thymine dimer with respect to the thymine-thymine (6-4) adduct is rationalized here by comparing the reactive ability of the seam of intersections leading respectively to the cyclobutane thymine dimer and the oxetane precursor of the thymine-thymine (6-4) adduct at the CASPT2 level of theory.

Multireference Electron Correlation Methods: Journeys along Potential Energy Surfaces

Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way and, therefore, can yield accurate and predictive results even when single-reference methods or multiconfigurational self-consistent field theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces. This includes the optimization of molecular geometries, such as equilibrium geometries and conical intersections and on-the-fly photodynamics simulations, both of which depend heavily on the ability of the method to properly explore the potential energy surface. Because such applications require nuclear gradients and derivative couplings, the availability of analytical nuclear gradients greatly enhances the scope of quantum chemical methods. This review focuses on the developments and advances made in the past two decades. A detailed account of the analytical nuclear gradient and derivative coupling theories is presented. Emphasis is given to the software infrastructure that allows one to make use of these methods. Notable applications of multireference electron correlation methods to chemistry, including geometry optimizations and on-the-fly dynamics, are summarized at the end followed by a discussion of future prospects.

Insights into the Mechanism of Nonadiabatic Photodissociation from Product Vibrational Distributions. The Remarkable Case of Phenol

The fate of a photoexcited molecule is often strongly influenced by electronic degeneracies, such as conical intersections, which break the Born-Oppenheimer separation of electronic and nuclear motion. Detailed information concerning internal energy redistribution in a nonadiabatic process can be extracted from the product state distribution of a photofragment in photodissociation. Here, we focus on the nonadiabatic photodissociation of phenol and discuss the internal excitation of the phenoxyl fragment using both symmetry analysis and wave packet dynamics. It is shown that unique and general selection rules exist, which can be attributed to the geometric phase in the adiabatic representation. Further, our results provide a reinterpretation of the experimental data, shedding light on the impact of conical intersections on the product state distribution.

Spin-flip methods in quantum chemistry

This Perspective discusses salient features of the spin-flip approach to strong correlation and describes different methods that sprung from this idea. The spin-flip treatment exploits the different physics of low-spin and high-spin states and is based on the observation that correlation is small for same-spin electrons. By using a well-behaved high-spin state as a reference, one can access problematic low-spin states by deploying the same formal tools as in the excited-state treatments (i.e., linear response, propagator, or equation-of-motion theories). The Perspective reviews applications of this strategy within wave function and density functional theory frameworks as well as the extensions for molecular properties and spectroscopy. The utility of spin-flip methods is illustrated by examples. Limitations and proposed future directions are also discussed.

Machine learning enables long time scale molecular photodynamics simulations

Photo-induced processes are fundamental in nature but accurate simulations of their dynamics are seriously limited by the cost of the underlying quantum chemical calculations, hampering their application for long time scales. Here we introduce a method based on machine learning to overcome this bottleneck and enable accurate photodynamics on nanosecond time scales, which are otherwise out of reach with contemporary approaches. Instead of expensive quantum chemistry during molecular dynamics simulations, we use deep neural networks to learn the relationship between a molecular geometry and its high-dimensional electronic properties. As an example, the time evolution of the methylenimmonium cation for one nanosecond is used to demonstrate that machine learning algorithms can outperform standard excited-state molecular dynamics approaches in their computational efficiency while delivering the same accuracy.

Quantum Chemistry in the Age of Quantum Computing

Practical challenges in simulating quantum systems on classical computers have been widely recognized in the quantum physics and quantum chemistry communities over the past century. Although many approximation methods have been introduced, the complexity of quantum mechanics remains hard to appease. The advent of quantum computation brings new pathways to navigate this challenging and complex landscape. By manipulating quantum states of matter and taking advantage of their unique features such as superposition and entanglement, quantum computers promise to efficiently deliver accurate results for many important problems in quantum chemistry, such as the electronic structure of molecules. In the past two decades, significant advances have been made in developing algorithms and physical hardware for quantum computing, heralding a revolution in simulation of quantum systems. This Review provides an overview of the algorithms and results that are relevant for quantum chemistry. The intended audience is both quantum chemists who seek to learn more about quantum computing and quantum computing researchers who would like to explore applications in quantum chemistry.

Theoretical Investigation of Excited-State Intramolecular Proton Transfer and Photoisomerization of 2-(Iminomethyl)phenol

Time-dependent density functional theory and high-level ab initio calculations were performed to investigate the excited-state intramolecular proton transfer (ESIPT) and subsequent isomerization of 2-(iminomethyl)phenol (IMP). According to the results of the correlated theoretical methods, ESIPT is a barrierless process; subsequently, the isomerization (rotation of the torsion angle) of IMP also readily occurs. Fictitious intermediates are found due to an insufficient theoretical level. The molecular structure of the conical intersection (CI) during the isomerization process is optimized, and its branching plane is characterized. Both the gradient difference vector and the derivative coupling vector are significantly correlated to the C═O and H₂N-C antiparallel stretching coordinates, and the dynamic electron correlation effect is crucial to optimize the molecular structure of the real CI of IMP. The relaxation pathway from the CI in the S₀ state was examined; the dominant pathway proceeds to the trans-keto form of IMP. However, if the C═O and H₂N-C antiparallel stretching mode is sufficiently populated, then the reaction proceeds to the cis-keto form of IMP and can eventually recover to the cis-enol form.

The Simplest Possible Approach for Simulating S0- S1 Conical Intersections with DFT/TDDFT: Adding One Doubly Excited Configuration

A simple combination of density functional theory/time-dependent density functional theory (DFT/TDDFT) and configuration interaction is presented to fix the incorrect topology of the S₀- S₁ conical intersection (CI) and allow a description of bond making and bond breaking in photoinduced dynamics. The proposed TDDFT-1D method includes one lone optimized doubly excited configuration in addition to the DFT/TDDFT singly excited states within the context of a large configuration interaction Hamiltonian. Results for ethylene and stilbene are provided to demonstrate that this ansatz can yield physically meaningful potential energy surfaces near S₀- S₁ avoided crossings without changing the vertical excitation energies far from the relevant crossings. We also investigate the famous linear water example to show that the algorithm calculates the correct topology of the S₀- S₁ CI and yields the correct geometric phase.

Conical Intersections at the Nanoscale: Molecular Ideas for Materials

The ability to predict and describe nonradiative processes in molecules via the identification and characterization of conical intersections is one of the greatest recent successes of theoretical chemistry. Only recently, however, has this concept been extended to materials science, where nonradiative recombination limits the efficiencies of materials for various optoelectronic applications. In this review, we present recent advances in the theoretical study of conical intersections in semiconductor nanomaterials. After briefly introducing conical intersections, we argue that specific defects in materials can induce conical intersections between the ground and first excited electronic states, thus introducing pathways for nonradiative recombination. We present recent developments in theoretical methods, computational tools, and chemical intuition for the prediction of such defect-induced conical intersections. Through examples in various nanomaterials, we illustrate the significance of conical intersections for nanoscience. We also discuss challenges facing research in this area and opportunities for progress.

Excited state tracking during the relaxation of coordination compounds

The ability to locate minima on electronic excited states (ESs) potential energy surfaces both in the case of bright and dark states is crucial for a full understanding of photochemical reactions. This task has become a standard practice for small- to medium-sized organic chromophores thanks to the constant developments in the field of computational photochemistry. However, this remains a very challenging effort when it comes to the optimization of ESs of transition metal complexes (TMCs), not only due to the presence of several electronic ESs close in energy, but also due to the complex nature of the ESs involved. In this article, we present a simple yet powerful method to follow an ES of interest during a structural optimization in the case of TMCs, based on the use of a compact hole-particle representation of the electronic transition, namely the natural transition orbitals (NTOs). State tracking using NTOs is unambiguously accomplished by computing the mono-electronic wave function overlap between consecutive steps of the optimization. Here, we demonstrate that this simple but robust procedure works not only in the case of the cytosine but also in the case of the ES optimization of a ruthenium nitrosyl complex which is very problematic with standard approaches. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.

Dynamics of electronically excited states in the eumelanin building block 5,6-dihydroxyindole

We report the first excited state dynamics study of gas-phase 5,6-dihydroxyindole (5,6-DHI), a key building block of eumelanin pigments that are found throughout nature and serve as important photo-protective compounds. Time-resolved ion-yield measurements over the 241-296 nm ultraviolet photoexcitation region revealed non-adiabatic processes occurring on up to three distinct timescales. These reflect ultrafast (i.e. sub-picosecond) internal conversion within the excited state singlet manifold, and much longer-lived processes ranging from 10 ps to in excess of 1 ns. Our investigation paves the way for precisely targeted future studies of 5,6-DHI that exploit more differential measurement techniques. The work was facilitated by the use of soft laser-based thermal desorption to introduce 5,6-DHI samples into the gas phase. This approach, based on low-cost, readily available diode lasers, is straightforward, easily controllable and potentially applicable to a wide range of non-volatile molecular species.

Experimental observation of nonadiabatic bifurcation dynamics at resonances in the continuum

The surface crossing of bound and unbound electronic states in multidimensional space often gives rise to resonances in the continuum. This situation happens in the πσ*-mediated photodissociation reaction of 2-fluorothioanisole; optically-bright bound S1 (ππ*) vibrational states of 2-fluorothioanisole are strongly coupled to the optically-dark S2 (πσ*) state, which is repulsive along the S-CH3 elongation coordinate. It is revealed here that the reactive flux prepared at such resonances in the continuum bifurcates into two distinct reaction pathways with totally different dynamics in terms of energy disposal and nonadiabatic transition probability. This indicates that the reactive flux in the Franck-Condon region may either undergo nonadiabatic transition funneling through the conical intersection from the upper adiabat, or follow a low-lying adiabatic path, along which multiple dynamic saddle points may be located. Since 2-fluorothioanisole adopts a nonplanar geometry in the S1 minimum energy, the quasi-degenerate S1/S2 crossing seam in the nonplanar geometry, which lies well below the planar S1/S2 conical intersection, is likely responsible for the efficient vibronic coupling, especially in the low S1 internal energy region. As the excitation energy increases, bound-to-continuum coupling is facilitated with the aid of intramolecular vibrational redistribution, along many degrees of freedom spanning the large structural volume. This leads to the rapid domination of the continuum character of the reactive flux. This work reports direct and robust experimental observations of the nonadiabatic bifurcation dynamics of the reactive flux occurring at resonances in the continuum of polyatomic molecules.

Up to a Sign. The Insidious Effects of Energetically Inaccessible Conical Intersections on Unimolecular Reactions

It is now well established that conical intersections play an essential role in nonadiabatic radiationless decay where their double-cone topography causes them to act as efficient funnels channeling wave packets from the upper to the lower adiabatic state. Until recently, little attention was paid to the effect of conical intersections on dynamics on the lower state, particularly when the total energy involved is significantly below that of the conical intersection seam. This energetic deficiency is routinely used as a sufficient condition to exclude consideration of excited states in ground state dynamics. In this account, we show that, this energy criterion notwithstanding, energy inaccessible conical intersections can and do exert significant influence on lower state dynamics. The origin of this influence is the geometric phase, a signature property of conical intersections, which is the fact that the real-valued electronic wave function changes sign when transported along a loop containing a conical intersection, making the wave function double-valued. This geometric phase is permitted by an often neglected property of the real-valued adiabatic electronic wave function; namely, it is determined only up to an overall sign. Noting that in order to change sign a normalized, continuous function must go through zero, for loops of ever decreasing radii, demonstrating the need for an electronic degeneracy (intersection) to accompany the geometric phase. Since the total wave function must be single-valued a compensating geometry dependent phase needs to be included in the total electronic-nuclear wave function. This Account focuses on how this consequence of the geometric phase can modify nuclear dynamics energetically restricted to the lower state, including tunneling dynamics, in directly measurable ways, including significantly altering tunneling lifetimes, thus confounding the relation between measured lifetimes and barrier heights and widths, and/or completely changing product rotational distributions. Some progress has been made in understanding the origin of this effect. It has emerged that for a system where the lower adiabatic potential energy surface exhibits a topography comprised of two saddle points separated by a high energy conical intersection, the effect of the geometric phase can be quite significant. In this case topologically distinct paths through the two adiabatic saddle points may lead to interference. This was pointed out by Mead and Truhlar almost 50 years ago and denoted the Molecular Aharonov-Bohm effect. Still, the difficulty in anticipating a significant geometric phase effect in tunneling dynamics due to energetically inaccessible conical intersections leads to the attribute insidious that appears in the title of this Account. Since any theory is only as relevant as the prevalence of the systems it describes, we include in this Account examples of real systems where these effects can be observed. The accuracy of the reviewed calculations is high since we use fully quantum mechanical dynamics and construct the geometric phase using an accurate diabatic state fit of high quality ab initio data, energies, energy gradients, and interstate couplings. It remains for future work to establish the prevalence of this phenomenon and its deleterious effects on the conventional wisdom discussed in this work.

A discontinuous basis enables numerically exact solution of the Schrödinger equation around conical intersections in the adiabatic representation

Solving the vibrational Schrödinger equation in the neighborhood of conical intersections in the adiabatic representation is a challenge. At the intersection point, first- and second-derivative nonadiabatic coupling matrix elements become singular, with the singularity in the second-derivative coupling (diagonal Born-Oppenheimer correction) being non-integrable. These singularities result from discontinuities in the vibronic functions associated with the individual adiabatic states, and our group has recently argued that these divergent matrix elements cancel when discontinuous adiabatic vibronic functions sum to a continuous total nonadiabatic wave function. Here we describe the realization of this concept: a novel scheme for the numerically exact solution of the Schrödinger equation in the adiabatic representation. Our approach is based on a basis containing functions that are discontinuous at the intersection point. We demonstrate that the individual adiabatic nuclear wave functions are themselves discontinuous at the intersection point. This proves that discontinuous basis functions are essential to any tractable method that solves the Schrödinger equation around conical intersections in the adiabatic representation with high numerical precision. We establish that our method provides numerically exact results by comparison to reference calculations performed in the diabatic representation. In addition, we quantify the energetic error associated with constraining the density to be zero at the intersection point, a natural approximation. Prospects for extending the present treatment of a two-dimensional model to systems of higher dimensionality are discussed.

Using density based indexes to characterize excited states evolution

With the aim of offering new computational tools helping in the description of photochemical reactions and phenomena occurring at the excited state, we present in this work the capability of a density based index (Π) in locating decay channels from higher to lower excited states. The Π index, previously applied to disclose non-radiative decay channels from the first excited state to the ground state, is very simple in its formulation and can be evaluated, practically with no extra computational cost, and coupled to any quantum method able to provide excited states densities. Indeed, this index relies only on the knowledge of energetics and electron densities of the different electronic states involved in the decay. In the present work, we show the proficiency of the Π index in the general case of decay between excited states by applying it to two model systems well characterized both theoretically and experimentally. In both cases, this descriptor was successful in spotting the regions where excited states are more likely to decay, thus suggesting its potential interest for further application in the design of new compounds. © 2018 Wiley Periodicals, Inc.

Ultrafast imaging of laser-controlled non-adiabatic dynamics in NO2 from time-resolved photoelectron emission

Imaging and controlling the ultrafast conical intersection dynamics in NO₂ using the latest advances in attosecond and light-synthesizer technology.

Coupled double triplet state in singlet fission

The highly unusual state, 1(TT), is a coupled, double triplet state that has recently garnered significant attention. This multiexcitonic state can be formed by a quantum transition from a single-photon bright state in a variety of organic semiconducting materials. 1(TT)'s transient nature and similarity to independent triplets, however, has led to significant difficulties in characterization and prediction of its properties. Recent progress describing 1(TT) from theory and experiment are breaking through these difficulties, and have greatly advanced our comprehension of this state. Starting from the early description of 1(TT) in polyenes, this perspective discusses formation mechanisms, spectroscopic signatures, and the scope of intertriplet interactions. When employing singlet fission to generate charge carriers in a solar cell, 1(TT) has a central role. Due to the variety of coupling strengths between triplet states in 1(TT) amongst different chromophores, two different strategies are discussed to enable efficient charge carrier extraction. Continued growth in our understanding of 1(TT) may lead to control over complex quantum states for intriguing applications beyond high-efficiency, organic solar cells.