Accurate localized resolution of identity approach for linear-scaling hybrid density functionals and for many-body perturbation theory

© 2015 IOP Publishing Ltd and Deutsche Physikalische Gesellschaft.A key component in calculations of exchange and correlation energies is the Coulomb operator, which requires the evaluation of two-electron integrals. For localized basis sets, these four-center integrals are most efficiently evaluated with the resolution of identity (RI) technique, which expands basis-function products in an auxiliary basis. In this work we show the practical applicability of a localized RI-variant ('RI-LVL'), which expands products of basis functions only in the subset of those auxiliary basis functions which are located at the same atoms as the basis functions. We demonstrate the accuracy of RI-LVL for Hartree-Fock calculations, for the PBE0 hybrid density functional, as well as for RPA and MP2 perturbation theory. Molecular test sets used include the S22 set of weakly interacting molecules, the G3 test set, as well as the G2-1 and BH76 test sets, and heavy elements including titanium dioxide, copper and gold clusters. Our RI-LVL implementation paves the way for linear-scaling RI-based hybrid functional calculations for large systems and for all-electron many-body perturbation theory with significantly reduced computational and memory cost.

Many-body entanglement in decoherence processes

A pure state decoheres into a mixed state as it entangles with an environment. When an entangled two-mode system is embedded in a thermal environment, however, each mode may not be entangled with its environment by their simple linear interaction. We consider an exactly solvable model to study the dynamics of a total system, which is composed of an entangled two-mode system and a thermal environment. The Markovian interaction with the environment is concerned with an array of infinite number of beam splitters. It is shown that many-body entanglement of the system and the environment may play a crucial role in the process of disentangling the system.

Many-body theory of positron-atom interactions

A many-body theory approach is developed for the problem of positron-atom scattering and annihilation. Strong electron- positron correlations are included nonperturbatively through the calculation of the electron-positron vertex function. It corresponds to the sum of an infinite series of ladder diagrams, and describes the physical effect of virtual positronium formation. The vertex function is used to calculate the positron-atom correlation potential and nonlocal corrections to the electron-positron annihilation vertex. Numerically, we make use of B-spline basis sets, which ensures rapid convergence of the sums over intermediate states. We have also devised an extrapolation procedure that allows one to achieve convergence with respect to the number of intermediate- state orbital angular momenta included in the calculations. As a test, the present formalism is applied to positron scattering and annihilation on hydrogen, where it is exact. Our results agree with those of accurate variational calculations. We also examine in detail the properties of the large correlation corrections to the annihilation vertex.

Independent particle descriptions of tunneling using the many-body quantum transport approach

Currents across thin insulators are commonly taken as single electrons moving across classically forbidden regions; this independent particle picture is well-known to describe most tunneling phenomena. Examining quantum transport from a different perspective, i.e., by explicit treatment of electron-electron interactions, we evaluate different single particle approximations with specific application to tunneling in metal-molecule-metal junctions. We find maximizing the overlap of a Slater determinant composed of single-particle states to the many-body current-carrying state is more important than energy minimization for defining single-particle approximations in a system with open boundary conditions. Thus the most suitable single particle effective potential is not one commonly in use by electronic structure methods, such as the Hartree-Fock or Kohn-Sham approximations.

Control-limited perfect state transfer, quantum stochastic resonance, and many-body entangling gate in imperfect qubit registers

We propose a protocol for perfect quantum state transfer that is resilient to a broad class of realistic experimental imperfections, including noise sources that could be modeled either as independent Markovian baths or as certain forms of spatially correlated environments. We highlight interesting connections between the fidelity of state transfer and quantum stochastic resonance effects. The scheme is flexible enough to act as an effective entangling gate for the generation of genuine multipartite entanglement in a control-limited setting. Possible experimental implementations using superconducting qubits are also briefly discussed.

Many-body theory of gamma spectra from positron-atom annihilation

A many-body theory approach to the calculation of gamma spectra of positron annihilation on many-electron atoms is developed. We evaluate the first-order correlation correction to the annihilation vertex and perform numerical calculations for the noble gas atoms. Extrapolation with respect to the maximal orbital momentum of the intermediate electron and positron states is used to achieve convergence. The inclusion of correlation corrections improves agreement with experimental gamma spectra.

Many-body theory calculations of positron binding to negative ions

A many-body theory approach developed by the authors [Phys. Rev. A 70, 032720 (2004)] is applied to positron bound states and annihilation rates in atomic systems. Within the formalism, full account of virtual positronium (Ps) formation is made by summing the electron-positron ladder diagram series, thus enabling the theory to include all important many-body correlation effects in the positron problem. Numerical calculations have been performed for positron bound states with the hydrogen and halogen negative ions, also known as Ps hydride and Ps halides. The Ps binding energies of 1.118, 2.718, 2.245, 1.873 and 1.393 eV and annihilation rates of 2.544, 2.482, 1.984, 1.913 and 1.809 ns^{-1}, have been obtained for PsH, PsF, PsCl, PsBr and PsI, respectively.

Electronic properties of interfaces and defects from many-body perturbation theory: Recent developments and applications

We review some recent developments in many body perturbation theory (MBPT) calculations that have enabled the study of interfaces and defects. Starting from the theoretical basis of MBPT, Hedin's equations are presented, leading to the CW and CWI' approximations. We introduce the perturbative approach, that is the one most commonly used for obtaining quasiparticle (QP) energies. The practical strategy presented for dealing with the frequency dependence of the self energy operator is based on either plasmon-pole models (PPM) or the contour deformation technique, with the latter being more accurate. We also discuss the extrapolar method for reducing the number of unoccupied states which need to be included explicity in the calculations. The use of the PAW method in the framework of MBPT is also described. Finally, results which have been obtained using, MBPT for band offsets a interfaces and for defects presented, with companies on the main difficulties and cancels.

Schematic representation of the QP corrections (marked with ) to the band edges (E and E-v) and a defect level (F) for a Si/SiO2 interface (Si and O atoms are represented in blue and red, respectively, in the ball and stick model) with an oxygen vacancy leading to a Si-Si bond (the Si atoms involved in this bond are colored light blue).

Electronic properties of zircon and hafnon from many-body perturbation theory

The electronic properties of zircon and hafnon, two wide-gap high-kappa materials, are investigated using many-body perturbation theory (MBPT) combined with the Wannier interpolation technique. For both materials, the calculated band structures differ from those obtained within density-functional theory and MBPT by (i) a slight displacement of the highest valence-band maximum from the Gamma point and (ii) an opening of the indirect band gap to 7.6 and 8.0 eV for zircon and hafnon, respectively. The introduction of vertex corrections in the many-body self-energy does not modify the results except for a global rigid shift of the many-body corrections.

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 , Phys. Rev. B 36, 6497 (1987)] and extend the conclusion to LiF and Ar, a wide-gap insulator and a noble-gas solid. (c) 2006 American Institute of Physics.

Emergent thermodynamics in a quenched quantum many-body system

We study the statistics of the work done, the fluctuation relations and the irreversible entropy production in a quantum many-body system subject to the sudden quench of a control parameter. By treating the quench as a thermodynamic transformation we show that the emergence of irreversibility in the nonequilibrium dynamics of closed many-body quantum systems can be accurately characterized. We demonstrate our ideas by considering a transverse quantum Ising model that is taken out of equilibrium by the instantaneous switching of the transverse field.

Entanglement Replication in Driven Dissipative Many-Body systems

We study the dissipative dynamics of two independent arrays of many-body systems, locally driven by a common entangled field. We showthat in the steady state the entanglement of the driving field is reproduced in an arbitrarily large series of inter-array entangled pairs over all distances. Local nonclassical driving thus realizes a scale-free entanglement replication and long-distance entanglement distribution mechanism that has immediate bearing on the implementation of quantum communication networks. 

Assessing the Nonequilibrium Thermodynamics in a Quenched Quantum Many-Body System via Single Projective Measurements

We analyze the nature of the statistics of the work done on or by a quantum many-body system brought out of equilibrium. We show that, for the sudden quench and for an initial state that commutes with the initial Hamiltonian, it is possible to retrieve the whole nonequilibrium thermodynamics via single projective measurements of observables. We highlight, in a physically clear way, the qualitative implications for the statistics of work coming from considering processes described by operators that either commute or do not commute with the unperturbed Hamiltonian of a given system. We consider a quantum many-body system and derive an expression that allows us to give a physical interpretation, for a thermal initial state, to all of the cumulants of the work in the case of quenched operators commuting with the unperturbed Hamiltonian. In the commuting case, the observables that we need to measure have an intuitive physical meaning. Conversely, in the noncommuting case, we show that, although it is possible to operate fully within the single-measurement framework irrespectively of the size of the quench, some difficulties are faced in providing a clear-cut physical interpretation to the cumulants. This circumstance makes the study of the physics of the system nontrivial and highlights the nonintuitive phenomenology of the emergence of thermodynamics from the fully quantum microscopic description. We illustrate our ideas with the example of the Ising model in a transverse field showing the interesting behavior of the high-order statistical moments of the work distribution for a generic thermal state and linking them to the critical nature of the model itself.