Pith. sign in

REVIEW 4 major objections 5 minor 20 references

Stability of bound states in multi-component DFT in absolute coordinate systems

T0 review · 4 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read The paper claims that a fully translation-invariant electron–nuclear Hamiltonian can support localized bound states, with stable windows at very high density (rs below 1.4) and very low density (rs above 10), while the intermediate regime…

desk verdict A genuinely new question undermined by an ad hoc rescaling of exchange that controls the claimed stability windows. read the letter →

arxiv 2506.07990 v1 pith:QEZOKI5S submitted 2025-06-09 cond-mat.mtrl-sci

Bander Linjawi This is my paper
classification cond-mat.mtrl-sci PACS 71.15.Mb31.15.E
keywords multi-componentdensityfunctionaltheoryelectron-nuclearcorrelationboundstateshomogeneouselectrongasGalileancoordinatesWigner-Seitzradiusexchangescalingsizeconsistency
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

Multi-component density functional theory usually breaks translational symmetry by treating nuclei as fixed external potentials. This paper asks whether that crutch is necessary, by testing a Hamiltonian written in absolute Galilean coordinates that has no external potential at all. The author compares the energy of two homogeneous Fermi gases (electrons and protons at equal density) with the energy of a localized trial state built from Gaussian Slater profiles, minimizing over the decay lengths with a BFGS optimizer. The central claim is that the localized state wins at high density (rs below 1.4) and at low density (rs above 10), with a middle window where the gas wins and electron–nuclear correlation is predicted to be decisive. If true, this supports constructing multi-component DFT functionals in Galilean coordinates without pinning nuclei, and it makes the intermediate density regime the place to look for correlation-driven bound-state physics.

What carries the argument

The load-bearing object is the energy-difference functional ΔE(α,β,rs,N) built from a localized trial state |Ψ⟩=|Ψ⟩FS|Ψ⟩S. Here |Ψ⟩FS is the product of electron and proton Fermi seas and |Ψ⟩S is a Slater determinant of Gaussians, giving localized density fluctuations around the homogeneous densities; α and β are the inverse decay lengths of the electronic and nuclear fluctuations. The machinery works because the Hartree terms cancel exactly for homogeneous densities, the Coulomb integrals of the Gaussians are known in closed form (for instance 20π²/α⁵ for like-charge overlap and 32π²(α²+3αβ+β²)/(α²β²(α+β)³) for cross terms), and exchange is treated in LDA form. The decisive step is a hand-set scaling of the localized-state exchange from its natural $N^{{4/3}}$ behavior to N-linear, which the author justifies by matching the homogeneous state's scaling and by references to N-dependent normalizations in density functional theory. Correlation is appended through an electron–proton functional in the center-of-mass coordinate s=(r+R)/2.

What would settle it

Recompute ΔE(α,β,rs,N) with the localized-state exchange kept at its orbital-resolved $N^{{4/3}}$ scaling, or with an exact exchange calculation, and check whether the sign of the energy difference still favors the localized state for rs<1.4 and rs>10. If the windows disappear, they are artifacts of the scaling choice rather than a property of the translationally invariant Hamiltonian.

Watch

Extended reading notes

Core claim

On the paper's own terms, the discovery is that a translationally invariant electron–nuclear Hamiltonian admits bound states relative to the homogeneous gas over two density windows. Starting from a product of electron and nuclear Fermi seas, the localized state is a Slater determinant of paired Gaussians with decay lengths α and β, producing density fluctuations δn(r)=n0(Nα³/(8π)$e^{{-n0^{1/3}}$αr}-1) and the analogous nuclear form. The energy difference ΔE(α,β,rs,N) contains kinetic, Hartree, exchange, and correlation pieces; the Hartree contribution is quadratic in N, the kinetic and exchange contributions are made linear in N, and the relation n0=3/(4πrs³) converts densities into the Wigner-Seitz radius rs. Minimizing over α and β at N=1000 yields stability for rs<1.4 and rs>10, and the optimum is rs-independent with equal electronic and nuclear profiles (α=β), at a localized energy of -0.0439 Ha per particle. The paper also reports that adding an electron–proton Colle-Salvetti correlation functional favors the homogeneous gas, leaving the intermediate regime as the key open problem.

Load-bearing premise

The result hinges on a hand-set rescaling of the localized state's exchange energy from its natural $N^{{4/3}}$ growth to linear-in-N growth so that it matches the homogeneous state; this rescaling is introduced for consistency, not derived, and the stability windows would move or vanish if the correct scaling is different.

Editorial extensions

If this is right

  • If the stability windows are correct, a translationally invariant multi-component Hamiltonian can have symmetry-broken localized ground states without any external potential, so DFT treatments of electron–nuclear systems need not pin nuclei to a fixed frame.
  • The high-density window (rs<1.4) means bound-state formation is possible even where the homogeneous gas kinetic energy is large, which is the opposite of the usual intuition that high densities favor delocalization.
  • The low-density window (rs>10) suggests that in dilute electron–proton systems, for example in warm dense matter regimes, localized electron–nuclear clusters can lower the energy below the gas.
  • The intermediate regime, roughly 1.4<rs<10, is predicted to prefer the homogeneous gas at the product-state level, with electron–nuclear correlation named as the factor most likely to change that ordering.
  • The framework yields explicit closed-form energy differences that can be reused as a starting point for more realistic trial states or for perturbing around the free-gas minimum.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Beyond the paper, the stability windows are only as strong as the ad hoc N-linear exchange rescaling; leaving the localized exchange at its orbital-resolved N^{4/3} scaling could move or close the windows, so the windows should be re-checked with an exact-exchange or orbital-resolved calculation.
  • The observed rs-independence of the optimal profiles with α=β suggests a hidden symmetry of the Gaussian product ansatz rather than a physical mechanism; breaking the shape symmetry between electrons and nuclei may open or shift the stable regions.
  • The hydrogen-atom correlation estimate (-0.0892 Ha total, -0.0446 Ha per particle) is comparable to the energy gap at the free-gas minimum, implying that a correlated trial state could stabilize the intermediate density regime; this is a testable extension within the same formalism.
  • A direct check in uniform-density quantum plasma simulations could look for spontaneous localization of protons and electrons at mean densities corresponding to rs<1.4 or rs>10, providing an independent test not available in the variational model.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. The manuscript asks whether a translationally invariant electron-nuclear Hamiltonian in absolute (Galilean) coordinates can support localized bound states relative to homogeneous electron and nuclear gases. A homogeneous charge-neutral system is perturbed by Gaussian-like localized density fluctuations, and the energy difference is computed from kinetic, Hartree, LDA exchange, and optionally a Colle-Salvetti-type electron-proton correlation functional, with decay exponents alpha and beta optimized by BFGS. The paper reports stability windows below r_s=1.4 and above r_s=10, an unstable intermediate regime around the free-gas optimum r_s=2.412, and discusses the merits of Galilean coordinates for constructing multi-component density functionals.

Significance. If the main result were correct, this would be a useful minimal demonstration that a translationally invariant multi-component Hamiltonian can have symmetry-broken localized states, and it would support the practical use of absolute-coordinate frame DFT. The paper is transparent in stating several limitations, including the heuristic connection between the two states and the possible overestimation of gas-state correlation. However, the central quantitative claim is not established: it relies on an explicitly ad hoc rescaling of the localized exchange energy from N^(4/3) to N in Eq. (15), and the density fluctuations in Eq. (4) do not integrate to zero in infinite space without an unstated finite-volume convention. In addition, the trial wavefunction in Eq. (3) appears inconsistent with the densities and kinetic energies used in the energy evaluation. These issues affect the sign and location of the reported stability windows, so the paper as written cannot support its headline conclusion.

major comments (4)
  1. [II, Eqs. (10)-(15)] The localized-state exchange energy is derived in Eq. (10) as V_X[n] = -(3/4)^4 (3/(8 pi^2))^{1/3} n_0^{1/3} alpha N^{4/3}, but Eq. (15) replaces it by an N-linear form, with the text stating that the scaling is fixed 'to give the same scaling as the homogeneous state.' No derivation or physical argument is supplied for this replacement. Exchange from a Slater determinant is a two-body interaction and should scale as N(N-1)/2 at fixed orbital shape or as N^{4/3} in LDA on a localized density; a strictly linear N scaling is the signature of a one-body term. At N=1000 the adopted rescaling changes the localized exchange by a factor N^{1/3} ~ 10, and since the stability windows are decided by the balance between Delta T ~ -1.105/r_s^2 and Delta V_X ~ 0.916/r_s in Eq. (17), the reported boundaries r_s<1.4 and r_s>10 are not robust. Please derive the correct scaling for the actual trial wavefunction or remove the quantitative stability claim.
  2. [II, Eq. (3)] The trial state in Eq. (3) is written as a product over all electron-nucleus pairs of exp[-(n_0^{1/3} alpha r_i + n_0^{1/3} beta R_I)^2]. This is not a product of independent electron and nuclear orbitals; for a single electron-nucleus pair, integrating over R leaves a uniform electron density, which contradicts the localized density in Eq. (4). The normalization prefactor and the kinetic energy in Eq. (12) appear to assume independent Gaussian orbitals, so the energy evaluation is not consistent with the stated wavefunction. Please specify how the densities in Eq. (4) are obtained from Eq. (3) and derive the kinetic energy from the actual wavefunction.
  3. [II, Eq. (4)] The density fluctuations delta n and delta m are stated to integrate to zero, but over infinite space the term -n_0 integral d^3 r diverges. The equality holds only in a finite normalization volume V = N/n_0, which is not specified in the manuscript. Because Eq. (7) subtracts the homogeneous-state energy and relies on cancellation of Hartree terms between the uniform backgrounds, the missing volume convention affects the energy difference and must be stated. As written, the derivation is not valid in infinite space.
  4. [III, Eq. (23) and Conclusion] At the free-gas optimum r_s = 2.412, the included electron-proton correlation functional gives -0.012 Ha for the localized state and -0.0462 Ha for the gas, i.e., correlation strongly favors the gas at the density where the uncorrelated model is closest to stability. The conclusion then states that the two states are 'only heuristically connected.' The paper does not compute the correlation contribution across the full r_s range, so the claimed stability windows are properties of the exchange-rescaled, correlation-neglected model rather than of the Hamiltonian. The abstract's 'Regions of stability' should be qualified accordingly.
minor comments (5)
  1. [I] The word 'nucelar' in the introduction is a typo for 'nuclear.'
  2. [II] The text refers to the 'Wingerseitz radius'; this should be the 'Wigner-Seitz radius.'
  3. [III] The sentence introducing Eq. (23) contains the typo 'a, b, care taken' and should read 'a, b, c are taken.'
  4. [IV] The stability-window statement in the abstract should carry the same caveat as the conclusion, namely 'in the regime investigated,' since the result depends on the chosen trial ansatz and on N=1000.
  5. [Figure 1] Figure 1 is referenced in the text but no figure or caption appears in the manuscript; please include the actual figure with labeled axes and a clear caption.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the variational calculation is self-contained, though the ad hoc exchange rescaling is a robustness concern, not a circular step.

full rationale

The paper's central claim—stability windows below rs=1.4 and above rs=10—is obtained by minimizing the energy difference in Eq. (17) with respect to the variational decay lengths alpha and beta. That optimization is a genuine energy minimization, not a fit to the reported stability boundaries. The Hamiltonian, the Slater-product ansatz, and the explicit evaluations of Hartree, exchange, and kinetic terms in Eqs. (8)-(12) are all stated in the paper, and the subsequent minimization is arithmetic from those expressions. The one passage that could raise concern is Section II, after Eq. (13), where the localized exchange is rescaled: 'In this study, we fix the scaling to give the same scaling as the homogeneous state.' This changes the N-dependence of the localized exchange from N^(4/3) to N, and the location of the stability windows is sensitive to that choice. However, this is an openly stated modeling assumption rather than a circular reduction: the windows are not defined as the assumption itself, nor are the reported rs values inserted into the model. The rescaling is not fitted to the target result, and it is not justified by a self-citation; the cited works on N-dependent normalizations are external. Similarly, the electron-proton correlation functional used in Section III is imported from Refs. [18,19], which are independent of this author, and the conclusion drawn from it (the gas state remains more stable at rs=2.412) is not used to construct the stability windows. The calculation is therefore self-contained in the sense required for circularity analysis. The arbitrary exchange rescaling is a legitimate scientific weakness that affects the robustness of the prediction, but it belongs under correctness or assumption-risk, not circularity.

Assumptions & free parameters 6 free parameters · 6 assumptions · 0 invented entities

The model introduces no new physical entities. Its free parameters are the two variational decay widths and the imported electron-proton correlation constants. The most important assumption is the ad hoc rescaling of the exchange scaling, which is load-bearing for the stability claim.

free parameters (6)
  • alpha (electronic decay exponent) = approx. 1.352 (dimensionless, at variational optimum)
    Variational parameter controlling the width of the localized electron density; minimized with BFGS in Eq. (17); stability windows shift with its value.
  • beta (nuclear decay exponent) = approx. 1.352 (dimensionless, at variational optimum)
    Variational parameter controlling the width of the localized nuclear density; at equal densities the optimum matches alpha.
  • N (particle number) = 1000
    Chosen to approach the thermodynamic limit; energy differences scale as N or N^2 and the stability windows depend on this choice.
  • a (epc correlation parameter) = 2.35
    Imported from Yang et al. 2017 where it was fitted to accurate proton densities; used in Eq. (23) for the correlation energy comparison.
  • b (epc correlation parameter) = 2.4
    Imported from Yang et al. 2017; appears in the denominator of Eq. (23).
  • c (epc correlation parameter) = 3.2
    Imported from Yang et al. 2017; appears in the denominator of Eq. (23).
assumptions (6)
  • domain assumption The full electron-nuclear Hamiltonian is translationally invariant and its ground state may be represented by a symmetry-broken product state.
    Invoked in Section II; the paper notes this is a constrained treatment and that interactions between the gases and the fluctuation are neglected.
  • domain assumption Electron and nuclear densities are equal (m0 = n0, Z = 1) to maintain charge neutrality.
    Set in Section II after Eq. (4); restricts the model to neutral equal-density systems.
  • ad hoc to paper Density fluctuations integrate to zero over the normalization volume.
    Required by Eq. (2) but inconsistent with the unbounded expressions in Eq. (4); a finite box is assumed without being stated.
  • ad hoc to paper Exchange can be approximated by LDA and the localized-state exchange scaling is reset from N^(4/3) to N.
    The rescaling is asserted in Section II and is not derived; it changes the energy balance.
  • domain assumption Correlation energy is evaluated with the epc functional using parameters from Ref. [18].
    Used in Section III; those parameters were fitted in prior work.
  • standard math BFGS minimization over alpha and beta yields the variational optimum for each rs.
    Numerical optimization described in Section II; no guarantee of a global minimum is stated.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Stability of bound states in multi-component DFT in absolute coordinate systems." pith.science (2026). https://pith.science/paper/QEZOKI5S

@misc{pith2026250607990,
  author       = {Pith},
  title        = {Pith review of: Stability of bound states in multi-component DFT in absolute coordinate systems},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/QEZOKI5S}},
  note         = {Machine review of arXiv:2506.07990}
}
read the original abstract

Homogeneous electron and nuclear gases are transformed to a localized trial density in absolute coordinates of the multi-component hamiltonian to determine the stability of forming bound states. Regions of stability were found both at the high density and low density regimes, where electron-nuclear correlations could play a critical role in the intermediate density regime. The use of Galilean coordinates is motivated for its use in density functional theory to develop kinetic and potential density functionals, from which suitable coordinate transformations to capture electron-nuclear correlations are applied.

Figures

Figures reproduced from arXiv: 2506.07990 by the authors.

Figure 1
Figure 1. FIG. 1. Energies of the bound state and the homogeneous [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

20 extracted references · 20 canonical work pages

  1. [1]

    Non-born–oppenheimer density functional theory of molecular systems

    Joseph F Capitani, Roman F Nalewajski, and Robert G Parr. Non-born–oppenheimer density functional theory of molecular systems. The Journal of Chemical Physics, 76(1):568–573, 1982

    work page 1982

  2. [2]

    The 2021 room-temperature superconductivity roadmap

    Boeri Lilia, Richard Hennig, Peter Hirschfeld, Gianni Profeta, Antonio Sanna, Eva Zurek, Warren E Pickett, Maximilian Amsler, Ranga Dias, Mikhail I Eremets, et al. The 2021 room-temperature superconductivity roadmap. Journal of Physics: Condensed Matter, 34(18):183002, 2022

    work page 2021

  3. [3]

    Unraveling electronic correlations in warm dense quantum plasmas

    Tobias Dornheim, Tilo D¨ oppner, Panagiotis Tolias, Max- imilian P B¨ ohme, Luke B Fletcher, Thomas Gawne, Frank R Graziani, Dominik Kraus, Michael J MacDon- ald, Zhandos A Moldabekov, et al. Unraveling electronic correlations in warm dense quantum plasmas. Nature Communications, 16(1):1–11, 2025

    work page 2025

  4. [4]

    Electron- proton correlation for hydrogen tunneling systems

    Michael V Pak and Sharon Hammes-Schiffer. Electron- proton correlation for hydrogen tunneling systems. Phys- ical review letters, 92(10):103002, 2004

    work page 2004

  5. [5]

    Development of electron-proton density functionals¡? format?¿ for multicomponent density functional theory

    Arindam Chakraborty, Michael V Pak, and Sharon Hammes-Schiffer. Development of electron-proton density functionals¡? format?¿ for multicomponent density functional theory. Physical review letters , 101(15):153001, 2008

    work page 2008

  6. [6]

    Multicomponent density-functional theory for electrons and nuclei

    Thomas Kreibich, Robert van Leeuwen, and EKU Gross. Multicomponent density-functional theory for electrons and nuclei. Physical Review A, 78(2):022501, 2008

    work page 2008

  7. [7]

    Multicomponent density- functional theory for electrons and nuclei

    T Kreibich and EKU Gross. Multicomponent density- functional theory for electrons and nuclei. Physical Re- view Letters, 86(14):2984, 2001

    work page 2001

  8. [8]

    Nuclear–electronic orbital approach to quantiza- tion of protons in periodic electronic structure calcula- tions

    Jianhang Xu, Ruiyi Zhou, Zhen Tao, Christopher Mal- bon, Volker Blum, Sharon Hammes-Schiffer, and Yosuke Kanai. Nuclear–electronic orbital approach to quantiza- tion of protons in periodic electronic structure calcula- tions. The Journal of Chemical Physics, 156(22), 2022

    work page 2022

Show all 20 references
  1. [9]

    Multicomponent density-functional the- ory for time-dependent systems

    O Butriy, H Ebadi, PL de Boeij, R van Leeuwen, and EKU Gross. Multicomponent density-functional the- ory for time-dependent systems. Physical Review A, 76(5):052514, 2007

    2007

  2. [10]

    First-principles approach to the electron-phonon interaction

    Robert van Leeuwen. First-principles approach to the electron-phonon interaction. Physical Review B, 69(11):115110, 2004

    2004

  3. [11]

    P. Ziesche. Pair–density functional theory. Phys. Lett. A, 195:213–220, 1994

    1994

  4. [12]

    The pair-density functional of the kinetic energy and its simple scaling property

    Mel Levy and Paul Ziesche. The pair-density functional of the kinetic energy and its simple scaling property. J. Chem. Phys., 115:9110–9112, 2001

    2001

  5. [13]

    E[n0+δn, m0+δm] = ⟨HF S⊗IS+IF S⊗HS⟩−⟨HF ⊗IS⟩ = ⟨IF S⊗ HS⟩

  6. [14]

    Introduction to quantum mechanics

    David J Griffiths and Darrell F Schroeter. Introduction to quantum mechanics. Cambridge university press, 2018

    2018

  7. [15]

    Approxi- mate normalizations for approximate density functionals

    Adam Clay, Kiril Datchev, Wenlan Miao, Adam Wasser- man, Kimberly J Daas, and Kieron Burke. Approxi- mate normalizations for approximate density functionals. arXiv preprint arXiv:2504.03845, 2025

    2025

  8. [16]

    Relevance of coordinate and particle-number scaling in density- functional theory

    Eduardo Fabiano and Lucian A Constantin. Relevance of coordinate and particle-number scaling in density- functional theory. Physical Review A—Atomic, Molec- ular, and Optical Physics, 87(1):012511, 2013

    2013

  9. [17]

    Approximate calcula- tion of the correlation energy for the closed shells

    Renato Colle and Oriano Salvetti. Approximate calcula- tion of the correlation energy for the closed shells. The- oretica chimica acta, 37:329–334, 1975

    1975

  10. [18]

    Development of a practical multicomponent density functional for electron- proton correlation to produce accurate proton densities

    Yang Yang, Kurt R Brorsen, Tanner Culpitt, Michael V Pak, and Sharon Hammes-Schiffer. Development of a practical multicomponent density functional for electron- proton correlation to produce accurate proton densities. The Journal of Chemical Physics, 147(11):114113, 2017

    2017

  11. [19]

    [18] as 2.35, 2.4, and 3.2 respectively

    functional without gradient terms since they evaluate to 0 for Fermi Seas VC[n, m] = − Z d3s n(s)m(s) a − b n 1 2 (s)m 1 2 (s) + c n(s)m(s) (23) where a, b, care taken from Ref. [18] as 2.35, 2.4, and 3.2 respectively. The coordinate s is the center of mass coordinate s = R′ C...

  12. [20]

    Multicomponent density functional theory: Including the density gradient in the electron-proton correlation functional for hydrogen and deuterium

    Zhen Tao, Yang Yang, and Sharon Hammes-Schiffer. Multicomponent density functional theory: Including the density gradient in the electron-proton correlation functional for hydrogen and deuterium. The Journal of Chemical Physics, 151(12):124102, 2019

    2019

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.