The correlation energy functional within the GW-RPA approximation: exact forms, approximate forms and challenges

In principle, the Luttinger-Ward Green's function formalism allows one to compute simultaneously the total energy and the quasiparticle band structure of a many-body electronic system from first principles. We present approximate and exact expressions for the correlation energy within the GW-RPA approximation that are more amenable to computation and allow for developing efficient approximations to the self-energy operator and correlation energy. The exact form is a sum over differences between plasmon and interband energies. The approximate forms are based on summing over screened interband transitions. We also demonstrate that blind extremization of such functionals leads to unphysical results: imposing physical constraints on the allowed solutions (Green's functions) is necessary. Finally, we present some relevant numerical results for atomic systems.Comment: 3 figures and 3 tables, under review at Physical Review

A multigrid scheme for 3D Monge-Amp\`ere equations

The elliptic Monge-Amp\`ere equation is a fully nonlinear partial differential equation which has been the focus of increasing attention from the scientific computing community. Fast three dimensional solvers are needed, for example in medical image registration but are not yet available. We build fast solvers for smooth solutions in three dimensions using a nonlinear full-approximation storage multigrid method. Starting from a second-order accurate centered finite difference approximation, we present a nonlinear Gauss-Seidel iterative method which has a mechanism for selecting the convex solution of the equation. The iterative method is used as an effective smoother, combined with the full-approximation storage multigrid method. Numerical experiments are provided to validate the accuracy of the finite difference scheme and illustrate the computational efficiency of the proposed multigrid solver.Comment: 18 pages, 1 figure, 7 tables, 41 references. Accepted by International Journal of Computer Mathematics (published online: 21 Nov 2016

Efficient Algorithm for Asymptotics-Based Configuration-Interaction Methods and Electronic Structure of Transition Metal Atoms

Asymptotics-based configuration-interaction (CI) methods [G. Friesecke and B. D. Goddard, Multiscale Model. Simul. 7, 1876 (2009)] are a class of CI methods for atoms which reproduce, at fixed finite subspace dimension, the exact Schr\"odinger eigenstates in the limit of fixed electron number and large nuclear charge. Here we develop, implement, and apply to 3d transition metal atoms an efficient and accurate algorithm for asymptotics-based CI. Efficiency gains come from exact (symbolic) decomposition of the CI space into irreducible symmetry subspaces at essentially linear computational cost in the number of radial subshells with fixed angular momentum, use of reduced density matrices in order to avoid having to store wavefunctions, and use of Slater-type orbitals (STO's). The required Coulomb integrals for STO's are evaluated in closed form, with the help of Hankel matrices, Fourier analysis, and residue calculus. Applications to 3d transition metal atoms are in good agreement with experimental data. In particular we reproduce the anomalous magnetic moment and orbital filling of Chromium in the otherwise regular series Ca, Sc, Ti, V, Cr.Comment: 14 pages, 1 figur

Exchange-assisted tunneling in the classical limit

The exchange interaction and correlations may produce a power-law decay instead of the usual exponential decrease of the wave function under potential barrier. The exchange-assisted tunneling vanishes in the classical limit, however, the dependence on the Planck constant h is different from that for a conventional single-particle tunneling

A theoretical study of the C- 4So_3/2 and 2Do_{3/2,5/2} bound states and C ground configuration: fine and hyperfine structures, isotope shifts and transition probabilities

This work is an ab initio study of the 2p3 4So_3/2, and 2Do_{3/2,5/2} states of C- and 2p2 3P_{0,1,2}, 1D_2, and 1S_0 states of neutral carbon. We use the multi-configuration Hartree-Fock approach, focusing on the accuracy of the wave function itself. We obtain all C- detachment thresholds, including correlation effects to about 0.5%. Isotope shifts and hyperfine structures are calculated. The achieved accuracy of the latter is of the order of 0.1 MHz. Intra-configuration transition probabilities are also estimated.Comment: 15 pages, 2 figures, 12 table