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

Non-Hermitian Rayleigh-Schroedinger Perturbation Theory

We devise a non-Hermitian Rayleigh-Schroedinger perturbation theory for the single- and the multireference case to tackle both the many-body problem and the decay problem encountered, for example, in the study of electronic resonances in molecules. A complex absorbing potential (CAP) is employed to facilitate a treatment of resonance states that is similar to the well-established bound-state techniques. For the perturbative approach, the full CAP-Schroedinger Hamiltonian, in suitable representation, is partitioned according to the Epstein-Nesbet scheme. The equations we derive in the framework of the single-reference perturbation theory turn out to be identical to those obtained by a time-dependent treatment in Wigner-Weisskopf theory. The multireference perturbation theory is studied for a model problem and is shown to be an efficient and accurate method. Algorithmic aspects of the integration of the perturbation theories into existing ab initio programs are discussed, and the simplicity of their implementation is elucidated.Comment: 10 pages, 1 figure, RevTeX4, submitted to Physical Review

Quasiparticle band structure of infinite hydrogen fluoride and hydrogen chloride chains

We study the quasiparticle band structure of isolated, infinite HF and HCl bent (zigzag) chains and examine the effect of the crystal field on the energy levels of the constituent monomers. The chains are one of the simplest but realistic models of the corresponding three-dimensional crystalline solids. To describe the isolated monomers and the chains, we set out from the Hartree-Fock approximation, harnessing the advanced Green's function methods "local molecular orbital algebraic diagrammatic construction" (ADC) scheme and "local crystal orbital ADC" (CO-ADC) in a strict second order approximation, ADC(2,2) and CO-ADC(2,2), respectively, to account for electron correlations. The configuration space of the periodic correlation calculations is found to converge rapidly only requiring nearest-neighbor contributions to be regarded. Although electron correlations cause a pronounced shift of the quasiparticle band structure of the chains with respect to the Hartree-Fock result, the bandwidth essentially remains unaltered in contrast to, e.g., covalently bound compounds.Comment: 11 pages, 6 figures, 6 tables, RevTeX4, corrected typoe

Ab initio Green's function formalism for band structures

Using the Green's function formalism, an ab initio theory for band structures of crystals is derived starting from the Hartree-Fock approximation. It is based on the algebraic diagrammatic construction scheme for the self-energy which is formulated for crystal orbitals (CO-ADC). In this approach, the poles of the Green's function are determined by solving a suitable Hermitian eigenvalue problem. The method is not only applicable to the outer valence and conduction bands, it is also stable for inner valence bands where strong electron correlations are effective. The key to the proposed scheme is to evaluate the self-energy in terms of Wannier orbitals before transforming it to a crystal momentum representation. Exploiting the fact that electron correlations are mainly local, one can truncate the lattice summations by an appropriate configuration selection scheme. This yields a flat configuration space; i.e., its size scales only linearly with the number of atoms per unit cell for large systems and, under certain conditions, the computational effort to determine band structures also scales linearly. As a first application of the new formalism, a lithium fluoride crystal has been chosen. A minimal basis set description is studied, and a satisfactory agreement with previous theoretical and experimental results for the fundamental band gap and the width of the F 2p valence band complex is obtained.Comment: 20 pages, 3 figures, 1 table, RevTeX4, new section on lithium fluorid