Basis set construction for molecular electronic structure theory: Natural orbital and Gauss-Slater basis for smooth pseudpotentials

A simple yet general method for constructing basis sets for molecular electronic structure calculations is presented. These basis sets consist of atomic natural orbitals from a multi-configurational self-consistent field calculation supplemented with primitive functions, chosen such that the asymptotics are appropriate for the potential of the system. Primitives are optimized for the homonuclear diatomic molecule to produce a balanced basis set. Two general features that facilitate this basis construction are demonstrated. First, weak coupling exists between the optimal exponents of primitives with different angular momenta. Second, the optimal primitive exponents for a chosen system depend weakly on the particular level of theory employed for optimization. The explicit case considered here is a basis set appropriate for the Burkatzki-Filippi-Dolg pseudopotentials. Since these pseudopotentials are finite at nuclei and have a Coulomb tail, the recently proposed Gauss-Slater functions are the appropriate primitives. Double- and triple-zeta bases are developed for elements hydrogen through argon. These new bases offer significant gains over the corresponding Burkatzki-Filippi-Dolg bases at various levels of theory. Using a Gaussian expansion of the basis functions, these bases can be employed in any electronic structure method. Quantum Monte Carlo provides an added benefit: expansions are unnecessary since the integrals are evaluated numerically.Comment: 9 pages, 7 figure

Towards a quantum-chemical description of crystalline insulators: A Wannier-function-based Hartree-Fock study of Li2O and Na2O

A recently proposed approach for performing electronic-structure calculations on crystalline insulators in terms of localized orthogonal orbitals is applied to the oxides of lithium and sodium, Li2O and Na2O. Cohesive energies, lattice constants and bulk moduli of the aforementioned systems are determined at the Hartree-Fock level, and the corresponding values are shown to be in excellent agreement with the values obtained by a traditional Bloch-orbital-based Hartree-Fock approach. The present Wannier-function-based approach is expected to be advantageous in the treatment of electron-correlation effects in an infinite solid by conventional quantum-chemical methods.Comment: 15 Pages, RevTex, 3 postscript figures (included), to appear in the Journal of Chemical Physics, May 15, 199

Ab initio many-body calculations on infinite carbon and boron-nitrogen chains

In this paper we report first-principles calculations on the ground-state electronic structure of two infinite one-dimensional systems: (a) a chain of carbon atoms and (b) a chain of alternating boron and nitrogen atoms. Meanfield results were obtained using the restricted Hartree-Fock approach, while the many-body effects were taken into account by second-order M{\o}ller-Plesset perturbation theory and the coupled-cluster approach. The calculations were performed using 6-31$G^{**}$ basis sets, including the d-type polarization functions. Both at the Hartree-Fock (HF) and the correlated levels we find that the infinite carbon chain exhibits bond alternation with alternating single and triple bonds, while the boron-nitrogen chain exhibits equidistant bonds. In addition, we also performed density-functional-theory-based local density approximation (LDA) calculations on the infinite carbon chain using the same basis set. Our LDA results, in contradiction to our HF and correlated results, predict a very small bond alternation. Based upon our LDA results for the carbon chain, which are in agreement with an earlier LDA calculation calculation [ E.J. Bylaska, J.H. Weare, and R. Kawai, Phys. Rev. B 58, R7488 (1998).], we conclude that the LDA significantly underestimates Peierls distortion. This emphasizes that the inclusion of many-particle effects is very important for the correct description of Peierls distortion in one-dimensional systems.Comment: 3 figures (included). To appear in Phys. Rev.

Cohesive energies of cubic III-V semiconductors

Cohesive energies for twelve cubic III-V semiconductors with zincblende structure have been determined using an ab-initio scheme. Correlation contributions, in particular, have been evaluated using the coupled-cluster approach with single and double excitations (CCSD). This was done by means of increments obtained for localized bond orbitals and for pairs and triples of such bonds. Combining these results with corresponding Hartree-Fock data, we recover about 92 \% of the experimental cohesive energies.Comment: 16 pages, 1 figure, late

A priori Wannier functions from modified Hartree-Fock and Kohn-Sham equations

The Hartree-Fock equations are modified to directly yield Wannier functions following a proposal of Shukla et al. [Chem. Phys. Lett. 262, 213-218 (1996)]. This approach circumvents the a posteriori application of the Wannier transformation to Bloch functions. I give a novel and rigorous derivation of the relevant equations by introducing an orthogonalizing potential to ensure the orthogonality among the resulting functions. The properties of these, so-called a priori Wannier functions, are analyzed and the relation of the modified Hartree-Fock equations to the conventional, Bloch-function-based equations is elucidated. It is pointed out that the modified equations offer a different route to maximally localized Wannier functions. Their computational solution is found to involve an effort that is comparable to the effort for the solution of the conventional equations. Above all, I show how a priori Wannier functions can be obtained by a modification of the Kohn-Sham equations of density-functional theory.Comment: 7 pages, RevTeX4, revise

A Hartree-Fock ab initio band-structure calculation employing Wannier-type orbitals

An ab initio Wannier-function-based approach to electronic ground-state calculations for crystalline solids is outlined. In the framework of the linear combination of atomic orbitals method the infinite character of the solid is rigorously taken into account. The Hartree-Fock ground-state energy, cohesive energy, lattice constant and bulk modulus are calculated in a fully ab initio manner as it is demonstrated for sodium chloride, NaCl, using basis sets close to the Hartree-Fock limit. It is demonstrated that the Hartree-Fock band-structure can easily be recovered with the current approach and agrees with the one obtained from a more conventional Bloch-orbital-based calculation. It is argued that the advantage of the present approach lies in its capability to include electron correlation effects for crystalline insulators by means of well-established quantum chemical procedures.Comment: 15 Pages, LaTex, 1 postscript figure (included), to appear in Chem. Phys. Letters (1998

Comparison of Smooth Hartree-Fock Pseudopotentials

The accuracy of two widely used scalar-relativistic Hartree-Fock pseudopotentials, the Trail-Needs-Dirac-Fock (TNDF) and the Burkatzki-Filippi-Dolg (BFD) pseudopotentials, is assessed. The performance of the pseudopotentials is tested for a chemically representative set of 34 first-row molecules. All comparisons are made at the Hartree-Fock level of theory, and both sets of pseudopotentials give good results. The all-electron equilibrium geometries, molecular dissociation energies, and zero-point vibrational energies are reproduced a little more accurately by the TNDF pseudopotentials than the BFD ones.Comment: 9 pages, 3 figures, 1 tabl

Ab initio Wannier-function-based correlated calculations of Born effective charges of crystalline Li2_{2}O and LiCl

In this paper we have used our recently developed ab initio Wannier-function-based methodology to perform extensive Hartree-Fock and correlated calculations on Li$_{2}$O and LiCl to compute their Born effective charges. Results thus obtained are in very good agreement with the experiments. In particular, for the case of Li$_{2}$O, we resolve a controversy originating in the experiment of Osaka and Shindo {[}Solid State Commun. 51 (1984) 421] who had predicted the effective charge of Li ions to be in the range 0.58--0.61, a value much smaller compared to its nominal value of unity, thereby, suggesting that the bonding in the material could be partially covalent. We demonstrate that effective charge computed by Osaka and Shindo is the Szigeti charge, and once the Born charge is computed, it is in excellent agreement with our computed value. Mulliken population analysis of Li$_{2}$O also confirms ionic nature of the bonding in the substance.Comment: 11 pages, 1 figure. To appear in Phys. Rev. B (Feb 2008