Excitonic optical spectrum of semiconductors obtained by time-dependent density-functional theory with the exact-exchange kernel

Adiabatic connection interaction strength interpolation method made accurate for the uniform electron gas

The adiabatic connection interaction strength interpolation (ISI)-like method provides a high-level expression for the correlation energy, being, in principle, exact not only in the weak-interaction limit, where it recovers the second-order Görling-Levy perturbation term, but also in the strong-interaction limit that is described by the strictly correlated electron approach. In this work, we construct a genISI functional made accurate for the uniform electron gas, a solid-state physics paradigm that is a very difficult test for ISI-like correlation functionals. We assess the genISI functional for various jellium spheres with the number of electrons Z ≤ 912 and for the non-relativistic noble atoms with Z ≤ 290. For the jellium clusters, the genISI is remarkably accurate, while for the noble atoms, it shows a good performance, similar to other ISI-like methods. Then, the genISI functional can open the path using the ISI-like method in solid-state calculations.

Exciton dispersion and exciton-phonon interaction in solids by time-dependent density functional theory

Understanding, predicting, and ultimately controlling exciton band structure and exciton dynamics are central to diverse chemical and materials problems. Here, we have developed a first-principles method to determine exciton dispersion and exciton-phonon interaction in semiconducting and insulating solids based on time-dependent density functional theory. The first-principles method is formulated in planewave bases and pseudopotentials and can be used to compute exciton band structures, exciton charge density, ionic forces, the non-adiabatic coupling matrix between excitonic states, and the exciton-phonon coupling matrix. Based on the spinor formulation, the method enables self-consistent noncollinear calculations to capture spin-orbital coupling. Hybrid exchange-correlation functionals are incorporated to deal with long-range electron-hole interactions in solids. A sub-Hilbert space approximation is introduced to reduce the computational cost without loss of accuracy. For validations, we have applied the method to compute the exciton band structure and exciton-phonon coupling strength in transition metal dichalcogenide monolayers; both agree very well with the previous GW-Bethe-Salpeter equation and experimental results. This development paves the way for accurate determinations of exciton dynamics in a wide range of solid-state materials.

An efficient adaptive variational quantum solver of the Schrödinger equation based on reduced density matrices

Recently, adaptive variational quantum algorithms, e.g., Adaptive Derivative-Assembled Pseudo-Trotter-Variational Quantum Eigensolver (ADAPT-VQE) and Iterative Qubit-Excitation Based-Variational Quantum Eigensolver (IQEB-VQE), have been proposed to optimize the circuit depth, while a huge number of additional measurements make these algorithms highly inefficient. In this work, we reformulate the ADAPT-VQE with reduced density matrices (RDMs) to avoid additional measurement overhead. With Valdemoro's reconstruction of the three-electron RDM, we present a revised ADAPT-VQE algorithm, termed ADAPT-V, without any additional measurements but at the cost of increasing variational parameters compared to the ADAPT-VQE. Furthermore, we present an ADAPT-Vx algorithm by prescreening the anti-Hermitian operator pool with this RDM-based scheme. ADAPT-Vx requires almost the same variational parameters as ADAPT-VQE but a significantly reduced number of gradient evaluations. Numerical benchmark calculations for small molecules demonstrate that ADAPT-V and ADAPT-Vx provide an accurate description of the ground- and excited-state potential energy curves. In addition, to minimize the quantum resource demand, we generalize this RDM-based scheme to circuit-efficient IQEB-VQE algorithm and achieve significant measurement reduction.

Two-level iterative solver for linear response time-dependent density functional theory with plane wave basis set

We present a two-level iterative solver for linear response time-dependent density functional theory by combining two forms of the Casida equation in the Kohn-Sham orbital representation and in Hutter's formulation. This two-level iterative solver has been implemented with the plane wave pseudopotential method for excited-state simulations of molecular and low-dimensional solid materials. Numerical studies with the Davidson algorithm demonstrate that this two-level iterative solver yields excited-state properties for molecules (benzene C₆H₆ and fullerene C₆₀) and low-dimensional semiconductors [two-dimensional molybdenum disulfide MoS₂ monolayer and rutile titanium dioxide TiO₂(110) surface] with significantly reduced computational cost and storage requirement compared with standard iterative algorithms. We apply our approach to investigate the photoinduced charge separation of methanol molecules adsorption on the rutile TiO₂(110) surface from the exciton perspective and validate that the photogenerated hole can be captured by methanol molecules.

Molecular excitations from meta-generalized gradient approximations in the Kohn-Sham scheme

Meta-Generalized Gradient Approximations (meta-GGAs) can, in principle, include spatial and temporal nonlocality in time-dependent density functional theory at a much lower computational cost than functionals that use exact exchange. We here test whether a meta-GGA that has recently been developed with a focus on capturing nonlocal response properties and the particle number discontinuity can realize such features in practice. To this end, we extended the frequency-dependent Sternheimer formalism to the meta-GGA case. Using the Krieger-Li-Iafrate (KLI) approximation, we calculate the optical response for the selected paradigm molecular systems and compare the meta-GGA Kohn-Sham response to the one found with exact exchange and conventional (semi-)local functionals. We find that the new meta-GGA captures important properties of the nonlocal exchange response. The KLI approximation, however, emerges as a limiting factor in the evaluation of charge-transfer excitations.

Modified Interaction-Strength Interpolation Method as an Important Step toward Self-Consistent Calculations

The modified point charge plus continuum (mPC) model [ConstantinL. A.; Phys. Rev. B2019, 99, 085117] solves the important failures of the original counterpart, namely, the divergences when the reduced gradient of the density is large, such as in the tail of the density and in quasi-dimensional density regimes. The mPC allows us to define a modified interaction-strength interpolation (mISI) method inheriting these good features, which are important steps toward the full self-consistent treatment. Here, we provide an assessment of mISI for molecular systems (i.e., considering thermochemistry properties, correlation energies, vertical ionization potentials, and several noncovalent interactions), harmonium atoms, and functional derivatives in the strong-interaction limit. For all our tests, mISI provides a systematic improvement over the original ISI method. Semilocal approximations of the second-order Görling–Levy (GL2) perturbation theory are also considered in the mISI method, showing considerable worsening of the results. Possible further development of mISI is briefly discussed.

Oxygen K-edge X-ray Absorption Spectra

We review oxygen K-edge X-ray absorption spectra of both molecules and solids. We start with an overview of the main experimental aspects of oxygen K-edge X-ray absorption measurements including X-ray sources, monochromators, and detection schemes. Many recent oxygen K-edge studies combine X-ray absorption with time and spatially resolved measurements and/or operando conditions. The main theoretical and conceptual approximations for the simulation of oxygen K-edges are discussed in the Theory section. We subsequently discuss oxygen atoms and ions, binary molecules, water, and larger molecules containing oxygen, including biomolecular systems. The largest part of the review deals with the experimental results for solid oxides, starting from s- and p-electron oxides. Examples of theoretical simulations for these oxides are introduced in order to show how accurate a DFT description can be in the case of s and p electron overlap. We discuss the general analysis of the 3d transition metal oxides including discussions of the crystal field effect and the effects and trends in oxidation state and covalency. In addition to the general concepts, we give a systematic overview of the oxygen K-edges element by element, for the s-, p-, d-, and f-electron systems.

Strong Quantum Coherence between Fermi Liquid Mahan Excitons

In modulation doped quantum wells, the excitons are formed as a result of the interactions of the charged holes with the electrons at the Fermi edge in the conduction band, leading to the so-called "Mahan excitons." The binding energy of Mahan excitons is expected to be greatly reduced and any quantum coherence destroyed as a result of the screening and electron-electron interactions. Surprisingly, we observe strong quantum coherence between the heavy hole and light hole excitons. Such correlations are revealed by the dominating cross-diagonal peaks in both one-quantum and two-quantum two-dimensional Fourier transform spectra. Theoretical simulations based on the optical Bloch equations where many-body effects are included phenomenologically reproduce well the experimental spectra. Time-dependent density functional theory calculations provide insight into the underlying physics and attribute the observed strong quantum coherence to a significantly reduced screening length and collective excitations of the many-electron system.

Ab initio study on optical properties of glycine sodium nitrate: a novel semiorganic nonlinear optical crystal

The electronic structure, and linear and nonlinear optical susceptibilities of crystalline glycine-sodium nitrate (GSN) has been studied using the full potential linear augmented plane wave method within density-functional theory.

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.

Optics of semiconductors from meta-generalized-gradient-approximation-based time-dependent density-functional theory

We calculate the optical spectra of silicon, germanium, and zinc blende semiconductors in the adiabatic time-dependent density-functional formalism, making use of kinetic energy density-dependent [meta-generalized-gradient-approximation (GGA)] exchange-correlation functionals. We find excellent agreement between theory and experiment. The success of the theory on this notoriously difficult problem is traced to the fact that the exchange-correlation kernel of meta-GGA supports a singularity of the form α/q(2) (where q is the wave vector and α is a constant), whereas previously employed approximations (e.g., local-density and generalized gradient approximations) do not. Thus, the use of the adiabatic meta-GGA opens a new path for handling the extreme nonlocality of the time-dependent exchange-correlation potential in solid-state systems.

Real-time electron dynamics with exact-exchange time-dependent density-functional theory

The exact exchange potential in time-dependent density-functional theory is defined as an orbital functional through the time-dependent optimized effective potential (TDOEP) method. We numerically solve the TDOEP integral equation for the real-time nonlinear intersubband electron dynamics in a semiconductor quantum well with two occupied subbands. It is found that memory effects become significant in the vicinity of intersubband resonances.

Density functionals from many-body perturbation theory: The band gap for semiconductors and insulators

Theoretically the Kohn-Sham band gap differs from the exact quasiparticle energy gap by the derivative discontinuity of the exchange-correlation functional. In practice for semiconductors and insulators the band gap calculated within any local or semilocal density approximations underestimates severely the experimental energy gap. On the other hand, calculations with an "exact" exchange potential derived from many-body perturbation theory via the optimized effective potential suggest that improving the exchange-correlation potential approximation can yield a reasonable agreement between the Kohn-Sham band gap and the experimental gap. The results in this work show that this is not the case. In fact, we add to the exact exchange the correlation that corresponds to the dynamical (random phase approximation) screening in the GW approximation. This accurate exchange-correlation potential provides band structures similar to the local density approximation with the corresponding derivative discontinuity that contributes 30%-50% to the energy gap. Our self-consistent results confirm substantially the results for Si and other semiconductors obtained perturbatively [R. W. Godby et al., Phys. Rev. B 36, 6497 (1987)] and extend the conclusion to LiF and Ar, a wide-gap insulator and a noble-gas solid.

Beyond time-dependent exact exchange: The need for long-range correlation

In the description of the interaction between electrons beyond the classical Hartree picture, bare exchange often yields a leading contribution. Here we discuss its effect on optical spectra of solids, comparing three different frameworks: time-dependent Hartree-Fock, a recently introduced combined density-functional and Green's function approaches applied to the bare exchange self-energy, and time-dependent exact exchange within time-dependent density-functional theory (TD-EXX). We show that these three approximations give rise to identical excitonic effects in solids; these effects are drastically overestimated for semiconductors. They are partially compensated by the usual overestimation of the quasiparticle band gap within Hartree-Fock. The physics that lacks in these approaches can be formulated as screening. We show that the introduction of screening in TD-EXX indeed leads to a formulation that is equivalent to previously proposed functionals derived from many-body perturbation theory. It can be simulated by reducing the long-range part of the Coulomb interaction: this produces absorption spectra of semiconductors in good agreement with experiment.

Away from generalized gradient approximation: Orbital-dependent exchange-correlation functionals

The local-density approximation of density functional theory (DFT) is remarkably accurate, for instance, for geometries and frequencies, and the generalized gradient approximations have also made bond energies quite reliable. Sometimes, however, one meets with failure in individual cases. One of the possible routes towards better functionals would be the incorporation of orbital dependence (which is an implicit density dependency) in the functionals. We discuss this approach both for energies and for response properties. One possibility is the use of the Hartree-Fock-type exchange energy expression as orbital-dependent functional. We will argue that in spite of the increasing popularity of this approach, it does not offer any advantage over Hartree-Fock for energies. We will advocate not to apply the separation of exchange and correlation, which is so ingrained in quantum chemistry, but to model both simultaneously. For response properties the energies and shapes of the virtual orbitals are crucial. We will discuss the benefits that Kohn-Sham potentials can offer which are derived from either an orbital-dependent energy functional, including the exact-exchange functional, or which can be obtained directly as orbital-dependent functional. We highlight the similarity of the Hartree-Fock and Kohn-Sham occupied orbitals and orbital energies, and the essentially different meanings the virtual orbitals and orbital energies have in these two models. We will show that these differences are beneficial for DFT in the case of localized excitations (in a small molecule or in a fragment), but are detrimental for charge-transfer excitations. Again, orbital dependency, in this case in the exchange-correlation kernel, offers a solution.

Orbital- and state-dependent functionals in density-functional theory

Shortcomings of present density-functional methods are considered. Kohn-Sham and time-dependent density-functional methods using orbital- and state-dependent functionals for exchange-correlation energies, potentials, and kernels are discussed as possible remedy for some of these shortcomings. A view on the Kohn-Sham formalism is presented which differs somewhat from the one conventionally taken. The crucial step of constructing local multiplicative exchange-correlation potentials in Kohn-Sham methods based on orbital- and state-dependent functionals is discussed. The description of open-shell systems via a symmetrized Kohn-Sham formalism employing state-dependent exchange-correlation functionals is elucidated. The generalized adiabatic connection Kohn-Sham approach for the self-consistent treatment of excited states within a density-functional framework is considered. In the latter approach orbital- and state-dependent exchange-correlation functionals occur in a density-functional framework which is no longer based on the Hohenberg-Kohn theorem but on a more general relation between electron densities and local multiplicative potentials.

TIME-DEPENDENT DENSITY FUNCTIONAL THEORY

Time-dependent density functional theory (TDDFT) can be viewed as an exact reformulation of time-dependent quantum mechanics, where the fundamental variable is no longer the many-body wave function but the density. This time-dependent density is determined by solving an auxiliary set of noninteracting Schrodinger equations, the Kohn-Sham equations. The nontrivial part of the many-body interaction is contained in the so-called exchange-correlation potential, for which reasonably good approximations exist. Within TDDFT two regimes can be distinguished: (a) If the external time-dependent potential is "small," the complete numerical solution of the time-dependent Kohn-Sham equations can be avoided by the use of linear response theory. This is the case, e.g., for the calculation of photoabsorption spectra. (b) For a "strong" external potential, a full solution of the time-dependent Kohn-Sham equations is in order. This situation is encountered, for instance, when matter interacts with intense laser fields. In this review we give an overview of TDDFT from its theoretical foundations to several applications both in the linear and in the nonlinear regime.

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

A robust and efficient frequency dependent and nonlocal exchange correlation f(xc)(r,r(');omega) is derived by imposing time-dependent density-functional theory (TDDFT) to reproduce the many-body diagrammatic expansion of the Bethe-Salpeter polarization function. As an illustration, we compute the optical spectra of LiF, SiO2, and diamond and the finite momentum transfer energy-loss spectrum of LiF. The TDDFT results reproduce extremely well the excitonic effects embodied in the Bethe-Salpeter approach, both for strongly bound and resonant excitons. We provide a working expression for f(xc) that is fast to evaluate and easy to implement.

Parameter-free calculation of response functions in time-dependent density-functional theory

We have established and implemented a fully ab initio method which allows one to calculate optical absorption spectra, including excitonic effects, without solving the cumbersome Bethe-Salpeter equation, but obtaining results of the same precision. This breakthrough has been achieved in the framework of time-dependent density-functional theory, using new exchange-correlation kernels f(xc) that are free of any empirical parameter. We show that the same excitonic effects in the optical spectra can be reproduced through different f(xc)'s, ranging from frequency-dependent ones to a static one, by varying the kernel's spatial degrees of freedom. This indicates that the key quantity is not f(xc), but f(xc) combined with a response function. We present results for the optical absorption of bulk Si and SiC in good agreement with experiment, almost indistinguishable from those of the Bethe-Salpeter approach.