Non-empirical hyper-generalized-gradient functionals constructed from the Lieb-Oxford bound

A simple and completely general representation of the exact exchange-correlation functional of density-functional theory is derived from the universal Lieb-Oxford bound, which holds for any Coulomb-interacting system. This representation leads to an alternative point of view on popular hybrid functionals, providing a rationale for why they work and how they can be constructed. A similar representation of the exact correlation functional allows to construct fully non-empirical hyper-generalized-gradient approximations (HGGAs), radically departing from established paradigms of functional construction. Numerical tests of these HGGAs for atomic and molecular correlation energies and molecular atomization energies show that even simple HGGAs match or outperform state-of-the-art correlation functionals currently used in solid-state physics and quantum chemistry.Comment: v2: Major revison. Added information on relation to the gradient expansion and to local hybrids, improved discussion of size consistency and of performance relative to other functional

How tight is the Lieb-Oxford bound?

Density-functional theory requires ever better exchange-correlation (xc) functionals for the ever more precise description of many-body effects on electronic structure. Universal constraints on the xc energy are important ingredients in the construction of improved functionals. Here we investigate one such universal property of xc functionals: the Lieb-Oxford lower bound on the exchange-correlation energy, $E_{xc}[n] \ge -C \int d^3r n^{4/3}$, where $C\leq C_{LO}=1.68$. To this end, we perform a survey of available exact or near-exact data on xc energies of atoms, ions, molecules, solids, and some model Hamiltonians (the electron liquid, Hooke's atom and the Hubbard model). All physically realistic density distributions investigated are consistent with the tighter limit $C \leq 1$. For large classes of systems one can obtain class-specific (but not fully universal) similar bounds. The Lieb-Oxford bound with $C_{LO}=1.68$ is a key ingredient in the construction of modern xc functionals, and a substantial change in the prefactor $C$ will have consequences for the performance of these functionals.Comment: 10 pages, 3 figure

Preparation of monotectic alloys having a controlled microstructure by directional solidification under dopant-induced interface breakdown

Monotectic alloys having aligned spherical particles of rods of the minor component dispersed in a matrix of the major component are prepared by forming a melt containing predetermined amounts of the major and minor components of a chosen monotectic system, providing in the melt a dopant capable of breaking down the liquid solid interface for the chosen alloy, and directionally solidfying the melt at a selected temperature gradient and a selected rate of movement of the liquid-solid interface (growth rate). Shaping of the minor component into spheres or rods and the spacing between them are controlled by the amount of dopant and the temperature gradient and growth rate values. Specific alloy systems include Al Bi, Al Pb and Zn Bi, using a transition element such as iron

Work functions of self-assembled monolayers on metal surfaces

Using first-principles calculations we show that the work function of noble metals can be decreased or increased by up to 2 eV upon the adsorption of self-assembled monolayers of organic molecules. We identify the contributions to these changes for several (fluorinated) thiolate molecules adsorbed on Ag(111), Au(111) and Pt(111) surfaces. The work function of the clean metal surfaces increases in this order, but adsorption of the monolayers reverses the order completely. Bonds between the thiolate molecules and the metal surfaces generate an interface dipole, whose size is a function of the metal, but it is relatively independent of the molecules. The molecular and bond dipoles can then be added to determine the overall work function.Comment: 5 pages, 2 figure

Ground-state degeneracies leave recognizable topological scars in the one-particle density

In Kohn-Sham density functional theory (KS-DFT) a fictitious system of non-interacting particles is constructed having the same ground-state (GS) density as the physical system of interest. A fundamental open question in DFT concerns the ability of an exact KS calculation to spot and characterize the GS degeneracies in the physical system. In this article we provide theoretical evidence suggesting that the GS density, as a function of position on a 2D manifold of parameters affecting the external potential, is "topologically scarred" in a distinct way by degeneracies. These scars are sufficiently detailed to enable determination of the positions of degeneracies and even the associated Berry phases. We conclude that an exact KS calculation can spot and characterize the degeneracies of the physical system

Adiabatic connection at negative coupling strengths

The adiabatic connection of density functional theory (DFT) for electronic systems is generalized here to negative values of the coupling strength $\alpha$ (with {\em attractive} electrons). In the extreme limit $\alpha\to-\infty$ a simple physical solution is presented and its implications for DFT (as well as its limitations) are discussed. For two-electron systems (a case in which the present solution can be calculated exactly), we find that an interpolation between the limit $\alpha\to-\infty$ and the opposite limit of infinitely strong repulsion ($\alpha\to+\infty$) yields a rather accurate estimate of the second-order correlation energy E\cor\glt[\rho] for several different densities $\rho$, without using virtual orbitals. The same procedure is also applied to the Be isoelectronic series, analyzing the effects of near-degeneracy.Comment: 9 pages, submitted to PR

Density-density functionals and effective potentials in many-body electronic structure calculations

We demonstrate the existence of different density-density functionals designed to retain selected properties of the many-body ground state in a non-interacting solution starting from the standard density functional theory ground state. We focus on diffusion quantum Monte Carlo applications that require trial wave functions with optimal Fermion nodes. The theory is extensible and can be used to understand current practices in several electronic structure methods within a generalized density functional framework. The theory justifies and stimulates the search of optimal empirical density functionals and effective potentials for accurate calculations of the properties of real materials, but also cautions on the limits of their applicability. The concepts are tested and validated with a near-analytic model.Comment: five figure

Exact results for a charged, harmonically trapped quantum gas at arbitrary temperature and magnetic field strength

An analytical expression for the first-order density matrix of a charged, two-dimensional, harmonically confined quantum gas, in the presence of a constant magnetic field is derived. In contrast to previous results available in the literature, our expressions are exact for any temperature and magnetic field strength. We also present a novel factorization of the Bloch density matrix in the form of a simple product with a clean separation of the zero-field and field-dependent parts. This factorization provides an alternative way of analytically investigating the effects of the magnetic field on the system, and also permits the extension of our analysis to other dimensions, and/or anisotropic confinement.Comment: To appear in Phys. Rev.

Kohn-Sham Exchange Potential for a Metallic Surface

The behavior of the surface barrier that forms at the metal-vacuum interface is important for several fields of surface science. Within the Density Functional Theory framework, this surface barrier has two non-trivial components: exchange and correlation. Exact results are provided for the exchange component, for a jellium metal-vacuum interface, in a slab geometry. The Kohn-Sham exact-exchange potential $V_{x}(z)$ has been generated by using the Optimized Effective Potential method, through an accurate numerical solution, imposing the correct boundary condition. It has been proved analytically, and confirmed numerically, that $V_{x}(z\to \infty)\to - e^{2}/z$; this conclusion is not affected by the inclusion of correlation effects. Also, the exact-exchange potential develops a shoulder-like structure close to the interface, on the vacuum side. The issue of the classical image potential is discussed.Comment: Phys. Rev. Lett. (to appear

Low energy electron scattering from DNA and RNA bases: shape resonances and radiation damage

Calculations are carried out to determine elastic scattering cross sections and resonance energies for low energy electron impact on uracil and on each of the DNA bases (thymine, cytosine, adenine, guanine), for isolated molecules in their equilibrium geometry. Our calculations are compared with available theory and experiment. We also attempt to correlate this information with experimental dissociation patterns through an analysis of the temporary anion structures that are formed by electron capture in shape resonances.Comment: 20 pages, 12 figures, submitted to J. Chem. Phy