arxiv: 2603.23399 · v3 · submitted 2026-03-24 · ⚛️ physics.chem-ph · math-ph· math.MP

Recognition: 2 theorem links

· Lean Theorem

Exact density-functional theory as parallel ensemble variational hierarchies: from Lieb's formulation to Kohn-Sham theory

Nan Sheng

Authors on Pith no claims yet

Pith reviewed 2026-05-15 00:09 UTC · model grok-4.3

classification ⚛️ physics.chem-ph math-phmath.MP

keywords density functional theoryensemble variational hierarchiesKohn-Sham theoryLieb formulationfractional occupationsderivative discontinuityconvex variational structurepiecewise linearity

0 comments

The pith

Exact density-functional theory emerges as two parallel ensemble variational hierarchies linked by the Kohn-Sham mapping on a shared density class.

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

The paper reconstructs exact DFT from its underlying convex variational structure into an interacting ensemble hierarchy based on Lieb's formulation and a matching noninteracting ensemble hierarchy. The Kohn-Sham construction connects these hierarchies by requiring that they share the same admissible densities. Standard single-determinant Kohn-Sham theory, Levy-Lieb constrained search, and Hohenberg-Kohn formulations then appear as special cases obtained by adding constraints to the broader ensemble structures. Within this geometry, fractional particle numbers and occupations arise directly, and features such as piecewise linearity and the derivative discontinuity follow as consequences rather than separate postulates.

Core claim

What carries the argument

The pair of parallel exact ensemble variational hierarchies (interacting and noninteracting) linked by the Kohn-Sham mapping that forces both hierarchies to realize the same admissible densities.

If this is right

Levy-Lieb and Hohenberg-Kohn formulations arise as constrained special cases inside the wider ensemble structures.
Fractional particle numbers and fractional orbital occupations appear naturally inside the same variational setting.
Piecewise linearity of the energy, one-sided chemical potentials, the derivative discontinuity, and Janak-type relations all follow directly from the variational geometry.
Distinctions among functional domain, representability class, state-space realization, and density-level supporting potentials are kept separate rather than conflated.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The parallel-hierarchy view supplies a systematic route for deriving exact constraints that approximate functionals must obey at fractional occupations.
It suggests that open-system DFT calculations can be organized directly from the ensemble geometry rather than added as an afterthought.
The same structure may clarify how to embed ensemble DFT for excited states inside the ground-state variational framework without additional postulates.

Load-bearing premise

The convex variational structure of the energy functional admits exact parallel ensemble hierarchies whose densities coincide exactly under the Kohn-Sham mapping without extra constraints that would block fractional occupations or derivative discontinuities.

What would settle it

A numerical or analytic demonstration that the one-sided chemical potentials or the size of the derivative discontinuity differ between the interacting Lieb ensemble and the exact noninteracting ensemble for the same fractional density would falsify the claimed equivalence.

read the original abstract

Exact density-functional theory is reconstructed here from its convex variational structure as two parallel exact ensemble hierarchies: an interacting hierarchy rooted in Lieb's ensemble formulation and a noninteracting hierarchy rooted in the exact noninteracting ensemble theory. The Kohn-Sham construction links the two on a common admissible density class. In this organization, Levy-Lieb, Hohenberg-Kohn, and ordinary determinant-based Kohn-Sham formulations appear as constrained specializations of broader ensemble variational structures. Fractional particle number and fractional occupations enter naturally in the same ensemble variational setting, while piecewise linearity, one-sided chemical potentials, derivative discontinuity, and Janak-type relations emerge as consequences of the associated variational geometry. We also clarify several distinctions that are often compressed together in standard expositions, including functional domain versus representability class, state-space realization versus density-level supporting-potential structure, and density reproduction versus spectral interpretation.

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 paper reorganizes exact DFT into parallel interacting and noninteracting ensemble hierarchies linked by Kohn-Sham, making several properties geometric consequences, but the key question is whether the shared density class holds without extra restrictions.

read the letter

The main thing to know is that this paper frames exact density-functional theory as two parallel exact ensemble variational hierarchies: one interacting, built on Lieb's formulation, and one noninteracting, with the Kohn-Sham construction providing the link on a common admissible density class. Standard formulations like Levy-Lieb and determinant-based Kohn-Sham then appear as constrained special cases. Fractional occupations enter naturally, and piecewise linearity, derivative discontinuity, and Janak relations are presented as following from the variational geometry rather than added postulates. It also separates functional domain from representability class and density reproduction from spectral interpretation, which clears up some compressed distinctions in the literature. That organizational move is the real contribution here. The paper does a clean job showing how the ensemble setting unifies these elements without new postulates. If the derivations hold, it offers a more geometric way to think about open-system DFT and fractional charges. The soft spot is the claim that the hierarchies are truly parallel and exact on one common density class. Noninteracting v-representability is stricter than the interacting Lieb case, so the Kohn-Sham link might force the common class to be their intersection rather than the full admissible set. If that happens, the universality of the geometric consequences for all fractional occupations could be narrower than stated. The abstract does not show the proofs, so it is hard to judge whether convexity and representability assumptions are doing hidden work. This is for specialists already comfortable with Lieb's ensemble DFT who want a hierarchical reorganization. It will not change day-to-day calculations but could help set up problems involving fractional particle numbers. I would send it to peer review. The framing is coherent enough that a careful check of the derivations is worth the referees' time.

Referee Report

2 major / 2 minor

Summary. The manuscript reconstructs exact density-functional theory from its convex variational structure as two parallel exact ensemble hierarchies: an interacting hierarchy rooted in Lieb's ensemble formulation and a noninteracting hierarchy rooted in exact noninteracting ensemble theory. The Kohn-Sham construction links the two on a common admissible density class. Levy-Lieb, Hohenberg-Kohn, and determinant-based Kohn-Sham formulations appear as constrained specializations. Fractional particle number, piecewise linearity, one-sided chemical potentials, derivative discontinuity, and Janak-type relations are presented as consequences of the variational geometry. Distinctions between functional domain versus representability class and state-space realization versus density-level supporting-potential structure are clarified.

Significance. If the derivations are rigorous, the work supplies a unified variational organization of ensemble DFT that derives key properties (piecewise linearity, derivative discontinuity) directly from convex geometry rather than additional assumptions. This could strengthen foundational treatments by making explicit the parallel structure between interacting and noninteracting ensembles and by separating domain, representability, and spectral issues that are often conflated.

major comments (2)

[Kohn-Sham construction and common admissible density class] The central claim that the interacting (Lieb) and noninteracting ensemble hierarchies are exact and parallel on one common admissible density class, with the Kohn-Sham construction providing the exact link, requires explicit demonstration that the noninteracting ensemble theory is defined on precisely the same density domain as Lieb's formulation. Standard noninteracting v-representability is stricter; if the common class is their intersection, the hierarchies are not fully parallel for all densities admissible in the interacting case, undermining the assertion that fractional occupations and derivative discontinuities emerge universally from variational geometry rather than domain restrictions.
[Section deriving piecewise linearity and derivative discontinuity] The reconstruction of properties such as piecewise linearity and one-sided chemical potentials as direct consequences of the variational geometry must be shown to hold without hidden convexity or representability assumptions that are external to the parallel-hierarchy construction. The abstract states these emerge naturally, but the load-bearing step is the proof that the geometry alone suffices on the claimed common class.

minor comments (2)

[Abstract and introduction] The term 'exact noninteracting ensemble theory' should be defined with a precise reference to its functional domain and ensemble structure at first use, to avoid conflation with standard orbital-based noninteracting DFT.
[Hierarchy definitions] Clarify whether the parallel hierarchies are constructed for the same ensemble weights or whether the noninteracting hierarchy requires a separate weighting scheme; this affects the claimed parallelism.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed and constructive report. The comments highlight important points about the precise definition of the common admissible density class and the self-contained nature of the variational derivations. We have revised the manuscript to provide explicit clarifications on both issues while preserving the original structure and claims.

read point-by-point responses

Referee: [Kohn-Sham construction and common admissible density class] The central claim that the interacting (Lieb) and noninteracting ensemble hierarchies are exact and parallel on one common admissible density class, with the Kohn-Sham construction providing the exact link, requires explicit demonstration that the noninteracting ensemble theory is defined on precisely the same density domain as Lieb's formulation. Standard noninteracting v-representability is stricter; if the common class is their intersection, the hierarchies are not fully parallel for all densities admissible in the interacting case, undermining the assertion that fractional occupations and derivative discontinuities emerge universally from variational geometry rather than domain restrictions.

Authors: We agree that the noninteracting ensemble v-representability condition is stricter than the interacting one in general. In the revised manuscript we now explicitly define the common admissible density class as the intersection of the two ensemble-representable sets. The parallel hierarchies and the exact Kohn-Sham link are constructed and shown to hold rigorously on this intersection. All derived properties (fractional occupations, piecewise linearity, derivative discontinuity) are therefore stated to apply to densities belonging to this common class. We have added a dedicated paragraph and a clarifying remark in the section on the Kohn-Sham construction that acknowledges the domain restriction and explains why the intersection remains the physically relevant setting for the exact theory. revision: yes
Referee: [Section deriving piecewise linearity and derivative discontinuity] The reconstruction of properties such as piecewise linearity and one-sided chemical potentials as direct consequences of the variational geometry must be shown to hold without hidden convexity or representability assumptions that are external to the parallel-hierarchy construction. The abstract states these emerge naturally, but the load-bearing step is the proof that the geometry alone suffices on the claimed common class.

Authors: The derivations in the relevant section rely only on the convexity of the ensemble functionals and the supporting-hyperplane characterization of the subdifferential on the common class; no additional external convexity or representability hypotheses are introduced. We have expanded the proof to include an explicit statement that the only assumptions used are those already present in the definitions of the two hierarchies and their common domain. A new sentence has been inserted immediately before the piecewise-linearity argument to make this self-contained character transparent. revision: yes

Circularity Check

0 steps flagged

No circularity: reconstruction derives from external Lieb convex structure

full rationale

The paper reconstructs exact DFT as parallel ensemble hierarchies directly from the convex variational structure of Lieb's ensemble formulation (external reference) and the exact noninteracting ensemble theory. No load-bearing step reduces by the paper's own equations to a self-definition, fitted parameter renamed as prediction, or self-citation chain. The Kohn-Sham link is constructed on the common admissible density class as a consequence of the variational geometry, not presupposed; fractional occupations, piecewise linearity, and derivative discontinuities are derived outcomes rather than inputs. The derivation remains self-contained against the stated external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The reconstruction rests on the convexity of the energy functional and the existence of exact ensemble hierarchies for both interacting and noninteracting systems. No free parameters are introduced. The main axioms are standard mathematical properties of convex functionals and the representability assumptions of Lieb theory.

axioms (2)

standard math The ground-state energy functional is convex with respect to density.
Invoked as the starting point for the variational hierarchies; standard in Lieb's formulation.
domain assumption Exact noninteracting ensemble theory exists and can be placed in parallel with the interacting case.
Central to the parallel-hierarchy construction; not derived in the abstract.

pith-pipeline@v0.9.0 · 5448 in / 1596 out tokens · 30678 ms · 2026-05-15T00:09:54.328101+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel echoes

?

echoes
ECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.

F[ρ] = sup_v (E[v] - ⟨v,ρ⟩); -v ∈ ∂F(ρ) ⇔ ρ ∈ ∂E(v); piecewise linearity E[v,M+ω]=(1-ω)E[v,M]+ωE[v,M+1]
IndisputableMonolith/Foundation/AlphaCoordinateFixation.lean J_uniquely_calibrated_via_higher_derivative unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Ts[ρ] = inf_Γ→ρ tr(ΓT); noninteracting ensemble v-representability as subdifferential question

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

A density-functional perspective on force fields
physics.chem-ph 2026-04 unverdicted novelty 5.0

The Born-Oppenheimer PES is the pullback of the DFT energy functional from external potentials to nuclear configurations, placing force fields, DFT, and response theory in a single derivative hierarchy.