Recognition: 3 theorem links
· Lean TheoremAb initio evidence for a framework-preserving spin-polarized high-DOS state in D-type carbon schwarzite C136
Pith reviewed 2026-05-08 18:58 UTC · model grok-4.3
The pith
Spin-polarized relaxation of D-type carbon schwarzite C136 preserves the framework while maintaining a high density of states at the Fermi level.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Ab initio calculations on C136 show that fixed-cell spin-polarized relaxation reduces energy by approximately 0.213 eV per 136-atom cell with an RMS atomic displacement of 0.098 Angstrom and maximum displacement of 0.200 Angstrom. The resulting geometry supports a magnetic state with total energy -2490.35442340 Ry, total magnetization 10.63 muB/cell, absolute magnetization 12.94 muB/cell, and a robust high density of states at the Fermi level that remains essentially unchanged when the k-grid is refined from 3x3x3 to 4x4x4.
What carries the argument
The fixed-cell spin-polarized ionic relaxation combined with subsequent self-consistent-field and spin-resolved density-of-states calculations performed on the high-symmetry D-type schwarzite framework.
If this is right
- The moderate distortion preserves the negative-curvature carbon network and its high density of states near the Fermi level.
- C136 behaves as a parent phase near coupled spin-lattice and high-DOS electronic instabilities.
- The high-DOS character is robust against k-point refinement between 3x3x3 and 4x4x4 grids.
- The results motivate searches for stabilized, distorted, doped, or intercalated descendants of D-type carbon schwarzites.
Where Pith is reading between the lines
- Other negative-curvature carbon networks may exhibit similar responses to spin polarization.
- Doping or intercalation could be used to stabilize the observed magnetic or high-DOS features for further study.
- Higher-level methods or variable-cell relaxations would provide a direct test of whether the reported magnetization and DOS persist.
Load-bearing premise
That the fixed-cell spin-polarized relaxation and follow-on SCF calculation capture the physical ground state without requiring full variable-cell optimization or validation against experiment or higher-level theory.
What would settle it
A full variable-cell optimization that produces atomic displacements much larger than 0.2 Angstrom, introduces bonds shorter than 1.2 Angstrom, or yields a lower-energy non-magnetic state would falsify the claim of framework preservation and high-DOS retention.
Figures
read the original abstract
Negative Gaussian curvature provides an unusual route for designing electronic structure in extended sp2 carbon networks. Here I report ab initio density-functional calculations on the D-type carbon schwarzite C136, focusing on the response of the ideal high-symmetry framework to spin polarization and fixed-cell ionic distortion. A partial spin-polarized fixed-cell relaxation lowers the total energy by approximately 0.213 eV per 136-atom cell over six completed ionic steps. The distortion remains moderate: the RMS atomic displacement is approximately 0.098 Angstrom, the maximum atomic displacement is approximately 0.200 Angstrom, the RMS C-C bond-length change for the 170 reference bonds shorter than 1.80 Angstrom is only approximately 0.0107 Angstrom, and no unphysically short C-C contacts below 1.20 Angstrom are found. A separate clean from-scratch spin-polarized SCF calculation on the last saved distorted geometry converges successfully to a magnetic state with total energy -2490.35442340 Ry, total magnetization 10.63 muB/cell, and absolute magnetization 12.94 muB/cell. Spin-resolved DOS calculations further show that the distorted geometry retains a high density of states near the Fermi level. A 3x3x3 diagnostic DOS gives N(EF) approximately 42.84 states/eV/cell, while a 4x4x4 validation DOS gives N(EF) approximately 42.85 states/eV/cell, demonstrating that the high-DOS character is robust with respect to this k-point refinement. These results support the interpretation of C136 as a negative-curvature carbon parent phase near coupled spin-lattice and high-DOS electronic instabilities. Superconductivity is not established here; rather, the results motivate a search for stabilized, distorted, doped, or intercalated descendants of D-type carbon schwarzites.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports ab initio density-functional calculations on D-type carbon schwarzite C136, showing that a spin-polarized fixed-cell ionic relaxation over six steps lowers the total energy by ~0.213 eV per 136-atom cell with moderate distortions (RMS displacement 0.098 Å, RMS bond-length change 0.0107 Å) while preserving a high density of states at the Fermi level (~42.85 states/eV/cell across 3x3x3 and 4x4x4 k-grids) and yielding a magnetization of 10.63 μB/cell; this is interpreted as evidence that C136 is a parent phase near coupled spin-lattice and high-DOS electronic instabilities.
Significance. If the reported state is robust, the work supplies concrete numerical evidence that negative-curvature sp2 carbon frameworks can sustain spin polarization and high DOS without framework collapse, providing a computational foundation for exploring schwarzite-derived materials with coupled magnetic and electronic instabilities. The k-grid convergence check on N(EF) and the small, physically reasonable distortions are strengths of the presented data.
major comments (2)
- [Abstract] Abstract and relaxation procedure: The central claim of proximity to 'coupled spin-lattice ... instabilities' rests on a fixed-cell relaxation that by construction excludes volume, shear, or lattice-strain modes; the modest 0.213 eV energy drop and preserved high DOS could change under full variable-cell optimization, which is required to substantiate the interpretation.
- Computational methodology: The manuscript provides no information on the exchange-correlation functional, pseudopotentials, plane-wave cutoff, or convergence thresholds, rendering the quantitative outputs (energy difference, magnetization, DOS values) impossible to reproduce or assess for systematic errors.
minor comments (1)
- [Abstract] The abstract reports several quantities with 'approximately'; supplying the exact computed values together with any estimated uncertainties would improve precision and allow direct comparison with future work.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and the constructive comments. We address each of the major comments below and will make the necessary revisions to improve the paper.
read point-by-point responses
-
Referee: [Abstract] Abstract and relaxation procedure: The central claim of proximity to 'coupled spin-lattice ... instabilities' rests on a fixed-cell relaxation that by construction excludes volume, shear, or lattice-strain modes; the modest 0.213 eV energy drop and preserved high DOS could change under full variable-cell optimization, which is required to substantiate the interpretation.
Authors: The referee correctly notes that our fixed-cell relaxation excludes volume and shear modes. Our intent was to demonstrate that a high-DOS spin-polarized state can be sustained under small, framework-preserving ionic distortions without causing the structure to collapse, as the title specifies 'framework-preserving'. The modest energy gain and unchanged high DOS support the idea that C136 is near such instabilities even when lattice parameters are held fixed. However, we acknowledge that a variable-cell calculation is needed for a complete assessment. In the revision, we will add a paragraph discussing this limitation and include results from a variable-cell relaxation to check the robustness of the high-DOS state. revision: yes
-
Referee: [—] Computational methodology: The manuscript provides no information on the exchange-correlation functional, pseudopotentials, plane-wave cutoff, or convergence thresholds, rendering the quantitative outputs (energy difference, magnetization, DOS values) impossible to reproduce or assess for systematic errors.
Authors: We regret the oversight in not including the computational details. We will revise the manuscript to include a section on the computational methodology, specifying the exchange-correlation functional, pseudopotentials, plane-wave cutoff, and all convergence thresholds used in the calculations. revision: yes
Circularity Check
No significant circularity; results are direct ab initio outputs
full rationale
The paper presents numerical results from density-functional calculations: an energy lowering of 0.213 eV/cell after six fixed-cell ionic steps, RMS displacement of 0.098 Å, total magnetization of 10.63 μB/cell, and N(EF) ≈ 42.85 states/eV/cell from spin-resolved DOS on the relaxed geometry. These quantities are computed outputs, not quantities defined in terms of themselves or fitted to match a target prediction. No equations, ansatzes, or self-citations are used to derive the central interpretation; the claim that C136 is near coupled instabilities follows from the independent simulation results rather than reducing to input assumptions by construction. The fixed-cell constraint is a methodological choice whose consequences are stated explicitly, not a hidden definitional loop.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Density functional theory accurately captures the relative energetics and electronic structure of spin-polarized carbon schwarzites.
Forward citations
Cited by 1 Pith paper
-
Hole-Doping Suppresses Competing Magnetism in High-DOS C136 Carbon Schwarzite: A Computational Route Toward Superconductivity in Negative-Curvature Carbon Networks
Hole doping in C136 schwarzite monotonically suppresses magnetism from ~11 to 4.76 Bohr magnetons per cell while preserving high DOS of ~44.7 states/eV/cell at the Fermi level.
Reference graph
Works this paper leans on
-
[1]
A. L. Mackay and H. Terrones, “Diamond from graphite,”Nature352, 762 (1991). doi:10.1038/352762a0
-
[2]
The geometry of hypothetical curved graphite structures,
H. Terrones and A. L. Mackay, “The geometry of hypothetical curved graphite structures,” Carbon30, 1251–1260 (1992). doi:10.1016/0008-6223(92)90066-6
-
[3]
Negative-curvature fullerene analog of C60,
D. Vanderbilt and J. Tersoff, “Negative-curvature fullerene analog of C60,”Physical Review Letters68, 511–513 (1992). doi:10.1103/PhysRevLett.68.511
-
[4]
Energetics of negatively curved graphitic carbon,
T. Lenosky, X. Gonze, M. Teter, and V. Elser, “Energetics of negatively curved graphitic carbon,” Nature355, 333–335 (1992). doi:10.1038/355333a0
-
[5]
Negatively curved graphitic sheet model of amorphous carbon,
S. J. Townsend, T. J. Lenosky, D. A. Muller, C. S. Nichols, and V. Elser, “Negatively curved graphitic sheet model of amorphous carbon,”Physical Review Letters69, 921–924 (1992). doi:10.1103/PhysRevLett.69.921
-
[6]
Generalized gradient approximation made simple
J. P. Perdew, K. Burke, and M. Ernzerhof, “Generalized gradient approximation made simple,” Physical Review Letters77, 3865–3868 (1996). doi:10.1103/PhysRevLett.77.3865. 8
-
[7]
doi:10.1088/0953-8984/21/39/395502 , url =
P. Giannozziet al., “QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials,”Journal of Physics: Condensed Matter21, 395502 (2009). doi:10.1088/0953-8984/21/39/395502
-
[8]
Advanced capabilities for materials modelling with Quantum
P. Giannozziet al., “Advanced capabilities for materials modelling with Quantum ESPRESSO,” Journal of Physics: Condensed Matter29, 465901 (2017). doi:10.1088/1361-648X/aa8f79
-
[9]
Pseudopotentials periodic table: From H to Pu,
A. Dal Corso, “Pseudopotentials periodic table: From H to Pu,”Computational Materials Science95, 337–350 (2014). doi:10.1016/j.commatsci.2014.07.043
-
[10]
H. J. Monkhorst and J. D. Pack, “Special points for Brillouin-zone integrations,”Physical Review B13, 5188–5192 (1976). doi:10.1103/PhysRevB.13.5188
-
[11]
Thermal Contraction and Disordering of the Al(110) Surface,
N. Marzari, D. Vanderbilt, A. De Vita, and M. C. Payne, “Thermal contraction and disordering of the Al(110) surface,”Physical Review Letters82, 3296–3299 (1999). doi:10.1103/PhysRevLett.82.3296
-
[12]
Collective electron ferromagnetism,
E. C. Stoner, “Collective electron ferromagnetism,”Proceedings of the Royal Society of London A165, 372–414 (1938)
1938
-
[13]
Ferromagnetism in single-band Hubbard models with a partially flat band,
A. Mielke, “Ferromagnetism in single-band Hubbard models with a partially flat band,”Physical Review Letters82, 4312–4315 (1999). doi:10.1103/PhysRevLett.82.4312
-
[14]
Ferromagnetism in Hubbard models,
H. Tasaki, “Ferromagnetism in Hubbard models,”Physical Review Letters75, 4678–4681 (1995). doi:10.1103/PhysRevLett.75.4678
-
[15]
R. E. Peierls,Quantum Theory of Solids. Oxford University Press, Oxford (1955)
1955
-
[16]
Gruner,Density Waves in Solids
G. Gruner,Density Waves in Solids. Addison-Wesley, Reading, MA (1994)
1994
-
[17]
First principles phonon calculations in materials science
A. Togo and I. Tanaka, “First principles phonon calculations in materials science,”Scripta Materialia108, 1–5 (2015). doi:10.1016/j.scriptamat.2015.07.021. 9
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.