Search

The doublet-spin nature of radical emitters is advantageous for applications in organic light-emitting diodes, as it avoids the formation of triplet excitons that limit the electroluminescence efficiency of non-radical emitters. However, radicals generally show low optical absorption and photoluminescence yields. Here we explain the poor optical properties of radicals based on alternant hydrocarbons, and establish design rules to increase the absorption and luminescence yields for donor-acceptor-type radicals. We show that non-alternant systems are necessary to lift the degeneracy of the lowest energy orbital excitations; moreover, intensity borrowing from an intense high-lying transition by the low-energy charge-transfer excitation enhances the oscillator strength of the emitter. We apply these rules to design tris(2,4,6-trichlorophenyl)methyl-pyridoindolyl derivatives with a high photoluminescence quantum yield (>90%). Organic light-emitting diodes based on these molecules showed a pure-red emission with an over 12% external quantum efficiency. These insights may be beneficial for the rational design and discovery of highly luminescent doublet emitters.

The discovery of molecules with tailored optoelectronic properties such as specific frequency and intensity of absorption or emission is a major challenge in creating next-generation organic light-emitting diodes (OLEDs) and photovoltaics. This raises the question: how can we predict a potential chemical structure from these properties? Approaches that attempt to tackle this inverse design problem include virtual screening, active machine learning and genetic algorithms. However, these approaches rely on a molecular database or many electronic structure calculations, and significant computational savings could be achieved if there was prior knowledge of (i) whether the optoelectronic properties of a parent molecule could easily be improved and (ii) what morphing operations on a parent molecule could improve these properties. In this perspective we address both of these challenges from first principles. We firstly adapt the Thomas-Reiche-Kuhn sum rule to organic chromophores and show how this indicates how easily the absorption and emission of a molecule can be improved. We then show how by combining electronic structure theory and intensity borrowing perturbation theory we can predict whether or not the proposed morphing operations will achieve the desired spectral alteration, and thereby derive widely-applicable design rules. We go on to provide proof-of-concept illustrations of this approach to optimizing the visible absorption of acenes and the emission of radical OLEDs. We believe this approach can be integrated into genetic algorithms by biasing morphing operations in favour of those which are likely to be successful, leading to faster molecular discovery and greener chemistry., Comment: The following article has been accepted by The Journal of Chemical Physics. After it is published, it will be found at https://aip.scitation.org/doi/10.1063/5.0082311

Since the demonstration in 2018 that organic radicals can be used to make highly efficient organic light-emitting diodes there has been an explosion of interest in their capabilities and many experimental and computational studies of their performance. Here we take a theoretical view and describe the electronic structure of radicals from an algebraic perspective. By rediscovering and adapting historic investigations of organic radicals we show how many experimentally useful properties can be determined without synthesis or computation, but simply from knowledge of the molecular structure and in particular whether or not the radical is an alternant hydrocarbon. We explain these results in the context of modern organic light-emitting design in order to inform future investigations.

Efficient charge photogeneration in conjugated polymers typically requires the presence of a second component to act as electron acceptor. Here, we report a novel low band-gap conjugated polymer with a donor/orthogonal acceptor motif: poly-2,6-(4,4-dihexadecyl-4H-cyclopenta [2,1-b:3,4-b′]dithiophene)-alt-2,6-spiro [cyclopenta[2,1-b:3,4-b′]dithiophene-4,9′-fluorene]-2′,7′-dicarbonitrile, referred to as PCPDT-sFCN. The role of the orthogonal acceptor is to spatially isolate the LUMO from the HOMO, allowing for negligible exchange energy between electrons in these orbitals and minimising the energy gap between singlet and triplet charge transfer states. We employ ultrafast and microsecond transient absorption spectroscopy to demonstrate that, even in the absence of a separate electron acceptor, PCPDT-sFCN shows efficient charge photogeneration in both pristine solution and film. This efficient charge generation is a result of an isoenergetic singlet/triplet charge transfer state equilibrium acting as a reservoir for charge carrier formation. Furthermore, clear evidence of enhanced triplet populations, which form in less than 1 ps, is observed. Using group theory, we show that this ultrafast triplet formation is due to highly efficient, quantum mechanically allowed intersystem crossing between the bright, initially photoexcited local singlet state and the triplet charge transfer state. Remarkably, the free charges that form via the charge transfer state are extraordinarily long-lived with millisecond lifetimes, possibly due to the stabilisation imparted by the spatial separation of PCPDT-sFCN's donor and orthogonal acceptor motifs. The efficient generation of long-lived charge carriers in a pristine polymer paves the way for single-material applications such as organic photovoltaics and photodetectors., The spatial separation of PCPDT-sFCN's donor and orthogonal acceptor motifs allows efficient photogeneration of extraordinarily long-lived charge carriers in the pristine polymer, providing an important step towards single-material optoelectronics.

Recent synthetic studies on the organic molecules tetracene and pentacene have found certain dimers and oligomers to exhibit an intense absorption in the visible region of the spectrum which is not present in the monomer or many previously-studied dimers. In this article we combine experimental synthesis with electronic structure theory and spectral computation to show that this absorption arises from an otherwise dark charge-transfer excitation 'borrowing intensity' from an intense UV excitation. Further, by characterizing the role of relevant monomer molecular orbitals, we arrive at a design principle that allows us to predict the presence or absence of an additional absorption based on the bonding geometry of the dimer. We find this rule correctly explains the spectra of a wide range of acene derivatives and solves an unexplained structure-spectrum phenomenon first observed seventy years ago. These results pave the way for the design of highly absorbent chromophores with applications ranging from photovoltaics to liquid crystals., 29 pages, 63 pages with SI, 5 figures

Zombie States are a recently introduced formalism to describe coupled coherent Fermionic states which address the Fermionic sign problem in a computationally tractable manner. Previously it has been shown that Zombie States with fractional occupations of spin-orbitals obeyed the correct Fermionic creation and annihilation algebra and presented results for real-time evolution [Dmitrii V. Shalashilin, J. Chem. Phys. 148, 194109 (2018)]. In this work we extend and build on this formalism by developing efficient algorithms for evaluating the Hamiltonian and other operators between Zombie States and address their normalization. We also show how imaginary time propagation can be used to find the ground state of a system. We also present a biasing method, for setting up a basis set of random Zombie States, that allow much smaller basis sizes to be used while still accurately describing the electronic structure Hamiltonian and its ground state and describe a technique of wave function "cleaning" which removes the contributions of configurations with the wrong number of electrons, improving the accuracy further. We also show how low-lying excited states can be calculated efficiently using a Gram-Schmidt orthogonalization procedure.The proposed algorithm of imaginary time propagation on a biased random grids of Zombie States may present an alternative to existing Quantum Monte Carlo methods., Comment: 45 pages, 9 figures, supplementary material 3 pages, submitted to the Journal of Chemical Physics

Organic light-emitting diodes (OLEDs)1–5, quantum-dot-based LEDs6–10, perovskite-based LEDs11–13 and micro-LEDs14,15 have been championed to fabricate lightweight and flexible units for next-generation displays and active lighting. Although there are already some high-end commercial products based on OLEDs, costs must decrease whilst maintaining high operational efficiencies for the technology to realise wider impact. Here we demonstrate efficient action of radical-based OLEDs16, whose emission originates from a spin doublet, rather than a singlet or triplet exciton. While the emission process is still spin-allowed in these OLEDs, the efficiency limitations imposed by triplet excitons are circumvented for doublets. Using a luminescent radical emitter, we demonstrate an OLED with maximum external quantum efficiency of 27 per cent at a wavelength of 710 nanometres—the highest reported value for deep-red and infrared LEDs. For a standard closed-shell organic semiconductor, holes and electrons occupy the highest occupied and lowest unoccupied molecular orbitals (HOMOs and LUMOs), respectively, and recombine to form singlet or triplet excitons. Radical emitters have a singly occupied molecular orbital (SOMO) in the ground state, giving an overall spin-1/2 doublet. If—as expected on energetic grounds—both electrons and holes occupy this SOMO level, recombination returns the system to the ground state, giving no light emission. However, in our very efficient OLEDs, we achieve selective hole injection into the HOMO and electron injection to the SOMO to form the fluorescent doublet excited state with near-unity internal quantum efficiency. Organic light-emitting devices containing radical emitters can achieve an efficiency of 27 per cent at deep-red and infrared wavelengths based on the excitation of spin doublets, rather than singlet or triplet states.

Thermal quantum time-correlation functions are of fundamental importance in quantum dynamics, allowing experimentally-measurable properties such as reaction rates, diffusion constants and vibrational spectra to be computed from first principles. Since the exact quantum solution scales exponentially with system size, there has been considerable effort in formulating reliable linear-scaling methods involving exact quantum statistics and approximate quantum dynamics modelled with classical-like trajectories. Here we review recent progress in the field with the development of methods including Centroid Molecular Dynamics (CMD), Ring Polymer Molecular Dynamics (RPMD) and Thermostatted RPMD (TRPMD). We show how these methods have recently been obtained from `Matsubara dynamics', a form of semiclassical dynamics which conserves the quantum Boltzmann distribution. We also rederive t->0+ quantum transition-state theory (QTST) in the Matsubara dynamics formalism showing that Matsubara-TST, like RPMD-TST, is equivalent to QTST. We end by surveying areas for future progress., 25 pages, 8 figures. Submitted as a New View article to Molecular Physics on 11th January 2017

Carbene-metal-amides (CMAs) are a promising family of donor–bridge–acceptor molecular charge-transfer (CT) emitters for organic light-emitting diodes. A universal approach is demonstrated to tune the energy of their CT emission. A blueshift of up to 210 meV is achievable in solid state via dilution in a polar host matrix. The origin of this shift has two components: constraint of thermally-activated triplet diffusion, and electrostatic interactions between guest and polar host. This allows the emission of mid-green CMA archetypes to be tuned to sky blue without chemical modifications. Monte-Carlo simulations based on a Marcus-type transfer integral successfully reproduce the concentration- and temperature-dependent triplet diffusion process, revealing a substantial shift in the ensemble density of states in polar hosts. In gold-bridged CMAs, this shift does not lead to a significant change in luminescence lifetime, thermal activation energy, reorganization energy, or intersystem crossing rate. These discoveries offer new insight into coupling between the singlet and triplet manifolds in CMA materials, revealing a dominant interaction between states of CT character. The same approach is employed using materials which have been chemically modified to alter the energy of their CT state directly, shifting the emission of sky-blue chromophores into the practical blue range.

The energies of molecular excited states arise as solutions to the electronic Schr\"{o}dinger equation and are often compared to experiment. At the same time, nuclear quantum motion is known to be important and to induce a red-shift of excited state energies. However, it is thus far unclear whether incorporating nuclear quantum motion in molecular excited state calculations leads to a systematic improvement of their predictive accuracy, making further investigation necessary. Here we present such an investigation by employing two first-principles methods for capturing the effect of quantum fluctuations on excited state energies, which we apply to the Thiel set of organic molecules. We show that accounting for zero-point motion leads to much improved agreement with experiment, compared to `static' calculations which only account for electronic effects, and the magnitude of the red-shift can become as large as 1.36 eV. Moreover, we show that the effect of nuclear quantum motion on excited state energies largely depends on the molecular size, with smaller molecules exhibiting larger red-shifts. Our methodology also makes it possible to analyze the contribution of individual vibrational normal modes to the red-shift of excited state energies, and in several molecules we identify a limited number of modes dominating this effect. Overall, our study provides a foundation for systematically quantifying the shift of excited state energies due to nuclear quantum motion, and for understanding this effect at a microscopic level.

International audience; Singlet fission in organic semiconductors causes a singlet exciton to decay into a pair of triplet excitons and holds potential for increasing the efficiency of photovoltaic devices. In this combined experimental and theoretical study, we reveal that a covalent dimer of the organic semiconductor tetracene undergoes activated singlet fission by qualitatively different mechanisms depending on the solvent environment. We show that intramolecular vibrations are an integral part of this mechanism, giving rise to mixing between charge transfer and triplet pair excitations. Either coherent or incoherent singlet fission can occur, depending on the transient solvent-induced energetic proximity between the states, giving rise to complex variation of the singlet fission mechanism and time scale in the different environments. Our results suggest a more general principle for controlling the efficiency of photochemical reactions by utilizing transient interactions to tune the energetics of reactant and product states and switch between incoherent and coherent dynamic

We investigate the mechanisms of condensed phase proton-coupled electron transfer (PCET) using Mapping-Variable Ring Polymer Molecular Dynamics (MV-RPMD), a recently developed method that employs an ensemble of classical trajectories to simulate nonadiabatic excited state dynamics. Here, we construct a series of system-bath model Hamiltonians for the PCET, where four localized electron-proton states are coupled to a thermal bath via a single solvent mode, and we employ MV-RPMD to simulate state population dynamics. Specifically, for each model, we identify the dominant PCET mechanism, and by comparing against rate theory calculations, we verify that our simulations correctly distinguish between concerted PCET, where the electron and proton transfer together, and sequential PCET, where either the electron or the proton transfers first. This work represents a first application of MV-RPMD to multi-level condensed phase systems; we introduce a modified MV-RPMD expression that is derived using a symmetric rather than asymmetric Trotter discretization scheme and an initialization protocol that uses a recently derived population estimator to constrain trajectories to a dividing surface. We also demonstrate that, as expected, the PCET mechanisms predicted by our simulations are robust to an arbitrary choice of the initial dividing surface.

We extend the Mixed Quantum-Classical Initial Value Representation (MQC-IVR), a semiclassical method for computing real-time correlation functions, to electronically nonadiabatic systems using the Meyer-Miller-Stock-Thoss (MMST) Hamiltonian to treat electronic and nuclear degrees of freedom (dofs) within a consistent dynamic framework. We introduce an efficient symplectic integration scheme, the MInt algorithm, for numerical time-evolution of the nuclear and electronic phase space variables as well as the Monodromy matrix, under the non-separable MMST Hamiltonian. We then calculate the probability of transmission through a curve-crossing in model two-level systems and show that in the quantum limit MQC-IVR is in good agreement with the exact quantum results, whereas in the classical limit the method yields results in keeping with mean-field approaches like the Linearized Semiclassical IVR. Finally, exploiting the ability of MQC-IVR to quantize different dofs to different extents, we present a detailed study of the extents to which quantizing the nuclear and electronic dofs improves numerical convergence properties without significant loss of accuracy., Comment: 16 pages, 6 figures. An earlier version of this manuscript has been submitted to JCP

Name: Timothy HeleInstitutions: Cornell University, USA; Cambridge University, UKInvited article: Thermal quantum time-correlation functions from classical-like dynamicsFinding out how atoms and el...

In a previous article [T. J. H. Hele and S. C. Althorpe, J. Chem. Phys. 138, 084108 (2013)], we showed that the t → 0+ limit of ring-polymer molecular dynamics (RPMD) rate-theory is also the t → 0+ limit of a new type of quantum flux-side time-correlation function, in which the dividing surfaces are invariant to imaginary-time translation; in other words, that RPMD transition-state theory (RMPD-TST) is a t → 0+ quantum transition-state theory (QTST). Recently, Jang and Voth [J. Chem. Phys. 144, 084110 (2016)] rederived this quantum t → 0+ limit and claimed that it gives instead the centroid-density approximation. Here we show that the t → 0+ limit derived by Jang and Voth is in fact RPMD-TST.

We derive an exact quantum propagator for nonadiabatic dynamics in multi-state systems using the mapping variable representation, where classical-like Cartesian variables are used to represent both continuous nuclear degrees of freedom and discrete electronic states. The resulting expression is a Moyal series that, when suitably approximated, can allow for the use of classical dynamics to efficiently model large systems. We demonstrate that different truncations of the exact propagator lead to existing approximate semiclassical and mixed quantum-classical methods and we derive an associated error term for each method. Furthermore, by combining the imaginary-time path-integral representation of the Boltzmann operator with the exact propagator, we obtain an analytic expression for thermal quantum real-time correlation functions. These results provide a rigorous theoretical foundation for the development of accurate and efficient classical-like dynamics to compute observables such as electron transfer reaction rates in complex quantized systems., Comment: 22 pages, 2 figures. Submitted to the "Reaction Rate Theory" Faraday Discussion (i.e. conference) to be held in Cambridge, UK in September 2016

Clathrate hydrates hold considerable promise as safe and economical materials for hydrogen storage. Here we present a quantum mechanical study of H$_2$ and D$_2$ diffusion through a hexagonal face shared by two large cages of clathrate hydrates over a wide range of temperatures. Path integral molecular dynamics simulations are used to compute the free-energy profiles for the diffusion of H$_2$ and D$_2$ as a function of temperature. Ring polymer molecular dynamics rate theory, incorporating both exact quantum statistics and approximate quantum dynamical effects, is utilized in the calculations of the H$_2$ and D$_2$ diffusion rates in a broad temperature interval. We find that the shape of the quantum free-energy profiles and their height relative to the classical free energy barriers at a given temperature, as well as the rate of diffusion, are profoundly affected by competing quantum effects: above 25 K, zero-point energy (ZPE) perpendicular to the reaction path for diffusion between cavities decreases the quantum rate compared to the classical rate, whereas at lower temperatures tunneling outcompetes the ZPE and as result the quantum rate is greater than the classical rate., Comment: 9 pages, 4 figures

We apply Thermostatted Ring Polymer Molecular Dynamics (TRPMD), a recently-proposed approximate quantum dynamics method, to the computation of thermal reaction rates. Its short-time Transition-State Theory (TST) limit is identical to rigorous Quantum Transition-State Theory, and we find that its long-time limit is independent of the location of the dividing surface. TRPMD rate theory is then applied to one-dimensional model systems, the atom-diatom bimolecular reactions H+H$_2$, D+MuH and F+H$_2$, and the prototypical polyatomic reaction H+CH$_4$. Above the crossover temperature, the TRPMD rate is virtually invariant to the strength of the friction applied to the internal ring-polymer normal modes, and beneath the crossover temperature the TRPMD rate generally decreases with increasing friction, in agreement with the predictions of Kramers theory. We therefore find that TRPMD is approximately equal to, or less accurate than, Ring Polymer Molecular Dynamics (RPMD) for symmetric reactions, and for certain asymmetric systems and friction parameters closer to the quantum result, providing a basis for further assessment of the accuracy of this method., 13 pages, 7 figures, 2 tables. Minor alterations and corrections to the previous version in light of referees' comments

We recently obtained a quantum-Boltzmann-conserving classical dynamics by making a single change to the derivation of the “Classical Wigner” approximation. Here, we show that the further approximation of this “Matsubara dynamics” gives rise to two popular heuristic methods for treating quantum Boltzmann time-correlation functions: centroid molecular dynamics (CMD) and ring-polymer molecular dynamics (RPMD). We show that CMD is a mean-field approximation to Matsubara dynamics, obtained by discarding (classical) fluctuations around the centroid, and that RPMD is the result of discarding a term in the Matsubara Liouvillian which shifts the frequencies of these fluctuations. These findings are consistent with previous numerical results and give explicit formulae for the terms that CMD and RPMD leave out.

We show that a single change in the derivation of the linearized semiclassical-initial value representation (LSC-IVR or classical Wigner approximation) results in a classical dynamics which conserves the quantum Boltzmann distribution. We rederive the (standard) LSC-IVR approach by writing the (exact) quantum time-correlation function in terms of the normal modes of a free ring-polymer (i.e. a discrete imaginary-time Feynman path), taking the limit that the number of polymer beads $N \to \infty$, such that the lowest normal-mode frequencies take their Matsubara values. The change we propose is to truncate the quantum Liouvillian, not explicitly in powers of $\hbar^2$ at $\hbar^0$ (which gives back the standard LSC-IVR approximation), but in the normal-mode derivatives corresponding to the lowest Matsubara frequencies. The resulting Matsubara dynamics is inherently classical (since all terms $\mathcal{O}\left(\hbar^{2}\right)$ disappear from the Matsubara Liouvillian in the limit $N \to \infty$), and conserves the quantum Boltzmann distribution because the Matsubara Hamiltonian is symmetric with respect to imaginary-time translation. Numerical tests show that the Matsubara approximation to the quantum time-correlation function converges with respect to the number of modes, and gives better agreement than LSC- IVR with the exact quantum result. Matsubara dynamics is too computationally expensive to be applied to complex systems, but its further approximation may lead to practical methods.

It was shown recently that there exists a true quantum transition-state theory (QTST) corresponding to the t->0+ limit of a (new form of) quantum flux-side time-correlation function. Remarkably, this QTST is identical to ring-polymer molecular dynamics (RPMD) TST. Here we provide evidence which suggests very strongly that this QTST (= RPMD-TST) is unique, in the sense that the t->0+ limit of any other flux-side time-correlation function gives either non-positive-definite quantum statistics or zero. We introduce a generalized flux-side time-correlation function which includes all other (known) flux-side time-correlation functions as special limiting cases. We find that the only non-zero t->0+ limit of this function that contains positive-definite quantum statistics is RPMD-TST., Comment: 10 pages, 1 figure. Typographical errors corrected, references updated and added

In Part I [J. Chem. Phys. 138, 084108 (2013)] we derived a quantum transition-state theory by taking the t->0+ (short-time) limit of a new form of quantum flux-side time-correlation function containing a ring-polymer dividing surface. This t->0+ limit appears to be unique in giving positive-definite Boltzmann statistics, and is identical to ring-polymer molecular dynamics (RPMD) TST. Here, we show that quantum TST (i.e. RPMD-TST) is exact if there is no recrossing (by the real-time quantum dynamics) of the ring-polymer dividing surface, nor of any surface orthogonal to it in the space describing fluctuations in the polymer-bead positions along the reaction coordinate. In practice, this means that RPMD-TST gives a good approximation to the exact quantum rate for direct reactions, provided the temperature is not too far below the cross-over to deep tunnelling. We derive these results by comparing the long-time limit of the ring-polymer flux-side time-correlation function with that of a hybrid flux-side time-correlation function (containing a ring-polymer flux operator and a Miller-Schwarz-Tromp side function), and by representing the resulting ring-polymer momentum integrals as hypercubes. Together with Part I, the results of this article validate a large number of RPMD calculations of reaction rates., Comment: 14 pages, 4 figures. Argument and wording clarified, typographical errors corrected and references added

Surprisingly, there exists a quantum flux-side time-correlation function which has a non-zero short-time (t->0+) limit, and thus yields a rigorous quantum generalization of classical transition-state theory (TST). In this Part I of two articles, we introduce the new time-correlation function, and derive its short-time limit. The new ingredient is a generalized Kubo transform which allows the flux and side dividing surfaces to be the same function of path-integral space. Choosing this common dividing surface to be a single point gives a short-time limit which is identical to an expression introduced on heuristic grounds by Wigner in 1932, but which does not give positive-definite quantum statistics, causing it to fail while still in the shallow-tunnelling regime. Choosing the dividing surface to be invariant to imaginary-time translation gives, uniquely, a short-time limit that gives the correct positive- definite quantum statistics at all temperatures, and which is identical to ring-polymer molecular dynamics (RPMD) TST. We find that the RPMD-TST rate is not a strict upper bound to the exact quantum rate, but a good approximation to one if real-time coherence effects are small. Part II will show that the RPMD-TST rate is equal to the exact quantum rate in the absence of recrossing., Comment: 14 pages, 5 figures, 1 table