Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Physical Review

Cumulant Expansion of Localized-Electron Model for Antiferromagnetic Insulators

Tadashi Arai

Published Jan 10, 1968• 6 citations
Publisher

Abstract

A many-body perturbation theory is developed for the calculation of the antiferromagnetic ground state and lower excited states of one-electron atoms in the insulating state. It is assumed that the interatomic distances are moderately small and that the interactions are well described by exchange coupling, but that distortion of atoms in the lattice is still negligible. The method is essentially the Heitler-London approach: A set of Slater determinants is constructed from nonorthogonal atomic orbitals with all possible spin arrangements. The unperturbed ground state $|0〉$ has an alternating spin configuration, while excited states are generated one by one from $|0〉$ by interchanges of spins and by projecting out the lower excited states. The resulting orthonormal states introduce the concept of quasielectron (or quasihole) states and define creation and destruction operators of those quasiparticles. Although the quasiparticles are correlated with each other through the overlap integrals in the projected wave functions, the energy matrix elements can be decomposed into linked clusters of localized electrons and holes. It is found that the clusters thus generated are cumulants. The expression for the energy matrix is then transformed into the form of an effective Hamiltonian consisting of cumulants and of creation and destruction operators of quasifermions, and the concept of the exchange interaction is defined in terms of those operators. The new formulation not only facilitates the use of the familiar many-body perturbation theory, but also eliminates difficulties in handling spin variables in the Heisenberg theory of magnetism. If only the two-body Coulomb and exchange diagrams are retained and the higher-order diagrams are neglected, however, the effective Hamiltonian takes a form reminiscent of the anisotropic exchange Hamiltonian, ${H}_{\mathrm{spin}}=\ensuremath{-}{\ensuremath{\Sigma}}_{k,h}{J}_{\mathrm{kh}}[{{S}_{k}}^{z}{{S}_{h}}^{z}+\ensuremath{\lambda}({{S}_{k}}^{x}{{S}_{h}}^{x}+{{S}_{k}}^{y}{{S}_{h}}^{y})]$, which is used in the generalized Heisenberg theory and which includes the Ising model in its limit $\ensuremath{\lambda}\ensuremath{\rightarrow}0$. Finally, it is shown that excitations accompanying electron transfers are readily included in the present method as additional diagrams by an extension of the prescriptions for the generation of excited spin states. This yields, in its first approximation, the intra-atomic interaction and the electron-transfer interaction, leading to the Hamiltonian assumed by Hubbard and Anderson.

DOI
10.1103/physrev.165.706

Related papers

Study of the excited states of the quantum antiferromagnets

Mohammad G. Merdan

We investigate the quantum dynamics of the spins on different Heisenberg antiferromagnetic spin lattice systems. Firstly, we applied the coupled-cluster method to the spin-1/2 antiferromagnetic \emph{XXZ} model on a square lattice by employing an approximation which contains two-body long-range correlations and high-order four-body local correlations. Improvement is found for the ground-state energy, sublattice magnetization, and the critical anisotropy when comparing with the approximation including the two-body correlations alone. We also obtain the full excitation spectrum which is in good agreement with the quantum Monte Carlo results and the high-order spin-wave theory.Secondly, we study the longitudinal excitations of quantum antiferromagnets on a triangular lattice by a recently proposed microscopic many-body approach based on magnon-density waves. We calculate the full longitudinal excitation spectra of the antiferromagnetic Heisenberg model for a general spin quantum number in the isotropic limit. Similar to the square lattice model, we find that, at the center of the first hexagonal Brillouin zone $\Gamma(\mathbf q=0)$ and at the magnetic ordering wavevectors $\pm[\mathbf Q= (4\pi/3,0)]$, the excitation spectra become gapless in the thermodynamic limit, due to the slow, logarithmic divergence of the structure factor. However, these longitudinal modes on two-dimensional models may be considered as quasi-gapped, as any finite-size effect or small anisotropy will induce a large energy gap, when compared with the counterpart of the transverse spin-wave excitations. We have also investigated the excited states of the quasi-one-dimensional quantum antiferromagnets on hexagonal lattices, including the longitudinal modes based on the magnon-density waves. A model Hamiltonian with a uniaxial single-ion anisotropy is first studied by a spin-wave theory based on the one-boson method; the ground state thus obtained is employed for the study of the longitudinal modes. The full energy spectra of both the transverse modes (i.e., magnons) and the longitudinal modes are obtained as functions of the nearest-neighbor coupling and the anisotropy constants. We have found two longitudinal modes due to the non-collinear nature of the triangular antiferromagnetic order, similar to that of the phenomenological field theory approach by Affleck. The excitation energy gaps due to the anisotropy and the energy gaps of the longitudinal modes without anisotropy are then investigated. We then compares our results for the longitudinal energy gaps at the magnetic wavevectors with the experimental results for several antiferromagnetic compounds with both integer and non-integer spin quantum numbers, and we find good agreements after the higher-order contributions are included in our calculations.

•Nov 10, 2013•0 citations

The Many Guises of Quasiparticles in Correlated Electron Systems

Ivan Pasqua

Many strongly correlated electron systems fall outside the conventional quasiparticle paradigm. Their low-energy fermionic excitations need not be adiabatically connected to the bare electrons or appear as poles of the physical electron single-particle Green's function. This thesis investigates to what extent an effective quasiparticle description can nevertheless survive, and which emergent degrees of freedom provide the appropriate low-energy structure. First, the interaction-driven evolution from a quantum spin Hall insulator to a Mott insulator is studied using the dynamical cluster approximation. Dispersive zeros of the single-particle Green's function emerge near the Mott transition and retain a nontrivial topological character. When symmetry breaking is allowed, a nontopological excitonic insulator is found to intrude between quantum spin Hall and Mott insulators. Encoding both poles and zeros within a low-energy quasiparticle Hamiltonian provides a continuous one-body description of the three phases. On the Mott side, its gap evolution supports interpreting the soft excitonic mode as a bound state of quasiparticle and quasihole excitations tied to the Green's function zeros, rather than to the widely separated Hubbard bands. Second, the intrinsic anomalous Hall conductivity of a topological metal is analyzed within a multiband extension of Landau Fermi-liquid theory. Residual interactions among quasiparticles at the Fermi surface dress the dynamic current vertex and generate corrections to the Berry-curvature contribution obtained from the quasiparticle bands alone. Our result supports recent claims that the correct expressions for topological observables include vertex corrections besides the topological invariants built just upon the Green's function. It also demonstrates that such corrections are naturally accounted for by Landau Fermi-liquid theory. Extracting the relevant ingredients for an effective description of material-relevant correlated systems requires computational tools that remain reliable at low temperatures. To this end, a deterministic strong-coupling impurity solver for dynamical mean-field theory is discussed. By providing static, dynamical, and thermodynamic observables without stochastic sampling, it offers a promising route toward multiorbital and spin-orbit coupled problems where conventional Monte Carlo methods currently struggle. Finally, the auxiliary quasiparticle framework underlying the ghost-Gutzwiller Ansatz is investigated. On a model for correlated quantum spin Hall insulators, we show that the topology of the interacting state is encoded in the auxiliary one-body Hamiltonian. The resulting band-structure reveals topological Hubbard bands whose topological character can be tuned through a finite magnetization. When applied to the single-band t-J model, the same framework yields a fractionalized Fermi liquid (FL*). In this phase, dispersive neutral spinons, with vanishing physical-electron spectral weight, coexist with conventional quasiparticles forming a small hole-like Fermi surface, thus violating Luttinger's theorem. The corresponding temperature--doping phase diagram contains a low-doping FL*, a d-wave superconducting dome, and an overdoped conventional Fermi liquid, thereby reproducing key qualitative features of cuprate phenomenology.

•Sep 14, 2026•0 citations

Many-body perturbation theory: The GW approximation

Christoph Friedrich, Arno Schindlmayr

In this lecture we present many-body perturbation theory as a method to determine quasiparticle excitations in solids, especially electronic band structures, accurately from first principles.The main ingredient is the electronic self-energy that, in principle, contains all many-body exchange and correlation effects beyond the Hartree potential.As its exact mathematical expression is unknown, approximations must be used in practical calculations.The approximation is obtained using a systematic algebraic approach on the basis of Green function techniques.It constitutes an expansion of the self-energy up to linear order in the screened Coulomb potential, which describes the interaction between the quasiparticles and includes dynamic screening through the creation of exchange-correlation holes around the bare particles.The implementation of the approximation relies on a perturbative treatment starting from density functional theory.Besides a detailed mathematical discussion we focus on the underlying physical concepts and show some illustrative applications.

•Jan 1, 2006•19 citations

The Study of Many-Electron Systems

Yu Ting Zhou

Various methods and approximation schemes are used to study many-electron interacting systems. Two important many-particle models, the Anderson model and the Hubbard model, and their electromagnetic properties have been investigated in many parameter regimes, and applied to physical systems. An Anderson single-impurity model Hamiltonian based calculation of the magnetic susceptibility is performed for YbN in the presence of crystal fields using an alteration of the Non-Crossing Approximation proposed by Zwicknagl et.al., incorporating parameters obtained from ab initio band structure calculations. It yields good agreement with experimental data. For the Anderson lattice model, a variational scheme which uses specific many-electron wavefunctions as basis is applied to both one- and two-dimensional systems represented by symmetric Anderson lattice Hamiltonians. Without much computational effort, the ground state energy is well approximated, especially in strong-coupling limit. Some electronic properties are examined using the variational ground state wavefunction. The one-dimensional Hubbard model has been solved exactly for small-size clusters by diagonalizing the Hamiltonian in the basis of many-electron Bloch states. The results for the energy spectrum and eigenfunctions of the ground state and low-lying excited states are presented. Also, mean field calculations of the two-dimensional single-band Hubbard model and Cu-O lattice model (three-band Hubbard model) are carried out for various physical quantities including the energy, occupation probability, staggered magnetization, momentum distribution Fermi surface and density of states, by using a projection operator formalism. To develop a systematic approach to solving many-electron problems, the many-particle partition function for the free electron gas system is explored using a cumulant expansion scheme. Starting from the ground state, the partition function can be approximated to any order in terms of excitation energy. Its application to interacting systems such as the Anderson model and the Hubbard model is briefly discussed.

•Feb 1, 1991•0 citations

Many-electron theory: Density functional approach generalized to treat spin eigenfunctions and relation to spinless low-order density matrices

N. H. March

Many-electron theory is treated here with much attention focused on the synthesis of the density functional and density matrix theory. This synthesis shows that the exchange–correlation potential, normally obtained as the functional derivative with respect to the electron density of an assumed approximate exchange–correlation energy, can be written exactly in terms of the “near-diagonal” behavior of the fully interacting first-order density matrix and the diagonal of the second-order density matrix, thereby bypassing the functional derivative. Also involved, since the exchange–correlation potential is an “artifact” of reducing the many-electron Schrödinger equation to one-body Slater–Kohn–Sham equations, are the one-body orbitals derived from these latter equations, not just the density thereby obtained. This situation has then motivated a discussion of what constraints are placed on the low-order density matrices of atoms and molecules, described by an assumed spin-independent Hamiltonian, by the requirement that these matrices are constructed from an N-electron wave function which is an eigenfunction of the spin operators S2 and Sz. Some briefer discussion is added on the optimally convergent determinantal expansion of such an N-electron wave function, the relation between Brueckner and Löwdin natural orbitals being commented on as a by-product. In this latter context, it is noted that the exchange-only potential of density functional theory can be represented as a sum of two line integrals, each of which is path-dependent, the sum being, of course, path-independent. But in addition to knowledge of Slater–Kohn–Sham orbitals appearing as in the Harbola–Sahni exchange-only potential (which is one of the path-dependent terms), there is a kinetic correction. This involves Löwdin natural orbitals and their occupation numbers, but only now in the weak coupling limit of the electron–electron interaction, as emphasized in the kinetic “correction” to the Harbola–Sahni formula exhibited by Levy and the writer. Some directions for further work are briefly set out. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 74: 163–176, 1999

International Journal of Quantum Chemistry•Jan 1, 1999•2 citations

Many-Body Perturbation Theory for Quantum Crystals. I

Heinz Horner

The possibility of a perturbation treatment for quantum crystals is investigated. It is found that the usual form of the unperturbed Hamiltonian, being a single-particle operator, is not suitable for this situation, since it does not produce correlations necessary in a theory of quantum crystals. These correlations prevent the unperturbed ground state or low excited states from containing multiply occupied lattice sites. An unperturbed Hamiltonian is proposed which generates these correlations but is no longer a single-particle operator. The rules for a diagrammatic representation of a perturbation expansion starting from this unperturbed Hamiltonian are given. They differ from the usual rules for Goldstone diagrams. In particular, the linked-cluster theorem is not valid. Adding and subtracting pairs of mutually cancelling diagrams, in a manner similar to that in which Pauli-principle-violating diagrams are treated in the standard perturbation theory for fermions, we find new rules which, at least at low temperatures, are identical to those for spinless fermions. This result is essentially independent of the statistics of the actual lattice particles. In the case where the lattice particles are actually fermions with spin \textonehalf{}, e.g., $^{3}\mathrm{He}$, we find an expansion formally identical to the problem of spinless fermions interacting among themselves and with a set of spins \textonehalf{} localized at lattice sites.

Physical Review A•Jun 1, 1970•58 citations
A many-body perturbation theory is developed for the calculation of the antiferromagnetic ground state and lower excited states of one-electron atoms in the insulating state. It is assumed that the interatomic distances are moderately small and that the interactions are well described by exchange coupling, but that distortion of atoms in the lattice is still negligible. The method is essentially the Heitler-London approach: A set of Slater determinants is constructed from nonorthogonal atomic orbitals with all possible spin arrangements. The unperturbed ground state ∣0〉|0〉∣0〉 has an alternating spin configuration, while excited states are generated one by one from ∣0〉|0〉∣0〉 by interchanges of spins and by projecting out the lower excited states. The resulting orthonormal states introduce the concept of quasielectron (or quasihole) states and define creation and destruction operators of those quasiparticles. Although the quasiparticles are correlated with each other through the overlap integrals in the projected wave functions, the energy matrix elements can be decomposed into linked clusters of localized electrons and holes. It is found that the clusters thus generated are cumulants. The expression for the energy matrix is then transformed into the form of an effective Hamiltonian consisting of cumulants and of creation and destruction operators of quasifermions, and the concept of the exchange interaction is defined in terms of those operators. The new formulation not only facilitates the use of the familiar many-body perturbation theory, but also eliminates difficulties in handling spin variables in the Heisenberg theory of magnetism. If only the two-body Coulomb and exchange diagrams are retained and the higher-order diagrams are neglected, however, the effective Hamiltonian takes a form reminiscent of the anisotropic exchange Hamiltonian, Hspin=\ensuremath−\ensuremathΣk,hJkh[SkzShz+\ensuremathλ(SkxShx+SkyShy)]{H}_{\mathrm{spin}}=\ensuremath{-}{\ensuremath{\Sigma}}_{k,h}{J}_{\mathrm{kh}}[{{S}_{k}}^{z}{{S}_{h}}^{z}+\ensuremath{\lambda}({{S}_{k}}^{x}{{S}_{h}}^{x}+{{S}_{k}}^{y}{{S}_{h}}^{y})]Hspin​=\ensuremath−\ensuremathΣk,h​Jkh​[Sk​zSh​z+\ensuremathλ(Sk​xSh​x+Sk​ySh​y)], which is used in the generalized Heisenberg theory and which includes the Ising model in its limit \ensuremathλ\ensuremath→0\ensuremath{\lambda}\ensuremath{\rightarrow}0\ensuremathλ\ensuremath→0. Finally, it is shown that excitations accompanying electron transfers are readily included in the present method as additional diagrams by an extension of the prescriptions for the generation of excited spin states. This yields, in its first approximation, the intra-atomic interaction and the electron-transfer interaction, leading to the Hamiltonian assumed by Hubbard and Anderson.
We investigate the quantum dynamics of the spins on different Heisenberg antiferromagnetic spin lattice systems. Firstly, we applied the coupled-cluster method to the spin-1/2 antiferromagnetic \emph{XXZ} model on a square lattice by employing an approximation which contains two-body long-range correlations and high-order four-body local correlations. Improvement is found for the ground-state energy, sublattice magnetization, and the critical anisotropy when comparing with the approximation including the two-body correlations alone. We also obtain the full excitation spectrum which is in good agreement with the quantum Monte Carlo results and the high-order spin-wave theory.Secondly, we study the longitudinal excitations of quantum antiferromagnets on a triangular lattice by a recently proposed microscopic many-body approach based on magnon-density waves. We calculate the full longitudinal excitation spectra of the antiferromagnetic Heisenberg model for a general spin quantum number in the isotropic limit. Similar to the square lattice model, we find that, at the center of the first hexagonal Brillouin zone Γ(q=0)\Gamma(\mathbf q=0)Γ(q=0) and at the magnetic ordering wavevectors ±[Q=(4π/3,0)]\pm[\mathbf Q= (4\pi/3,0)]±[Q=(4π/3,0)], the excitation spectra become gapless in the thermodynamic limit, due to the slow, logarithmic divergence of the structure factor. However, these longitudinal modes on two-dimensional models may be considered as quasi-gapped, as any finite-size effect or small anisotropy will induce a large energy gap, when compared with the counterpart of the transverse spin-wave excitations. We have also investigated the excited states of the quasi-one-dimensional quantum antiferromagnets on hexagonal lattices, including the longitudinal modes based on the magnon-density waves. A model Hamiltonian with a uniaxial single-ion anisotropy is first studied by a spin-wave theory based on the one-boson method; the ground state thus obtained is employed for the study of the longitudinal modes. The full energy spectra of both the transverse modes (i.e., magnons) and the longitudinal modes are obtained as functions of the nearest-neighbor coupling and the anisotropy constants. We have found two longitudinal modes due to the non-collinear nature of the triangular antiferromagnetic order, similar to that of the phenomenological field theory approach by Affleck. The excitation energy gaps due to the anisotropy and the energy gaps of the longitudinal modes without anisotropy are then investigated. We then compares our results for the longitudinal energy gaps at the magnetic wavevectors with the experimental results for several antiferromagnetic compounds with both integer and non-integer spin quantum numbers, and we find good agreements after the higher-order contributions are included in our calculations.
Many strongly correlated electron systems fall outside the conventional quasiparticle paradigm. Their low-energy fermionic excitations need not be adiabatically connected to the bare electrons or appear as poles of the physical electron single-particle Green's function. This thesis investigates to what extent an effective quasiparticle description can nevertheless survive, and which emergent degrees of freedom provide the appropriate low-energy structure. First, the interaction-driven evolution from a quantum spin Hall insulator to a Mott insulator is studied using the dynamical cluster approximation. Dispersive zeros of the single-particle Green's function emerge near the Mott transition and retain a nontrivial topological character. When symmetry breaking is allowed, a nontopological excitonic insulator is found to intrude between quantum spin Hall and Mott insulators. Encoding both poles and zeros within a low-energy quasiparticle Hamiltonian provides a continuous one-body description of the three phases. On the Mott side, its gap evolution supports interpreting the soft excitonic mode as a bound state of quasiparticle and quasihole excitations tied to the Green's function zeros, rather than to the widely separated Hubbard bands. Second, the intrinsic anomalous Hall conductivity of a topological metal is analyzed within a multiband extension of Landau Fermi-liquid theory. Residual interactions among quasiparticles at the Fermi surface dress the dynamic current vertex and generate corrections to the Berry-curvature contribution obtained from the quasiparticle bands alone. Our result supports recent claims that the correct expressions for topological observables include vertex corrections besides the topological invariants built just upon the Green's function. It also demonstrates that such corrections are naturally accounted for by Landau Fermi-liquid theory. Extracting the relevant ingredients for an effective description of material-relevant correlated systems requires computational tools that remain reliable at low temperatures. To this end, a deterministic strong-coupling impurity solver for dynamical mean-field theory is discussed. By providing static, dynamical, and thermodynamic observables without stochastic sampling, it offers a promising route toward multiorbital and spin-orbit coupled problems where conventional Monte Carlo methods currently struggle. Finally, the auxiliary quasiparticle framework underlying the ghost-Gutzwiller Ansatz is investigated. On a model for correlated quantum spin Hall insulators, we show that the topology of the interacting state is encoded in the auxiliary one-body Hamiltonian. The resulting band-structure reveals topological Hubbard bands whose topological character can be tuned through a finite magnetization. When applied to the single-band t-J model, the same framework yields a fractionalized Fermi liquid (FL*). In this phase, dispersive neutral spinons, with vanishing physical-electron spectral weight, coexist with conventional quasiparticles forming a small hole-like Fermi surface, thus violating Luttinger's theorem. The corresponding temperature--doping phase diagram contains a low-doping FL*, a d-wave superconducting dome, and an overdoped conventional Fermi liquid, thereby reproducing key qualitative features of cuprate phenomenology.
In this lecture we present many-body perturbation theory as a method to determine quasiparticle excitations in solids, especially electronic band structures, accurately from first principles.The main ingredient is the electronic self-energy that, in principle, contains all many-body exchange and correlation effects beyond the Hartree potential.As its exact mathematical expression is unknown, approximations must be used in practical calculations.The approximation is obtained using a systematic algebraic approach on the basis of Green function techniques.It constitutes an expansion of the self-energy up to linear order in the screened Coulomb potential, which describes the interaction between the quasiparticles and includes dynamic screening through the creation of exchange-correlation holes around the bare particles.The implementation of the approximation relies on a perturbative treatment starting from density functional theory.Besides a detailed mathematical discussion we focus on the underlying physical concepts and show some illustrative applications.
Various methods and approximation schemes are used to study many-electron interacting systems. Two important many-particle models, the Anderson model and the Hubbard model, and their electromagnetic properties have been investigated in many parameter regimes, and applied to physical systems. An Anderson single-impurity model Hamiltonian based calculation of the magnetic susceptibility is performed for YbN in the presence of crystal fields using an alteration of the Non-Crossing Approximation proposed by Zwicknagl et.al., incorporating parameters obtained from ab initio band structure calculations. It yields good agreement with experimental data. For the Anderson lattice model, a variational scheme which uses specific many-electron wavefunctions as basis is applied to both one- and two-dimensional systems represented by symmetric Anderson lattice Hamiltonians. Without much computational effort, the ground state energy is well approximated, especially in strong-coupling limit. Some electronic properties are examined using the variational ground state wavefunction. The one-dimensional Hubbard model has been solved exactly for small-size clusters by diagonalizing the Hamiltonian in the basis of many-electron Bloch states. The results for the energy spectrum and eigenfunctions of the ground state and low-lying excited states are presented. Also, mean field calculations of the two-dimensional single-band Hubbard model and Cu-O lattice model (three-band Hubbard model) are carried out for various physical quantities including the energy, occupation probability, staggered magnetization, momentum distribution Fermi surface and density of states, by using a projection operator formalism. To develop a systematic approach to solving many-electron problems, the many-particle partition function for the free electron gas system is explored using a cumulant expansion scheme. Starting from the ground state, the partition function can be approximated to any order in terms of excitation energy. Its application to interacting systems such as the Anderson model and the Hubbard model is briefly discussed.
Many-electron theory is treated here with much attention focused on the synthesis of the density functional and density matrix theory. This synthesis shows that the exchange–correlation potential, normally obtained as the functional derivative with respect to the electron density of an assumed approximate exchange–correlation energy, can be written exactly in terms of the “near-diagonal” behavior of the fully interacting first-order density matrix and the diagonal of the second-order density matrix, thereby bypassing the functional derivative. Also involved, since the exchange–correlation potential is an “artifact” of reducing the many-electron Schrödinger equation to one-body Slater–Kohn–Sham equations, are the one-body orbitals derived from these latter equations, not just the density thereby obtained. This situation has then motivated a discussion of what constraints are placed on the low-order density matrices of atoms and molecules, described by an assumed spin-independent Hamiltonian, by the requirement that these matrices are constructed from an N-electron wave function which is an eigenfunction of the spin operators S2 and Sz. Some briefer discussion is added on the optimally convergent determinantal expansion of such an N-electron wave function, the relation between Brueckner and Löwdin natural orbitals being commented on as a by-product. In this latter context, it is noted that the exchange-only potential of density functional theory can be represented as a sum of two line integrals, each of which is path-dependent, the sum being, of course, path-independent. But in addition to knowledge of Slater–Kohn–Sham orbitals appearing as in the Harbola–Sahni exchange-only potential (which is one of the path-dependent terms), there is a kinetic correction. This involves Löwdin natural orbitals and their occupation numbers, but only now in the weak coupling limit of the electron–electron interaction, as emphasized in the kinetic “correction” to the Harbola–Sahni formula exhibited by Levy and the writer. Some directions for further work are briefly set out. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 74: 163–176, 1999
The possibility of a perturbation treatment for quantum crystals is investigated. It is found that the usual form of the unperturbed Hamiltonian, being a single-particle operator, is not suitable for this situation, since it does not produce correlations necessary in a theory of quantum crystals. These correlations prevent the unperturbed ground state or low excited states from containing multiply occupied lattice sites. An unperturbed Hamiltonian is proposed which generates these correlations but is no longer a single-particle operator. The rules for a diagrammatic representation of a perturbation expansion starting from this unperturbed Hamiltonian are given. They differ from the usual rules for Goldstone diagrams. In particular, the linked-cluster theorem is not valid. Adding and subtracting pairs of mutually cancelling diagrams, in a manner similar to that in which Pauli-principle-violating diagrams are treated in the standard perturbation theory for fermions, we find new rules which, at least at low temperatures, are identical to those for spinless fermions. This result is essentially independent of the statistics of the actual lattice particles. In the case where the lattice particles are actually fermions with spin \textonehalf{}, e.g., 3He^{3}\mathrm{He}3He, we find an expansion formally identical to the problem of spinless fermions interacting among themselves and with a set of spins \textonehalf{} localized at lattice sites.