{"id":"fc64cb2c-22dc-4d91-a685-a9e0280c0219","arxiv_id":"2607.21525","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"On a square lattice with coexisting bond and chiral-flux charge order, s-wave superconductivity develops sequential topological phases with Chern numbers +2 and -2.","lead":"This paper computes what happens when a square-lattice superconductor also has two kinds of charge-order modulation: one that changes bond strengths and one that adds chiral loop currents. It finds two topological phases carrying Chern numbers +2 and -2, with quantized thermal Hall plateaus as a signature.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Predicted C=-2 phase and 'clear fingerprint' rely on fixed Δ=0.2 in a rigid CBO+CFP coexistence; no self-consistency check, and the paper's own tiny gap (ΔE≈0.005t) undercuts the experimental claim.","rationale":"The reader's weakest assumption identifies the same load-bearing concern: the topological phase diagram is computed for a rigid coexistence of CBO, CFP, and uniform s-wave pairing. The paper provides no energetic or self-consistent justification for this coexistence, and its own discussion of CFP-induced pairing suppression and the very small gap in the C=-2 phase make the concern concrete. If a self-consistent calculation were to show that Δ_self vanishes before the C=+2→C=-2 transition, the central claim would fail for realistic conditions. This is a correctness risk, not merely a scope limitation, because the abstract and summary assert a 'clear experimental fingerprint' and a 'pristine platform' for engineering topological superconductivity. The internal consistency of the numerical BdG calculation (Chern numbers matching edge spectra, band-resolved Chern numbers summing correctly) gives the model result independent mathematical support, but it does not address physical stability. I also note the analytic boundary λ_CFP=sqrt(Δ²+μ²/4) appears to disagree with the numerical transition at (μ=2, λCFP=0.5025) by roughly a factor of two; this is a secondary internal inconsistency that should be corrected, but it does not by itself invalidate the existence of the topological phases. Therefore the reader's CONDITIONAL verdict is appropriate; our read does not change it. The authors should be asked to add a self-consistency analysis or, failing that, to substantially temper the experimental fingerprint claim and clarify the model's artificial status.","tokens_in":14061,"tokens_out":18462,"duration_ms":188733,"concrete_test":"Perform self-consistent BdG calculations on the same square lattice with a microscopic attractive interaction (e.g., V Σ_i n_{i↑}n_{i↓}) that generates the s-wave pairing and with terms generating the CBO and CFP orders. Fix μ=0.2, λCBO=0.2, and solve for the equilibrium Δ_self and Chern number as λCFP is increased from 0 to 3. Check whether Δ_self remains finite at λCFP≈1.0 (C=+2) and λCFP≈2.8 (C=-2). If Δ_self collapses to zero before the second transition, the C=-2 phase is not a ground state of the model.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central sequence trivial→C=+2→C=-2 is obtained from the BdG Hamiltonian (Eqs. 2-5) with λCBO, λCFP, and Δ=0.2 treated as rigid, independent inputs. No microscopic term or free-energy minimization selects these orders or fixes their relative phase; the paper even states that CFP 'suppresses the superconducting pairing' and reports bulk gaps of ΔE=0.005 and 0.009 for λCFP=2.8 and 3.0 (Fig. 3d, bottom). This is not a small quantitative detail: if a self-consistent treatment were performed (e.g., from a Hubbard-type interaction), Δ would be a function of λCFP and could vanish before the C=-2 boundary is reached, eliminating that phase. The 'clear experimental fingerprint' is also undermined by the text's own statement that these small gaps suppress the quantized plateau even at very low T. Thus the central claim is load-bearing on the unexamined assumption that the coexisting CDW+SC state is energetically stable, not merely an artificial construction.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a square-lattice BdG Hamiltonian with uniform s-wave pairing coexisting with a 2x2 charge bond order (CBO) and chiral flux phase (CFP). It reports two topological superconducting phases with Chern numbers C=+2 and C=-2, separated by two phase transitions: a C=0 to C=+2 transition via quadratic band touching at Γ, and a C=+2 to C=-2 transition via Dirac cones along X-M. The authors compute band-resolved Chern numbers, edge spectra under partial open boundary conditions, and the thermal Hall conductivity, which at low T is quantized as κxy = Cκ0. The mechanism is attributed to normal-state band reconstruction from the combined real and imaginary hopping modulations.","tokens_in":14400,"tokens_out":5005,"duration_ms":50003,"significance":"If the predictions hold, the paper establishes a minimal square-lattice model in which the coexistence of real and imaginary bond modulations produces sequential Chern-number transitions in a superconductor, with chiral Majorana edge modes and quantized thermal Hall conductance as observable signatures. The central results are not obtained by fitting: Chern numbers and phase boundaries are computed directly from the stated Hamiltonian, and the thermal Hall plateaus follow from the standard identity κxy = Cκ0. The k·p symmetry argument explaining the λCBO-independence of the first transition is elegant. However, the manuscript provides no self-consistency check for the assumed coexistence of rigid CDW order and superconductivity, and the analytical boundary for the first transition appears inconsistent with the numerical phase diagram. The reported bulk gaps for the C=-2 phase are extremely small, which undermines the proposed experimental fingerprint. No code or data is provided to aid reproducibility.","major_comments":[{"comment":"The text claims the C=0→C=2 boundary is given by λCFP = sqrt(Δ² + μ²/4), independent of λCBO. This closed-form expression is not derived in Appendix A; the appendix only shows that the first-order effective Hamiltonian vanishes, not that this specific condition holds. Moreover, the formula disagrees with the numerics: for μ=2 and Δ=0.2 it gives λCFP≈1.02, while Fig. 2(a) and Fig. 5(a) place the transition at λCFP≈0.5025 (point A). Since the horizontal boundary in Fig. 2(b) is only shown for λCBO=0.2, this discrepancy needs to be resolved. Either derive the correct boundary condition or remove the unsubstantiated analytical expression.","section":"Sec. III (analytical solution at Γ) and Appendix A"},{"comment":"The C=-2 phase is obtained with Δ fixed at 0.2 while λCFP is varied. The paper itself notes that CFP suppresses superconducting pairing and reports gaps of only 0.005t and 0.009t at λCFP=2.8 and 3.0. A self-consistent treatment in which Δ depends on the CDW order could drive Δ to zero before the C=2→C=-2 boundary, eliminating that phase. The assumption that the rigid 2x2 CBO+CFP pattern coexists with uniform s-wave pairing across the entire parameter range is not energetically justified. The authors should at least estimate the mean-field free-energy competition or show that the conclusions survive for a range of Δ values that includes possible self-consistent values.","section":"Sec. II (Eqs. 3-4) and Sec. III (Fig. 3d)"},{"comment":"The abstract calls the quantized thermal Hall plateaus a 'clear experimental fingerprint'. This is contradicted by the text's own statement that the extremely small gaps for λCFP=2.8 and 3.0 'suppress the quantized plateau even at very low temperatures.' For the C=-2 phase, the proposed fingerprint is therefore not clear in any experimentally accessible temperature window. Please temper the claim or identify a parameter regime with a larger bulk gap.","section":"Abstract and Sec. III (Fig. 3d)"},{"comment":"The statement that 'CFP alone does not induce topology' is not demonstrated. The paper only shows that at λCBO=0 the gap closes at the would-be transition point and does not reopen; no phase diagram or Chern-number computation is presented for λCBO=0 across the (μ,λCFP) plane. Since CFP breaks time-reversal symmetry, the BdG bands could in principle acquire a nonzero Chern number without a transition at λCBO=0. This claim should be supported explicitly or qualified.","section":"Sec. III (first paragraph on CFP-only case) and Abstract"}],"minor_comments":[{"comment":"Typos: 'sufficient' should be 'sufficient'; 'effective' appears with non-ASCII ligature. Also, 'kagome' is inconsistently capitalized (e.g., 'Kagome' vs 'kagome').","section":"Throughout"},{"comment":"Reference [32] is cited as a pnictide example ('Na2Ti2Pn2O'), but the listed reference concerns the Kitaev material Na2Co2TeO6. This appears to be the wrong reference.","section":"Ref. [32]"},{"comment":"The symmetry operator P is defined only in Appendix A (Eq. A8) but is used in the main text before it is introduced. Move the definition to Sec. II or add a forward reference.","section":"Sec. III and Appendix A"},{"comment":"The notation ⟨αα'⟩ for the bond sum is not defined. Clarify that the sum is over bonds between sublattice sites within the enlarged 2x2 unit cell and specify the position vectors rα.","section":"Eq. (5)"},{"comment":"In Fig. 3(c), the inset resolving the edge modes should indicate the relevant kx range and the line colors corresponding to the two edges. The caption should also state the value of Ny used for the strip.","section":"Fig. 3"},{"comment":"The concept of a 'composite band' is introduced when discussing Figs. 2(c) and 2(d) but is not defined in Sec. II. Please define it explicitly, e.g., as a set of bands that remain degenerate at high-symmetry points and must be treated together for the Chern-number assignment.","section":"Sec. III (Composite band)"}],"recommendation":"major_revision","confidential_remarks":"The paper's central numerical results appear internally consistent, and the k·p symmetry argument is a nice contribution. However, the analytical boundary claim needs correction or removal, and the rigidity of the CDW+SC coexistence is a genuine physical concern. I would encourage the editor to request a self-consistency estimate or a clear parameter scan for Δ. Also, the manuscript relies heavily on self-citations (Refs. 44, 46, 47); the novelty relative to these prior square-lattice DDW studies should be stated more sharply. No code/data is provided, which is particularly relevant given the very small C=-2 gaps."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Punchline: this is a straightforward model study of a square-lattice s-wave superconductor with coexisting real bond order (CBO) and imaginary hopping (CFP). The new piece is the sequence trivial → C=+2 → C=-2 as the CFP strength grows, with the first boundary independent of CBO strength. The numerics look internally consistent: Chern numbers match the edge spectra, and the thermal Hall plateaus follow from the standard identity κxy=Cκ0. That part deserves credit.\n\nThe soft spots are real. The analytic boundary in Sec. III, λ_CFP=sqrt(Δ²+μ²/4), is stated without derivation, and it disagrees with Fig. 2(a) by about a factor of two for the point (μ=2, λCBO=0.2). The appendix does the k·p projection but never produces that formula. Second, the 'clear experimental fingerprint' in the abstract is undercut by the paper's own numbers: the bulk gaps of 0.005t and 0.009t in the C=-2 phase suppress the quantized thermal Hall plateau even at very low T, as the text admits. Third, the model keeps Δ=0.2 fixed while noting that CFP suppresses pairing; there is no self-consistency check on whether the coexisting CDW+SC state is energetically stable. That is a fair worry, though not fatal for a model study.\n\nThe stress-test concern about self-consistency is legitimate but speculative. The phase diagram is computed directly from the stated Hamiltonian, no parameter was fitted to get the central result, and the Chern numbers are what they are for this rigid-order model. The factor-of-two mismatch is more concrete and should be fixed in revision.\n\nWho this is for: people working on topological superconductivity from charge order, especially in kagome and square-lattice contexts. It is a moderate-significance model result, not a new framework. I would send it to a serious referee, but I would expect a revision that fixes the boundary formula, provides numerical details or artifacts, and tempers the experimental claim to match the small-gap region.","headline":"Clean model study with a genuinely new C=±2 mechanism, but the factor-of-two mismatch in the stated analytic boundary and the tiny bulk gaps in the C=-2 phase need attention before I'd take the experimental claim seriously.","tokens_in":14872,"tokens_out":2495,"would_cite":false,"duration_ms":26475,"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":"In a square-lattice superconductor, a real bond modulation plus a time-reversal-breaking chiral flux drives two successive topological phase transitions, ending in Chern numbers +2 and −2.","keywords":["topological superconductivity","Chern number","chiral flux phase","charge bond order","square lattice","thermal Hall conductivity","Majorana edge modes","charge density wave"],"falsifier":"Run a self-consistent mean-field calculation in which the pairing amplitude and the CBO/CFP order parameters are all determined by the same interaction, and check whether any region of parameter space has all three nonzero and falls inside the predicted C=+2 and C=−2 regions; alternatively, measure thermal Hall conductivity in a candidate square-lattice material with coexisting bond order and loop currents and look for the predicted plateau at κ0 or 2κ0.","tokens_in":13927,"feed_emoji":"🌀","tokens_out":6951,"duration_ms":57013,"temperature":0.7,"pith_summary":"The paper tries to establish that a square-lattice superconductor can be made topological not by either charge order alone, but only by their combination. When a real bond modulation (CBO) and an imaginary, time-reversal-breaking hopping modulation (CFP) coexist with uniform s-wave pairing, increasing the CFP strength drives the system through two topological phase transitions: from trivial to Chern number C=+2, then to C=−2. The first transition is a quadratic band touching at Γ whose critical flux is exactly sqrt(Δ²+μ²/4), independent of CBO strength; the second comes from Dirac cones along X–M. If correct, this makes the square lattice a clean, frustration-free host for chiral topological superconductivity, with chiral Majorana edge modes and quantized thermal Hall plateaus as observable signatures.","feed_headline":"Square-lattice superconductor reaches Chern +2 then −2 via charge order","feed_subtitle":"Real and imaginary bond modulations flip the Chern number twice; thermal Hall plateaus would reveal it.","key_machinery":"The engine is the 2×2 CDW pattern of Eq. (3) and Fig. 1: ξij=±1 on inter-/intra-cell bonds (CBO, a 2D analog of the Su–Schrieffer–Heeger modulation) and ηij=±1 along loop-current arrows (CFP, imaginary hopping). With uniform s-wave pairing Δ, the Bogoliubov–de Gennes Hamiltonian is in class D. The argument rides on the symmetry operator P=σx sx τ0: at Γ it commutes with the zeroth-order Hamiltonian but anticommutes with the linear k·p velocities, so the first-order effective Hamiltonian vanishes and the gap closes quadratically, forcing ΔC=±2. Along X–M that constraint is lifted, so linear Dirac cones govern the second transition. The persistent degeneracy at M means the second and third ban","core_discovery":"The central claim is that a square-lattice s-wave superconductor endowed with a composite 2×2 charge density wave — a real bond modulation (CBO) plus an imaginary hopping modulation (CFP) that breaks time-reversal symmetry — hosts two topologically nontrivial superconducting phases with Chern numbers C=+2 and C=−2. The paper maps the phase diagram in chemical potential, CBO strength, and CFP strength, and shows the two transitions are controlled by distinct symmetries: at Γ a P-protected quadratic band touching yields ΔC=±2 and a phase boundary independent of λCBO, while along X–M eight Dirac cones (four inequivalent) yield the second transition. It verifies the bulk–edge correspondence by v","pith_inferences":["If the C=−2 phase's tiny gap (≈0.005t) is generic rather than a parameter accident, the thermal Hall plateau for that phase may be practically unobservable unless pairing or band structure is engineered to widen it.","The mechanism likely generalizes: on any bipartite lattice whose imaginary flux vanishes at a symmetry point where a real bond order splits bands, the same composite-band Chern sequence could appear.","A self-consistent treatment in which Δ, λCBO, and λCFP all emerge from the same interaction would test whether the predicted Chern regions survive; the paper fixes them by hand.","Materials with coexisting bond order and loop currents, such as certain kagome superconductors, might exhibit a similar sequence, though the square-lattice version is cleaner to test."],"forward_implications":["The square lattice, previously thought immune to CFP-induced topology, becomes a platform for topological superconductivity once CBO and CFP coexist.","The first phase boundary is fixed analytically at λCFP = sqrt(Δ²+μ²/4) and is independent of λCBO, so this specific transition is a robust quantitative prediction.","Each nontrivial phase supports |C| chiral Majorana edge modes, directly visible in the strip-geometry spectra.","The thermal Hall conductivity at low temperature is quantized to C·κ0, so a plateau at κ0 or 2κ0 is a direct readout of the Chern number.","ARPES and STM can see the coexistence: a zero-energy DOS peak from CBO, a V-shaped DOS from CFP, and peak splitting when both are strong."],"fun_headline_variants":["Square-lattice superconductor flips Chern +2 to −2 via charge order","Two topological phases from real and imaginary bond modulations","Chern ±2 phases emerge from charge order in square-lattice superconductor","Thermal Hall plateaus track Chern flips in chiral CDW superconductor","Sequential topological superconductivity from composite charge order"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The paper fixes the s-wave pairing amplitude and the CDW pattern by hand and never checks whether the orders would actually coexist self-consistently; if the charge order suppresses or expels superconductivity, or if CBO and CFP do not share the same unit cell and phase relation, the predicted Chern phases may not exist.","fun_headline_variants_meta":{"raw":{"variants":["Square-lattice superconductor flips Chern +2 to −2 via charge order","Two topological phases from real and imaginary bond modulations","Chern ±2 phases emerge from charge order in square-lattice superconductor","Thermal Hall plateaus track Chern flips in chiral CDW superconductor","Sequential topological superconductivity from composite charge order"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000723,"raw_usage":{"total_tokens":3057,"prompt_tokens":698,"completion_tokens":2359,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":442,"completion_tokens_details":{"reasoning_tokens":2275}},"tokens_in":442,"tokens_out":2359,"duration_ms":15460,"temperature":1.0,"reasoning_tokens":2275,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T07:12:05.375730+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a self-consistent mean-field calculation in which the pairing amplitude and the CBO/CFP order parameters are all determined by the same interaction, and check whether any region of parameter space has all three nonzero and falls inside the predicted C=+2 and C=−2 regions; alternatively, measure thermal Hall conductivity in a candidate square-lattice material with coexisting bond order and loop currents and look for the predicted plateau at κ0 or 2κ0.","supporting_citations":[],"review_version":1}