Exchange and correlation as a functional of the local density of states

A functional $E_{xc}[\rho(\r,\epsilon)]$ is presented, in which the exchange and correlation energy of an electron gas depends on the local density of occupied states. A simple local parametrization scheme is proposed, entirely from first principles, based on the decomposition of the exchange-correlation hole in scattering states of different relative energies. In its practical Kohn-Sham-like form, the single-electron orbitals become the independent variables, and an explicit formula for the functional derivative is obtained.Comment: 5 pages. Expanded version. Will appear in Phys. Rev.

Spin Resolution of the Electron-Gas Correlation Energy: Positive same-spin contribution

The negative correlation energy per particle of a uniform electron gas of density parameter $r_s$ and spin polarization $\zeta$ is well known, but its spin resolution into up-down, up-up, and down-down contributions is not. Widely-used estimates are incorrect, and hamper the development of reliable density functionals and pair distribution functions. For the spin resolution, we present interpolations between high- and low-density limits that agree with available Quantum Monte Carlo data. In the low-density limit for $\zeta = 0$, we find that the same-spin correlation energy is unexpectedly positive, and we explain why. We also estimate the up and down contributions to the kinetic energy of correlation.Comment: new version, to appear in PRB Rapid Communicatio

Density-functional calculation of ionization energies of current-carrying atomic states

Current-density-functional theory is used to calculate ionization energies of current-carrying atomic states. A perturbative approximation to full current-density-functional theory is implemented for the first time, and found to be numerically feasible. Different parametrizations for the current-dependence of the density functional are critically compared. Orbital currents in open-shell atoms turn out to produce a small shift in the ionization energies. We find that modern density functionals have reached an accuracy at which small current-related terms appearing in open-shell configurations are not negligible anymore compared to the remaining difference to experiment.Comment: 7 pages, 2 tables, accepted by Phys. Rev.

Comparative study of density functional theories of the exchange-correlation hole and energy in silicon

We present a detailed study of the exchange-correlation hole and exchange-correlation energy per particle in the Si crystal as calculated by the Variational Monte Carlo method and predicted by various density functional models. Nonlocal density averaging methods prove to be successful in correcting severe errors in the local density approximation (LDA) at low densities where the density changes dramatically over the correlation length of the LDA hole, but fail to provide systematic improvements at higher densities where the effects of density inhomogeneity are more subtle. Exchange and correlation considered separately show a sensitivity to the nonlocal semiconductor crystal environment, particularly within the Si bond, which is not predicted by the nonlocal approaches based on density averaging. The exchange hole is well described by a bonding orbital picture, while the correlation hole has a significant component due to the polarization of the nearby bonds, which partially screens out the anisotropy in the exchange hole.Comment: 16 pages, 5 figures, RevTeX, added conten

Carbon clusters near the crossover to fullerene stability

The thermodynamic stability of structural isomers of $\mathrm{C}_{24}$, $\mathrm{C}_{26}$, $\mathrm{C}_{28}$ and $\mathrm{C}_{32}$, including fullerenes, is studied using density functional and quantum Monte Carlo methods. The energetic ordering of the different isomers depends sensitively on the treatment of electron correlation. Fixed-node diffusion quantum Monte Carlo calculations predict that a $\mathrm{C}_{24}$ isomer is the smallest stable graphitic fragment and that the smallest stable fullerenes are the $\mathrm{C}_{26}$ and $\mathrm{C}_{28}$ clusters with $\mathrm{C}_{2v}$ and $\mathrm{T}_{d}$ symmetry, respectively. These results support proposals that a $\mathrm{C}_{28}$ solid could be synthesized by cluster deposition.Comment: 4 pages, includes 4 figures. For additional graphics, online paper and related information see http://www.tcm.phy.cam.ac.uk/~prck

Nucleation of a sodium droplet on C60

We investigate theoretically the progressive coating of C60 by several sodium atoms. Density functional calculations using a nonlocal functional are performed for NaC60 and Na2C60 in various configurations. These data are used to construct an empirical atomistic model in order to treat larger sizes in a statistical and dynamical context. Fluctuating charges are incorporated to account for charge transfer between sodium and carbon atoms. By performing systematic global optimization in the size range 1<=n<=30, we find that Na_nC60 is homogeneously coated at small sizes, and that a growing droplet is formed above n=>8. The separate effects of single ionization and thermalization are also considered, as well as the changes due to a strong external electric field. The present results are discussed in the light of various experimental data.Comment: 17 pages, 10 figure

Fatty Acid Methyl Esters as Biosolvents of Epoxy Resins: A Physicochemical Study

The C8 to C18 fatty acid methyl esters (FAME) have been compared as solvents for two epoxy resin pre-polymers, bisphenol A diglycidyl ether (DGEBA) and triglycidyl paminophenol ether (TGPA). It was found that the solubilization limits vary according to the ester and that methyl caprylate is the best solvent of both resins. To explain these solubility performances, physical and chemical properties of FAME were studied, such as the Hansen parameters, viscosity, binary diffusion coefficient and vaporization enthalpy. Determination of the physicochemical parameters of FAME was carried out by laboratory experimentations and by calculation from bibliographic data. The Hansen parameters of FAME and epoxy resins pre-polymers were theoretically and experimentally determined. The FAME chain length showed a long dependence on the binary diffusion parameters and kinematic viscosity, which are mass and momentum transport properties. Moreover, the vaporization enthalpy of these compounds was directly correlated with the solubilization limits

Ab Initio Evidence for the Formation of Impurity d(3z^2-r^2) Holes in Doped La_{2-x}Sr_xCuO_4

Using the spin unrestricted Becke-3-Lee-Yang-Parr density functional, we computed the electronic structure of explicitly doped La_{2-x}Sr_xCuO_4 (x = 0.125, 0.25, and 0.5). At each doping level, an impurity hole band is formed within the undoped insulating gap. This band is well-localized to CuO_6 octahedra adjacent to the Sr impurities. The nature of the impurity hole is A_{1g} in symmetry, formed primarily from the z^2 orbital on the Cu and p_z orbitals on the apical O's. There is a strong triplet coupling of this hole with the intrinsic B_{1g} Cu x^2-y^2/O1 p_{sigma} hole on the same site. Optimization of the c coordinate of the apical O's in the doped CuO_6 octahedron lead to an asymmetric anti-Jahn-Teller distortion of the O2 atoms toward the central Cu. In particular, the O2 atom between the Cu and Sr is displaced 0.26 A while the O2 atom between the Cu and La is displaced 0.10 A. Contrary to expectations, investigation of a 0.1 A enhanced Jahn-Teller distortion of this octahedron does not force formation of an x^2-y^2 hole, but instead leads to migration of the z^2 hole to the four other CuO_6 octahedra surrounding the Sr impurity. This latter observation offers a simple explanation for the bifurcation of the Sr-O2 distance revealed in x-ray absorption fine structure data.Comment: Submitted to Phys. Rev. B. See http://www.firstprinciples.com for more informatio

Electronic structure of CaCuO2 from the B3LYP hybrid density functional

Published versio

Crystal Structures and Electronic Properties of Haloform-Intercalated C60

Using density functional methods we calculated structural and electronic properties of bulk chloroform and bromoform intercalated C60, C60 2CHX3 (X=Cl,Br). Both compounds are narrow band insulator materials with a gap between valence and conduction bands larger than 1 eV. The calculated widths of the valence and conduction bands are 0.4-0.6 eV and 0.3-0.4 eV, respectively. The orbitals of the haloform molecules overlap with the $\pi$ orbitals of the fullerene molecules and the p-type orbitals of halogen atoms significantly contribute to the valence and conduction bands of C60 2CHX3. Charging with electrons and holes turns the systems to metals. Contrary to expectation, 10 to 20 % of the charge is on the haloform molecules and is thus not completely localized on the fullerene molecules. Calculations on different crystal structures of C60 2CHCl3 and C60 2CHBr3 revealed that the density of states at the Fermi energy are sensitive to the orientation of the haloform and C60 molecules. At a charging of three holes, which corresponds to the superconducting phase of pure C60 and C60 2CHX3, the calculated density of states (DOS) at the Fermi energy increases in the sequence DOS(C60) < DOS(C60 2CHCl3) < DOS(C60 2CHBr3).Comment: 11 pages, 7 figures, 4 table