Nuclear dynamics of singlet exciton fission in pentacene single crystals

Nonadiabatic molecular dynamics with subsystem density functional theory: application to crystalline pentacene

In this work, we report the development and assessment of the nonadiabatic molecular dynamics approach with the electronic structure calculations based on the linearly scaling subsystem density functional method. The approach is implemented in an open-source embedded Quantum Espresso/Libra software specially designed for nonadiabatic dynamics simulations in extended systems. As proof of the applicability of this method to large condensed-matter systems, we examine the dynamics of nonradiative relaxation of excess excitation energy in pentacene crystals with the simulation supercells containing more than 600 atoms. We find that increased structural disorder observed in larger supercell models induces larger nonadiabatic couplings of electronic states and accelerates the relaxation dynamics of excited states. We conduct a comparative analysis of several quantum-classical trajectory surface hopping schemes, including two new methods proposed in this work (revised decoherence-induced surface hopping and instantaneous decoherence at frustrated hops). Most of the tested schemes suggest fast energy relaxation occurring with the timescales in the 0.7-2.0 ps range, but they significantly overestimate the ground state recovery rates. Only the modified simplified decay of mixing approach yields a notably slower relaxation timescales of 8-14 ps, with a significantly inhibited ground state recovery.

Phonon-Driven Femtosecond Dynamics of Excitons in Crystalline Pentacene from First Principles

Nonradiative exciton relaxation processes are critical for energy transduction and transport in optoelectronic materials, but how these processes are connected to the underlying crystal structure and the associated electron, exciton, and phonon band structures, as well as the interactions of all these particles, is challenging to understand. Here, we present a first-principles study of exciton-phonon relaxation pathways in pentacene, a paradigmatic molecular crystal and optoelectronic semiconductor. We compute the momentum- and band-resolved exciton-phonon interactions, and use them to analyze key scattering channels. We find that both exciton intraband scattering and interband scattering to parity-forbidden dark states occur on the same ∼100  fs timescale as a direct consequence of the longitudinal-transverse splitting of the bright exciton band. Consequently, exciton-phonon scattering exists as a dominant nonradiative relaxation channel in pentacene. We further show how the propagation of an exciton wave packet is connected with crystal anisotropy, which gives rise to the longitudinal-transverse exciton splitting and concomitant anisotropic exciton and phonon dispersions. Our results provide a framework for understanding the role of exciton-phonon interactions in exciton nonradiative lifetimes in molecular crystals and beyond.

Optical Absorption Properties in Pentacene/Tetracene Solid Solutions

Modifying the optical and electronic properties of crystalline organic thin films is of great interest for improving the performance of modern organic semiconductor devices. Therein, the statistical mixing of molecules to form a solid solution provides an opportunity to fine-tune optical and electronic properties. Unfortunately, the diversity of intermolecular interactions renders mixed organic crystals highly complex, and a holistic picture is still lacking. Here, we report a study of the optical absorption properties in solid solutions of pentacene and tetracene, two prototypical organic semiconductors. In the mixtures, the optical properties can be continuously modified by statistical mixing at the molecular level. Comparison with time-dependent density functional theory calculations on occupationally disordered clusters unravels the electronic origin of the low energy optical transitions. The disorder partially relaxes the selection rules, leading to additional optical transitions that manifest as optical broadening. Furthermore, the contribution of diabatic charge-transfer states is modified in the mixtures, reducing the observed splitting in the 0-0 vibronic transition. Additional comparisons with other blended systems generalize our results and indicate that changes in the polarizability of the molecular environment in organic thin-film blends induce shifts in the absorption spectrum.

Computational Approaches for Organic Semiconductors: From Chemical and Physical Understanding to Predicting New Materials

While a complete understanding of organic semiconductor (OSC) design principles remains elusive, computational methods─ranging from techniques based in classical and quantum mechanics to more recent data-enabled models─can complement experimental observations and provide deep physicochemical insights into OSC structure-processing-property relationships, offering new capabilities for in silico OSC discovery and design. In this Review, we trace the evolution of these computational methods and their application to OSCs, beginning with early quantum-chemical methods to investigate resonance in benzene and building to recent machine-learning (ML) techniques and their application to ever more sophisticated OSC scientific and engineering challenges. Along the way, we highlight the limitations of the methods and how sophisticated physical and mathematical frameworks have been created to overcome those limitations. We illustrate applications of these methods to a range of specific challenges in OSCs derived from π-conjugated polymers and molecules, including predicting charge-carrier transport, modeling chain conformations and bulk morphology, estimating thermomechanical properties, and describing phonons and thermal transport, to name a few. Through these examples, we demonstrate how advances in computational methods accelerate the deployment of OSCsin wide-ranging technologies, such as organic photovoltaics (OPVs), organic light-emitting diodes (OLEDs), organic thermoelectrics, organic batteries, and organic (bio)sensors. We conclude by providing an outlook for the future development of computational techniques to discover and assess the properties of high-performing OSCs with greater accuracy.

Phonon-Induced Localization of Excitons in Molecular Crystals from First Principles

The spatial extent of excitons in molecular systems underpins their photophysics and utility for optoelectronic applications. Phonons are reported to lead to both exciton localization and delocalization. However, a microscopic understanding of phonon-induced (de)localization is lacking, in particular, how localized states form, the role of specific vibrations, and the relative importance of quantum and thermal nuclear fluctuations. Here, we present a first-principles study of these phenomena in solid pentacene, a prototypical molecular crystal, capturing the formation of bound excitons, exciton-phonon coupling to all orders, and phonon anharmonicity, using density functional theory, the ab initio GW-Bethe-Salpeter equation approach, finite-difference, and path integral techniques. We find that for pentacene zero-point nuclear motion causes uniformly strong localization, with thermal motion providing additional localization only for Wannier-Mott-like excitons. Anharmonic effects drive temperature-dependent localization, and, while such effects prevent the emergence of highly delocalized excitons, we explore the conditions under which these might be realized.

Direct Observation of Ultrafast Lattice Distortions during Exciton-Polaron Formation in Lead Halide Perovskite Nanocrystals

The microscopic origin of slow hot-carrier cooling in lead halide perovskites remains debated and has direct implications for applications. Slow hot-carrier cooling of several picoseconds has been attributed to either polaron formation or a hot-phonon bottleneck effect at high excited carrier densities (>10¹⁸ cm^-3). These effects cannot be unambiguously disentangled with optical experiments alone. However, they can be distinguished by direct observations of ultrafast lattice dynamics, as these effects are expected to create qualitatively distinct fingerprints. To this end, we employ femtosecond electron diffraction and directly measure the sub-picosecond lattice dynamics of weakly confined CsPbBr₃ nanocrystals following above-gap photoexcitation. While we do not observe signatures of a hot-phonon bottleneck lasting several picoseconds, the data reveal a light-induced structural distortion appearing on a time scale varying between 380 and 1200 fs depending on the excitation fluence. We attribute these dynamics to the effect of exciton-polarons on the lattice and the slower dynamics at high fluences to slower sub-picosecond hot-carrier cooling, which slows down the establishment of the exciton-polaron population. Further analysis and simulations show that the distortion is consistent with motions of the [PbBr₃]^- octahedral ionic cage, and closest agreement with the data is obtained for Pb-Br bond lengthening. Our work demonstrates how direct studies of lattice dynamics on the sub-picosecond time scale can discriminate between competing scenarios proposed in the literature to explain the origin of slow hot-carrier cooling in lead halide perovskites.

Spectroscopic analysis of vibrational coupling in multi-molecular excited states

Multi-molecular excited states accompanied by intra- and inter-molecular geometric relaxation are commonly encountered in optical and electrooptical studies and applications of organic semiconductors as, for example, excimers or charge transfer states. Understanding the dynamics of these states is crucial to improve organic devices such as light emitting diodes and solar cells. Their full microscopic description, however, demands sophisticated tools such as ab initio quantum chemical calculations which come at the expense of high computational costs and are prone to errors by assumptions as well as iterative algorithmic procedures. Hence, the analysis of spectroscopic data is often conducted at a phenomenological level only. Here, we present a toolkit to analyze temperature dependent luminescence data and gain first insights into the relevant microscopic parameters of the molecular system at hand. By means of a Franck-Condon based approach considering a single effective inter-molecular vibrational mode and different potentials for the ground and excited state we are able to explain the luminescence spectra of such multi-molecular states. We demonstrate that by applying certain reasonable simplifications the luminescence of charge transfer states as well as excimers can be satisfactorily reproduced for temperatures ranging from cryogenics to above room temperature. We present a semi-classical and a quantum-mechanical description of our model and, for both cases, demonstrate its applicability by analyzing the temperature dependent luminescence of the amorphous donor-acceptor heterojunction tetraphenyldibenzoperiflanthene:C₆₀ as well as polycrystalline zinc-phthalocyanine to reproduce the luminescence spectra and extract relevant system parameters such as the excimer binding energy.

From a one-mode to a multi-mode understanding of conical intersection mediated ultrafast organic photochemical reactions

This review discusses how ultrafast organic photochemical reactions are controlled by conical intersections, highlighting that decay to the ground-state at multiple points of the intersection space results in their multi-mode character.

Nanoscale coherent phonon spectroscopy

Coherent phonon spectroscopy can provide microscopic insight into ultrafast lattice dynamics and its coupling to other degrees of freedom under nonequilibrium conditions. Ultrafast optical spectroscopy is a well-established method to study coherent phonons, but the diffraction limit has hampered observing their local dynamics directly. Here, we demonstrate nanoscale coherent phonon spectroscopy using ultrafast laser-induced scanning tunneling microscopy in a plasmonic junction. Coherent phonons are locally excited in ultrathin zinc oxide films by the tightly confined plasmonic field and are probed via the photoinduced tunneling current through an electronic resonance of the zinc oxide film. Concurrently performed tip-enhanced Raman spectroscopy allows us to identify the involved phonon modes. In contrast to the Raman spectra, the phonon dynamics observed in coherent phonon spectroscopy exhibit strong nanoscale spatial variations that are correlated with the distribution of the electronic local density of states resolved by scanning tunneling spectroscopy.

Direct observation of ultrafast singlet exciton fission in three dimensions

We present quantitative ultrafast interferometric pump-probe microscopy capable of tracking of photoexcitations with sub-10 nm spatial precision in three dimensions with 15 fs temporal resolution, through retrieval of the full transient photoinduced complex refractive index. We use this methodology to study the spatiotemporal dynamics of the quantum coherent photophysical process of ultrafast singlet exciton fission. Measurements on microcrystalline pentacene films grown on glass (SiO₂) and boron nitride (hBN) reveal a 25 nm, 70 fs expansion of the joint-density-of-states along the crystal a,c-axes accompanied by a 6 nm, 115 fs change in the exciton density along the crystal b-axis. We propose that photogenerated singlet excitons expand along the direction of maximal orbital π-overlap in the crystal a,c-plane to form correlated triplet pairs, which subsequently electronically decouples into free triplets along the crystal b-axis due to molecular sliding motion of neighbouring pentacene molecules. Our methodology lays the foundation for the study of three dimensional transport on ultrafast timescales.

A New Frontier in Exciton Transport: Transient Delocalization

Efficient exciton transport is crucial to the application of organic semiconductors (OSCs) in light-harvesting devices. While the physics of exciton transport in highly disordered media is well-explored, the description of transport in structurally and energetically ordered OSCs is less established, despite such materials being favorable for devices. In this Perspective we describe and highlight recent research pointing toward a highly efficient exciton transport mechanism which occurs in ordered OSCs, transient delocalization. Here, exciton-phonon couplings play a critical role in allowing localized exciton states to temporarily access higher-energy delocalized states whereupon they move large distances. The mechanism shows great promise for facilitating long-range exciton transport and may allow for improved device efficiencies and new device architectures. However, many fundamental questions on transient delocalization remain to be answered. These questions and suggested next steps are summarized.

Chemical Modifications Suppress Anharmonic Effects in the Lattice Dynamics of Organic Semiconductors

The lattice dynamics of organic semiconductors has a significant role in determining their electronic and mechanical properties. A common technique to control these macroscopic properties is to chemically modify the molecular structure. These modifications are known to change the molecular packing, but their effect on the lattice dynamics is relatively unexplored. Therefore, we investigate how chemical modifications to a core [1]benzothieno[3,2-b]benzothiophene (BTBT) semiconducting crystal affect the evolution of the crystal structural dynamics with temperature. Our study combines temperature-dependent polarization-orientation (PO) low-frequency Raman measurements with first-principles calculations and single-crystal X-ray diffraction measurements. We show that chemical modifications can indeed suppress specific expressions of vibrational anharmonicity in the lattice dynamics. Specifically, we detect in BTBT a gradual change in the PO Raman response with temperature, indicating a unique anharmonic expression. This anharmonic expression is suppressed in all examined chemically modified crystals (ditBu-BTBT and diC8-BTBT, diPh-BTBT, and DNTT). In addition, we observe solid–solid phase transitions in the alkyl-modified BTBTs. Our findings indicate that π-conjugated chemical modifications are the most effective in suppressing these anharmonic effects.

Evaluating the Electronic Structure of Coexisting Excitonic and Multiexcitonic States in Periodic Systems: Significance for Singlet Fission

Singlet fission (SF) in organic molecular solids is an example of a process that is challenging to describe with the most common electronic structure approaches. It involves optically bright singlet excited states delocalized over many molecules, which could be efficiently treated by density functional theory, and multiexcitonic localized states that have to be studied with wavefunction methods, usually with small clusters considering their expensive computational costs. In this work, we propose a methodology to combine multiconfigurational wavefunction calculations with reduced Hamiltonian to investigate the electronic structure of large clusters or fully periodic systems. The method is applied to the prototypical SF materials tetracene and pentacene. The results allow one to study how states of different natures (excitonic, charge-transfer, and multiexcitonic) coexist and are contaminated by their couplings in large or periodic systems. Novel insights are therefore possible. For example, because the excitonic bands are relatively broad with respect to the multiexcitonic states, there are limited regions of the crystal momentum space where the transition between the two is more likely.