{"id":"4659a750-709f-4b48-a116-5e55a4bdbf30","arxiv_id":"2608.09415","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"The Chern number of a 2D p_x+ip_y superconductor equals -n(Gamma) + 2 n(X) - n(M), a combination of ARPES-measurable electron occupations at three high-symmetry points.","lead":"This paper proposes a way to read off the topological invariant of a p_x+ip_y superconductor from the electron density measured at three specific momenta using ARPES. If the formula holds in real materials, it would give experimentalists a relatively simple bulk probe of topological superconductivity.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (8) cannot distinguish p_x+ip_y from p_x-ip_y: the ARPES occupation at Γ, X, M depends only on |Δ_k|, while the Chern number changes sign under Δ→Δ^*, so the marker is not a self-contained measurement of the invariant.","rationale":"At T=0 the formula Eq. (8) is exact for the fixed chirality ζ=π/2: the Chern number of the two-band BdG model is 1/2(sgn ξ_Γ - 2 sgn ξ_X + sgn ξ_M), which equals Eq. (8). Thus the reader's concern about additivity of Berry curvature, while valid as a criticism of the paper's derivation, is not the weakest point of the central claim. The load-bearing issue is that the observable used in Eq. (8) is the electron occupation, which is even under Δ→Δ^*, whereas the Chern number is odd. A p_x+ip_y and a p_x-ip_y superconductor with identical band parameters have identical occupations at all k (since |Δ_k| is the same) and hence identical markers, but opposite Chern numbers. The paper's geometric construction in Eq. (4) contains sin ζ and would change sign, so the claimed equivalence between Eq. (4) and Eq. (8) is not chirality-covariant. This does not invalidate the model calculation for the specific ζ=π/2 Hamiltonian, but it undermines the advertised ARPES-based measurement of the topological invariant: the measured density cannot determine the sign of the invariant without additional phase-sensitive information. A concrete numerical test with ζ=-π/2 settles the question. The verdict remains conditional: the central formula is correct for the stated model, but the experimental interpretation needs a caveat about chirality, and the derivation should clarify the ζ-dependence.","tokens_in":20262,"tokens_out":23896,"duration_ms":244605,"concrete_test":"Take t=1, Δ=0.2, μ=-2t. (i) Compute the BCS occupations at Γ, X, M from Eq. (7) and evaluate Eq. (8); the marker is -1. (ii) Compute the Chern number of the same Hamiltonian but with Δ_k = 2Δ(sin k_x - i sin k_y) (ζ=-π/2) on a 200×200 lattice using the Fukui-Hatsugai-Suzuki discretized Berry curvature (Ref. [33]). Since |Δ_k| is unchanged, the occupations and hence Eq. (8) again give -1, while the true Chern number is +1. Repeating (ii) with ζ=+π/2 gives C=-1 and confirms the marker tracks the chirality only through the externally fixed ζ. This isolates the chirality-blindness of Eq. (8).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central formula Eq. (8), C = -<n_Γ> + 2<n_X> - <n_M>, is built from the BCS occupation <c†_k c_k> = 1/2(1 - ξ_k/E_k), which depends on ξ_k and |Δ_k| only; the phase of Δ_k enters nowhere. At the four gap nodes Γ, X, Y, M the occupations are identical for the pairing Δ_k = 2Δ(sin k_x + i sin k_y) (ζ=+π/2) and its complex conjugate Δ_k^* = 2Δ(sin k_x - i sin k_y) (ζ=-π/2). But the Chern number of the BdG Hamiltonian is odd under complex conjugation of the pairing: the d-vector winding reverses, so C changes sign in every topological phase. Therefore Eq. (8) assigns the same quantized value to the two chiralities, even though their Chern numbers are opposite. The paper's own geometric sum Eq. (4) contains the factor sin ζ and does change sign when ζ→-ζ, so Eq. (8) is only accidentally equivalent to Eq. (4) for one sign of the chirality. Because ARPES measures the phase-insensitive occupation, the proposed marker cannot by itself fix the sign of the topological invariant; the sign must be supplied by prior knowledge of the pairing chirality. This directly weakens the abstract's claim of an ARPES-measurable topological invariant for the p_x+ip_y superconductor.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a local topological marker for the square-lattice p_x + ip_y superconductor, defined as C = -⟨n_Γ⟩ + 2⟨n_X⟩ - ⟨n_M⟩, where each ⟨n_k⟩ is the BCS electron occupation at the time-reversal-invariant points Γ, X, and M. The authors show that this combination reproduces the Chern-number phase diagram, extend it to finite temperature and finite quasiparticle lifetime, and provide a real-space correlation-function representation. The occupation identity in Eq. (8) is exact within the BCS model, and the paper contains several nontrivial checks. However, the proposed marker depends only on |Δ_k| and ξ_k, so it cannot distinguish the two chiralities p_x + ip_y and p_x - ip_y; this limitation is not stated and conflicts with the abstract's claim of an ARPES-measurable topological invariant.","tokens_in":20611,"tokens_out":10132,"duration_ms":119530,"significance":"If the central claim were fully valid, the paper would provide a remarkably simple ARPES-accessible formula for the Chern number of a chiral superconductor from local momentum-resolved densities alone, together with an exact real-space expression and a finite-temperature stability statement. The exact identity in Eq. (8) is a genuine and useful result within the fixed-chirality model, and the finite-temperature and lifetime calculations are consistent. The main limitation is that the marker is phase-insensitive: it is identical for the two opposite chiralities, so it cannot by itself determine the sign of the topological invariant. The manuscript therefore needs a clear statement of what the marker does and does not measure, and an analytic justification of the disk additivity used in the geometric derivation.","major_comments":[{"comment":"The central marker is phase-insensitive and cannot distinguish the two chiralities of the pairing. Since ⟨c†_k c_k⟩ = 1/2(1 - ξ_k/E_k) depends only on ξ_k and |Δ_k|, the occupations at Γ, X, and M are identical for Δ_k = 2Δ(sin k_x + i sin k_y) and for its complex conjugate, whereas the BdG Chern number changes sign under ζ → -ζ because the d-vector winding reverses. Consequently, Eq. (8) assigns the same quantized value to p_x + ip_y and p_x - ip_y in every topological phase. The paper should state explicitly that the marker gives the chirality-fixed Chern number (with the chirality supplied by other information) or the absolute value |C|, rather than claiming that the local electron density alone measures the topological invariant. This is a load-bearing issue for the abstract's main claim, though it does not invalidate Eq. (8) as an exact identity within a fixed-chirality convention.","section":"Eqs. (7)-(8) and abstract"},{"comment":"The disk regularization is tuned to produce the desired half-quantized jumps. The radius is chosen as R² = m²/(2|Δ|² sin ζ), which is not positive definite when sin ζ < 0, so the geometric derivation is not defined for the opposite chirality. In addition, the additivity of the integrated Berry curvatures around X and Y is asserted and checked numerically, but no analytic proof is given. Because Eq. (8) is exact without any Berry-curvature integration, the disk construction is not needed for the main formula; however, the text presents Eqs. (4), (8), (13), and (14) as equivalent definitions, and that equivalence is not established for general ζ. Please supply an analytic derivation of the additivity or explicitly label the disk construction as a heuristic illustration rather than a proof.","section":"Sec. I, Eqs. (25)-(26) and Eq. (4)"}],"minor_comments":[{"comment":"There are typographical errors, including 'superfuid' for 'superfluid' and the inconsistent rendering of 'p_x + ip_y' in the title and opening lines.","section":"Main text after Eq. (2)"},{"comment":"The symbol ΔS is introduced in Eq. (4) without a definition, and it is later equated to C in Eq. (8); please define the notation clearly before first use.","section":"Eq. (4)"},{"comment":"The caption states that the marker is shown for different temperatures but does not identify which curve corresponds to which temperature nor the values of Δ and the line color legend; please add these details.","section":"Fig. 2"},{"comment":"The sentence 'The factor 1/2 ensures that the jump will always be an integer' is not an explanation; the factor follows from the chosen normalization of the integrated Berry curvature, and this should be stated explicitly.","section":"Sec. I after Eq. (26)"},{"comment":"The 'topological superfluid density' n_s = sin ζ |Δ|²/m² is introduced by analogy with the Karplus-Luttinger velocity; while interesting, this identification is not derived from a response function and should be clearly labeled as a phenomenological analogy.","section":"Sec. II, Eq. (19)"}],"recommendation":"major_revision","confidential_remarks":"The exact occupation identity is a solid core and deserves publication after the chirality limitation is addressed. The paper's novelty relative to Ref. [14] and related work would be clearer if the authors explicitly separated the exact identity from the heuristic geometric justification. A statement about the sign ambiguity of the ARPES-accessible marker is essential before acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I'll get straight to it. The genuinely new thing here is the exact identity C = -n_Γ + 2 n_X - n_M for the square-lattice p_x+ip_y superconductor: the BCS occupation at those three time-reversal-invariant momenta reproduces the Chern number, with a clean finite-temperature generalization. The Sec. V analysis of lifetime and frequency-window robustness is careful and correct. The real-space correlation-function form of the marker is a nice extra. That part of the paper is solid.\n\nThe soft spot that matters is the chirality blindness. The occupations depend only on |Δ_k|, so Eq. (8) assigns the same value to p_x-ip_y, whose Chern number is opposite. The paper's own geometric sum Eq. (4) has the factor sin ζ and flips sign under ζ→-ζ, but the density formula has no such factor. So the marker is the Chern number only if you already know the pairing chirality. ARPES alone cannot measure the invariant's sign. The abstract's phrasing 'may be measured via ARPES' is stronger than what the relation delivers. The authors should either state the marker is the invariant up to the chirality sign or find a second observable that fixes it.\n\nTwo smaller issues. The regularization radius in Eqs. (25)-(26) is tuned to get the desired 1/2 factors, and the additivity of the local Berry curvature is checked numerically, not derived. Because Eq. (8) is an exact algebraic identity, this does not sink the main result; it just makes the geometric narrative less convincing. The 'topological superfluid density' in Eq. (19) is an analogy, not a derivation, and reads as speculative. No code or data is shipped, though the numerics are simple.\n\nThe paper is for people working on local topological markers and on ARPES probes of paired states; they will find the exact identity useful. It deserves a serious referee. The sign ambiguity is a revision-worthy flaw, not a desk-reject. I would send it to review, asking the authors to address the chirality caveat head-on.","headline":"A clean exact density-to-Chern marker for p_x+ip_y, but it is blind to the sign of the pairing chirality, so ARPES alone does not measure the invariant.","tokens_in":21133,"tokens_out":5993,"would_cite":true,"duration_ms":60620,"reading_group":"yes","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 the Chern number of the square-lattice $p_x+ip_y$ superconductor equals the local marker $C=-\\langle n_{\\Gamma}\\rangle+2\\langle n_{X}\\rangle-\\langle n_{M}\\rangle$, a combination of electron densities measurable with…","keywords":["p_x+ip_y superconductor","Chern number","local topological marker","electron density","ARPES","Berry curvature","topological phase transition","Van Hove singularity"],"falsifier":"Compute the full Brillouin-zone Chern number and the marker $-\\langle n_{\\Gamma}\\rangle+2\\langle n_{X}\\rangle-\\langle n_{M}\\rangle$ in a model with an added next-nearest-neighbor hopping $t'$ large enough to shift the gapless points off the four special momenta; any mismatch between the two values away from the topological transitions would falsify the central claim. An ARPES measurement of the three densities in a candidate chiral superconductor that returns a non-integer combination would likewise rule out the marker as a bulk invariant.","tokens_in":20079,"feed_emoji":"🧲","tokens_out":9796,"duration_ms":94232,"temperature":0.7,"pith_summary":"The paper proposes a way to extract the topological invariant of a two-dimensional $p_x+ip_y$ superconductor from local, momentum-resolved electron densities instead of from edge transport. Its central formula is $C=-\\langle n_{\\Gamma}\\rangle+2\\langle n_{X}\\rangle-\\langle n_{M}\\rangle$, where each $\\langle n_{\\mathbf{k}}\\rangle$ is the occupation of a high-symmetry electron state, a quantity accessible to angle-resolved photoemission spectroscopy (ARPES). The paper argues that this combination equals the Chern number, takes the quantized values $-1,0,+1$, and remains a plateau at finite temperature within the BCS model. It derives the marker from two one-dimensional invariants on the Bloch sphere, relates it to a Pfaffian-like sign formula, and links it to a topological superfluid density. If the claim is right, a genuinely global property of the superconducting state becomes a small set of local density measurements.","feed_headline":"Chern number of p_x+ip_y superconductor equals three ARPES densities","feed_subtitle":"Local electron occupations at Γ, X, M reproduce the quantized bulk invariant and survive finite temperature.","key_machinery":"The central object is the local topological marker obtained by adding Berry curvatures integrated over small disks centered at the four special points $\\Gamma$, $X$, $Y$, $M$; near a transition the disk radius is set by $R^2=m^2/(2|\\Delta|^2\\sin\\zeta)$ and each disk contributes $\\pm\\frac12\\,\\mathrm{sgn}(\\cdot)$ to the sum. The paper then rewrites the marker in the Anderson pseudo-spin representation, where $\\langle S_z\\rangle=2\\langle n_{\\mathbf{k}}\\rangle-1$, and shows that the sum is exactly the sum of two 1D invariants $C_x$ and $C_y$ obtained by walking on the Bloch sphere from $\\Gamma$ to $X$ and from $X$ to $M$. That rewriting is what converts the geometrical invariant into a density formula, and it also gives the Pfaffian-type expression $C=|C_y|-|C_x|$ as a difference of products of signs of the kinetic energies.","core_discovery":"On the paper's own terms, the central claim is that every value of the chemical potential in the square-lattice model with pairing $\\Delta_{\\mathbf{k}}=2\\Delta(\\sin k_x+e^{i\\zeta}\\sin k_y)$ is captured by the local marker $C=-\\langle n_{\\Gamma}\\rangle+2\\langle n_{X}\\rangle-\\langle n_{M}\\rangle$, with $\\langle n_{Y}\\rangle=\\langle n_{X}\\rangle$ by symmetry and $\\langle n_{\\mathbf{k}}\\rangle=\\int d\\omega\\,n_F(\\omega)A(\\mathbf{k},\\omega)$ the electron occupation obtained from the spectral function. At the points $\\Gamma$, $X$, $Y$, and $M$, where the pairing function vanishes, the occupation is $\\langle n_{\\mathbf{k}}\\rangle=\\frac12-\\frac12\\,\\mathrm{sgn}(\\xi_{\\mathbf{k}})\\tanh(\\beta E_{\\mathbf{k}}/2)$, so the marker reduces to signs of the kinetic energy and changes at $\\mu=-4t$, $\\mu=0$, and $\\mu=+4t$. The paper identifies $C$ with the sum $C_x+C_y$ of two 1D $\\mathbb{Z}$ invariants defined on the Bloch sphere, and shows that at half-filling the transition is a protected semimetal characterized by $\\langle n_{X}\\rangle=1/2$.","pith_inferences":["Beyond the paper: if the marker remains equal to the Chern number under perturbations that keep the four special points gapless, the same three-density combination could serve as a fast bulk diagnostic for chiral superconductivity in candidate materials.","Beyond the paper: the construction suggests a general recipe for other lattice models and pairing symmetries: whenever the pairing vanishes at a small set of high-symmetry points, the Chern number may be read from a weighted sum of occupations at those points, with weights fixed by kinetic-energy signs.","Beyond the paper: the protected value $\\langle n_X\\rangle=1/2$ at half-filling gives a momentum-space analogue of a half-filled Landau level, and it would be worth testing whether this semimetal description survives interactions beyond BCS mean field."],"forward_implications":["Within BCS theory, the topological invariant of a $p_x+ip_y$ superconductor can in principle be obtained from bulk ARPES measurements of electron occupation at $\\Gamma$, $X$, and $M$, without measuring edge transport.","The plateau in $C$ survives finite temperature; the practical upper limit is the superconducting critical temperature, not the thermal smearing of the density step.","A finite experimental resolution in momentum space does not destroy the marker: integrating the density over a small disk around $X$ or $Y$ preserves the quantized value while turning the sharp derivative into a logarithmic Van Hove profile.","The real-space version of the marker, written with correlation functions $\\langle c_i^\\dagger c_j\\rangle$, is expressed in a form the paper states remains usable when weak disorder is present."],"supporting_citations":[{"why":"Introduces the chiral p+ip paired state and the weak- and strong-paired phases that the paper uses as its physical context and reference model.","marker":"[9]"},{"why":"Relates Green's function properties at zero frequency to topological phase transitions, supporting the use of local responses at gap-closing points.","marker":"[13]"},{"why":"Supplies the sphere-and-poles method that rewrites global topological invariants as Berry phases at two points, which the paper generalizes from Chern insulators to two-dimensional superconductors.","marker":"[14]"},{"why":"Defines the one-dimensional Kitaev p-wave chain whose $\\mathbb{Z}_2$/Pfaffian invariants and Majorana edge modes underlie the 1D invariants $C_x$ and $C_y$.","marker":"[26]"},{"why":"Provides the Anderson pseudo-spin representation of BCS states, the object through which local electron density is connected to the topological marker.","marker":"[29]"},{"why":"Gives the numerical method used to compute Berry curvature and Chern numbers on the discretized Brillouin zone, the benchmark against which Eq. (8) is checked.","marker":"[33]"},{"why":"Connects the many-body Chern number to Luttinger's theorem through Streda's formula, motivating the search for a density-based topological marker.","marker":"[36]"}],"fun_headline_variants":["Three ARPES densities yield Chern number in p_x+ip_y superconductor","Chern number from Γ,X,M electron densities","Local marker: Chern number from three spectral densities","p_x+ip_y Chern number encoded in three ARPES occupations","Topological invariant from electron density at symmetry points"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the integrated Berry curvatures on small disks around $\\Gamma$, $X$, $Y$, and $M$ add up to the global Chern number; the paper checks this numerically for every $\\mu$ but gives no analytic proof of the additivity, so if the additivity failed away from the transitions the marker would deviate from the true invariant.","fun_headline_variants_meta":{"raw":{"variants":["Three ARPES densities yield Chern number in p_x+ip_y superconductor","Chern number from Γ,X,M electron densities","Local marker: Chern number from three spectral densities","p_x+ip_y Chern number encoded in three ARPES occupations","Topological invariant from electron density at symmetry points"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000841,"raw_usage":{"total_tokens":3707,"prompt_tokens":1028,"completion_tokens":2679,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":644,"completion_tokens_details":{"reasoning_tokens":2595}},"tokens_in":644,"tokens_out":2679,"duration_ms":21117,"temperature":1.0,"reasoning_tokens":2595,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T17:44:27.437871+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the full Brillouin-zone Chern number and the marker $-\\langle n_{\\Gamma}\\rangle+2\\langle n_{X}\\rangle-\\langle n_{M}\\rangle$ in a model with an added next-nearest-neighbor hopping $t'$ large enough to shift the gapless points off the four special momenta; any mismatch between the two values away from the topological transitions would falsify the central claim. An ARPES measurement of the three densities in a candidate chiral superconductor that returns a non-integer combination would likewise rule out the marker as a bulk invariant.","supporting_citations":[{"cited_title":"Le Hur, Interacting topological quantum aspects with light and geometrical functions, Physics Reports1104, pp.1-42 (2025)","cited_arxiv_id":null,"evidence_quote":"Supplies the sphere-and-poles method that rewrites global topological invariants as Berry phases at two points, which the paper generalizes from Chern insulators to two-dimensional superconductors."},{"cited_title":"Kitaev, Unpaired Majorana fermions in quantum wires, Uspekhi Fizicheskikh Nauk, (2001), Volume171, supplement 10, Pages 131-136","cited_arxiv_id":null,"evidence_quote":"Defines the one-dimensional Kitaev p-wave chain whose $\\mathbb{Z}_2$/Pfaffian invariants and Majorana edge modes underlie the 1D invariants $C_x$ and $C_y$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Anderson pseudo-spin representation of BCS states, the object through which local electron density is connected to the topological marker."},{"cited_title":"Fukui, Y","cited_arxiv_id":null,"evidence_quote":"Gives the numerical method used to compute Berry curvature and Chern numbers on the discretized Brillouin zone, the benchmark against which Eq. (8) is checked."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Connects the many-body Chern number to Luttinger's theorem through Streda's formula, motivating the search for a density-based topological marker."}],"review_version":1}