Diverging exchange force and form of the exact density matrix functional
Abstract
For translationally invariant one-band lattice models, we exploit the ab initio knowledge of the natural orbitals to simplify reduced density matrix functional theory (RDMFT). Striking underlying features are discovered: First, within each symmetry sector, the interaction functional F depends only on the natural occupation numbers n. The respective sets P1N and E1N of pure and ensemble N-representable one-matrices coincide. Second, and most importantly, the exact functional is strongly shaped by the geometry of the polytope E1N P1N , described by linear constraints D(j)(n)≥ 0. For smaller systems, it follows as F[n]=Σi,i' Vi,i' D(i)(n)D(i')(n). This generalizes to systems of arbitrary size by replacing each D(i) by a linear combination of \D(j)(n)\ and adding a non-analytical term involving the interaction V. Third, the gradient dF/dn is shown to diverge on the boundary ∂E1N, suggesting that the fermionic exchange symmetry manifests itself within RDMFT in the form of an "exchange force". All findings hold for systems with non-fixed particle number as well and V can be any p-particle interaction. As an illustration, we derive the exact functional for the Hubbard square.
Turn this paper into a lesson
ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.