The ab initio Bethe-Salpeter equation (BSE) approach, an established method for the study of excitons in materials, is typically solved in a limit where only static screening from electrons is captured. Here, we generalize this framework to include dynamical screening from phonons at lowest order in the electron-phonon interaction. We apply this generalized BSE approach to a series of inorganic lead halide perovskites, CsPbX_{3}, with X=Cl, Br, and I. We find that inclusion of screening from phonons significantly reduces the computed exciton binding energies of these systems. By deriving a simple expression for phonon screening effects, we reveal general trends for their importance in semiconductors and insulators, based on a hydrogenic exciton model. We demonstrate that the magnitude of the phonon screening correction in isotropic materials can be reliably predicted using four material specific parameters: the reduced effective mass, static and optical dielectric constants, and frequency of the most strongly coupled longitudinal-optical phonon mode. This framework helps to elucidate the importance of phonon screening and its relation to excitonic properties in a broad class of semiconductors.

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.

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.

The dependence of ab initio many-body perturbation theory within the GW approximation on the eigensystem used in calculating quasiparticle corrections limits this method's predictive power. Here, we investigate the accuracy of the recently developed Wannier-localized optimally tuned screened range-separated hybrid (WOT-SRSH) functional as a generalized Kohn-Sham starting point for single-shot GW (G0W0) calculations for a range of semiconductors and insulators. Comparison to calculations based on well-established functionals, namely, PBE, PBE0, and HSE, as well as to self-consistent GW schemes and to experiment, shows that band gaps computed via G0W0@WOT-SRSH have a level of precision and accuracy that is comparable to that of more advanced methods such as quasiparticle self-consistent GW and eigenvalue self-consistent GW. We also find that G0W0@WOT-SRSH improves the description of states deeper in the valence band manifold. Finally, we show that G0W0@WOT-SRSH significantly reduces the sensitivity of computed band gaps to ambiguities in the underlying WOT-SRSH tuning procedure.

Accurate prediction of fundamental band gaps of crystalline solid-state systems entirely within density functional theory is a long-standing challenge. Here, we present a simple and inexpensive method that achieves this by means of nonempirical optimal tuning of the parameters of a screened range-separated hybrid functional. The tuning involves the enforcement of an ansatz that generalizes the ionization potential theorem to the removal of an electron from an occupied state described by a localized Wannier function in a modestly sized supercell calculation. The method is benchmarked against experiment for a set of systems ranging from narrow band-gap semiconductors to large band-gap insulators, spanning a range of fundamental band gaps from 0.2 to 14.2 electronvolts (eV), and is found to yield quantitative accuracy across the board, with a mean absolute error of ∼0.1 eV and a maximal error of ∼0.2 eV.

The intrinsic weak and highly nonlocal dielectric screening of two-dimensional materials is well-known to lead to high sensitivity of their optoelectronic properties to environment. Less studied theoretically is the role of free carriers in those properties. Here, we use ab initio GW and Bethe-Salpeter equation calculations, with a rigorous treatment of dynamical screening and local-field effects, to study the doping dependence of the quasiparticle and optical properties of a monolayer transition-metal dichalcogenide, 2H MoTe₂. We predict a quasiparticle band gap renormalization of several hundreds of meV for experimentally attainable carrier densities and a similarly sizable decrease in the exciton binding energy. This results in an almost constant excitation energy for the lowest-energy exciton resonance with an increasing doping density. Using a newly developed and generally applicable plasmon-pole model and a self-consistent solution of the Bethe-Salpeter equation, we reveal the importance of accurately capturing both dynamical and local-field effects to understand detailed photoluminescence measurements.

Atomic spin centers in 2D materials are a highly anticipated building block for quantum technologies. Here, we demonstrate the creation of an effective spin-1/2 system via the atomically controlled generation of magnetic carbon radical ions (CRIs) in synthetic two-dimensional transition metal dichalcogenides. Hydrogenated carbon impurities located at chalcogen sites introduced by chemical doping are activated with atomic precision by hydrogen depassivation using a scanning probe tip. In its anionic state, the carbon impurity is computed to have a magnetic moment of 1 μ_B resulting from an unpaired electron populating a spin-polarized in-gap orbital. We show that the CRI defect states couple to a small number of local vibrational modes. The vibronic coupling strength critically depends on the spin state and differs for monolayer and bilayer WS₂. The carbon radical ion is a surface-bound atomic defect that can be selectively introduced, features a well-understood vibronic spectrum, and is charge state controlled.

We present band structure and optical absorption spectra obtained from density functional theory (DFT) and linear response time-dependent DFT (TDDFT) calculations using a screened range-separated hybrid (SRSH) functional, including spin-orbit coupling, for seven prototypical semiconductors. The results are compared to those obtained from highly converged many-body perturbation theory calculations using the GW approximation and the GW plus Bethe-Salpeter equation (GW-BSE) approaches. We use a single empirical parameter for our SRSH calculations, fit such that the SRSH band gap reproduces the GW band gap at the Γ point. We then find that ground-state generalized Kohn-Sham SRSH eigenvalues accurately reproduce the band structure obtained from GW calculations, typically to within 0.1-0.2 eV, and optical absorption spectra obtained using TDDFT with the SRSH functional agree well with those of GW-BSE, with a mean deviation of 0.03 and 0.11 eV for the location of the first and second absorption peaks, respectively, at a fraction of the computational cost.

Despite the weak, van der Waals interlayer coupling, photoinduced charge transfer vertically across atomically thin interfaces can occur within surprisingly fast, sub-50 fs time scales. An early theoretical understanding of charge transfer is based on a noninteracting picture, neglecting excitonic effects that dominate optical properties of such materials. We employ an ab initio many-body perturbation theory approach, which explicitly accounts for the excitons and phonons in the heterostructure. Our large-scale first-principles calculations directly probe the role of exciton-phonon coupling in the charge dynamics of the WS₂/MoS₂ heterobilayer. We find that the exciton-phonon interaction induced relaxation time of photoexcited excitons at the K valley of MoS₂ and WS₂ is 67 and 15 fs at 300 K, respectively, which sets a lower bound to the intralayer-to-interlayer exciton transfer time and is consistent with experiment reports. We further show that electron-hole correlations facilitate novel transfer pathways that are otherwise inaccessible to noninteracting electrons and holes.