A density-matrix derivation of the Hartree--Fock equations in a nonorthogonal atomic-orbital basis

Thomas Kj{\ae}rgaard

arxiv: 2605.03761 · v2 · submitted 2026-05-05 · 🪐 quant-ph

A density-matrix derivation of the Hartree--Fock equations in a nonorthogonal atomic-orbital basis

Thomas Kj{\ae}rgaard This is my paper

Pith reviewed 2026-05-12 03:27 UTC · model grok-4.3

classification 🪐 quant-ph

keywords Hartree-Fockdensity matrixnonorthogonal basisatomic orbitalssecond quantizationresponse theoryexponential parametrization

0 comments

The pith

The standard Hartree-Fock stationarity condition follows from the exponential parametrization of the one-particle density matrix in a nonorthogonal atomic-orbital basis.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper derives the Hartree-Fock equations in a nonorthogonal atomic-orbital basis by starting from the second-quantized density-matrix formalism. It shows that the usual stationarity condition arises directly once the one-particle density matrix receives an exponential parametrization. This creates a compact connection between elementary Hartree-Fock theory and the density-matrix tools employed in modern response calculations and linear-scaling methods. A sympathetic reader sees the derivation as a natural unification that avoids separate orthogonality assumptions.

Core claim

By parametrizing the one-particle density matrix exponentially inside the second-quantized atomic-orbital formalism, the Hartree-Fock energy stationarity condition is recovered as the requirement that the energy be stationary with respect to variations in that matrix, reproducing the standard AO Hartree-Fock equations without additional basis transformations.

What carries the argument

Exponential parametrization of the one-particle density matrix in the second-quantized atomic-orbital density-matrix formalism.

Load-bearing premise

The second-quantized atomic-orbital density-matrix formalism developed in the cited prior work is a valid and complete starting point for deriving the Hartree-Fock equations.

What would settle it

A direct algebraic check in which the exponential parametrization is substituted into the energy functional and the resulting stationarity condition fails to match the known AO Hartree-Fock equations would falsify the derivation.

read the original abstract

We present a pedagogical derivation of the Hartree--Fock equations using the second-quantization atomic-orbital density-matrix formalism developed by Kj{\ae}rgaard, J{\o}rgensen, Olsen, Coriani, and Helgaker for AO-based response theory. The purpose is to introduce an alternative derivation of the Hartree--Fock equation, showing that the standard AO Hartree--Fock stationarity condition follows naturally from the exponential parametrization of the one-particle density matrix in a nonorthogonal AO basis. This route provides a compact bridge between elementary Hartree--Fock theory and the density-matrix machinery used in modern response theory and linear-scaling formulations.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This is a clean pedagogical re-derivation of the known AO Hartree-Fock stationarity condition from the exponential density-matrix parametrization, with no new physics or results.

read the letter

Hi, the core point is that this paper walks through how the standard nonorthogonal AO Hartree-Fock equations drop out when you vary the energy with respect to the parameters in the exponential map of the one-particle density matrix. It uses the second-quantized AO density-matrix setup from the lead author's earlier papers with Jørgensen, Olsen, Coriani, and Helgaker. Nothing here is new in substance—the stationarity condition with the Fock matrix projected by the overlap is already standard—but the route is presented as a compact bridge to the machinery used in response theory and linear-scaling work. That connection is the main thing it does well: it shows the HF equations emerging naturally once the parametrization and constraints are in place, without extra machinery. For someone already comfortable with that formalism, this gives a tidy way to see the link back to elementary HF. The logic appears consistent on the high-level outline, with no internal contradictions or unsupported jumps flagged in the stress test. The soft spots are straightforward. The derivation rests completely on the prior AO density-matrix formalism, so the paper is not self-contained for readers who have not worked through those earlier references; that makes it less useful as a standalone teaching piece. It also supplies no numerical verification or new applications, which is fine for a pure derivation but keeps the scope narrow. No free parameters or invented entities are involved, and the citation pattern is honest about building on the group's own prior work. This is for quantum chemists who already use AO-based response methods and want to see the HF foundation laid out in that language. A reader looking for alternative derivations or clearer connections between basic theory and modern implementations could find it helpful. It is not essential reading for most, but the thinking is clear and the presentation is direct. I would send it to peer review. A referee can verify the explicit steps in the full manuscript and suggest minor clarifications on accessibility, but the work is solid enough on its own terms to merit that step rather than a desk reject.

Referee Report

0 major / 2 minor

Summary. The manuscript presents a pedagogical derivation of the Hartree-Fock equations in a nonorthogonal atomic-orbital basis. It employs the second-quantized AO density-matrix formalism of Kjærgaard et al. to define an exponential parametrization of the one-particle density matrix, then demonstrates that the standard AO Hartree-Fock stationarity condition (involving the Fock matrix projected onto the overlap) follows directly from requiring the energy to be stationary with respect to the parameters of that map, subject to the trace and idempotency constraints in the nonorthogonal basis.

Significance. If the steps are fully explicit and free of gaps, the work supplies a compact bridge between elementary Hartree-Fock theory and the density-matrix machinery already used in response theory and linear-scaling formulations. The derivation contains no free parameters or ad-hoc entities and re-uses an established prior formalism, which is a strength for pedagogical clarity and for showing how the conventional AO stationarity condition emerges naturally from the exponential map.

minor comments (2)

The abstract states that the stationarity condition 'follows naturally,' but the manuscript should explicitly flag the points at which the prior second-quantized AO formalism is invoked (e.g., the definition of the density-matrix parametrization and the trace/idempotency constraints) so that readers can trace the logic without consulting the cited references.
Notation for the overlap matrix S and the projected Fock matrix should be introduced with a brief reminder of their action on the nonorthogonal AO basis vectors, even if standard, to keep the pedagogical flow self-contained.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript and for recommending minor revision. No specific major comments were raised in the report, so we have no individual points to address point-by-point. We will ensure the revised version makes all derivation steps fully explicit with no gaps, consistent with the referee's significance remarks.

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The paper derives the known AO Hartree-Fock stationarity condition from the exponential parametrization of the one-particle density matrix using the cited second-quantized AO density-matrix formalism. This is an application of an external foundational formalism (developed for response theory) to obtain a standard result, without any self-definitional loop, fitted inputs renamed as predictions, or load-bearing self-citation that reduces the central claim to unverified inputs. The derivation chain connects parametrization, trace/idempotency constraints, and energy variation to the Fock matrix projection in a nonorthogonal basis, providing independent content as a pedagogical bridge. No steps exhibit the enumerated circular patterns.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The derivation starts from the second-quantized AO density-matrix formalism published earlier by the same group; no new free parameters or invented entities are introduced.

axioms (1)

domain assumption The second-quantized atomic-orbital density-matrix formalism is a valid representation of the one-particle density matrix in a nonorthogonal basis
Invoked as the starting point for the exponential parametrization

pith-pipeline@v0.9.0 · 5410 in / 1085 out tokens · 30635 ms · 2026-05-12T03:27:20.721828+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Reference graph

Works this paper leans on

12 extracted references · 12 canonical work pages

[1]

The density matrixD(equivalently, ∆) char- acterizes the Hartree–Fock state

work page
[2]

The idempotency conditionDSD=Dcon- strains the state to be a single Slater determi- nant

work page
[3]

The exponential parametrization ˜D = eiκSDe−iSκgenerates all allowed variations while preserving idempotency and particle number

work page
[4]

Stationarity of the energy with respect to κ yields the Hartree–Fock equationFDS− SDF= 0

work page
[5]

Hartree–Fock and Kohn–Sham time- dependent response theory in a second-quantization atomic-orbital formalism suitable for linear scaling,

T. Kjærgaard, P. Jørgensen, J. Olsen, S. Coriani, and T. Helgaker, “Hartree–Fock and Kohn–Sham time- dependent response theory in a second-quantization atomic-orbital formalism suitable for linear scaling,” J. Chem. Phys.129, 054106 (2008)

work page 2008
[6]

Hartree–Fock and Kohn–Sham atomic- orbital based time-dependent response theory,

H. Larsen, P. Jørgensen, J. Olsen, and T. Hel- gaker, “Hartree–Fock and Kohn–Sham atomic- orbital based time-dependent response theory,” J. Chem. Phys.113, 8908 (2000)

work page 2000
[7]

New developments in molecular orbital theory,

C. C. J. Roothaan, “New developments in molecular orbital theory,” Rev. Mod. Phys.23, 69–89 (1951)

work page 1951
[8]

The molecular orbital theory of chemical valency. VIII. A method of calculating ionization potentials,

G. G. Hall, “The molecular orbital theory of chemical valency. VIII. A method of calculating ionization potentials,” Proc. R. Soc. Lond. A205, 541–552 (1951)

work page 1951
[9]

Löwdin, “Quantum theory of many-particle sys- tems

P.-O. Löwdin, “Quantum theory of many-particle sys- tems. I. Physical interpretations by means of density matrices, natural spin-orbitals, and convergence prob- lems in the method of configurational interaction,” Phys. Rev.97, 1474–1489 (1955)

work page 1955
[10]

Some recent advances in density ma- trix theory,

R. McWeeny, “Some recent advances in density ma- trix theory,” Rev. Mod. Phys.32, 335–369 (1960)

work page 1960
[11]

Helgaker, P

T. Helgaker, P. Jørgensen, and J. Olsen,Molecu- lar Electronic-Structure Theory(Wiley, Chichester, 2000)

work page 2000
[12]

Szabo and N

A. Szabo and N. S. Ostlund,Modern Quantum Chem- istry: Introduction to Advanced Electronic Structure Theory(Dover, Mineola, 1996)

work page 1996

[1] [1]

The density matrixD(equivalently, ∆) char- acterizes the Hartree–Fock state

work page

[2] [2]

The idempotency conditionDSD=Dcon- strains the state to be a single Slater determi- nant

work page

[3] [3]

The exponential parametrization ˜D = eiκSDe−iSκgenerates all allowed variations while preserving idempotency and particle number

work page

[4] [4]

Stationarity of the energy with respect to κ yields the Hartree–Fock equationFDS− SDF= 0

work page

[5] [5]

Hartree–Fock and Kohn–Sham time- dependent response theory in a second-quantization atomic-orbital formalism suitable for linear scaling,

T. Kjærgaard, P. Jørgensen, J. Olsen, S. Coriani, and T. Helgaker, “Hartree–Fock and Kohn–Sham time- dependent response theory in a second-quantization atomic-orbital formalism suitable for linear scaling,” J. Chem. Phys.129, 054106 (2008)

work page 2008

[6] [6]

Hartree–Fock and Kohn–Sham atomic- orbital based time-dependent response theory,

H. Larsen, P. Jørgensen, J. Olsen, and T. Hel- gaker, “Hartree–Fock and Kohn–Sham atomic- orbital based time-dependent response theory,” J. Chem. Phys.113, 8908 (2000)

work page 2000

[7] [7]

New developments in molecular orbital theory,

C. C. J. Roothaan, “New developments in molecular orbital theory,” Rev. Mod. Phys.23, 69–89 (1951)

work page 1951

[8] [8]

The molecular orbital theory of chemical valency. VIII. A method of calculating ionization potentials,

G. G. Hall, “The molecular orbital theory of chemical valency. VIII. A method of calculating ionization potentials,” Proc. R. Soc. Lond. A205, 541–552 (1951)

work page 1951

[9] [9]

Löwdin, “Quantum theory of many-particle sys- tems

P.-O. Löwdin, “Quantum theory of many-particle sys- tems. I. Physical interpretations by means of density matrices, natural spin-orbitals, and convergence prob- lems in the method of configurational interaction,” Phys. Rev.97, 1474–1489 (1955)

work page 1955

[10] [10]

Some recent advances in density ma- trix theory,

R. McWeeny, “Some recent advances in density ma- trix theory,” Rev. Mod. Phys.32, 335–369 (1960)

work page 1960

[11] [11]

Helgaker, P

T. Helgaker, P. Jørgensen, and J. Olsen,Molecu- lar Electronic-Structure Theory(Wiley, Chichester, 2000)

work page 2000

[12] [12]

Szabo and N

A. Szabo and N. S. Ostlund,Modern Quantum Chem- istry: Introduction to Advanced Electronic Structure Theory(Dover, Mineola, 1996)

work page 1996