{"id":"7aa3277d-f6fa-43c3-ab01-4e6eb2111096","arxiv_id":"2608.26078","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"For gapless free fermions with a single massless Dirac cone, the modular commutator J satisfies 3J/pi = c_minus + 1/2, where c_minus is the chiral central charge, providing an entanglement-based measurement of the parity anomaly.","lead":"This paper computes a quantum-information quantity called the modular commutator for gapless free-fermion models, and finds it equals half a Chern number at a massless Dirac point. A generalist might read it because the result connects an entanglement measure to the famous parity anomaly, which is otherwise only seen in transport or partition-function arguments.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The proof's reflection-oddness of the low-energy modular-commutator contribution is not justified for a lattice with a single Dirac node, because spatial reflection maps the gapless cone to the gapped partner valley; half-quantization therefore rests on an unproven symmetry assumption.","rationale":"The paper's numerical observation (Fig. 2) is clear, and the robustness checks in Figs. 6-7 are valuable, so I do not dispute the empirical claim. The problem is that the analytical support for the central theorem is structurally incomplete: Eq. (4) is derived from a symmetry that the lattice model does not possess at the single-node critical point. Appendix D is honest about being heuristic, but the abstract's 'prove' overstates the status. The high-energy continuity flagged by the reader is a real assumption, but it is the less dangerous one: a gapped high-energy sector should be smooth in m. The more fundamental gap is the oddness of Phi_l, since the parity transformation used to establish it maps the gapless cone to the gapped partner valley and therefore does not act within the low-energy sector. I therefore agree only partially with the reader's weakest_assumption. My recommendation is unchanged conditional acceptance: the result is likely correct and interesting, but the proof should be labeled as a symmetry-based conjecture or supplemented by a symmetry-independent derivation (e.g., from the eta-invariant), as the conclusion already suggests.","tokens_in":10837,"tokens_out":16286,"duration_ms":181513,"concrete_test":"Perform the same modular-commutator calculation on a modified Haldane model in which an inversion-breaking term (e.g., a staggered next-nearest-neighbor hopping or a position-dependent phase) shifts the single gapless Dirac node away from the K/K' high-symmetry momenta to a generic momentum that is not fixed by any spatial reflection. Extrapolate 3J/pi at the critical point for the same tripartition geometry. If the thermodynamic limit differs from 1/2, the reflection-oddness assumption is exposed as the load-bearing step; if it remains 1/2, the half-quantization is more robust than the symmetry proof suggests, and the proof should be replaced by a symmetry-independent argument.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The proof of Eq. (4) hinges on the claim that the low-energy contribution Phi_l(m) is odd in the mass m, derived from an 'emergent reflection symmetry' psi(q_x,q_y) -> sigma_x psi(q_x,-q_y). In the lattice Haldane model at the single-node critical point, the gapless node sits at K, while the parity operation q_y -> -q_y maps the Brillouin-zone momentum K to K'=-K, where the partner node is gapped with a fixed mass of order 2M_c. The continuum parity transformation therefore does not preserve the low-energy subspace used to define Phi_l: it sends low-energy states around K to high-energy states around K'. Consequently, the oddness relation Phi_l(-m) = -Phi_l(m) is not a consequence of a microscopic or even an emergent low-energy symmetry of the model; it is an additional assumption. If this oddness fails, Eq. (6) only fixes the sum J(delta)+J(-delta), and gives no reason for J(0) to equal (C_+ + C_-)/2; the half-quantization would then depend on detailed high-energy physics rather than following from the gapped values. The paper itself labels the argument 'heuristic' in Sec. III.B, yet the abstract states 'We prove...' — this gap is the load-bearing step. The high-energy continuity assumption flagged by the reader is real, but the unverified reflection-oddness of Phi_l is the more fundamental problem.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the modular commutator J of free-fermion Haldane-type models tuned to gapless critical points. It reports numerically that for a single massless Dirac cone, 3J/π converges to a half-integer value (C_+ + C_-)/2, with power-law finite-size corrections, and to c_- + 1/2 when protected chiral edge modes coexist with the bulk node. The authors interpret the half-quantization as an information-theoretic manifestation of the parity anomaly, with the massive partner node playing the role of a Pauli-Villars regulator. They also report that quadratic nodes and Fermi surfaces do not exhibit robust half-quantization. The central derivation is a symmetry/heuristic argument based on separating high- and low-energy contributions.","tokens_in":11191,"tokens_out":10465,"duration_ms":99456,"significance":"If established, the half-quantization of the modular commutator would be a valuable extension of this entanglement diagnostic beyond gapped phases, and its connection to the parity anomaly is conceptually appealing. The numerical part is clean: exact diagonalization of free-fermion correlation matrices, checks of the entanglement area law and Markov-property violation, robustness under tripartition deformation and parameter tuning, and a useful comparison between Dirac, quadratic, and Fermi-surface gaplessness. However, the proof of the central formula is not rigorous as written, and the finite-size numerics alone do not constitute a derivation. The gap between the abstract's 'We prove' and the text's 'heuristic argument' is a substantive issue.","major_comments":[{"comment":"The proof of Eq. (4) hinges on the oddness relation Φ_l(-m) = -Φ_l(m), justified by an 'emergent reflection symmetry' ψ(q_x,q_y) -> σ_x ψ(q_x,-q_y) of the Dirac Hamiltonian. On the lattice, however, the gapless node sits at K while the parity-related point is K', where the partner node has a mass of order 2M_c. The proposed transformation does not preserve the low-energy subspace used to define Φ_l: it maps low-energy states near K to high-energy states near K'. Thus the oddness is an additional assumption, not a consequence of a microscopic or emergent low-energy symmetry. Without it, Eq. (6) only fixes J(δ)+J(-δ), and J(0) is not constrained to equal (C_+ + C_-)/2. The statement in the abstract that the half-quantization is proved is therefore not supported by the argument given.","section":"Sec. III.B, Eq. (6)"},{"comment":"The continuity of the high-energy contribution Φ_h(m) used in Eq. (6) is not established. The argument that Q_>(m) → 0 relies on assumed correlation functions G(r) ~ e^{-mr}/r^2 at low energy and G(r) ~ e^{-Λ r} at high energy, together with a separation of scales 1/R << Λ << Δ_UV. These scaling forms are heuristics; a large excitation energy does not by itself make the matrix elements in Eq. (D6) small, and the modular commutator is a nonlocal functional of the reduced density matrices rather than a simple integral of the single-particle Green function. The step 'Q_>(m) → 0' therefore needs a controlled derivation before Eq. (4) can be called a proof.","section":"Appendix D, Eqs. (D6)-(D9)"},{"comment":"The manuscript simultaneously says 'We prove the half-quantization' (abstract) and describes the derivation as a 'heuristic symmetry argument' (Sec. III.B), while Sec. VI lists a first-principle demonstration as an open problem. This mismatch is not cosmetic: the central claim of the paper is the derivation, not only the numerical observation. The authors should either supply a rigorous derivation of Φ_l oddness and Φ_h continuity, or explicitly reframe the result as a numerical conjecture supported by a heuristic argument.","section":"Abstract, Sec. III.B, Sec. VI"}],"minor_comments":[{"comment":"The notation Φ_l(0) is used both for the value at m = 0 (set to zero by oddness) and for the one-sided limits Φ_l(0±) = ±1/2; please clarify these distinct uses.","section":"Sec. III.B"},{"comment":"The power-law decay of the finite-size correction is stated but the fitted exponent and the system sizes used in the fit are not given; adding them would strengthen the claim.","section":"Fig. 2(d)"},{"comment":"The claim that the conditional mutual information 'saturates' at the critical point should be quantified by specifying the plateau value and the range of subsystem sizes over which saturation is observed.","section":"Appendix A, Fig. 5(b)"},{"comment":"The interlayer coupling term written as 't_c 1' should be typeset as t_c times the 2x2 identity matrix to avoid ambiguity.","section":"Eq. (7)"}],"recommendation":"major_revision","confidential_remarks":"The numerical results are likely correct and would be a useful contribution if framed appropriately, but the proof gap in the central claim and the overstatement in the abstract need to be addressed. I would welcome a revision that either completes the derivation or clearly labels the half-quantization as a numerically supported conjecture."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: new numerical result, useful framing, overclaimed proof. The half-quantized modular commutator at a single Dirac node, 3J/pi = C_- + 1/2, is cleanly established by exact diagonalization of free-fermion correlation matrices, with power-law convergence and robustness to tripartition geometry, velocity, and anisotropy. That is a real addition to the modular-commutator literature, which has so far been gapped-phase-only. The Pauli-Villars/parity-anomaly interpretation is a nice conceptual hook, and the extension to 3J/pi = c_- + 1/2 when chiral edge modes coexist is a natural and well-illustrated step. The loss of quantization for quadratic nodes and Fermi surfaces sharpens the claim that Dirac nodes are special.\n\nThe soft spots are real, and the stress-test note lands. The derivation in Sec. III.B and Appendix D does not close. The reflection-oddness of the low-energy contribution Phi_l(m) assumes parity maps the massless cone to itself. In the lattice model, parity maps K to K', where the partner cone is gapped with a nonzero mass, so the low-energy subspace is not preserved. That step needs a real justification, or the argument needs to be labeled a conjecture supported by numerics. The energy-scale separation and the continuity of the high-energy piece are also heuristic, as the paper itself admits in Sec. III.B. The abstract's \"We prove\" overshoots what the paper actually delivers. No code or data is shipped, which is a minor but avoidable omission for a numerical paper; the quadratic-node section is also thinner than the Dirac case.\n\nWhere the paper is honest, it is honest. The main text flags the symmetry argument as heuristic, and the appendices lay out the assumptions clearly enough for a careful reader to see the gap. That is the right scientific posture, but the abstract should match it.\n\nWho this is for: anyone working on modular commutators, entanglement bootstrap, gapless SPT phases, or the parity anomaly in lattice models. I would send it to a serious referee. With the proof claim softened or a genuine derivation supplied, this becomes a solid PRL-type result. As it stands, conditional acceptance is the right call.","headline":"Solid numerical discovery; the proof claim outruns the argument, but the paper deserves a serious referee.","tokens_in":11703,"tokens_out":1499,"would_cite":true,"duration_ms":17376,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"For a gapless 2D electron system with a single massless Dirac cone, the modular commutator — an entanglement probe of topological data — converges to a half-quantized value $3J/\\pi = 1/2$, revealing the parity anomaly from the…","keywords":["modular commutator","parity anomaly","massless Dirac cone","chiral central charge","half-quantization","entanglement","gapless topological phases","Haldane model"],"falsifier":"Compute $3J/\\pi$ at criticality for a lattice model with a single massless Dirac node but two different high-energy regularizations—say, Haldane next-nearest-neighbor hopping phase $\\phi$ replaced by $\\pi-\\phi$ or an added next-nearest-neighbor hopping term that keeps the node but changes the massive partner—and check whether the thermodynamic-limit value remains exactly $1/2$. If the value moves off $1/2$, the assumed continuity of the high-energy contribution is violated and the central claim is falsified.","tokens_in":10611,"feed_emoji":"⚛️","tokens_out":8810,"duration_ms":67653,"temperature":0.7,"pith_summary":"The paper asks what the modular commutator $J(A,B,C)=i\\langle[K_{AB},K_{BC}]\\rangle$ — an entanglement quantity that reads off the chiral central charge of a gapped two-dimensional state — measures when the bulk becomes gapless. For a system with a single massless Dirac cone it finds that $J$ remains sharp, converging to half the gapped value, $3J/\\pi=\\frac12$, with corrections that decay as a power law in subsystem size. The paper interprets this half-quantization as the entanglement shadow of the parity anomaly: the massive partner of the massless cone acts as a Pauli-Villars regulator that breaks parity, and $J$ detects that breaking. When a protected chiral edge mode coexists with the bulk node, the value becomes $3J/\\pi=c_-+\\frac12$, so $J$ measures total chirality. This matters because it extends a wavefunction-only diagnostic of topological data to gapless states, turning the parity anomaly into an information-theoretic observable.","feed_headline":"At one Dirac cone, an entanglement number locks to one-half","feed_subtitle":"A wavefunction-only probe reads off the parity anomaly in gapless 2D matter, converging to 1/2 with power-law corrections.","key_machinery":"The load-bearing object is the modular commutator $J(A,B,C)\\equiv i\\langle[K_{AB},K_{BC}]\\rangle$ with $K_X=-\\log\\rho_X$, evaluated on Gaussian free-fermion states through the single-particle modular Hamiltonian $K_X=\\log[(1-G_X)/G_X]$. The proof of half-quantization combines two pieces: an energy-scale decomposition of $J$ into high- and low-energy parts (Eq. (5)), and an emergent reflection symmetry of the massless Dirac cone, implemented by $q_y\\to -q_y$, $\\psi(q_x,q_y)\\to\\sigma_x\\psi(q_x,-q_y)$, which makes the low-energy contribution odd in the mass $m$. The high-energy part is assumed continuous across the gap-closing transition, so averaging $J(\\pm\\delta)$ isolates the jump and yields the half-integer. The gapped massive cone at the other valley acts as the parity-breaking Pauli-Villars regulator whose effect $J$ detects.","core_discovery":"The paper's central claim is that for a gapless 2D free-fermion ground state with one massless Dirac cone, the modular commutator is well defined and satisfies $3J/\\pi=(C_++C_-)/2=C_-+\\frac12$ in the thermodynamic limit, where $C_\\pm$ are the Chern numbers of the two gapped phases adjacent to the transition. The result is a half-integer, the finite-size corrections are power law rather than exponential, and the value is robust to deforming the tripartition, tuning the Dirac velocity, and making the cone anisotropic. The half-quantized value is carried by the gapped massive partner of the massless node, which plays the role of the Pauli-Villars regulator; this is why the paper reads $J$ as an information-theoretic measurement of the parity anomaly. When chiral edge modes coexist with the bulk Dirac node, the formula becomes $3J/\\pi=c_-+\\frac12$. Quadratic nodes and Fermi surfaces do not show this robust half-quantization: the quadratic-node value depends on curvature, and the Fermi-surface value depends on the tripartition.","pith_inferences":["Because the half-quantization proof uses only emergent reflection symmetry and an energy-scale separation, I would expect $3J/\\pi=\\frac12$ to survive in interacting Dirac critical points even though the Gaussian free-fermion numerics no longer apply; the paper leaves this as an open direction.","A direct modular-Hamiltonian derivation would likely connect the half-quantized value to the Atiyah-Patodi-Singer eta invariant of the Dirac operator, making $J$ a wavefunction-only realization of the parity anomaly rather than an analogy; the paper says this derivation remains open.","A practical diagnostic suggests itself: in numerical studies of candidate gapless Dirac materials, a converged $3J/\\pi=\\frac12$ would point to a single parity-odd Dirac node, while a curvature-dependent value would indicate a quadratic node and a tripartition-dependent value would indicate a Fermi surface."],"forward_implications":["At a Chern transition between $C_+$ and $C_-$, the modular commutator of the gapless ground state converges to $(C_++C_-)/2$, so entanglement alone can detect what happens to a topological invariant as the gap closes.","When protected chiral edge modes coexist with a bulk Dirac node, $3J/\\pi=c_-+\\frac12$; the modular commutator then measures the system's total chirality rather than separating bulk and edge contributions.","The related FSV wavefunction formula for Hall conductance also returns the half-quantized value at the Dirac node, so the half-integer is not specific to the modular commutator but appears in other single-wavefunction probes.","Quadratic nodes and Fermi surfaces destroy the robust quantization: the quadratic-node value depends on the node's curvature, while the Fermi-surface value depends on the shape of the tripartition, making $J$ a distinguishing diagnostic of the kind of gaplessness.","At the Dirac point the conditional mutual information saturates rather than decaying exponentially, showing that the area law of a single Dirac cone is a different type from the one assumed by the original entanglement-bootstrap derivation of the gapped formula."],"supporting_citations":[{"why":"Establishes that for gapped 2D states $3J/\\pi=c_-$; this is the relation the paper generalizes to gapless Dirac criticality.","marker":"[7]"},{"why":"Defines the modular commutator framework and its entanglement-bootstrap derivation, supplying the object $J$ being studied.","marker":"[8]"},{"why":"Shows how to obtain reduced density matrices for free fermions from the two-point correlation matrix, the technical tool for all numerical results.","marker":"[15]"},{"why":"Provides the single-particle modular Hamiltonian formula $K_X=\\log((1-G_X)/G_X)$ used to evaluate $J$ in Gaussian states.","marker":"[16]"},{"why":"Introduces the Haldane honeycomb model realizing a lattice Dirac fermion with a parity anomaly and a massive partner acting as the Pauli-Villars regulator.","marker":"[17]"},{"why":"Gives the Widom formula that the Fermi-surface entanglement entropy violates the area law, used to contrast the Dirac-node case.","marker":"[19]"},{"why":"Shows point-node 2D quantum critical points retain an area law, supporting use of the modular commutator at the Dirac critical point.","marker":"[22]"},{"why":"Provides the FSV wavefunction-only formula for Hall conductance which the paper shows is also half-quantized at the Dirac node.","marker":"[34]"}],"fun_headline_variants":["Half-quantized modular commutator exposes parity anomaly in gapless Dirac matter","At a massless Dirac node, modular commutator locks to one-half","Modular commutator hits 1/2 for gapless Dirac cone, revealing parity anomaly","Gapless Dirac state yields half-integer modular commutator, a parity anomaly probe","Half-quantization of modular commutator in gapless Dirac systems encodes parity anomaly"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes that an energy-scale cutoff cleanly separates the modular commutator's response into a low-energy part, whose correlation function is $G(r)\\sim e^{-mr}/r^2$, and a high-energy part, whose correlation function is $G(r)\\sim e^{-\\Lambda r}$ and whose contribution $Q_>(m)$ vanishes, so the high-energy contribution is continuous across the transition.","fun_headline_variants_meta":{"raw":{"variants":["Half-quantized modular commutator exposes parity anomaly in gapless Dirac matter","At a massless Dirac node, modular commutator locks to one-half","Modular commutator hits 1/2 for gapless Dirac cone, revealing parity anomaly","Gapless Dirac state yields half-integer modular commutator, a parity anomaly probe","Half-quantization of modular commutator in gapless Dirac systems encodes parity anomaly"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000855,"raw_usage":{"total_tokens":3828,"prompt_tokens":1175,"completion_tokens":2653,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":791,"completion_tokens_details":{"reasoning_tokens":2548}},"tokens_in":791,"tokens_out":2653,"duration_ms":17528,"temperature":1.0,"reasoning_tokens":2548,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-27T18:42:55.273360+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute $3J/\\pi$ at criticality for a lattice model with a single massless Dirac node but two different high-energy regularizations—say, Haldane next-nearest-neighbor hopping phase $\\phi$ replaced by $\\pi-\\phi$ or an added next-nearest-neighbor hopping term that keeps the node but changes the massive partner—and check whether the thermodynamic-limit value remains exactly $1/2$. If the value moves off $1/2$, the assumed continuity of the high-energy contribution is violated and the central claim is falsified.","supporting_citations":[{"cited_title":"Shi and I","cited_arxiv_id":null,"evidence_quote":"Establishes that for gapped 2D states $3J/\\pi=c_-$; this is the relation the paper generalizes to gapless Dirac criticality."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows how to obtain reduced density matrices for free fermions from the two-point correlation matrix, the technical tool for all numerical results."},{"cited_title":"Swingle, Entanglement entropy and the Fermi surface, Phys","cited_arxiv_id":null,"evidence_quote":"Shows point-node 2D quantum critical points retain an area law, supporting use of the modular commutator at the Dirac critical point."}],"review_version":1}