Publications

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.

Heterostructures of layered transition metal dichalcogenides (TMDs) host long-lived, tunable excitons, making them intriguing candidates for material-based quantum information applications. Light absorption in these systems induces a plethora of optically excited states that hybridize both interlayer and intralayer characteristics, providing a distinctive starting point for their relaxation processes, in which the interplay between generated electron-hole pairs and their scattering with phonons play a key role. We present a first-principles theoretical approach to compute phonon-induced exciton decomposition due to rapid occupation of electron-hole pairs with finite momentum and opposite spin. Using the MoSe2/WSe2 heterostructure as a case study, we observe a reduction in the optical activity of bright states upon phonon scattering already in the first few femtoseconds after a photoexcitation, driving exciton interlayer delocalization and subsequent variations in the exciton spin. Our results reveal an unexpected and previously unexplored starting point for exciton relaxation dynamics, suggesting increased availability for coherent interactions and nonradiative processes through ultrafast changes in exciton momentum, spatial, and spin properties upon light excitation.

We theoretically study the contribution of quantum effects to the exciton diffusion coefficient in atomically thin crystals. It is related to the weak localization caused by the interference of excitonic wavefunctions on the trajectories with closed loops. Due to the weak inelasticity of the exciton-phonon interaction, the effect is present even if the excitons are scattered by long-wavelength acoustic phonons. We consider exciton interaction with longitudinal acoustic phonons with linear dispersion and flexural phonons with quadratic dispersion. We identify the regimes where the weak localization effect can be particularly pronounced. We also briefly address the role of free charge carriers in the exciton quantum transport and, within the self-consistent theory of localization, the weak localization effects beyond the lowest order.

Excited-state processes at organic-inorganic interfaces consisting of molecular crystals are essential in energy conversion applications. While advances in experimental methods allow direct observation and detection of exciton transfer across such junctions, a detailed understanding of the underlying excitonic properties due to crystal packing and interface structure is still largely lacking. In this work, we use many-body perturbation theory to study structure-property relations of excitons in molecular crystals upon adsorption on a gold surface. We explore the case of the experimentally-studied octyl perylene diimide (C8-PDI) as a prototypical system, and use the GW and Bethe-Salpeter equation (BSE) approach to quantify the change in quasiparticle and exciton properties due to intermolecular and substrate screening. Our findings provide a close inspection of both local and environmental structural effects dominating the excitation energies and the exciton binding and nature, as well as their modulation upon the metal-organic interface composition.

We determine the phase diagram of excitons in a symmetric transition-metal dichalcogenide 3-layer heterostructure. First principles calculations reveal interlayer exciton states of a symmetric quadrupole, from which higher energy asymmetric dipole states are composed. We find quantum phase transitions between a repulsive quadrupolar and an attractive staggered dipolar lattice phases, driven by a competition between interactions and single exciton energies. The different internal quantum state of excitons in each phase is a striking example of a system where single-particle and interacting many-body states are coupled.

Control of impurity concentrations in semiconducting materials is essential to device technology. Because of their intrinsic confinement, the properties of two-dimensional semiconductors such as transition metal dichalcogenides (TMDs) are more sensitive to defects than traditional bulk materials. The technological adoption of TMDs is dependent on the mitigation of deleterious defects and guided incorporation of functional foreign atoms. The first step toward impurity control is the identification of defects and assessment of their electronic properties. Here, we present a comprehensive study of point defects in monolayer tungsten disulfide (WS2) grown by chemical vapor deposition using scanning tunneling microscopy/spectroscopy, CO-tip noncontact atomic force microscopy, Kelvin probe force spectroscopy, density functional theory, and tight-binding calculations. We observe four different substitutional defects: chromium (Cr-W) and molybdenum (Mo-W) at a tungsten site, oxygen at sulfur sites in both top and bottom layers (O-S top/bottom), and two negatively charged defects (CD type I and CD type II). Their electronic fingerprints unambiguously corroborate the defect assignment and reveal the presence or absence of in-gap defect states. Cr-W forms three deep unoccupied defect states, two of which arise from spin-orbit splitting. The formation of such localized trap states for Cr-W differs from the Mo-W case and can be explained by their different d shell energetics and local strain, which we directly measured. Utilizing a tight-binding model the electronic spectra of the isolectronic substitutions O-S and Cr-W are mimicked in the limit of a zero hopping term and infinite on-site energy at a S and W site, respectively. The abundant CDs are negatively charged, which leads to a significant band bending around the defect and a local increase of the contact potential difference. In addition, CD-rich domains larger than 100 nm are observed, causing a work function increase of 1.1 V. While most defects are electronically isolated, we also observed hybrid states formed between Cr-W dimers. The important role of charge localization, spin-orbit coupling, and strain for the formation of deep defect states observed at substitutional defects in WS2 as reported here will guide future efforts of targeted defect engineering and doping of TMDs.

We study the effect of point-defect chalcogen vacancies on the optical properties of monolayer transition metal dichalcogenides using ab initio GW and Bethe-Salpeter equation calculations. We find that chalcogen vacancies introduce unoccupied in-gap states and occupied resonant defect states within the quasiparticle continuum of the valence band. These defect states give rise to a number of strongly bound defect excitons and hybridize with excitons of the pristine system, reducing the valley-selective circular dichroism. Our results suggest a pathway to tune spin-valley polarization and other optical properties through defect engineering.

We develop a new first-principles approach to predict and understand rates of singlet fission with an ab initio Green's-function formalism based on many-body perturbation theory. Starting with singlet and triplet excitons computed from a GW plus Bethe-Salpeter equation approach, we calculate the exciton-biexciton coupling to lowest order in the Coulomb interaction, assuming a final state consisting of two noninteracting spin-correlated triplets with finite center-of-mass momentum. For crystalline pentacene, symmetries dictate that the only purely Coulombic fission decay process from a bright singlet state requires a final state consisting of two inequivalent nearly degenerate triplets of nonzero, equal and opposite, center-of-mass momenta. For such a process, we predict a singlet lifetime of 30-70 fs, in very good agreement with experimental data, indicating that this process can dominate singlet fission in crystalline pentacene. Our approach is general and provides a framework for predicting and understanding multiexciton interactions in solids.

Charge migration for electron transfer via the polypeptide matrix of proteins is a key process in biological energy conversion and signaling systems. It is sensitive to the sequence of amino acids composing the protein and, therefore, offers a tool for chemical control of charge transport across biomaterial-based devices. We designed a series of linear oligoalanine peptides with a single tryptophan substitution that acts as a "dopant," introducing an energy level closer to the electrodes' Fermi level than that of the alanine homopeptide. We investigated the solid-state electron transport (ETp) across a selfassembled monolayer of these peptides between gold contacts. The single tryptophan "doping" markedly increased the conductance of the peptide chain, especially when its location in the sequence is close to the electrodes. Combining inelastic tunneling spectroscopy, UV photoelectron spectroscopy, electronic structure calculations by advanced density-functional theory, and dc current-voltage analysis, the role of tryptophan in ETp is rationalized by charge tunneling across a heterogeneous energy barrier, via electronic states of alanine and tryptophan, and by relatively efficient direct coupling of tryptophan to a Au electrode. These results reveal a controlled way of modulating the electrical properties of molecular junctions by tailormade "building block" peptides.

Many novel applications in bioelectronics rely on the interaction between biomolecules and electronically conducting substrates. However, crucial knowledge about the relation between electronic transport via peptides and their amino-acid composition is still absent. Here, we report results of electronic transport measurements via several homopeptides as a function of their structural properties and temperature. We demonstrate that the conduction through the peptide depends on its length and secondary structure as well as on the nature of the constituent amino acid and charge of its residue. We support our experimental observations with high-level electronic structure calculations and suggest off-resonance tunneling as the dominant conduction mechanism via extended peptides. Our findings indicate that both peptide composition and structure can affect the efficiency of electronic transport across peptides.

A self-consistent optimally tuned range-separated hybrid density functional (scOT-RSH) approach is developed. It can simultaneously predict accurate geometries, vibrational modes, and frontier orbital energies. This is achieved by optimizing the range-separation parameter, γ to both satisfy the ionization energy theorem and minimize interatomic forces. We benchmark our approach against an established hybrid functional, B3LYP, using the G2 test set. scOT-RSH greatly improves the accuracy of occupied frontier orbital energies, with a mean absolute error (MAE) of only 0.2 eV relative to experimental ionization energies compared to 2.96 eV with B3LYP. Geometries do not change significantly compared to those obtained from B3LYP, with a bond length MAE of 0.012 Å compared to 0.008 Å for B3LYP, and a 6.5% MAE for zero-point energies, slightly larger than that of B3LYP (3.1%). scOT-RSH represents a new paradigm in which accurate geometries and ionization energies can be predicted simultaneously from a single functional approach.

We present a method for obtaining outer-valence quasiparticle excitation energies from a density-functional-theory-based calculation, with an accuracy that is comparable to that of many-body perturbation theory within the GW approximation. The approach uses a range-separated hybrid density functional, with an asymptotically exact and short-range fractional Fock exchange. The functional contains two parameters, the range separation and the short-range Fock fraction. Both are determined nonempirically, per system, on the basis of the satisfaction of exact physical constraints for the ionization potential and frontier-orbital many-electron self-interaction, respectively. The accuracy of the method is demonstrated on four important benchmark organic molecules: perylene, pentacene, 3,4,9,10-perylene-tetracarboxylic-dianydride (PTCDA), and 1,4,5,8-naphthalene-tetracarboxylic-dianhydride (NTCDA). We envision that for the outer-valence excitation spectra of finite systems the approach could provide an inexpensive alternative to GW, opening the door to the study of presently out of reach large-scale systems.