{"id":"4ada99dc-088b-4710-abad-a61d572f9df5","arxiv_id":"2501.07935","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Adding a mass term to the ψ theory leaves a gapped bulk but generates a boundary flat band of Majorana zero modes, yielding ground state entropy proportional to the boundary area.","lead":"A free fermionic field theory with a broken subsystem symmetry stays gapped in the bulk, but develops a huge number of zero-energy modes when placed on a slab. The ground state degeneracy grows exponentially with the boundary area, producing residual entropy proportional to that area.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Case (ii) degeneracy is overstated: Eq. (3.27) has boundary modes only for 0 < k_y k_z/m < L_x, so the approximate zero-mode count is not the same as in Case (i).","rationale":"The reader's conditional verdict is appropriate. The central exact-zero-mode result for Case (i) appears internally consistent: the boundary condition (3.5) is a legitimate variational boundary condition, the zero-mode ansatz satisfies the equations of motion and boundary conditions, and the even number of Majorana zero modes avoids the 1+1-dimensional obstruction. The reader's weakest assumption about completeness of the boundary-condition classification is a reasonable caution, but it is not the most concrete defect in the manuscript. The sharper, checkable issue is the Case (ii) counting claim. Equation (3.27) has a positive solution only when 0 < k_y k_z / m < L_x, so the number of approximate boundary modes is at most half of the Case (i) count and depends on L_x and the ultraviolet cutoff. Since the paper's conclusion explicitly states that the degeneracy and gap are the same in both cases, this is a real overstatement. It does not destroy the area-law entropy claim for Case (i), so the verdict should remain conditional rather than reject; the authors should either correct the Case (ii) statements or restrict the conclusion to Case (i).","tokens_in":14418,"tokens_out":24746,"duration_ms":263487,"concrete_test":"Fix m, L_x, and the anti-periodic transverse momenta from Eq. (3.2). For every transverse pair (k_y, k_z), numerically solve tanh(ΛL_x) = Λ k_y k_z / m for Λ > 0 and count the solutions. Compare this count with the number of transverse momenta in Case (i). Repeat with m L_x smaller than the smallest positive k_y k_z; if the count drops to zero while Case (i) still has N modes, the claimed equality in Section 4 is false.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Section 3.2.2 and Section 4 claim that Case (ii) has the same vacuum degeneracy and gap as Case (i). This is contradicted by Eq. (3.27). Solving tanh(ΛL_x) = Λ k_y k_z / m for Λ > 0, a nonzero solution exists only when k_y k_z / m > 0 and k_y k_z < m L_x. For k_y k_z / m < 0 the right-hand side is negative while tanh(ΛL_x) is positive, so no boundary-localized mode exists; for positive product with k_y k_z > m L_x, the slope is too large and only Λ = 0 survives. Thus, even in the regime where the footnote's large-L_x assumption holds, Case (ii) has approximate zero modes for only the roughly half of transverse momenta with positive k_y k_z, not for all momenta as in Case (i). The degeneracy exponent is at most half of the Case (i) exponent, and can be reduced further if m L_x is not large enough. The gap statement is also wrong: when these approximate modes exist their frequency is exponentially small, not O(|m|), while in Case (i) the first excited state above the degenerate manifold has energy |m|. The exact zero modes of Case (i) and the area-law entropy are not invalidated by this, but the paper's summary for Case (ii) is an overstatement.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":14701,"tokens_out":6321,"duration_ms":62288,"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":[{"comment":"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).","section":"Sec. 3.2.2, Eq. (3.27)"},{"comment":"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.","section":"Sec. 3.1, Eq. (3.5)"},{"comment":"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.","section":"Sec. 3.2 (mode expansion)"}],"minor_comments":[{"comment":"The title contains an apparent typo: \"degenera cy\" should be \"degeneracy\".","section":"Title"},{"comment":"The phrase \"reminder term\" should be \"remainder term\".","section":"Sec. 3.2.2, Eq. (3.28)"},{"comment":"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.","section":"Sec. 3.2, counting N"},{"comment":"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.","section":"Sec. 3.2.2, discussion of (3.27)"}],"recommendation":"major_revision","confidential_remarks":"The core Case (i) result is interesting and appears correct, and the paper is within the scope of a hep-th journal. The main issue is the overstated Case (ii) claims, which are fixable by correcting the counting and gap statements. The missing completeness/self-adjointness discussion should also be addressed before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The one thing to know: the central result is solid. In Case (i) the paper exhibits exact Majorana zero modes for every transverse momentum pair, with normalizable wavefunctions, and gets residual entropy proportional to the boundary area. That is a genuinely new mechanism in free fermion field theory, distinct from gapless SSPT edge modes. The derivation is direct, self-contained, and has no fitted parameters; the variational argument for the boundary conditions and the explicit mode expansion are easy to follow. Credit where due: this is a real calculation, not a conjecture dressed up as one.\n\nThe soft spots are real but localized. The treatment of Case (ii) is overstated. Solving tanh(ΛL_x) = Λ k_y k_z / m shows a nonzero solution exists only for 0 < k_y k_z < m L_x. The paper's own Eq. (3.27) plus the positivity of tanh forces this, so the approximate zero modes exist only for about half of the transverse momenta, not all. The degeneracy exponent is therefore at most half that of Case (i), and the gap to the first excited state is exponentially small, not O(|m|). Section 4's claim that Case (ii) has the same degeneracy and gap as Case (i) is wrong and should be corrected. The paper also does not prove that the boundary conditions in (3.5) exhaust the self-adjoint possibilities, nor does it prove completeness of the mode expansion. That is a standard gap in physics papers, but it matters here because the claim of a huge degeneracy depends on having counted all modes. A short remark acknowledging this would be honest.\n\nThe Case (i) result survives all these objections. If I were refereeing, I would ask for revision: fix the Case (ii) counting, restate the gap correctly, and add a note on the boundary-condition completeness. This paper deserves a serious referee. It is aimed at readers working on fracton phases, subsystem symmetries, and boundary phenomena in free-field theories, and the Case (i) result is worth citing on its own.","headline":"The Case (i) boundary zero-mode result is real and clean; the paper oversells Case (ii).","tokens_in":15215,"tokens_out":4961,"would_cite":true,"duration_ms":48722,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["massive ψ theory","fermionic subsystem symmetry","boundary zero modes","vacuum degeneracy","residual entropy","area law","Majorana zero modes","gapped boundary"],"falsifier":"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.","tokens_in":14182,"feed_emoji":"⚛️","tokens_out":12059,"duration_ms":98227,"temperature":0.7,"pith_summary":"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.","feed_headline":"Boundary alone gives a fermion slab exponential vacuum degeneracy","feed_subtitle":"One zero mode per transverse momentum pair makes residual entropy scale with boundary area.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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|$."],"supporting_citations":[{"why":"Introduces the $\\psi$ theory and its fermionic subsystem symmetry, the model whose massive deformation this paper studies.","marker":"[6]"},{"why":"Provides the (1+1)-dimensional Majorana fermion analog and the odd-number zero-mode obstruction that makes the analogous boundary condition forbidden.","marker":"[26]"},{"why":"Earlier work on anomaly inflow for subsystem symmetries whose gapless boundary degrees of freedom this paper contrasts with its gapped boundary.","marker":"[21]"},{"why":"Previous construction of gapless edge modes in a related topologically massive theory, used as a comparison for why the present boundary degeneracy is different.","marker":"[22]"}],"fun_headline_variants":["Boundary-localized zero modes give exponential vacuum degeneracy","Zero modes at boundary yield area-proportional entropy","Boundary alone triggers exponential degeneracy in fermion slab","Majorana zero modes at boundary cause exponential degeneracy","Massive fermion slab with boundary gains exponential vacuum entropy"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Boundary-localized zero modes give exponential vacuum degeneracy","Zero modes at boundary yield area-proportional entropy","Boundary alone triggers exponential degeneracy in fermion slab","Majorana zero modes at boundary cause exponential degeneracy","Massive fermion slab with boundary gains exponential vacuum entropy"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001656,"raw_usage":{"total_tokens":6500,"prompt_tokens":793,"completion_tokens":5707,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":409,"completion_tokens_details":{"reasoning_tokens":5629}},"tokens_in":409,"tokens_out":5707,"duration_ms":35374,"temperature":1.0,"reasoning_tokens":5629,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:30:50.759195+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[],"review_version":1}