Bound excitons in time-dependent density-functional theory: optical and energy-loss spectra

Reformulation of Time-Dependent Density Functional Theory for Nonperturbative Dynamics: The Rabi Oscillation Problem Resolved

Rabi oscillations have long been thought to be out of reach in simulations using time-dependent density functional theory (TDDFT), a prominent symptom of the failure of the adiabatic approximation for nonperturbative dynamics. We present a reformulation of TDDFT which requires response quantities only, thus enabling an adiabatic approximation to predict such dynamics accurately because the functional is evaluated on a density close to the ground state, instead of on the fully nonperturbative density. Our reformulation applies to any real-time dynamics, redeeming TDDFT far from equilibrium. Examples of a resonantly-driven local excitation in a model He atom, and charge-transfer in the LiCN molecule are given.

Total Energy beyond GW: Exact Results and Guidelines for Approximations

The total energy and electron addition and removal spectra can, in principle, be obtained exactly from the one-body Green's function (GF). In practice, the GF is obtained from an approximate self-energy. In the framework of many-body perturbation theory, we derive different expressions that are based on an approximate self-energy, but that yield nevertheless, in principle, the exact exchange-correlation contribution to the total energy for any interaction strength. Response functions play a crucial role, which explains why, for example, ingredients from time-dependent density functional theory can be used to build these approximate self-energies. We show that the key requirement for obtaining exact results is the consistent combination of ingredients. Also when further approximations are made, as it is necessary in practice, this consistency remains the key to obtain good results. All findings are illustrated using the exactly solvable symmetric Hubbard dimer.

Non-adiabatic Exciton Dynamics in van der Waals Heterostructures

In this Perspective, we introduce a first-principles method that combines time-dependent density functional theory with non-adiabatic molecular dynamics (NAMD) to explore exciton dynamics in two-dimensional (2D) van der Waals (vdW) heterostructures. The theoretical foundation and computational efficiency of the method are discussed and compared with those of related methods (e.g., GW-BSE). Using three 2D vdW heterostructures as examples, we demonstrate that the proposed method can provide a reliable description of many-body electron-hole interactions crucial to exciton dynamics. With much lower computational costs than the GW-BSE method, the proposed method represents a particularly promising theoretical tool to probe exciton dynamics in solids. Moreover, we find that the NAMD simulations widely used in the literature cannot capture the excitonic effect in 2D materials and often yield incorrect results because they are formulated in a single-particle picture. The instances where the single-particle picture fails are pointed out and contrasted with the many-body simulation results.

Real-Time Exciton Dynamics with Time-Dependent Density-Functional Theory

Linear-response time-dependent density-functional theory (TDDFT) can describe excitonic features in the optical spectra of insulators and semiconductors, using exchange-correlation (xc) kernels behaving as -1/k^{2} to leading order. We show how excitons can be modeled in real-time TDDFT, using an xc vector potential constructed from approximate, long-range corrected xc kernels. We demonstrate, for various materials, that this real-time approach is consistent with frequency-dependent linear response, gives access to femtosecond exciton dynamics following short-pulse excitations, and can be extended with some caution into the nonlinear regime.

Materials Design and Optimization for Next-Generation Solar Cell and Light-Emitting Technologies

We review some of the most potent directions in the design of materials for next-generation solar cell and light-emitting technologies that go beyond traditional solid-state inorganic semiconductor-based devices, from both the experimental and computational standpoints. We focus on selected recent conceptual advances in tackling issues which are expected to significantly impact applied literature in the coming years. Specifically, we consider solution processability, design of dopant-free charge transport materials, two-dimensional conjugated polymeric semiconductors, and colloidal quantum dot assemblies in the fields of experimental synthesis, characterization, and device fabrication. Key modeling issues that we consider are calculations of optical properties and of effects of aggregation, including recent advances in methods beyond linear-response time-dependent density functional theory and recent insights into the effects of correlation when going beyond the single-particle ansatz as well as in the context of modeling of thermally activated fluorescence.

Advanced First-Principle Modeling of Relativistic Ruddlesden—Popper Strontium Iridates

In this review, we provide a survey of the application of advanced first-principle methods on the theoretical modeling and understanding of novel electronic, optical, and magnetic properties of the spin-orbit coupled Ruddlesden–Popper series of iridates Srn+1IrnO3n+1 (n = 1, 2, and ∞). After a brief description of the basic aspects of the adopted methods (noncollinear local spin density approximation plus an on-site Coulomb interaction (LSDA+U), constrained random phase approximation (cRPA), GW, and Bethe–Salpeter equation (BSE)), we present and discuss select results. We show that a detailed phase diagrams of the metal–insulator transition and magnetic phase transition can be constructed by inspecting the evolution of electronic and magnetic properties as a function of Hubbard U, spin–orbit coupling (SOC) strength, and dimensionality n, which provide clear evidence for the crucial role played by SOC and U in establishing a relativistic (Dirac) Mott–Hubbard insulating state in Sr2IrO4 and Sr3Ir2O7. To characterize the ground-state phases, we quantify the most relevant energy scales fully ab initio—crystal field energy, Hubbard U, and SOC constant of three compounds—and discuss the quasiparticle band structures in detail by comparing GW and LSDA+U data. We examine the different magnetic ground states of structurally similar n = 1 and n = 2 compounds and clarify that the origin of the in-plane canted antiferromagnetic (AFM) state of Sr2IrO4 arises from competition between isotropic exchange and Dzyaloshinskii–Moriya (DM) interactions whereas the collinear AFM state of Sr3Ir2O7 is due to strong interlayer magnetic coupling. Finally, we report the dimensionality controlled metal–insulator transition across the series by computing their optical transitions and conductivity spectra at the GW+BSE level from the the quasi two-dimensional insulating n = 1 and 2 phases to the three-dimensional metallic n=∞ phase.

Ultrafast charge ordering by self-amplified exciton-phonon dynamics in TiSe2

The origin of charge density waves (CDWs) in TiSe\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{2}$$\end{document}2 has long been debated, mainly due to the difficulties in identifying the timescales of the excitonic pairing and electron–phonon coupling (EPC). Without a time-resolved and microscopic mechanism, one has to assume simultaneous appearance of CDW and periodic lattice distortions (PLD). Here, we accomplish a complete separation of ultrafast exciton and PLD dynamics and unravel their interplay in our real-time time-dependent density functional theory simulations. We find that laser pulses knock off the exciton order and induce a homogeneous bonding–antibonding transition in the initial 20 fs, then the weakened electronic order triggers ionic movements antiparallel to the original PLD. The EPC comes into play after the initial 20 fs, and the two processes mutually amplify each other leading to a complete inversion of CDW ordering. The self-amplified dynamics reproduces the evolution of band structures in agreement with photoemission experiments. Hence we resolve the key processes in the initial dynamics of CDWs that help elucidate the underlying mechanism.

The physical origins of charge density waves in 1T-TiSe2 and their response to ultrafast excitation have long been a topic of theoretical and experimental debate. Here the authors present an ab initio theory that successfully captures the observed dynamics of charge density wave formation.

Excitons in two-dimensional atomic layer materials from time-dependent density functional theory: mono-layer and bi-layer hexagonal boron nitride and transition-metal dichalcogenides

Time-dependent density functional theory has been applied to the calculation of absorption spectra for two dimensional atomic layer materials: mono-layer and bi-layer hexagonal boron nitride and mono-layer transition metal dichalcogenides.

Design of auxiliary systems for spectroscopy

In this contribution, we advocate the possibility of designing auxiliary systems with effective potentials or kernels that target only the specific spectral properties of interest and are simpler than the self-energy of many-body perturbation theory or the exchange–correlation kernel of time-dependent density-functional theory.

Plasmonic performance of AuxAgyCu1-x-y alloys from many-body perturbation theory

We present a detailed appraisal of the optical and plasmonic properties of ordered alloys of the form Au_xAg_yCu_1-x-y, as predicted by means of first-principles many-body perturbation theory augmented by a semi-empirical Drude-Lorentz model. In benchmark simulations on elemental Au, Ag, and Cu, we find that the random-phase approximation (RPA) fails to accurately describe inter-band transitions when it is built upon semi-local approximate Kohn-Sham density-functional theory band-structures. We show that non-local electronic exchange-correlation interactions sufficient to correct this, particularly for the fully-filled, relatively narrow d-bands which contribute strongly throughout the low-energy spectral range (0-6 eV), may be modeled very expediently using band-stretching operators that imitate the effect of a perturbative [Formula: see text] self-energy correction incorporating quasiparticle (QP) mass renormalization. We thereby establish a convenient work-flow for carrying out approximated [Formula: see text] spectroscopic calculations on alloys and, in particular here, we have considered alloy concentrations down to 12.5% in [Formula: see text], including all possible crystallographic orderings of face-centred cubic type. We develop a pragmatic procedure for calculating the Drude plasmon frequency from first principles, including self-energy effects, as well as a semi-empirical scheme for interpolating the plasmon inverse lifetimes between stoichiometries. A distinctive M-shaped profile is observed in both quantities for binary alloys, in qualitative agreement with previous experimental findings. A range of optical and plasmonic figures of merit are discussed, and plotted for ordered [Formula: see text] at three representative solid-state laser wavelengths. On this basis, we predict that certain compositions may offer improved performance over elemental Au for particular application types. We predict that while the loss functions for both bulk and surface plasmons are typically diminished in strength through binary alloying, certain stoichiometric ratios may exhibit higher-quality (longer-lived) localized surface-plasmons and surface-plasmon polaritons, at technologically-relevant wavelengths, than those in elemental Au.

Many-body simulation of two-dimensional electronic spectroscopy of excitons and trions in monolayer transition metal dichalcogenides

Indications of coherently interacting excitons and trions in doped transition metal dichalcogenides have been measured as quantum beats in two-dimensional electronic spectroscopy, but the microscopic principles underlying the optical signals of exciton-trion coherence remain to be clarified. Here we present calculations of two-dimensional spectra of such monolayers based on a microscopic many-body formalism. We use a parameterized band structure and a static model dielectric function, although a full ab initio implementation of our formalism is possible in principle. Our simulated spectra are in excellent agreement with experiments, including the quantum beats, while revealing the interplay between excitons and trions in molybdenum- and tungsten-based transition metal dichalcogenides. Quantum beats are confirmed to unambiguously reflect the exciton-trion coherence time in molybdenum compounds, but are shown to provide a lower bound to the coherence time for tungsten analogues due to a destructive interference from coexisting singlet and triplet trions.

2D electronic spectroscopy found experimental indications of coherently interacting excitons and trions in doped transition metal dichalcogenides (TMDCs). Here, the authors perform simulations of 2D spectra of monolayer TMDCs based on a many-body formalism, allowing to relate exciton-trion coherence to quantum beats based on microscopic principles.

Finite-Field Approach to Solving the Bethe-Salpeter Equation

We present a method to compute optical spectra and exciton binding energies of molecules and solids based on the solution of the Bethe-Salpeter equation and the calculation of the screened Coulomb interaction in a finite field. The method does not require either the explicit evaluation of dielectric matrices or of virtual electronic states, and can be easily applied without resorting to the random phase approximation. In addition, it utilizes localized orbitals obtained from Bloch states using bisection techniques, thus greatly reducing the complexity of the calculation and enabling the efficient use of hybrid functionals to obtain single particle wave functions. We report exciton binding energies of several molecules and absorption spectra of condensed systems of unprecedented size, including water and ice samples with hundreds of atoms.

Performance of polarisation functionals for linear and nonlinear optical properties of bulk zinc chalcogenides ZnX (X = S, Se, and Te)

We calculated the frequency dependent macroscopic dielectric function and second-harmonic generation of cubic ZnS, ZnSe and ZnTe within time-dependent density-polarisation functional theory.

Fully Parameter-Free Calculation of Optical Spectra for Insulators, Semiconductors, and Metals from a Simple Polarization Functional

We present a fully parameter-free density-functional approach for the accurate description of optical absorption spectra of insulators, semiconductors, and metals. We show that this can be achieved within time-dependent current-density-functional theory using a simple dynamical polarization functional. We derive this functional from physical principles that govern optical spectra. Our method is truly predictive because not a single parameter is used. In particular, we do not use an ad hoc material-dependent broadening parameter to compare theory to experiment as is usually done. Our approach is numerically efficient; the cost equals that of a calculation within the random-phase approximation.

A systematic benchmark of the ab initio Bethe-Salpeter equation approach for low-lying optical excitations of small organic molecules

The predictive power of the ab initio Bethe-Salpeter equation (BSE) approach, rigorously based on many-body Green's function theory but incorporating information from density functional theory, has already been demonstrated for the optical gaps and spectra of solid-state systems. Interest in photoactive hybrid organic/inorganic systems has recently increased and so has the use of the BSE for computing neutral excitations of organic molecules. However, no systematic benchmarks of the BSE for neutral electronic excitations of organic molecules exist. Here, we study the performance of the BSE for the 28 small molecules in Thiel's widely used time-dependent density functional theory benchmark set [Schreiber et al., J. Chem. Phys. 128, 134110 (2008)]. We observe that the BSE produces results that depend critically on the mean-field starting point employed in the perturbative approach. We find that this starting point dependence is mainly introduced through the quasiparticle energies obtained at the intermediate GW step and that with a judicious choice of starting mean-field, singlet excitation energies obtained from BSE are in excellent quantitative agreement with higher-level wavefunction methods. The quality of the triplet excitations is slightly less satisfactory.

Estimating excitonic effects in the absorption spectra of solids: problems and insight from a guided iteration scheme

A major obstacle for computing optical spectra of solids is the lack of reliable approximations for capturing excitonic effects within time-dependent density functional theory. We show that the accurate prediction of strongly bound electron-hole pairs within this framework using simple approximations is still a challenge and that available promising results have to be revisited. Deriving a set of analytical formulas we analyze and explain the difficulties. We deduce an alternative approximation from an iterative scheme guided by previously available knowledge, significantly improving the description of exciton binding energies. Finally, we show how one can "read" exciton binding energies from spectra determined in the random phase approximation, without any further calculation.

Exciting: a full-potential all-electron package implementing density-functional theory and many-body perturbation theory

Linearized augmented planewave methods are known as the most precise numerical schemes for solving the Kohn-Sham equations of density-functional theory (DFT). In this review, we describe how this method is realized in the all-electron full-potential computer package, exciting. We emphasize the variety of different related basis sets, subsumed as (linearized) augmented planewave plus local orbital methods, discussing their pros and cons and we show that extremely high accuracy (microhartrees) can be achieved if the basis is chosen carefully. As the name of the code suggests, exciting is not restricted to ground-state calculations, but has a major focus on excited-state properties. It includes time-dependent DFT in the linear-response regime with various static and dynamical exchange-correlation kernels. These are preferably used to compute optical and electron-loss spectra for metals, molecules and semiconductors with weak electron-hole interactions. exciting makes use of many-body perturbation theory for charged and neutral excitations. To obtain the quasi-particle band structure, the GW approach is implemented in the single-shot approximation, known as G(0)W(0). Optical absorption spectra for valence and core excitations are handled by the solution of the Bethe-Salpeter equation, which allows for the description of strongly bound excitons. Besides these aspects concerning methodology, we demonstrate the broad range of possible applications by prototypical examples, comprising elastic properties, phonons, thermal-expansion coefficients, dielectric tensors and loss functions, magneto-optical Kerr effect, core-level spectra and more.

First-principles method of propagation of tightly bound excitons: verifying the exciton band structure of LiF with inelastic x-ray scattering

We propose a simple first-principles method to describe the propagation of tightly bound excitons. By viewing the exciton as a composite object (an effective Frenkel exciton in Wannier orbitals), we define an exciton kinetic kernel to encapsulate the exciton propagation and decay for all binding energies. Applied to prototypical LiF, our approach produces three exciton bands, which we verified quantitatively via inelastic x-ray scattering. The proposed real-space picture is computationally inexpensive and thus enables study of the full exciton dynamics, even in the presence of surfaces and impurity scattering. It also provides an intuitive understanding to facilitate practical exciton engineering in semiconductors, strongly correlated oxides, and their nanostructures.

A minimal model for excitons within time-dependent density-functional theory

The accurate description of the optical spectra of insulators and semiconductors remains an important challenge for time-dependent density-functional theory (TDDFT). Evidence has been given in the literature that TDDFT can produce bound as well as continuum excitons for specific systems, but there are still many unresolved basic questions concerning the role of dynamical exchange and correlation (xc). In particular, the roles of the long spatial range and the frequency dependence of the xc kernel f(xc) for excitonic binding are still not very well explored. We present a minimal model for excitons in TDDFT, consisting of two bands from a one-dimensional (1D) Kronig-Penney model and simple approximate xc kernels, providing an easily accessible model system for studying excitonic effects in TDDFT. For the 1D model system, it is found that adiabatic xc kernels can produce at most two bound excitons, confirming that the long spatial range of f(xc) is not a necessary condition. It is shown how the Wannier model, featuring an effective electron-hole interaction, emerges from TDDFT. The collective, many-body nature of excitons is explicitly demonstrated.

Van der Waals interactions in ionic and semiconductor solids

van der Waals (vdW) energy corrected density-functional theory [Phys. Rev. Lett. 102, 073005 (2009)] is applied to study the cohesive properties of ionic and semiconductor solids (C, Si, Ge, GaAs, NaCl, and MgO). The required polarizability and dispersion coefficients are calculated using the dielectric function obtained from time-dependent density-functional theory. Coefficients for "atoms in the solid" are then calculated from the Hirshfeld partitioning of the electron density. It is shown that the Clausius-Mossotti equation that relates the polarizability and the dielectric function is accurate even for covalently-bonded semiconductors. We find an overall improvement in the cohesive properties of Si, Ge, GaAs, NaCl, and MgO, when vdW interactions are included on top of the Perdew-Burke-Ernzerhof or Heyd-Scuseria-Ernzerhof functionals. The relevance of our findings for other solids is discussed.

Bootstrap approximation for the exchange-correlation kernel of time-dependent density-functional theory

A new parameter-free approximation for the exchange-correlation kernel f(xc) of time-dependent density-functional theory is proposed. This kernel is expressed as an algorithm in which the exact Dyson equation for the response, as well as an approximate expression for f(xc) in terms of the dielectric function, are solved together self-consistently, leading to a simple parameter-free kernel. We apply this to the calculation of optical spectra for various small band gap (Ge, Si, GaAs, AlN, TiO(2), SiC), large band gap (C, LiF, Ar, Ne), and magnetic (NiO) insulators. The calculated spectra are in very good agreement with the experiment for this diverse set of materials, highlighting the universal applicability of the new kernel.

Design of effective kernels for spectroscopy and molecular transport: time-dependent current-density-functional theory

Time-dependent current-density-functional theory (TDCDFT) provides an, in principle, exact scheme to calculate efficiently response functions for a very broad range of applications. However, the lack of approximations valid for a range of parameters met in experimental conditions has so far delayed its extensive use in inhomogeneous systems. On the other side, in many-body perturbation theory accurate approximations are available, but at a price of a higher computational cost. In the present work, the possibility of combining the advantages of both approaches is exploited. In this way, an exact equation for the exchange-correlation kernel of TDCDFT is obtained, which opens the way for a systematic improvement of the approximations adopted in practical applications. Finally, an approximate kernel for an efficient calculation of spectra of solids and molecular conductances is suggested and its validity is discussed.

Strong charge-transfer excitonic effects and the Bose-Einstein exciton condensate in graphane

Using first principles many-body theory methods (GW+Bethe-Salpeter equation) we demonstrate that the optical properties of graphane are dominated by localized charge-transfer excitations governed by enhanced electron correlations in a two-dimensional dielectric medium. Strong electron-hole interaction leads to the appearance of small radius bound excitons with spatially separated electron and hole, which are localized out of plane and in plane, respectively. The presence of such bound excitons opens the path towards an excitonic Bose-Einstein condensate in graphane that can be observed experimentally.

First-principles calculation of the electron dynamics in crystalline SiO(2)

We present a first-principles description for electron dynamics in crystalline SiO(2) induced by an optical field in both weak and intense regimes. We rely upon the time-dependent density-functional theory with the adiabatic local-density approximation, and a real-space and real-time method is employed to solve the time-dependent Kohn-Sham equation. The response calculation to a weak field provides us with information on the dielectric function, while the response to an intense field shows the optical dielectric breakdown. We discuss the critical threshold for the dielectric breakdown of crystalline SiO(2), in comparison with the results for diamond.

Excited state properties of liquid water

In this paper, we give an overview of the state of the art in calculations of the electronic band structure and absorption spectra of water. After an introduction to the main theoretical and computational schemes used, we present results for the electronic and optical excitations of water. We focus mainly on liquid water, but spectroscopic properties of ice and vapor phase are also described. The applicability and the accuracy of first-principles methods are discussed, and results are critically presented.

Optical saturation driven by exciton confinement in molecular chains: a time-dependent density-functional theory approach

We identify excitonic confinement in one-dimensional molecular chains (i.e., polyacetylene and H2) as the main driving force for the saturation of the chain polarizability as a function of the number of molecular units. This conclusion is based on first principles time-dependent density-functional theory calculations using a recently developed exchange-correlation kernel that accounts for excitonic effects. The failure of simple local and semilocal functionals is shown to be linked to the lack of memory effects, spatial ultranonlocality, and self-interaction corrections. These effects get smaller as the gap reduces, in which case such simple approximations do perform better.

Generalized Hedin's equations for quantum many-body systems with spin-dependent interactions

Hedin's equations for the electron self-energy and the vertex have originally been derived for a many-electron system with Coulomb interaction. In recent years, it has been increasingly recognized that spin interactions can play a major role in determining physical properties of systems such as nanoscale magnets or of interfaces and surfaces. We derive a generalized set of Hedin's equations for quantum many-body systems containing spin interactions, e.g., spin-orbit and spin-spin interactions. The corresponding spin-dependent GW approximation is constructed.

Time-dependent density functional theory: past, present, and future

Time-dependent density functional theory (TDDFT) is presently enjoying enormous popularity in quantum chemistry, as a useful tool for extracting electronic excited state energies. This article discusses how TDDFT is much broader in scope, and yields predictions for many more properties. We discuss some of the challenges involved in making accurate predictions for these properties.

Ab initio calculation of optical spectra of liquids: many-body effects in the electronic excitations of water

We present ab initio calculations of the excited state properties of liquid water in the framework of many-body Green's function formalism. Snapshots taken from molecular dynamics simulations are used as input geometries to calculate electronic and optical spectra, and the results are averaged over the different configurations. The optical absorption spectra with the inclusion of excitonic effects are calculated by solving the Bethe-Salpeter equation. The insensitivity of screening effects to a particular configuration make these calculations feasible. The resulting spectra, which are strongly modified by many-body effects, are in good agreement with experiments.