The self-consistent perturbation theory developed in earlier papers is extended to the open-shell case. Density matrices for both shells are calculated iteratively until first-order self-consistency is achieved. A numerical application indicates the importance of considering both shells separately when discussing the effects of polarization on the charge density and the spin density.
High-precision, large-basis-set calculations, in the local-density approximation (LDA) (using the allelectron, full-potential, linear combination of Gaussian orbitals, fitting-function technique), of the cohesive properties and electronic states (bare Kohn-Sham energies) of the isolated AB dilayer of graphite are reported. They show that the dilayer interplanar spacing (c axis) differs little from the value for ABABAB . crystalline graphite (0.7%%uo expansion relative to one calculation, 2.5%%uo contraction relative to another, 2% expansion relative to experiment). This result, which differs significantly from a preliminary report of strong c-axis contraction, is related to the weak interplanar binding. The intraplanar lattice spacing (a axis) is virtually identical with the crystalline value for both the graphite dilayer and monolayer. The interplanar binding energy (obtained directly via optimization of the monolayer ground state with the same techniques) is in excellent (perhaps fortuitous) agreement with the experimental value for the crystal, in contrast with crystalline calculations, which are too large (in magnitude) by 40 -100% or more. The dilayer cohesive energy agrees well with the crystalline value from an allelectron calculation. Both exceed the experimental value in magnitude by over 1 eV/atom, a problem already known to arise from inadequacies in the LDA treatment of the multiplet structure of the isolated C atom. The dilayer uniaxial compressibility is much larger than calculated for the crystal, apparently another manifestation of weak interplanar binding. Dilayer Kohn-Sham eigenvalues are largely consistent with those calculated self-consistently for the crystal using the same LDA model. Both differ substantially from the non-self-consistent band structure commonly used to parametrize graphite optical properties of interest in astrophysics. Calculated values of the dilayer work function are larger by about 0.6-0.7 eV than the crystalline experimental results. The dilayer density of states at the Fermi level is predicted to be much smaller than for the crystal, while the occupied bandwidth is in reasonable agreement with crystalline experimental results.
Limited configuration interaction methods suffer from size-extensivity errors. The origin and behavior of these errors is discussed and new versions of single and multireference corrections are presented. Accuracy of the new and various other size-extensivity corrections used in the literature is discussed and compared in a series of model calculations and calculations on small molecules. None of the commonly used multireference corrections restores the size extensivity of multireference configuration interaction calculations. Our correction behaves correctly for the special case of a reference state composed from all singly and doubly excited configurations. Formulas for size extensivity corrections in the variational-perturbation method are given and discussed.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.