Static Self-Energy and Effective Mass of the Homogeneous Electron Gas from Quantum Monte Carlo Calculations.

We discuss the methodology of quantum Monte Carlo calculations of the effective mass based on the static self-energy Σ(k,0). We then use variational Monte Carlo calculations of Σ(k,0) of the homogeneous electron gas at various densities to obtain results very close to perturbative G_{0}W_{0} calculations for values of the density parameter 1≤r_{s}≤10. The obtained values for the effective mass are close to diagrammatic Monte Carlo results and disagree with previous quantum Monte Carlo calculations based on a heuristic mapping of excitation energies to those of an ideal gas.

XPS Core-Level Chemical Shift by Ab Initio Many-Body Theory.

X-ray photoemission spectroscopy (XPS) provides direct information on atomic composition and stoichiometry by measuring core-electron binding energies. Moreover, from the shift of the binding energy, the so-called chemical shift, the precise chemical type of bonds can be inferred, which brings additional information on the local structure. In this work, we present a theoretical study of the chemical shift first by comparing different theories, from Hartree-Fock and density functional theory to many-body perturbation theory approaches like the GW approximation and its static version (COHSEX). The accuracy of each theory is assessed with respect to a carbon 1s chemical shift experimental benchmark measured on a set of gas-phase molecules. More importantly, by decomposing the chemical shift into different contributions according to terms in the total Hamiltonian, classical electrostatics is identified as the major contributor to the chemical shift, one order of magnitude larger than the correlation.

A formally exact one-frequency-only Bethe-Salpeter-like equation. Similarities and differences between GW+BSE and self-consistent RPA.

A formally exact Bethe-Salpeter-like equation for the linear-response function is introduced with a kernel which depends only on the one frequency of the applied field. This is in contrast with the standard Bethe-Salpeter equation (BSE) which involves multiple-frequency integrals over the kernel and response functions. From the one-frequency kernel, known approximations are straightforwardly recovered. However, the present formalism lends itself to more powerful approximations. This is demonstrated with the exact analytical solution of the Hubbard molecule. Similarities and differences of the GW + BSE approach with the self-consistent random-phase approximation are also discussed.

Helium Atom Excitations by the GW and Bethe-Salpeter Many-Body Formalism.

The helium atom is the simplest many-body electronic system provided by nature. The exact solution to the Schrödinger equation is known for helium ground and excited states, and it represents a benchmark for any many-body methodology. Here, we check the ab initio many-body GW approximation and the Bethe-Salpeter equation (BSE) against the exact solution for helium. Starting from the Hartree-Fock method, we show that the GW and the BSE yield impressively accurate results on excitation energies and oscillator strength, systematically improving the time-dependent Hartree-Fock method. These findings suggest that the accuracy of the BSE and GW approximations is not significantly limited by self-interaction and self-screening problems even in this few electron limit. We further discuss our results in comparison to those obtained by time-dependent density-functional theory.

Bethe-Salpeter equation insights into the photo-absorption function and exciton structure of chlorophyll a and b in light-harvesting complex II.

The photo-absorption process and the excitation of chlorophyll (Chl) is the primary and essential step of photosynthesis in green plants. By solving the Bethe-Salpeter equation (BSE) on top of the GW approximation within ab initio many-body perturbation theory, we calculate the photo-absorption function and the excitons structure of Chl a and b in their in vivo conformations as measured by X-ray diffraction in the light-harvesting complex (LHC) II. BSE optical absorption spectra are in good agreement with the experiment and we discuss residual discrepancies. The experimental evidence of multiple Chla forms in vivo is explained by BSE. The Chla and Chlb BSE exciton wavefunctions present important charge-transfer differences on the Soret band. Q excitons are almost identical, apart from charge (both electron and hole) localization on the Chlb C7 aldheide formyl group, absent on the Chla methyl C7, that is exactly the group where the two chlorophylls differ.

Momentum distribution of the homogeneous electron gas.

We calculate the off-diagonal density matrix of the homogeneous electron gas at zero temperature using unbiased reptation Monte Carlo calculations for various densities and extrapolate the momentum distribution and the kinetic and potential energies to the thermodynamic limit. Our results on the renormalization factor allow us to validate approximate G0W0 calculations concerning quasiparticle properties over a broad density region (1≤r(s)≲10) and show that, near the Fermi surface, vertex corrections and self-consistency aspects almost cancel each other out.

Momentum distribution and renormalization factor in sodium and the electron gas.

We present experimental and theoretical results on the momentum distribution and the quasiparticle renormalization factor in sodium. From an x-ray Compton-profile measurement of the valence-electron-momentum density, we derive its discontinuity at the Fermi wave vector. This yields an accurate measure of the renormalization factor that we compare with quantum Monte Carlo and G0W0 calculations performed both on crystalline sodium and on the homogeneous electron gas.

Giant excitonic absorption and emission in two-dimensional group-III nitrides.

Absorption and emission of pristine-like semiconducting monolayers of BN, AlN, GaN, and InN are systematically studied by ab-initio methods. We calculate the absorption spectra for in-plane and out-of-plane light polarization including quasiparticle and excitonic effects. Chemical trends with the cation of the absorption edge and the exciton binding are discussed in terms of the band structures. Exciton binding energies and localization radii are explained within the Rytova-Keldysh model for excitons in two dimensions. The strong excitonic effects are due to the interplay of low dimensionality, confinement effects, and reduced screening. We find exciton radiative lifetimes ranging from tenths of picoseconds (BN) to tenths of nanoseconds (InN) at room temperature, thus making 2D nitrides, especially InN, promising materials for light-emitting diodes and high-performance solar cells.

Transforming nonlocality into a frequency dependence: a shortcut to spectroscopy.

Measurable spectra are often derived from contractions of many-body Green's functions. One calculates hence more information than needed. Here we present and illustrate an in principle exact approach to construct effective potentials and kernels for the direct calculation of electronic spectra. In particular, a dynamical but local and real potential yields the spectral function needed to describe photoemission. We discuss for model solids the frequency dependence of this "photoemission potential" stemming from the nonlocality of the corresponding self-energy.

Can molecular projected density of states (PDOS) be systematically used in electronic conductance analysis?

Using benzenediamine and benzenedithiol molecular junctions as benchmarks, we investigate the widespread analysis of the quantum transport conductance in terms of the projected density of states (PDOS) onto molecular orbitals (MOs). We first consider two different methods for identifying the relevant MOs: (1) diagonalization of the Hamiltonian of the isolated molecule and (2) diagonalization of a submatrix of the junction Hamiltonian constructed by considering only basis elements localized on the molecule. We find that these two methods can lead to substantially different MOs and hence PDOS. Furthermore, within Method 1, the PDOS can differ depending on the isolated molecule chosen to represent the molecular junction (e.g., with or without dangling bonds); within Method 2, the PDOS depends on the chosen basis set. We show that these differences can be critical when the PDOS is used to provide a physical interpretation of the conductance (especially when its value is small, as it happens typically at zero bias). In this work, we propose a new approach in an attempt to reconcile the two traditional methods. Although some improvements were achieved, the main problems remain unsolved. Our results raise more general questions and doubts on a PDOS-based analysis of the conductance.

Coherent electronic transport through graphene constrictions: subwavelength regime and optical analogy.

Nanoelectronic devices smaller than the electron wavelength can be achieved in graphene with current lithography techniques. Here we show that the electronic quantum transport of graphene subwavelength nanodevices presents deep analogies with subwavelength optics. We introduce the concept of electronic diffraction barrier to represent the effect of constrictions and the rich transport phenomena of a variety of nanodevices. Results are presented for Bethe and Kirchhoff diffraction in graphene slits and Fabry-Perot interference oscillations in nanoribbons. The same concept applies to graphene quantum dots and gives new insight into recent experiments in these systems.

Ab initio GW many-body effects in graphene.

We present an ab initio numerical many-body GW calculation of the band plot in freestanding graphene. We consider the full ionic and electronic structure introducing e-e interaction and correlation effects via a self-energy containing non-Hermitian and dynamical terms. With respect to the density-functional theory local-density approximation, the Fermi velocity is renormalized with an increase of 17%, in better agreement with the experiment. Close to the Dirac point the linear dispersion is modified by the presence of a kink, as observed by angle-resolved photoemission spectroscopy. We demonstrate that the kink is due to low-energy pi-->pi* single-particle excitations and to the pi plasmon. The GW self-energy does not open the band gap.

Understanding correlations in vanadium dioxide from first principles.

Vanadium dioxide is a prototype material for the discussion of correlation effects in solids. First-principles density-functional theory does not describe the metal-insulator transition, whereas strongly correlated models reproduce the main features. Here we present a parameter-free GW calculation of VO2 and show that the correlation effects in the band structure of both the metallic and the insulating phases are correctly reproduced, provided that quasiparticle energies and wave functions are calculated self-consistently. Our calculations explain the satellite in the photoemission spectrum of the metal as due to a plasmon resonance in the energy-loss function and show that this feature disappears in the insulator.

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.

Signatures of short-range many-body effects in the dielectric function of silicon for finite momentum transfer.

We present an investigation of the dynamic structure factor and of the dielectric function epsilonM(Q,omega) of the prototypical semiconductor silicon for finite momentum transfer, combining inelastic x-ray scattering experiments and ab initio calculations. In contrast with optical spectra, for finite momentum transfer time-dependent density-functional theory in the adiabatic local-density approximation together with lifetime broadening describes the physics of valence excitations correctly. Major structures in the spectra, governed by short-range crystal and exchange-correlation local-field effects, are strongly influenced by a mixing of transitions of positive and negative energies, in striking difference to spectra for vanishing momentum transfer. This mixing gives rise to a pronounced Fano asymmetry.

Many-body perturbation theory using the density-functional concept: beyond the GW approximation.

We propose an alternative formulation of many-body perturbation theory that uses the density-functional concept. Instead of the usual four-point integral equation for the polarizability, we obtain a two-point one, which leads to excellent optical absorption and energy-loss spectra. The corresponding three-point vertex function and self-energy are then simply calculated via an integration, for any level of approximation. Moreover, we show the direct impact of this formulation on the time-dependent density-functional theory. Numerical results for the band gap of bulk silicon and solid argon illustrate corrections beyond the GW approximation for the self-energy.

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.

Excitonic effects in solids described by time-dependent density-functional theory.

Starting from the many-body Bethe-Salpeter equation we derive an exchange-correlation kernel f(xc) that reproduces excitonic effects in bulk materials within time-dependent density functional theory. The resulting f(xc) accounts for both self-energy corrections and the electron-hole interaction. It is static, nonlocal, and has a long-range Coulomb tail. Taking the example of bulk silicon, we show that the -alpha/q(2) divergency is crucial and can, in the case of continuum excitons, even be sufficient for reproducing the excitonic effects and yielding excellent agreement between the calculated and the experimental absorption spectrum.

Local field effects in the electron energy loss spectra of rutile TiO2.

We present an ab initio calculation of the electron energy loss spectrum of rutile TiO2 in the energy range of 0 to 60 eV, focusing our interest on the excitation from the titanium 3p semicore levels. The results are compared to our measurements. Local field effects turn out to be crucial at those energies, and their inclusion in the calculation yields excellent agreement between theory and experiment. We show how in rutile these effects induce an anisotropy in the otherwise isotropic transitions from quasispherical 3p semicore states to 3d states of almost cubic symmetry.

Ab Initio calculations of the anisotropic dielectric tensor of GaAs/AlAs superlattices.

The static dielectric properties of (001)(GaAs)(p)/(AlAs)(p) superlattices have been calculated as a function of their period p for 1< or = p < or =12, starting from density-functional theory. The interplay between quantum confinement and local field effects is shown to be crucial. For light polarized in the growth direction it leads to the otherwise surprising justification of the use of a classical effective medium theory, even for the smallest periods. Only the inclusion of both contributions allows in ab initio and in semiempirical calculations to reproduce the experimentally observed birefringence.