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

Nonadiabatic Simulation of Exciton Dynamics in Organic Semiconductors Using Neural Network-Based Frenkel Hamiltonian and Gradients

In this study, we present a multiscale method to simulate the propagation of Frenkel singlet excitons in organic semiconductors (OSCs). The approach uses neural network models to train a Frenkel-type Hamiltonian and its gradient, obtained by the long-range correction version of density functional tight-binding with self-consistent charges. Our models accurately predict site energies, excitonic couplings, and corresponding gradients, essential for the nonadiabatic molecular dynamics simulations. Combined with the fewest switches surface hopping algorithm, the method was applied to four representative OSCs: anthracene, pentacene, perylenediimide, and diindenoperylene. The simulated exciton diffusion constants align well with experimental and reported theoretical values and offer valuable insights into exciton dynamics in OSCs.

Influence of Electronic Correlations on Electron-Phonon Interactions of Molecular Systems with the GW and Coupled Cluster Methods

Electron-phonon interactions are of great importance to a variety of physical phenomena, and their accurate description is an important goal for first-principles calculations. Isolated examples of materials and molecular systems have emerged where electron-phonon coupling is enhanced over density functional theory (DFT) when using the Green's-function-based ab initio GW method, which provides a more accurate description of electronic correlations. It is, however, unclear how general this enhancement is and how employing high-end quantum chemistry methods, which further improve the description of electronic correlations, might further alter electron-phonon interactions over GW or DFT. Here, we address these questions by computing the renormalization of the highest occupied molecular orbital energies of Thiel's set of organic molecules by harmonic vibrations using DFT, GW, and equation-of-motion coupled-cluster calculations. We find that, depending on the amount of exact exchange included in the DFT starting point, GW can increase the magnitude of the electron-phonon coupling across Thiel's set of molecules by an average factor of 1.1-1.8 compared to the underlying DFT, while equation-of-motion coupled-cluster leads to an increase of 1.4-2. The electron-phonon coupling predicted with the ab initio GW method is generally in much closer agreement to coupled cluster values compared to DFT, establishing GW as a promising route for accurately computing electron-phonon phenomena in molecules and beyond at a much lower computational cost than higher-end quantum chemistry techniques.

Exciton Transport in the Nonfullerene Acceptor O-IDTBR from Nonadiabatic Molecular Dynamics

Theory, computation, and experiment have given strong evidence that charge carriers in organic molecular crystals form partially delocalized quantum objects that diffuse very efficiently via a mechanism termed transient delocalization. It is currently unclear how prevalent this mechanism is for exciton transport. Here we carry out simulation of singlet Frenkel excitons (FE) in a molecular organic semiconductor that belongs to the class of nonfullerene acceptors, O-IDTBR, using the recently introduced FE surface hopping nonadiabatic molecular dynamics method. We find that FE are, on average, localized on a single molecule in the crystal due to sizable reorganization energy and moderate excitonic couplings. Yet, our simulations suggest that the diffusion mechanism is more complex than simple local hopping; in addition to hopping, we observe frequent transient delocalization events where the exciton wave function expands over 10 or more molecules for a short period of time in response to thermal excitations within the excitonic band, followed by de-excitation and contraction onto a single molecule. The transient delocalization events lead to an increase in the diffusion constant by a factor of 3-4, depending on the crystallographic direction as compared to the situation where only local hopping events are considered. Intriguingly, O-IDTBR appears to be a moderately anisotropic 3D "conductor" for excitons but a highly anisotropic 2D conductor for electrons. Taken together with previous simulation results, two trends seem to emerge for molecular organic crystals: excitons tend to be more localized and slower than charge carriers due to higher internal reorganization energy, while exciton transport tends to be more isotropic than charge transport due to the weaker distance dependence of excitonic versus electronic coupling.

Phonon screening and dissociation of excitons at finite temperatures from first principles

The properties of excitons, or correlated electron-hole pairs, are of paramount importance to optoelectronic applications of materials. A central component of exciton physics is the electron-hole interaction, which is commonly treated as screened solely by electrons within a material. However, nuclear motion can screen this Coulomb interaction as well, with several recent studies developing model approaches for approximating the phonon screening of excitonic properties. While these model approaches tend to improve agreement with experiment, they rely on several approximations that restrict their applicability to a wide range of materials, and thus far they have neglected the effect of finite temperatures. Here, we develop a fully first-principles, parameter-free approach to compute the temperature-dependent effects of phonon screening within the ab initio [Formula: see text]-Bethe-Salpeter equation framework. We recover previously proposed models of phonon screening as well-defined limits of our general framework, and discuss their validity by comparing them against our first-principles results. We develop an efficient computational workflow and apply it to a diverse set of semiconductors, specifically AlN, CdS, GaN, MgO, and [Formula: see text]. We demonstrate under different physical scenarios how excitons may be screened by multiple polar optical or acoustic phonons, how their binding energies can exhibit strong temperature dependence, and the ultrafast timescales on which they dissociate into free electron-hole pairs.

Neural canonical transformations for vibrational spectra of molecules

The behavior of polyatomic molecules around their equilibrium positions can be regarded as that of quantum-coupled anharmonic oscillators. Solving the corresponding Schrödinger equations enables the interpretation or prediction of the experimental spectra of molecules. In this study, we developed a novel approach for solving the excited states of anharmonic vibrational systems. The normal coordinates of the molecules are transformed into new coordinates through a normalizing flow parameterized by a neural network. This facilitates the construction of a set of orthogonal many-body variational wavefunctions. This methodology has been validated on an exactly solvable 64-dimensional coupled harmonic oscillator, yielding numerical results with a relative error of 10-6. The neural canonical transformations are also applied to calculate the energy levels of two specific molecules, acetonitrile (CH3CN) and ethylene oxide (C2H4O). These molecules involve 12 and 15 vibrational modes, respectively. A key advantage of this approach is its flexibility concerning the potential energy surface, as it requires no specific form. Furthermore, this method can be readily implemented on large-scale distributed computing platforms, making it easy to extend to investigating complex vibrational structures.

Evolution of the pentacene exciton band width in pentacene-tetracene blends

Pentacene is one of the most investigated organic semiconductors. It is well known that the motion of excitons in pentacene and other organic semiconductors is determined by inter-molecular exciton coupling based on charge-transfer processes. In the present study, we demonstrate the impact of the admixture of tetracene, which has a larger band gap and interrupts the pentacene-pentacene interaction, on the exciton behavior in pentacene. Using a combination of optical absorption and electron energy-loss spectroscopy, we show that both the Davydov splitting and the exciton band width in pentacene strongly decrease with increasing tetracene concentration, while the decrease of the exciton band width is substantially larger.

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.

Energy-Conserving and Thermally Corrected Neglect of Back-Reaction Approximation Method for Nonadiabatic Molecular Dynamics

In this work, the energy-conserving and thermally corrected neglect of the back-reaction approximation approach for nonadiabatic molecular dynamics in extended atomistic systems is developed. The new approach introduces three key corrections to the original method: (1) it enforces the total energy conservation, (2) it introduces an explicit coupling of the system to its environment, and (3) it introduces a renormalization of nonadiabatic couplings to account for a difference between the instantaneous nuclear kinetic energy and the kinetic energy of guiding trajectories. In the new approach, an auxiliary kinetic energy variable is introduced as an independent dynamical variable. The new approach produces nonzero equilibrium populations, whereas the original neglect of the back-reaction approximation method does not. It yields population relaxation time scales that are favorably comparable to the reference values, and it introduces an explicit and controllable way of dissipating energy into a bath without an assumption of the bath being at equilibrium.

Tunable Interlayer Delocalization of Excitons in Layered Organic-Inorganic Halide Perovskites

Layered organic-inorganic halide perovskites exhibit remarkable structural and chemical diversity and hold great promise for optoelectronic devices. In these materials, excitons are thought to be strongly confined within the inorganic metal halide layers with interlayer coupling generally suppressed by the organic cations. Here, we present an in-depth study of the energy and spatial distribution of the lowest-energy excitons in layered organic-inorganic halide perovskites from first-principles many-body perturbation theory, within the GW approximation and the Bethe-Salpeter equation. We find that the quasiparticle band structures, linear absorption spectra, and exciton binding energies depend strongly on the distance and the alignment of adjacent metal halide perovskite layers. Furthermore, we show that exciton delocalization can be modulated by tuning the interlayer distance and alignment, both parameters determined by the chemical composition and size of the organic cations. Our calculations establish the general intuition needed to engineer excitonic properties in novel halide perovskite nanostructures.

Large Coherent States Formed from Disordered k-Regular Random Graphs

The present work is motivated by the need for robust, large-scale coherent states that can play possible roles as quantum resources. A challenge is that large, complex systems tend to be fragile. However, emergent phenomena in classical systems tend to become more robust with scale. Do these classical systems inspire ways to think about robust quantum networks? This question is studied by characterizing the complex quantum states produced by mapping interactions between a set of qubits from structure in graphs. We focus on maps based on k-regular random graphs where many edges were randomly deleted. We ask how many edge deletions can be tolerated. Surprisingly, it was found that the emergent coherent state characteristic of these graphs was robust to a substantial number of edge deletions. The analysis considers the possible role of the expander property of k-regular random graphs.

Directed exciton transport highways in organic semiconductors

Exciton bandwidths and exciton transport are difficult to control by material design. We showcase the intriguing excitonic properties in an organic semiconductor material with specifically tailored functional groups, in which extremely broad exciton bands in the near-infrared-visible part of the electromagnetic spectrum are observed by electron energy loss spectroscopy and theoretically explained by a close contact between tightly packing molecules and by their strong interactions. This is induced by the donor-acceptor type molecular structure and its resulting crystal packing, which induces a remarkable anisotropy that should lead to a strongly directed transport of excitons. The observations and detailed understanding of the results yield blueprints for the design of molecular structures in which similar molecular features might be used to further explore the tunability of excitonic bands and pave a way for organic materials with strongly enhanced transport and built-in control of the propagation direction.