Scrutinizing GW-Based Methods Using the Hubbard Dimer.

Using the simple (symmetric) Hubbard dimer, we analyze some important features of the GW approximation. We show that the problem of the existence of multiple quasiparticle solutions in the (perturbative) one-shot GW method and its partially self-consistent version is solved by full self-consistency. We also analyze the neutral excitation spectrum using the Bethe-Salpeter equation (BSE) formalism within the standard GW approximation and find, in particular, that 1) some neutral excitation energies become complex when the electron-electron interaction U increases, which can be traced back to the approximate nature of the GW quasiparticle energies; 2) the BSE formalism yields accurate correlation energies over a wide range of U when the trace (or plasmon) formula is employed; 3) the trace formula is sensitive to the occurrence of complex excitation energies (especially singlet), while the expression obtained from the adiabatic-connection fluctuation-dissipation theorem (ACFDT) is more stable (yet less accurate); 4) the trace formula has the correct behavior for weak (i.e., small U) interaction, unlike the ACFDT expression.

Photoemission Spectra from the Extended Koopman's Theorem, Revisited.

The Extended Koopman's Theorem (EKT) provides a straightforward way to compute charged excitations from any level of theory. In this work we make the link with the many-body effective energy theory (MEET) that we derived to calculate the spectral function, which is directly related to photoemission spectra. In particular, we show that at its lowest level of approximation the MEET removal and addition energies correspond to the so-called diagonal approximation of the EKT. Thanks to this link, the EKT and the MEET can benefit from mutual insight. In particular, one can readily extend the EKT to calculate the full spectral function, and choose a more optimal basis set for the MEET by solving the EKT secular equation. We illustrate these findings with the examples of the Hubbard dimer and bulk silicon.

Optical spectra of 2D monolayers from time-dependent density functional theory.

The optical spectra of two-dimensional (2D) periodic systems provide a challenge for time-dependent density-functional theory (TDDFT) because of the large excitonic effects in these materials. In this work we explore how accurately these spectra can be described within a pure Kohn-Sham time-dependent density-functional framework, i.e., a framework in which no theory beyond Kohn-Sham density-functional theory, such as GW, is required to correct the Kohn-Sham gap. To achieve this goal we adapted a recent approach we developed for the optical spectra of 3D systems [S. Cavo, J. A. Berger and P. Romaniello, Phys. Rev. B, 2020, 101, 115109] to those of 2D systems. Our approach relies on the link between the exchange-correlation kernel of TDDFT and the derivative discontinuity of ground-state density-functional theory, which guarantees a correct quasi-particle gap, and on a generalization of the polarization functional [J. A. Berger, Phys. Rev. Lett., 2015, 115, 137402], which describes the excitonic effects. We applied our approach to two prototypical 2D monolayers, h-BN and MoS2. We find that our protocol gives a qualitatively good description of the optical spectrum of h-BN, whereas improvements are needed for MoS2 to describe the intensity of the excitonic peaks.

Many-Body Effective Energy Theory: Photoemission at Strong Correlation.

In this work, we explore the performance of a recently derived many-body effective energy theory for the calculation of photoemission spectra in the regime of strong electron correlation. We apply the theory to paramagnetic MnO, FeO, CoO, and NiO, which are typical examples of strongly correlated materials and, therefore, a challenge for standard theories. We show that our method opens a correlation gap in all of the oxides studied without breaking the symmetry. Although the materials seem similar, we show that an analysis of the occupation numbers reveals that the nature of the gap is not the same for these materials. Overall, the results are very promising, although improvements are clearly required, since the band gap is overestimated for all of the systems studied. We indicate some possible strategies to further develop the theory.

Reduced density-matrix functional theory: Correlation and spectroscopy.

In this work, we explore the performance of approximations to electron correlation in reduced density-matrix functional theory (RDMFT) and of approximations to the observables calculated within this theory. Our analysis focuses on the calculation of total energies, occupation numbers, removal/addition energies, and spectral functions. We use the exactly solvable Hubbard dimer at 1/4 and 1/2 fillings as test systems. This allows us to analyze the underlying physics and to elucidate the origin of the observed trends. For comparison, we also report the results of the GW approximation, where the self-energy functional is approximated, but no further hypothesis is made concerning the approximations of the observables. In particular, we focus on the atomic limit, where the two sites of the dimer are pulled apart and electrons localize on either site with equal probability, unless a small perturbation is present: this is the regime of strong electron correlation. In this limit, using the Hubbard dimer at 1/2 filling with or without a spin-symmetry-broken ground state allows us to explore how degeneracies and spin-symmetry breaking are treated in RDMFT. We find that, within the used approximations, neither in RDMFT nor in GW, the signature of strong correlation is present, when looking at the removal/addition energies and spectral function from the spin-singlet ground state, whereas both give the exact result for the spin-symmetry broken case. Moreover, we show how the spectroscopic properties change from one spin structure to the other.

Towards time-dependent current-density-functional theory in the non-linear regime.

Time-Dependent Density-Functional Theory (TDDFT) is a well-established theoretical approach to describe and understand irradiation processes in clusters and molecules. However, within the so-called adiabatic local density approximation (ALDA) to the exchange-correlation (xc) potential, TDDFT can show insufficiencies, particularly in violently dynamical processes. This is because within ALDA the xc potential is instantaneous and is a local functional of the density, which means that this approximation neglects memory effects and long-range effects. A way to go beyond ALDA is to use Time-Dependent Current-Density-Functional Theory (TDCDFT), in which the basic quantity is the current density rather than the density as in TDDFT. This has been shown to offer an adequate account of dissipation in the linear domain when the Vignale-Kohn (VK) functional is used. Here, we go beyond the linear regime and we explore this formulation in the time domain. In this case, the equations become very involved making the computation out of reach; we hence propose an approximation to the VK functional which allows us to calculate the dynamics in real time and at the same time to keep most of the physics described by the VK functional. We apply this formulation to the calculation of the time-dependent dipole moment of Ca, Mg and Na2. Our results show trends similar to what was previously observed in model systems or within linear response. In the non-linear domain, our results show that relaxation times do not decrease with increasing deposited excitation energy, which sets some limitations to the practical use of TDCDFT in such a domain of excitations.

On transition rates in surface hopping.

Trajectory surface hopping (TSH) is one of the most widely used quantum-classical algorithms for nonadiabatic molecular dynamics. Despite its empirical effectiveness and popularity, a rigorous derivation of TSH as the classical limit of a combined quantum electron-nuclear dynamics is still missing. In this work, we aim to elucidate the theoretical basis for the widely used hopping rules. Naturally, we concentrate thereby on the formal aspects of the TSH. Using a Gaussian wave packet limit, we derive the transition rates governing the hopping process at a simple avoided level crossing. In this derivation, which gives insight into the physics underlying the hopping process, some essential features of the standard TSH algorithm are retrieved, namely (i) non-zero electronic transition rate ("hopping probability") at avoided crossings; (ii) rescaling of the nuclear velocities to conserve total energy; (iii) electronic transition rates linear in the nonadiabatic coupling vectors. The well-known Landau-Zener model is then used for illustration.

Nonlinear optical properties of Ni(Me(6)pzS(2))MX (M = Ni, Pd, Pt; X = Me(2)timdt, mnt).

In this work we investigate the second-order response of complexes of formula Ni(Me(6)pzS(2))MX (M = Ni, Pd, Pt; X = Me(2)timdt (monoanion of N,N'-disubstituted imidazolidine-2,4,5-trithione), mnt (maleonitriledithiolate)): by binding the porphyrazine to the metal-dithiolene, the electron asymmetry and pi-conjugation of the latter is increased, and its second-order response can result enhanced. By performing ab initio calculations of the ground-state and response properties of these compounds, we predict the molecular first hyperpolarizability, we elucidate its electronic origin, and we illustrate its dependence on the metal and the dithiolene ligand. Our study indicates that these complexes show a very high second-order response, comparable to that of organic "push-pull" materials, and that the appropriate metal-dithiolene combination can significantly enhance it.

Double excitations in finite systems.

Time-dependent density-functional theory (TDDFT) is widely used in the study of linear response properties of finite systems. However, there are difficulties in properly describing excited states, which have double- and higher-excitation characters, which are particularly important in molecules with an open-shell ground state. These states would be described if the exact TDDFT kernel were used; however, within the adiabatic approximation to the exchange-correlation (xc) kernel, the calculated excitation energies have a strict single-excitation character and are fewer than the real ones. A frequency-dependent xc kernel could create extra poles in the response function, which would describe states with a multiple-excitation character. We introduce a frequency-dependent xc kernel, which can reproduce, within TDDFT, double excitations in finite systems. In order to achieve this, we use the Bethe-Salpeter equation with a dynamically screened Coulomb interaction W(omega), which can describe these excitations, and from this we obtain the xc kernel. Using a two-electron model system, we show that the frequency dependence of W does indeed introduce the double excitations that are instead absent in any static approximation of the electron-hole screening.

The self-energy beyond GW: local and nonlocal vertex corrections.

It is commonly accepted that the GW approximation for the electron self-energy is successful for the description of the band structure of weakly to moderately correlated systems, whereas it will fail for strongly correlated materials. In the present work, we discuss two important aspects of this approximation: first, the "self-screening error," which is due to an incorrect treatment of induced exchange, and second, the atomic limit, in which, instead, correlation is directly responsible for the observed problem. Using the example of the removal of a particle from a box, we show that the self-screening error stems from the use of test charge-test charge screening and that it can be corrected by a two-point vertex contribution to the self-energy derived from time-dependent density functional theory (TDDFT). We explain why the addition of a particle, instead, requires the use of a different approximate vertex. This illustrates why the general vertex function, valid both for valence and conduction states, must be a three-point function. Moreover, we show that also the bad performance of GW in the atomic limit is due to the neglect of the vertex in the self-energy; in that case, the TDDFT-derived vertex correction is not sufficient in order to remove the error even qualitatively. We discuss the effects of the self-screening error as well as the atomic limit using GW for the exactly solvable two-site Hubbard model.

Relativistic two-component formulation of time-dependent current-density functional theory: application to the linear response of solids.

In this paper we derive the relativistic two-component formulation of time-dependent current-density-functional theory. To arrive at a two-component current-density formulation we apply a Foldy-Wouthuysen-type transformation to the time-dependent four-component Dirac-Kohn-Sham equations of relativistic density-functional theory. The two-component Hamiltonian is obtained as a regular expansion which is gauge invariant at each order of approximation, and to zeroth order it represents the time-dependent version of the relativistic zeroth order regular Hamiltonian obtained by van Lenthe et al., for the ground state [J. Chem. Phys.99, 4597 (1993)]. The corresponding zeroth order regular expression for the density is unchanged, whereas the current-density operator now comprises a paramagnetic, a diamagnetic, and a spin contribution, similar to the Gordon decomposition of the Dirac four current. The zeroth order current density is directly related to the mean velocity corresponding to the zeroth order Hamiltonian. These density and current density operators satisfy the continuity equation. This zeroth order approximation is therefore consistent and physically realistic. By combining this formalism with the formulation of the linear response of solids within time-dependent current-density functional theory [Romaniello and de Boeij, Phys. Rev. B71, 155108 (2005)], we derive a method that can treat orbital and spin contributions to the response in a unified way. The effect of spin-orbit coupling can now be taken into account. As first test we apply the method to calculate the relativistic effects in the linear response of several metals and nonmetals to a macroscopic electric field. Treatment of spin-orbit coupling yields visible changes in the spectra: a smooth onset of the interband transitions in the absorption spectrum of Au, a sharp onset with peak at about 0.46 eV in the absorption spectrum of W, and a low-frequency doublet structure in the absorption spectra of ZnTe, CdTe, and HgTe in agreement with experimental results.

The role of relativity in the optical response of gold within the time-dependent current-density-functional theory.

We included relativistic effects in the formulation of the time-dependent current-density-functional theory for the calculation of linear response properties of metals [P. Romaniello and P. L. de Boeij, Phys. Rev. B (to be published)]. We treat the dominant scalar-relativistic effects using the zeroth-order regular approximation in the ground-state density-functional theory calculations, as well as in the time-dependent response calculations. The results for the dielectric function of gold calculated in the spectral range of 0-10 eV are compared with experimental data reported in literature and recent ellipsometric measurements. As well known, relativistic effects strongly influence the color of gold. We find that the onset of interband transitions is shifted from around 3.5 eV, obtained in a nonrelativistic calculation, to around 1.9 eV when relativity is included. With the inclusion of the scalar-relativistic effects there is an overall improvement of both real and imaginary parts of the dielectric function over the nonrelativistic ones. Nevertheless some important features in the absorption spectrum are not well reproduced, but can be explained in terms of spin-orbit coupling effects. The remaining deviations are attributed to the underestimation of the interband gap (5d-6sp band gap) in the local-density approximation and to the use of the adiabatic local-density approximation in the response calculation.