Pith. sign in

REVIEW 3 major objections 4 minor 26 references

Exotic massive fermionic systems with huge vacuum degeneracy at boundaries

T0 review · 3 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash

Pith's one-line read The paper claims that a massive, otherwise trivially gapped free fermionic system develops a boundary vacuum degeneracy whose residual entropy is proportional to the area of the boundary.

desk verdict The Case (i) boundary zero-mode result is real and clean; the paper oversells Case (ii). read the letter →

arxiv 2501.07935 v1 pith:CZDUTCMK submitted 2025-01-14 hep-th cond-mat.str-el

classification hep-thcond-mat.str-el
keywords massiveψtheoryfermionicsubsystemsymmetryboundaryzeromodesvacuumdegeneracyresidualentropyarealawMajoranagapped
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

The paper studies the massive $\psi$ theory, a free non-relativistic fermionic model with two real fields $\psi_+$ and $\psi_-$, on a slab. The mass term explicitly breaks the fermionic subsystem symmetry and opens a bulk gap of size $|\mu|$. With the same-sign boundary condition $\gamma_y\gamma_z \psi_+ = \gamma_y\gamma_z \psi_-$ at both boundaries, the authors find exactly one Majorana zero mode for each transverse momentum pair, all with exactly zero frequency, so the slab has $2^{N/2}$ ground states with $N\sim L_y L_z/\varepsilon^2$ and residual entropy $(N/2)\log 2$ proportional to the boundary area. All nonzero modes have frequency at least $|\mu|$, so the system remains gapped despite the large degeneracy. With the opposite-sign boundary condition the zero modes acquire exponentially small frequencies, so the degeneracy becomes only approximate.

What carries the argument

The load-bearing object is the variational boundary condition $\gamma_y\gamma_z \psi_+ = \pm \gamma_y\gamma_z \psi_-$ at each boundary, chosen so that the boundary term in the variation of the action vanishes. The mode analysis superposes incident and reflected plane waves with evanescent boundary modes $\psi_\pm \sim A e^{-\Lambda x + ik_y y + ik_z z}$; the determinant condition for the boundary matrices selects $\Lambda = \mu/(k_y k_z)$ and forces the local dispersion $\omega_{\mathrm{local}}(\Lambda,k_y,k_z) = \sqrt{-\Lambda^2 k_y^2 k_z^2 + \mu^2}$ to vanish, producing the exact zero modes. The evenness of the Majorana zero-mode count under complex conjugation is what lets the same-sign boundary condition remain consistent.

What would settle it

A lattice discretization of the massive $\psi$ theory on a slab with boundary condition (i) could settle the claim: if numerical diagonalization finds no exact zero modes, or an odd number of Majorana zero modes, or a non-Hermitian Hamiltonian under this boundary condition, the central claim fails. A second check is to classify all self-adjoint boundary conditions for the continuum slab and test whether condition (i) is among them.

Watch

Extended reading notes

Core claim

The central claim is that the massive $\psi$ theory with boundary condition (i) develops boundary-localized Majorana zero modes: one for each transverse momentum pair $(k_y,k_z)$, with exactly vanishing frequency and profile decaying into the slab at rate $|\mu/(k_y k_z)|$. Because complex conjugation pairs these zero modes, their total number is even, so this boundary condition is consistent, unlike the analogous same-sign condition in the 1+1-dimensional Majorana fermion. The vacuum degeneracy is therefore $2^{N/2}$ with $N\sim L_yL_z/\varepsilon^2$, giving a residual entropy $(N/2)\log 2$ proportional to the boundary area. The bulk spectrum keeps a gap $|\mu|$, so the system is gapped and degenerate at the same time. For boundary condition (ii), the localized modes occur in the same number but with exponentially small frequencies, so the degeneracy is only approximate.

Load-bearing premise

The result rests on the same-sign boundary condition being a physically allowed way to close the slab; the paper does not prove that the listed boundary conditions are the only possible ones or that the mode expansion is complete, so if that boundary condition is inconsistent the zero modes and the area-law entropy would not occur.

Editorial extensions

If this is right

  • With boundary condition (i), the slab vacuum degeneracy is $2^{N/2}$ with $N\sim L_yL_z/\varepsilon^2$, so the residual entropy is $(N/2)\log 2$ and scales with boundary area.
  • In case (i) every nonzero frequency satisfies $|\omega|\geq |\mu|$, so the system remains gapped even with exact boundary zero modes.
  • The degeneracy is not protected by the fermionic subsystem symmetry, since the mass term explicitly breaks it; another, unidentified symmetry or mechanism must be responsible.
  • The same-sign boundary condition avoids the obstruction that forbids its analog in the 1+1-dimensional Majorana chain, because the Majorana zero modes here occur in even numbers.
  • In case (ii) the localized modes appear with exponentially small positive and negative frequencies, so the degeneracy is approximate and the first excited gap is exponentially small rather than of order $|\mu|$.

Reading between the lines

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

  • If the counting survives regularization, the model is a free-field example of boundary-dominated entropy without gapless edge modes, distinct from the anomaly-inflow mechanism usually invoked for subsystem-symmetry systems.
  • A direct extension is to put boundaries orthogonal to the $y$ or $z$ directions or to add corners; the same counting logic suggests corner contributions to the residual entropy, which the paper leaves open.
  • The unidentified protecting symmetry could be probed by perturbing the Hamiltonian with terms that break discrete translations or 90-degree rotations in the transverse directions; if the degeneracy lifts, that would point to the operative symmetry.
  • The low-temperature partition function of the slab should contain a boundary factor $\exp((N/2)\log 2)$ in case (i), a quantity that could be computed directly as a check.
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

3 major / 4 minor

Summary. The paper studies a (3+1)-dimensional free fermionic theory (the massive ψ theory) on a slab with two transverse directions compactified. The authors derive boundary conditions from the variational principle. For boundary condition (i), where ψ_+ and ψ_- are identified (up to a γ factor) at both boundaries, they find one exact zero-energy boundary-localized Majorana mode for each transverse momentum pair, leading to a vacuum degeneracy 2^{N/2} with N ~ L_y L_z / ε^2 and a residual entropy proportional to the boundary area, while all nonzero modes are gapped by |μ|. For boundary condition (ii), they find approximate zero modes with exponentially small frequencies and claim that the vacuum degeneracy and the gap are the same as in Case (i). The paper contrasts this behavior with SPT/SSPT phases and notes that the protecting symmetry, if any, is not yet identified.

Significance. If the Case (i) construction is a well-defined quantum theory, it provides a concrete free-field example of a gapped bulk that develops a huge, area-law vacuum degeneracy at a boundary, without the degeneracy being protected by the explicitly broken subsystem symmetry. The analytic solution is explicit, self-contained, and has no fitted parameters, which are clear strengths. However, the paper currently overstates the Case (ii) conclusion, and it omits a completeness/self-adjointness proof for the mode decomposition; both points need to be addressed before the general claim can be taken at face value.

major comments (3)
  1. [Sec. 3.2.2, Eq. (3.27)] The statement that the number of approximate zero modes in Case (ii) is the same as in Case (i) is not supported by Eq. (3.27). Since tanh(Λ L_x) is positive and bounded above by 1 for Λ > 0, a nonzero solution exists only when k_y k_z / μ > 0 and k_y k_z < μ L_x. Because k_y = (2n+1)π/L_y and k_z = (2m+1)π/L_z, only roughly half of the transverse momentum pairs satisfy k_y k_z > 0, and even those with k_y k_z ≥ μ L_x are excluded. Thus the number of approximate zero modes is at most about half of the Case (i) count, and can be further reduced if μ L_x is not large enough. Consequently, the claim in Sec. 4 that Case (ii) has the same vacuum degeneracy and the same gap as Case (i) is incorrect: the approximate modes have nonzero frequencies, so the exact ground state is unique, and the gap to the first excited state is exponentially small rather than O(|μ|). The counting and gap statements should be corrected, or the conclusion should be restricted to Case (i).
  2. [Sec. 3.1, Eq. (3.5)] The derivation shows that the boundary conditions (3.5) are sufficient to make the variational boundary term vanish, but it does not prove that they exhaust all possible self-adjoint boundary conditions for the Dirac-like operator, nor that condition (i) defines a well-defined Hamiltonian on the slab. This is load-bearing because the central result—one exact Majorana zero mode per transverse momentum pair—depends on the admissibility of condition (i). Please either prove that (3.5) are necessary (or at least that condition (i) yields a self-adjoint Hamiltonian with a complete mode expansion), or state this as an explicit assumption. Without such a proof, the physical realizability of the area-law residual entropy remains an open point.
  3. [Sec. 3.2 (mode expansion)] The paper lists bulk modes and boundary-localized modes but does not demonstrate that these modes form a complete set, that they satisfy the canonical anticommutation relations, or that they diagonalize the Hamiltonian. In particular, the vacuum degeneracy 2^{N/2} assumes that the zero modes are the only states at zero energy and that their number is exactly N. Please provide a completeness argument for the eigenfunctions of the differential operator with boundary condition (i), or cite a standard result that guarantees it. This is needed to turn the explicit mode solutions into a rigorous counting of the ground states.
minor comments (4)
  1. [Title] The title contains an apparent typo: "degenera cy" should be "degeneracy".
  2. [Sec. 3.2.2, Eq. (3.28)] The phrase "reminder term" should be "remainder term".
  3. [Sec. 3.2, counting N] The UV cutoff ε is introduced only in a footnote; because the area-law entropy N ~ L_y L_z / ε^2 depends on this cutoff, it should be stated prominently in the main text when the number of zero modes is counted.
  4. [Sec. 3.2.2, discussion of (3.27)] The condition for a nonzero solution of Eq. (3.27) is stated as "if and only if 1 < μ L_x / (k_y k_z)"; this implicitly assumes k_y k_z > 0. The sign condition should be stated explicitly to avoid confusion.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the boundary zero-mode spectrum is solved from the equations of motion with no fitted parameters or load-bearing self-citations.

full rationale

The paper's central result is an explicit free-field spectrum computation on a slab. The action, the boundary conditions (3.5), and the mode ansätze (3.9) and (3.16) fully determine the algebra; the zero-mode frequencies and the local mode conditions are obtained by solving determinant equations (3.18) and (3.27), not by assuming the conclusion. No parameter is fitted to the degeneracy, and no previous result is used to infer the number of boundary modes. The self-citations to the same author's earlier work ([6], [22]) introduce the model and motivate the gapped/boundary context, but the boundary calculation is restated independently and does not rely on those papers' conclusions. The paper's assertion that Case (ii) has the same number of approximate zero modes as Case (i) appears arithmetically inconsistent with Eq. (3.27), which admits nonzero Lambda only for 0 < k_y k_z/m < L_x; however, that is a potential correctness or overstatement issue, not a circularity. Thus the derivation is self-contained and the circularity score is 0.

Assumptions & free parameters 2 free parameters · 4 assumptions · 0 invented entities

The central claim rests on the free-field mode expansion and on the admissibility of the chosen boundary conditions. No fitted parameters or invented entities appear. The mass and cutoff are model inputs, and the main unproven input is the completeness and consistency of the boundary condition classification.

free parameters (2)
  • mass μ = nonzero real (not fitted)
    Model parameter in the action (2.2). The result holds for any μ ≠ 0; the gap scale and the condition for zero modes depend on μ, but no value is fitted to data.
  • UV cutoff / lattice spacing ε = implicit, not fixed
    Used to count boundary modes, giving N ~ L_y L_z / ε^2. The residual entropy coefficient is non-universal and depends on this cutoff.
assumptions (4)
  • standard math The γ_± matrices satisfy the Clifford algebra identities needed to derive the equations of motion (3.7) and the mode solutions.
    The mode expansion in Sec. 2 and the boundary-localized ansatz in Sec. 3.2 rely on standard gamma matrix algebra, including (γ_yγ_z)^2 = -1.
  • domain assumption The free-field mode expansion and canonical anticommutation relations extend to a manifold with boundaries without subtle consistency conditions beyond the vanishing of the boundary variation term (3.4).
    The paper quantizes by mode expansion on the slab without explicitly proving self-adjointness of the Hamiltonian or completeness of the mode set.
  • ad hoc to paper The boundary conditions (3.5), with ψ_+ = ± ψ_- at each boundary, are the possible and sufficient conditions to make the variational boundary term vanish.
    Four sign combinations are listed and two are treated. The paper does not prove that these exhaust all admissible boundary conditions or that the same-sign condition (i) is anomaly-free beyond the evenness of the Majorana zero mode count.
  • domain assumption Anti-periodic boundary conditions in the y and z directions are imposed for simplicity, ensuring k_y and k_z are never zero.
    This choice in Eq. (3.2) avoids a (0,0) mode and keeps the zero-mode counting finite; the result is expected to hold more generally but is only demonstrated here.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Exotic massive fermionic systems with huge vacuum degeneracy at boundaries." pith.science (2026). https://pith.science/paper/CZDUTCMK

@misc{pith2026250107935,
  author       = {Pith},
  title        = {Pith review of: Exotic massive fermionic systems with huge vacuum degeneracy at boundaries},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/CZDUTCMK}},
  note         = {Machine review of arXiv:2501.07935}
}
read the original abstract

We investigate a massive non-relativistic fermionic system exhibiting exotic features. When the mass parameter is set to zero, the system acquires the fermionic subsystem symmetry. Introducing the mass term explicitly breaks this symmetry, resulting in a trivially gapped system in the absence of boundaries. We demonstrate that with the introduction of boundaries, the system remains gapped, but it has a huge vacuum degeneracy. The residual entropy is proportional to the area of the boundary.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

26 extracted references · 4 canonical work pages

  1. [1]

    Quantum Glassiness,

    C. Chamon, “Quantum Glassiness,” Phys. Rev. Lett. 94 no. 4, (2005) 040402 , arXiv:cond-mat/0404182

  2. [2]

    Local stabilizer codes in three dimensions wit hout string logical operators,

    J. Haah, “Local stabilizer codes in three dimensions wit hout string logical operators,” Phys. Rev. A 83 no. 4, (2011) 042330 , arXiv:1101.1962 [quant-ph]

  3. [3]

    Fractons,

    R. M. Nandkishore and M. Hermele, “Fractons,” Ann. Rev. Condensed Matter Phys. 10 (2019) 295–313 , arXiv:1803.11196 [cond-mat.str-el]

  4. [4]

    Fracton Phases of Matter,

    M. Pretko, X. Chen, and Y . Y ou, “Fracton Phases of Matter, ” Int. J. Mod. Phys. A 35 no. 06, (2020) 2030003 , arXiv:2001.01722 [cond-mat.str-el]

  5. [5]

    Colloquium: Fracton matt er,

    A. Gromov and L. Radzihovsky, “Colloquium: Fracton matt er,” Rev. Mod. Phys. 96 no. 1, (2024) 011001 , arXiv:2211.05130 [cond-mat.str-el]

  6. [6]

    Supersymmetric quantum field theory with exotic symmetry in 3+1 dimensions and fermionic fracton phases,

    S. Y amaguchi, “Supersymmetric quantum field theory with exotic symmetry in 3+1 dimensions and fermionic fracton phases,” PTEP 2021 no. 6, (2021) 063B04 , arXiv:2102.04768 [hep-th]

  7. [7]

    Field Theories With a Vector Global Symmetr y,

    N. Seiberg, “Field Theories With a Vector Global Symmetr y,” SciPost Phys. 8 no. 4, (2020) 050 , arXiv:1909.10544 [cond-mat.str-el]

  8. [8]

    Exotic Symmetries, Duality, and Fractons in 2+1-Dimensional Quantum Field Theory,

    N. Seiberg and S.-H. Shao, “Exotic Symmetries, Duality, and Fractons in 2+1-Dimensional Quantum Field Theory,” arXiv:2003.10466 [cond-mat.str-el]

Show all 26 references
  1. [9]

    Exotic /u1D448(1) Symmetries, Duality, and Fractons in 3+1-Dimensional Quantum Field Theory,

    N. Seiberg and S.-H. Shao, “Exotic /u1D448(1) Symmetries, Duality, and Fractons in 3+1-Dimensional Quantum Field Theory,” SciPost Phys. 9 no. 4, (2020) 046 , arXiv:2004.00015 [cond-mat.str-el]

  2. [10]

    More Exotic Field Theories in 3+1 Dimensions,

    P . Gorantla, H. T. Lam, N. Seiberg, and S.-H. Shao, “More Exotic Field Theories in 3+1 Dimensions,” SciPost Phys. 9 (2020) 073, arXiv:2007.04904 [cond-mat.str-el]

  3. [11]

    Spontaneously broken supe rsymmetric fracton phases with fermionic subsystem symmetries,

    H. Katsura and Y . Nakayama, “Spontaneously broken supe rsymmetric fracton phases with fermionic subsystem symmetries,” JHEP 08 (2022) 072, arXiv:2204.01924 [hep-th] . 11

  4. [12]

    Scalar, fermionic and super symmetric field theories with subsystem symmetries in d + 1 dimensions,

    M. Honda and T. Nakanishi, “Scalar, fermionic and super symmetric field theories with subsystem symmetries in d + 1 dimensions,” JHEP 03 (2023) 188, arXiv:2212.13006 [hep-th]

  5. [13]

    Building fracton phases by majo rana manipulation,

    Y . Y ou and F. von Oppen, “Building fracton phases by majo rana manipulation,” Physical Review Research 1 no. 1, (2019) 013011

  6. [14]

    Jordan-Wigner Dualities for Tra nslation-Invariant Hamiltonians in Any Dimension: Emergent Fermions in Fracton Topological Or der,

    N. Tantivasadakarn, “Jordan-Wigner Dualities for Tra nslation-Invariant Hamiltonians in Any Dimension: Emergent Fermions in Fracton Topological Or der,” Phys. Rev. Res. 2 no. 2, (2020) 023353 , arXiv:2002.11345 [cond-mat.str-el]

  7. [15]

    Fractonic order and emergent fermionic ga uge theory,

    W. Shirley, “Fractonic order and emergent fermionic ga uge theory,” arXiv:2002.12026 [cond-mat.str-el]

  8. [16]

    Boson-fermion duali ty with subsystem symmetry,

    W. Cao, M. Y amazaki, and Y . Zheng, “Boson-fermion duali ty with subsystem symmetry,” Phys. Rev. B 106 no. 7, (2022) 075150 , arXiv:2206.02727 [cond-mat.str-el]

  9. [17]

    Anomalies and Fermio n Zero Modes on Strings and Domain Walls,

    C. G. Callan, Jr. and J. A. Harvey, “Anomalies and Fermio n Zero Modes on Strings and Domain Walls,” Nucl. Phys. B 250 (1985) 427–436

  10. [18]

    Subs ystem symmetry protected topological order,

    Y . Y ou, T. Devakul, F. J. Burnell, and S. L. Sondhi, “Subs ystem symmetry protected topological order,” Phys. Rev. B 98 no. 3, (2018) 035112 , arXiv:1803.02369 [cond-mat.str-el]

  11. [19]

    Classificatio n of subsystem symmetry-protected topological phases,

    T. Devakul, D. J. Williamson, and Y . Y ou, “Classificatio n of subsystem symmetry-protected topological phases,” Phys. Rev. B 98 no. 23, (2018) 235121 , arXiv:1808.05300 [cond-mat.str-el]

  12. [20]

    Strong planar subs ystem symmetry-protected topological phases and their dual fracton orders,

    T. Devakul, W. Shirley, and J. Wang, “Strong planar subs ystem symmetry-protected topological phases and their dual fracton orders,” Phys. Rev. Res. 2 no. 1, (2020) 012059 , arXiv:1910.01630 [cond-mat.str-el]

  13. [21]

    Anomaly inflow for subsystem symmetries,

    F. J. Burnell, T. Devakul, P . Gorantla, H. T. Lam, and S.- H. Shao, “Anomaly inflow for subsystem symmetries,” Phys. Rev. B 106 no. 8, (2022) 085113 , arXiv:2110.09529 [cond-mat.str-el]

  14. [22]

    Gapless edge modes in (4+1)-dimensiona l topologically massive tensor gauge theory and anomaly inflow for subsystem symmetry,

    S. Y amaguchi, “Gapless edge modes in (4+1)-dimensiona l topologically massive tensor gauge theory and anomaly inflow for subsystem symmetry,” PTEP 2022 no. 3, (2022) 033B08 , arXiv:2110.12861 [hep-th]

  15. [23]

    Anomaly of subsystem symmetries in exot ic and foliated BF theories,

    S. Shimamura, “Anomaly of subsystem symmetries in exot ic and foliated BF theories,” JHEP 06 (2024) 002, arXiv:2404.10601 [cond-mat.str-el] . 12

  16. [24]

    Anomaly inflow for CSS and fractonic lattice models and dualities via cluster state measurement,

    T. Okuda, A. Parayil Mana, and H. Sukeno, “Anomaly inflow for CSS and fractonic lattice models and dualities via cluster state measurement,” SciPost Phys. 17 no. 4, (2024) 113 , arXiv:2405.15853 [quant-ph]

  17. [25]

    Anomaly inflow for dipole symmetry and higher form foliated field theories,

    H. Ebisu, M. Honda, and T. Nakanishi, “Anomaly inflow for dipole symmetry and higher form foliated field theories,” JHEP 09 (2024) 061, arXiv:2406.04919 [cond-mat.str-el]

  18. [26]

    Unpaired Majorana fermions in quantum wire s,

    A. Kitaev, “Unpaired Majorana fermions in quantum wire s,” Phys. Usp. 44 no. 10S, (2001) 131–136 , arXiv:cond-mat/0010440. 13

Pith tools

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