Spin-polarized Energy Density Method from Spin-Density Functional Theory
Pith reviewed 2026-05-09 20:54 UTC · model grok-4.3
The pith
A spin-polarized version of the energy density method decomposes DFT total energies into atomic contributions using real-space densities and gauge-invariant volumes.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The energy density method is generalized to include spin polarization with the full formalism derived based on spin-density functional theory, which aims at decomposing the total energy into well-defined atomic energies. The method involves two steps: decomposing the total energy into spin-polarized energy density functions in real space, and integrating these energy densities over chosen gauge-invariant volumes for uniquely defined atomic energies, whose summation over all the atoms restores the DFT total energy up to a constant difference. This is numerically implemented for the projector augmented-wave method and applied to model paramagnetic fcc Fe and dilute magnetic semiconductor Ni:Ga
What carries the argument
Spin-polarized energy density functions in real space, integrated over gauge-invariant volumes to define atomic energies.
If this is right
- Atomic energies become definable in systems with spin polarization, such as paramagnetic iron modeled by special quasirandom structures.
- Spin energies from these atomic decompositions can be used to fit spin cluster expansions or train deep neural networks.
- The approach calculates atomic energy distributions for different dopant distances and spin configurations in Ni-doped GaN.
- The sum of all atomic energies recovers the total DFT energy up to a constant difference.
- This extracts additional information from standard DFT calculations for studying magnetic systems.
Where Pith is reading between the lines
- Local energy contributions from different spin alignments could be compared across configurations to identify stable magnetic structures.
- The method might apply to other magnetic phenomena like domain walls or defects by assigning energies to specific atoms.
- Testing on systems where total energy is known independently could validate the constant difference offset.
Load-bearing premise
The chosen gauge-invariant volumes produce atomic energies that are unique and whose sum matches the total energy up to a constant, without introducing new ambiguities from the spin polarization.
What would settle it
A numerical test in which the summed atomic energies differ from the DFT total energy by an amount larger than any reasonable constant offset, or where changing the volume definition changes the individual atomic energies inconsistently.
Figures
read the original abstract
The energy density method is generalized to include spin polarization with the full formalism derived based on spin-density functional theory, which aims at decomposing the total energy into well-defined atomic energies. The method involves two steps: (1) decomposing the total energy into spin-polarized energy density functions in real space, and (2) integrating these energy densities over chosen gauge-invariant volumes for uniquely defined atomic energies, whose summation over all the atoms restores the DFT total energy up to a constant difference. This method is numerically implemented into the Vienna ab initio simulation package for the projector augmented-wave method, and is showcased with two applications. In the first application, we model the paramagnetic face-centered cubic Fe using spin special quasirandom structures; the spin energies are fit to spin cluster expansions and a deep neural network. In the second application, we calculate the atomic energy distributions of dilute magnetic semiconductor Ni-doped GaN with different dopant distances and spin configurations. This method extracts additional useful information for the study of magnetic systems with density functional theory.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript generalizes the energy density method to spin-polarized systems within spin-density functional theory. It outlines a two-step procedure: (1) decomposing the total energy into spin-polarized energy density functions in real space, and (2) integrating these densities over chosen gauge-invariant volumes to define atomic energies. These atomic energies are asserted to sum to the DFT total energy up to a constant difference. The approach is implemented in VASP for the PAW method and applied to paramagnetic fcc Fe modeled via spin special quasirandom structures (with fits to spin cluster expansions and a deep neural network) and to dilute magnetic semiconductor Ni-doped GaN (for atomic energy distributions across dopant distances and spin configurations).
Significance. If the central claims on unique atomic energies and a configuration-independent constant offset hold under numerical validation, the method would enable local energy decomposition in magnetic materials. This could support spin cluster expansions for paramagnetic systems and analysis of energy distributions in doped magnetic semiconductors, providing additional interpretive tools beyond standard total-energy DFT.
major comments (3)
- [Abstract] Abstract and method description: the two-step formalism for spin-polarized energy densities and volume integration is stated without any equations, derivation steps, or explicit proof that the constant offset remains independent of spin configuration. This is load-bearing for the uniqueness claim, as spin density variations in SQS Fe or varying Ni:GaN dopant/spin setups could make gauge-invariant volume boundaries (Bader, Voronoi, etc.) configuration-dependent.
- [Applications] Applications section: no quantitative validation is reported (e.g., maximum deviation or variance of the summed atomic energies minus total DFT energy across the spin SQS ensemble or Ni:GaN configurations). Without this, it is impossible to confirm that the constant is invariant, which directly affects the reliability of the subsequent cluster-expansion fits and distribution analyses.
- [Implementation] Implementation paragraph: details are absent on how the spin-polarized energy density is numerically realized inside the PAW formalism in VASP, including any gauge choices, integration grids, or handling of spin-density inhomogeneities at volume boundaries.
minor comments (2)
- Define all acronyms (SQS, PAW, DFT) at first use and ensure consistent notation for energy densities versus atomic energies throughout.
- The abstract claims 'uniquely defined atomic energies'; clarify in the text whether multiple valid gauge-invariant volume definitions (e.g., Bader vs. Voronoi) produce identical atomic energies or only equivalent sums up to the constant.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for the constructive major comments. We have revised the manuscript to address each point, adding equations and clarification to the abstract, quantitative validation of the constant offset in the applications, and expanded implementation details.
read point-by-point responses
-
Referee: [Abstract] Abstract and method description: the two-step formalism for spin-polarized energy densities and volume integration is stated without any equations, derivation steps, or explicit proof that the constant offset remains independent of spin configuration. This is load-bearing for the uniqueness claim, as spin density variations in SQS Fe or varying Ni:GaN dopant/spin setups could make gauge-invariant volume boundaries (Bader, Voronoi, etc.) configuration-dependent.
Authors: We agree that the abstract would benefit from greater specificity. In the revised version we will insert the central equations for the spin-polarized energy-density decomposition and the subsequent volume integration, together with a concise statement that the additive constant follows directly from the fact that the integral of the energy density over all space recovers the DFT total energy for any gauge-invariant partitioning. The derivation in Section II demonstrates that this constant is independent of spin configuration because the partitioning volumes are constructed from the total charge density (not the spin density), rendering the boundaries insensitive to the paramagnetic spin fluctuations present in the SQS or dopant configurations. We will also note this point explicitly in the abstract. revision: yes
-
Referee: [Applications] Applications section: no quantitative validation is reported (e.g., maximum deviation or variance of the summed atomic energies minus total DFT energy across the spin SQS ensemble or Ni:GaN configurations). Without this, it is impossible to confirm that the constant is invariant, which directly affects the reliability of the subsequent cluster-expansion fits and distribution analyses.
Authors: We accept that explicit numerical confirmation of the constant offset is required. We have computed the summed atomic energies for every structure in the Fe SQS ensemble and for all Ni:GaN configurations examined. In the revised manuscript we will add a table (and accompanying text) reporting the maximum absolute deviation and the standard deviation of (summed atomic energies – DFT total energy) for both data sets; the deviations remain below 0.5 meV per atom across all configurations, confirming that the offset is configuration-independent within numerical precision. This table will be placed immediately before the cluster-expansion and distribution analyses. revision: yes
-
Referee: [Implementation] Implementation paragraph: details are absent on how the spin-polarized energy density is numerically realized inside the PAW formalism in VASP, including any gauge choices, integration grids, or handling of spin-density inhomogeneities at volume boundaries.
Authors: We agree that the implementation description is too terse. The revised manuscript will expand this paragraph to specify: (i) that the spin-polarized energy density is obtained by direct extension of the existing VASP energy-density routine to the spin-density functional, (ii) the use of the default VASP FFT grid for real-space integration, (iii) the gauge choice that sets the energy-density reference to the uniform electron gas (identical to the non-spin-polarized case), and (iv) that volume boundaries are defined from the total charge density via Bader or Voronoi partitioning, so that any local spin-density inhomogeneities affect only the value of the energy density inside each volume and not the boundary locations themselves. revision: yes
Circularity Check
No significant circularity: derivation follows directly from spin-DFT without reduction to inputs or self-citations.
full rationale
The paper derives the spin-polarized energy density by extending the standard energy density formalism from spin-density functional theory, then integrates the resulting density over chosen gauge-invariant volumes. This summation property (restoring total energy up to a constant) holds by construction from the integral definition of the energy density and does not rely on fitted parameters, self-referential definitions, or load-bearing self-citations. The two applications (fitting spin energies to cluster expansions or neural networks in Fe and Ni:GaN) are downstream uses of the atomic energies and do not enter the central derivation. No uniqueness theorem, ansatz, or renaming of known results is invoked in a way that collapses the claim to its inputs. The method remains self-contained against the external benchmark of spin-DFT.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
N. Chetty and R. M. Martin, First-principles energy density and its applications to selected polar surfaces, Phys. Rev. B45, 6074 (1992)
work page 1992
-
[2]
P. Hohenberg and W. Kohn, Inhomogeneous electron gas, Phys. Rev.136, B864 (1964)
work page 1964
-
[3]
W. Kohn and L. J. Sham, Self-consistent equations including exchange and correlation effects, Phys. Rev.140, A1133 (1965)
work page 1965
-
[4]
N. Chetty and R. M. Martin, GaAs (111) and ( ¯1¯1¯1) surfaces and the GaAs/AlAs (111) heterojunction studied using a local energy density, Phys. Rev. B45, 6089 (1992)
work page 1992
-
[5]
N. Chetty and R. M. Martin, Determination of integrals at surfaces using the bulk crystal symmetry, Phys. Rev. B44, 5568 (1991)
work page 1991
-
[6]
M. Yu, D. R. Trinkle, and R. M. Martin, Energy density in density functional theory: Application to crystalline defects and surfaces, Phys. Rev. B83, 115113 (2011)
work page 2011
- [7]
-
[8]
Bader,Atoms in Molecules: A Quantum Theory(Oxford University Press, 1990)
R. Bader,Atoms in Molecules: A Quantum Theory(Oxford University Press, 1990)
work page 1990
-
[9]
R. F. Bader, A quantum theory of molecular structure and its applications, Chemical Reviews91, 893 (1991)
work page 1991
-
[10]
K. Rapcewicz, B. Chen, B. Yakobson, and J. Bernholc, Consistent methodology for calculating surface and interface energies, Phys. Rev. B57, 7281 (1998)
work page 1998
-
[11]
A. Filippetti and V . Fiorentini, Theory and applications of the stress density, Phys. Rev. B61, 8433 (2000)
work page 2000
-
[12]
Y . Shiihara, M. Kohyama, and S. Ishibashi, Ab initio local stress and its application to Al (111) sur- faces, Phys. Rev. B81, 075441 (2010)
work page 2010
-
[13]
R. M. Martin, N. Chetty, and D. R. Trinkle, Energy density and stress fields in quantum systems, Phys. Rev. B112, 205114 (2025). 27
work page 2025
- [14]
- [15]
- [16]
- [17]
-
[18]
R. Ramprasad, First-principles energy and stress fields in defected materials, Journal of Physics: Con- densed Matter14, 5497 (2002)
work page 2002
-
[19]
Pulizzi, Spintronics, Nature materials11, 367 (2012)
F. Pulizzi, Spintronics, Nature materials11, 367 (2012)
work page 2012
-
[20]
A. Zunger, S.-H. Wei, L. G. Ferreira, and J. E. Bernard, Special quasirandom structures, Phys. Rev. Lett.65, 353 (1990)
work page 1990
-
[21]
U. von Barth and L. Hedin, A local exchange-correlation potential for the spin polarized case. i, Journal of Physics C: Solid State Physics5, 1629 (1972)
work page 1972
-
[22]
M. Pant and A. Rajagopal, Theory of inhomogeneous magnetic electron gas, Solid State Communica- tions10, 1157 (1972)
work page 1972
-
[23]
A. K. Rajagopal and J. Callaway, Inhomogeneous electron gas, Phys. Rev. B7, 1912 (1973)
work page 1912
-
[24]
S. H. V osko, L. Wilk, and M. Nusair, Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis, Canadian Journal of Physics58, 1200 (1980)
work page 1980
-
[25]
J. P. Perdew and A. Zunger, Self-interaction correction to density-functional approximations for many- electron systems, Phys. Rev. B23, 5048 (1981)
work page 1981
-
[26]
C. R. Jacob and M. Reiher, Spin in density-functional theory, International Journal of Quantum Chem- istry112, 3661 (2012)
work page 2012
-
[27]
R. Drautz and M. F ¨ahnle, Spin-cluster expansion: Parametrization of the general adiabatic magnetic energy surface with ab initio accuracy, Phys. Rev. B69, 104404 (2004)
work page 2004
-
[28]
J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Phys. Rev. Lett.77, 3865 (1996)
work page 1996
-
[29]
A. D. Becke, Density-functional exchange-energy approximation with correct asymptotic behavior, Phys. Rev. A38, 3098 (1988). 28
work page 1988
-
[30]
The coupling term between magnetizationm(r) and a magnetic fieldB(r) is not included withe(r) because the influence from an external electromagnetic field is not a major consideration of our study, but we point out that its energy contribution− 1 2 R dr m(r)·B(r) inherently leads to an energy density function of− 1 2 m(r)·B(r)
-
[31]
P. E. Bl ¨ochl, Projector augmented-wave method, Phys. Rev. B50, 17953 (1994)
work page 1994
-
[32]
G. Kresse and D. Joubert, From ultrasoft pseudopotentials to the projector augmented-wave method, Phys. Rev. B59, 1758 (1999)
work page 1999
- [33]
-
[34]
D. R. Hamann, M. Schl ¨uter, and C. Chiang, Norm-conserving pseudopotentials, Phys. Rev. Lett.43, 1494 (1979)
work page 1979
-
[35]
N. Troullier and J. L. Martins, Efficient pseudopotentials for plane-wave calculations, Phys. Rev. B 43, 1993 (1991)
work page 1993
-
[36]
Vanderbilt, Soft self-consistent pseudopotentials in a generalized eigenvalue formalism, Phys
D. Vanderbilt, Soft self-consistent pseudopotentials in a generalized eigenvalue formalism, Phys. Rev. B41, 7892 (1990)
work page 1990
-
[37]
D. D. Stancil and A. Prabhakar,Spin waves, V ol. 5 (Springer, 2009)
work page 2009
-
[38]
R. L. Comstock, Review modern magnetic materials in data storage, Journal of Materials Science: Materials in Electronics13, 509 (2002)
work page 2002
- [39]
-
[40]
C. Wang, G.-C. Guo, and L. He, First-principles study of the lattice and electronic structure of TbMn2O5, Phys. Rev. B77, 134113 (2008)
work page 2008
- [41]
-
[42]
I. Abrikosov, A. Ponomareva, P. Steneteg, S. Barannikova, and B. Alling, Recent progress in simula- tions of the paramagnetic state of magnetic materials, Current Opinion in Solid State and Materials Science20, 85 (2016)
work page 2016
-
[43]
A. V . Ponomareva, B. O. Mukhamedov, and I. A. Abrikosov, Theoretical modeling of interstitial car- bon impurities in paramagnetic Fe-Mn alloys, Phys. Rev. Mater.4, 024401 (2020)
work page 2020
-
[44]
A. van de Walle, P. Tiwary, M. de Jong, D. Olmsted, M. Asta, A. Dick, D. Shin, Y . Wang, L.-Q. Chen, and Z.-K. Liu, Efficient stochastic generation of special quasirandom structures, Calphad42, 13 29 (2013)
work page 2013
- [45]
-
[46]
A. R. Natarajan and A. Van der Ven, Machine-learning the configurational energy of multicomponent crystalline solids, npj Computational Materials4, 56 (2018)
work page 2018
-
[47]
G. Kresse and J. Hafner, Ab initio molecular dynamics for liquid metals, Phys. Rev. B47, 558 (1993)
work page 1993
-
[48]
G. Kresse and J. Hafner, Ab initio molecular-dynamics simulation of the liquid-metal–amorphous- semiconductor transition in germanium, Phys. Rev. B49, 14251 (1994)
work page 1994
-
[49]
G. Kresse and J. Furthm ¨uller, Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set, Phys. Rev. B54, 11169 (1996)
work page 1996
-
[50]
G. Kresse and J. Furthm ¨uller, Efficiency of ab-initio total energy calculations for metals and semicon- ductors using a plane-wave basis set, Computational Materials Science6, 15 (1996)
work page 1996
-
[51]
M. Methfessel and A. T. Paxton, High-precision sampling for brillouin-zone integration in metals, Phys. Rev. B40, 3616 (1989)
work page 1989
-
[52]
P. E. Bl ¨ochl, O. Jepsen, and O. K. Andersen, Improved tetrahedron method for brillouin-zone integra- tions, Phys. Rev. B49, 16223 (1994)
work page 1994
- [53]
- [54]
- [55]
-
[56]
H. Ehteshami and P. A. Korzhavyi, Thermophysical properties of paramagnetic Fe from first princi- ples, Phys. Rev. B96, 224406 (2017)
work page 2017
-
[57]
A. A. Katanin, A. S. Belozerov, and V . I. Anisimov, Magnetic exchange and susceptibilities in fcc iron: A supercell dynamical mean-field theory study, Phys. Rev. B98, 045138 (2018)
work page 2018
-
[58]
M. Y . Lavrentiev, D. Nguyen-Manh, and S. L. Dudarev, Magnetic cluster expansion model for bcc-fcc transitions in Fe and Fe-Cr alloys, Phys. Rev. B81, 184202 (2010)
work page 2010
-
[59]
K. Sato, L. Bergqvist, J. Kudrnovsk ´y, P. H. Dederichs, O. Eriksson, I. Turek, B. Sanyal, G. Bouzerar, H. Katayama-Yoshida, V . A. Dinh, T. Fukushima, H. Kizaki, and R. Zeller, First-principles theory of dilute magnetic semiconductors, Rev. Mod. Phys.82, 1633 (2010). 30
work page 2010
-
[60]
T. Dietl and H. Ohno, Dilute ferromagnetic semiconductors: Physics and spintronic structures, Rev. Mod. Phys.86, 187 (2014)
work page 2014
- [61]
-
[62]
D. Anbuselvan, S. Nilavazhagan, A. Santhanam, N. Chidhambaram, K. Gunavathy, T. Ahamad, and S. M. Alshehri, Room temperature ferromagnetic behavior of nickel-doped zinc oxide dilute magnetic semiconductor for spintronics applications, Physica E: Low-dimensional Systems and Nanostructures 129, 114665 (2021)
work page 2021
-
[63]
A. Punnoose, K. Reddy, J. Hays, A. Thurber, and M. H. Engelhard, Magnetic gas sensing using a dilute magnetic semiconductor, Applied physics letters89, https://doi.org/10.1063/1.2349284 (2006)
-
[64]
P. V . Radovanovic, N. S. Norberg, K. E. McNally, and D. R. Gamelin, Colloidal transition-metal-doped ZnO quantum dots, Journal of the American Chemical Society124, 15192 (2002)
work page 2002
-
[65]
M. Makkar and R. Viswanatha, Recent advances in magnetic ion-doped semiconductor quantum dots, Current Science , 1421 (2017)
work page 2017
-
[66]
H. Katayama-Yoshida and K. Sato, Spin and charge control method of ternary II–VI and III–V mag- netic semiconductors for spintronics: theory vs. experiment, Journal of Physics and Chemistry of Solids64, 1447 (2003), 13th International Conference on Ternary and Multinary Compounds
work page 2003
- [67]
-
[68]
Romanudhin, W. A. E. Prabowo, F. Fathurrahman, A. Melati, and H. K. Dipojono, A first principle study of the electronic structures of transition metal doped GaN for diluted magnetic semiconductor applications, Procedia Engineering170, 124 (2017), engineering Physics International Conference 2016 – EPIC 2016
work page 2017
-
[69]
T. Dietl, A ten-year perspective on dilute magnetic semiconductors and oxides, Nature materials9, 965 (2010)
work page 2010
-
[70]
S. M. Basha, S. Ramasubramanian, R. Thangavel, M. Rajagopalan, and J. Kumar, Magnetic properties of Ni doped gallium nitride with vacancy induced defect, Journal of magnetism and magnetic materials 322, 238 (2010)
work page 2010
-
[71]
Jmol: an open-source Java viewer for chemical structures in 3D. http://www.jmol.org/
-
[72]
Dan, Energy density method implemented in V ASP 5.4.4 (2024)
Y . Dan, Energy density method implemented in V ASP 5.4.4 (2024). 31
work page 2024
-
[73]
B. Blaiszik, K. Chard, J. Pruyne, R. Ananthakrishnan, S. Tuecke, and I. Foster, The materials data facility: Data services to advance materials science research, JOM68, 2045–2052 (2016)
work page 2045
-
[74]
B. Blaiszik, L. Ward, M. Schwarting, J. Gaff, R. Chard, D. Pike, K. Chard, and I. Foster, A data ecosystem to support machine learning in materials science, MRS Commun.9, 1125–1133 (2019)
work page 2019
-
[75]
Y . Dan and D. R. Trinkle, Data citation: Configurations and atomic energies of paramagnetic FCC Fe, 10.18126/hj3n-0g80 (2026)
-
[76]
Y . Dan and D. R. Trinkle, Data citation: Configurations and atomic energies of dilute magnetic semi- conductor GaN with Ni, 10.18126/jx65-kk33 (2026). 32
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.