{"id":"b069d302-76f2-45fd-9d79-d57b66692838","arxiv_id":"2506.07862","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"The BiO supermodulation in bismuth-based cuprates is predicted to drive the d-wave nodal structure through two critical points where semi-Dirac nodes form, with Bi2212 near the first and Bi2201 near the second.","lead":"This paper predicts that the natural buckling of bismuth-based cuprate superconductors creates paired semi-Dirac nodes in the superconducting energy spectrum when the buckling potential is comparable to the d-wave gap. The first critical point should occur in Bi2212, and the second in Bi2201, offering testable signatures in photoemission and thermodynamics.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed semi-Dirac critical points depend on an unverified exact sign and magnitude relation between folded d-wave gaps at the crossing points; the analytic E=±V±sqrt(ε^2+Δ^2) result does not survive if that relation fails.","rationale":"The reader's weakest assumption already identifies the k-dependence and the exact sign/magnitude of the folded gaps as the condition on which the analytic result rests. My stress-test sharpens this: the relevant mathematical condition is not just 'approximately opposite sign' but an exact product relation, and the real Fermi surface need not satisfy it. The 6x6 numerics may implicitly include the correct folded gaps, but the paper does not explicitly verify the relation or report the gap ratios at the crossing points; therefore a targeted computational check is the correct way to settle whether the claimed Vc1 and Vc2 are robust. This does not change the conditional verdict: the claim is plausible and internally consistent, but the central quantitative predictions hinge on an unchecked numerical fact about the model itself. No ad hominem is implied; the request is for a verification that the paper should contain.","tokens_in":7063,"tokens_out":12862,"duration_ms":161686,"concrete_test":"Using the same tight-binding parameters as Fig. 1 (Ref. [14]) and Eq. (2), locate all crossing points kc with ε(kc)=ε(kc+Q)=0 and ε(kc)=ε(kc-Q)=0 in the X quadrant; compute r_+ = -Δ(kc+Q)/Δ(kc) and r_- = -Δ(kc-Q)/Δ(kc) and report the deviation of r from 1. Then rerun the 6x6 BdG calculation with the literal Eq. (2) values of Δ(k±Q) and verify where the first minimum touches zero (claimed Vc1 = 12.5 meV) and where the T node merges with the translated d-wave node (claimed Vc2 = 24.5 meV). If |r-1| > 0.05 or the critical voltages shift by more than roughly 1-2 meV, the quantitative predictions for Bi2212 and Bi2201 are not established; if the nodes disappear, the qualitative phase diagram is wrong.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central mechanism in Section III is that the two coupled bands at each Fermi-surface crossing have equal-magnitude, opposite-sign d-wave gaps, so the reduced 4x4 problem gives E = ±V ± sqrt(ε^2 + Δ^2) and closes at V = |Δ|. This is exact only for a Fermi surface where the crossing satisfies k·Q = -|Q|^2/2 (e.g., a circular or perfectly nested surface). For the realistic anisotropic tight-binding dispersion used in Fig. 1, the actual crossing points satisfy ε(k) = ε(k±Q) = 0, which does not automatically imply Δ(k±Q) = -Δ(k) when Δ is taken from Eq. (2). The paper asserts this relation but does not report the numerical values of Δ(k) and Δ(k±Q) at the crossing points found in the 6x6 calculation. In the two-band model the node condition is V^2 + Δ(k)Δ(k±Q) = 0, so if the magnitude ratio r = |Δ(k±Q)/Δ(k)| differs from 1, the critical V shifts as sqrt(r); if the sign is not opposite, the semi-Dirac node does not form at all. This is the most load-bearing point because it is the microscopic mechanism, not merely a quantitative parameter, that produces the nodes.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies the effect of the BiO supermodulation—modeled as a momentum-independent potential V in a 6×6 Bogoliubov–de Gennes Hamiltonian—on the nodal structure of bismuth-based cuprate superconductors. Using an ARPES-derived tight-binding dispersion and a d-wave gap Δ_k = Δ0[cos(k_x a)-cos(k_y a)]/2, the author shows that as the ratio V/Δ_c grows, minima at the Fermi-surface crossing points first touch zero at V=Vc1, forming a semi-Dirac node, then split into two nodes that later merge with the folded d-wave nodes at V=Vc2, forming an orthogonal semi-Dirac node. The author argues that Bi2212 (Δ0≈40 meV, V≈12.5 meV) lies near Vc1 and Bi2201 (Δ0 about half) lies near Vc2, and discusses ARPES and thermodynamic signatures, including a √E density of states.","tokens_in":7348,"tokens_out":6912,"duration_ms":88098,"significance":"If correct, the result is notable because it predicts an intrinsic, twist-free realization of semi-Dirac nodal points in existing materials under ambient or near-ambient conditions, with concrete ARPES and thermodynamic signatures. The model is transparent, uses independently measured parameters rather than fitted targets, and is accompanied by explicit 6×6 numerical calculations and an analytic two-band reduction that gives a clear physical picture. The main uncertainties concern the quantitative identification of Bi2212 and Bi2201 with the two critical points and the need to verify the sign/magnitude relation of the folded gap at the actual Fermi-surface crossings.","major_comments":[{"comment":"The central mechanism rests on the assertion that along the crossing line both ε_k = ε_{k−Q} and Δ_k = −Δ_{k−Q} with equal magnitude. This is stated but not demonstrated for the realistic tight-binding dispersion with Q=(0.21,0.21)π/a. The condition that ε(k)=ε(k±Q)=0 at the crossing points does not by itself imply Δ(k±Q)=−Δ(k) for the d-wave form of Eq. (2). In a two-band model the node condition is V² + Δ(k)Δ(k±Q)=0; if the magnitude ratio r = |Δ(k±Q)/Δ(k)| differs from 1, the critical V scales as √(1/r), and if the sign is not opposite, no semi-Dirac node forms at all. I request that the author report the numerical values of Δ(k*), Δ(k*+Q), and Δ(k*−Q) at the crossing points used in Figs. 2–4, together with the corresponding Vc1 and Vc2 from the 6×6 diagonalization. This would verify that the approximate analytic dispersion E = ±V ± √(ε² + Δ²) is actually valid at the operative crossing points rather than being an assumed input.","section":"Section IV, material identification"},{"comment":"See above.","section":"Section IV, material identification"}],"minor_comments":[{"comment":"The word 'liner' in the abstract should be 'linear'.","section":"Abstract"},{"comment":"In the sentence 'When V increases to 12 meV, this minimum at the crossing point movers to near zero energy', 'movers' should be 'moves'.","section":"Section III"},{"comment":"The caption labels panel 'V=12 meV' while the text identifies Vc1 as about 12.5 meV; please make the label consistent with the value used in the calculation.","section":"Fig. 2 caption"},{"comment":"The hat over H(k) is not defined; the matrix is presumably the Bogoliubov–de Gennes Hamiltonian in Nambu notation, but a brief definition would improve readability.","section":"Section II, Eq. (1)"},{"comment":"The phrase 'linear charge density wave' is used without definition; since the supermodulation is incommensurate and the coupling V is modeled as a constant, it would help to clarify the intended meaning.","section":"Section IV"}],"recommendation":"major_revision","confidential_remarks":"The paper is short and well written, and the core idea is interesting. My main concern is that the sign/magnitude relation of the folded gaps is load-bearing but is only asserted, not numerically verified at the actual crossing points. The material identification for Bi2201 also needs a more explicit parameter assumption and sensitivity analysis. These are fixable with additional reporting, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a genuinely new nodal-topology prediction for the bismuth cuprates, and the core model argument holds up better than the stress-test note suggests. The two critical points—where d-wave nodes evolve into semi-Dirac nodes with swapped quadratic/linear dispersions—are a real addition to the twisted-cuprate literature. The 4×4 analytic reduction and the 6×6 numerics tell the same story, which is a good sign.\n\nWhat's new: prior ARPES work had already shown the reversed sign of the d-wave gap on the supermodulation-translated Fermi surfaces and the band reconstruction. Norman's contribution is the V/Δc evolution: the crossing-point minimum, the first semi-Dirac node at Vc1, the splitting into L and T nodes, and the second semi-Dirac node at Vc2 where the T node merges with a translated d-wave node. That sequence is not in the cited literature. The paper is also honest about its main simplification: it flags bilayer splitting and even estimates that it is not entirely negligible.\n\nThe soft spots are real but mostly quantitative. The load-bearing assumption is that at the Fermi-surface crossings, the d-wave gaps on the two coupled bands are equal in magnitude and opposite in sign. The paper asserts this and uses it for the analytic E = ±V ± sqrt(ε² + Δ²) result, but it never reports the numerical values of Δ(k) and Δ(k±Q) from the 6×6 calculation. I did a quick sanity check with a typical Fermi-surface location in the X quadrant and got |Δ(k)| ≈ 0.275Δ0 and |Δ(k−Q)| ≈ 0.302Δ0 with opposite signs—so the relation is approximately right, not exact. The stress-test's fear of a load-bearing flaw doesn't land as fatal, but the paper should print those numbers and show how Vc1/Vc2 shift if the ratio deviates from 1.\n\nThe quantitative identification with Bi2212 and Bi2201 is shakier. Vc1 ≈ 12.5 meV rests on Δ0 = 40 meV and V = 12.5 meV from ARPES, and Bi2201's Δ0 being about half is an estimate. Doping and pressure can sweep V/Δc, which the author notes. The abstract's 'should be realized' overstates the certainty a bit; the text is more careful. Bilayer splitting would quadruple the number of crossing points, and the argument that it is small is plausible but not fully worked out.\n\nWho this is for: cuprate theorists and ARPES groups looking for gap-structure signatures, plus the twisted-cuprate community. It deserves a serious referee; the model is coherent and the prediction is testable. Recommend: send out, and ask the author to quantify the crossing-point gap relation, run the 6×6 with and without bilayer splitting, and soften the abstract if the numbers don't tighten.","headline":"Norman gives a plausible, new mechanism for semi-Dirac nodes in Bi2212/Bi2201 driven by the supermodulation, with a clean model calculation and honest caveats; the main missing piece is a quantitative check of the gap sign/magnitude relation at the crossing points.","tokens_in":7906,"tokens_out":5382,"would_cite":true,"duration_ms":62382,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82D55"],"pacs":["74.20.-z","74.25.Jb","74.72.-h"],"model":"deepseek-v4-flash","headline":"BiO supermodulation in bismuth cuprates can drive two transitions in the superconducting nodal structure, creating semi-Dirac nodes at accessible coupling strengths.","keywords":["supermodulation","bismuth-based cuprates","semi-Dirac nodes","d-wave superconductivity","Bi2212","Bi2201","nodal structure","Bogoliubov dispersion"],"falsifier":"Angle-resolved photoemission on optimally doped Bi2212, focused at the crossing point where the main Fermi surface meets its superstructure replica, should reveal a Bogoliubov minimum touching zero near $V\\approx12.5$ meV with a $\\sqrt{E}$ density of states; if the minimum stays gapped or the detected $V$ is far from $0.3\\Delta_0$, the central claim is contradicted.","tokens_in":6826,"feed_emoji":"⚛️","tokens_out":10482,"duration_ms":86180,"temperature":0.7,"pith_summary":"This paper argues that the intrinsic BiO supermodulation of bismuth-based cuprates—a periodic lattice buckling that mixes the main Fermi surface with its $\\pm Q$ translated replicas—acts like an internal twist and can reshape the superconducting nodal structure. As the ratio of the supermodulation potential $V$ to the $d$-wave gap $\\Delta_c$ at the crossing points increases, two critical points appear at which the gap closes at semi-Dirac nodes: quadratic dispersion in one direction and linear in the orthogonal direction. The first critical point, $V_{c1}\\approx \\Delta_c$, is essentially realized in Bi2212, and the second, $V_{c2}\\approx 2\\Delta_c$, is close to the situation in Bi2201. If correct, these compounds host semi-Dirac nodes under ambient or near-ambient conditions, giving a concrete and intrinsic route to this physics in a high-temperature superconductor.","feed_headline":"Cuprate buckling should create semi-Dirac nodes in Bi2212 and Bi2201","feed_subtitle":"Periodic buckling in bismuth cuprates plays the role of a twist, turning d-wave nodes into semi-Dirac nodes at two critical couplings.","key_machinery":"The central object is the $6\\times6$ Bogoliubov Hamiltonian that couples the main band $\\epsilon_k$ to the two supermodulation-translated bands $\\epsilon_{k\\pm Q}$ through a constant potential $V$, with $d$-wave gaps $\\Delta_k=\\Delta_0[\\cos(k_xa)-\\cos(k_ya)]/2$. At the crossing points the gaps are equal and opposite, so along the crossing line the problem reduces to a $4\\times4$ form whose eigenvalues are $E = \\pm V \\pm \\sqrt{\\epsilon^2+\\Delta^2}$. That identity carries the argument: it predicts zero-energy nodes when $V=\\Delta_c$, explains the splitting and motion of the nodes with increasing $V$, and yields the second semi-Dirac node at $V=2\\Delta_c$.","core_discovery":"The central claim is that the Fermi-surface crossing points between the main band and the $\\pm Q$ supermodulation-translated bands are special because the $d$-wave order parameters there are equal in magnitude but opposite in sign. Treating the supermodulation as a constant potential $V$ and projecting the $6\\times6$ Bogoliubov problem onto the two crossing bands, the positive-energy dispersions along the crossing line become $E = \\pm V \\pm \\sqrt{\\epsilon^2+\\Delta^2}$. When $V=\\Delta_c\\approx 0.3\\Delta_0$, the minimum at the crossing point reaches zero and forms a semi-Dirac node with quadratic dispersion along the crossing line and linear dispersion perpendicular to it. For $V_{c1}<V<V_{c2}\\approx 2\\Delta_c$, this node splits into two nodes moving in orthogonal directions; at $V_{c2}$ they merge with a translated-band $d$-wave node to form a second semi-Dirac node with the quadratic and linear directions swapped, and beyond $V_{c2}$ the node lifts. With $\\Delta_0=40$ meV, the numerically located critical values are $V_{c1}\\approx12.5$ meV and $V_{c2}\\approx24.5$ meV.","pith_inferences":["If the supermodulation potential is momentum dependent rather than constant, the exact $E=\\pm V\\pm\\sqrt{\\epsilon^2+\\Delta^2}$ form is modified; the critical values would shift, and in some regions the nodes could become small pockets or open gaps instead of touching at zero.","The predicted $\\sqrt{E}$ density of states should be visible as a zero-bias anomaly in scanning tunneling spectroscopy centered on the crossing points, providing a spatially local probe of the scenario.","Because the supermodulation acts like a linear charge density wave, the mechanism is non-chiral; comparing its nodal evolution with the trilayer twisted case could distinguish intrinsic sign-change effects from twist-induced chirality."],"forward_implications":["Bi2212 at optimal doping, with $\\Delta_0\\approx40$ meV and $V\\approx12.5$ meV, sits essentially at $V_{c1}$; its crossing-point minimum should be a semi-Dirac node.","Bi2201, whose gap is roughly half that of Bi2212, sits near $V_{c2}$; its nodes should merge into the second type of semi-Dirac node.","Doping or $c$-axis pressure changes $\\Delta_0$ or $V$, sweeping the system through the critical points and changing the node count per quadrant in the sequence 3, 5, 7, 5, 3.","At a semi-Dirac node the density of states near the chemical potential scales as $\\sqrt{E}$, giving a thermodynamic signature of the transition.","The bilayer splitting, ignored in the main calculation, is small near the nodes (order $(0.3)^2$), so the nodal evolution survives as the dominant effect."],"supporting_citations":[{"why":"This work introduced the twisted-cuprate prediction that the paper argues the supermodulation mimics intrinsically.","marker":"[1]"},{"why":"These photoemission data show the d-wave gap changes sign at the crossing points, the condition that makes the nodal evolution possible.","marker":"[12]"},{"why":"This work supplies the tight-binding dispersion used for the Fermi surfaces and crossing-point locations.","marker":"[14]"},{"why":"This work sets the optimal-doping gap $\\Delta_0\\approx40$ meV for Bi2212, fixing the scale for $V_{c1}$.","marker":"[15]"},{"why":"This work determines the supermodulation potential $V\\approx12.5$ meV from avoided crossings, placing Bi2212 at $V_{c1}$.","marker":"[16]"},{"why":"This work confirms the selective hybridization between main and superstructure bands that supports the constant-$V$ model.","marker":"[17]"}],"fun_headline_variants":["Supermodulation creates semi-Dirac nodes in Bi2212 and Bi2201","BiO buckling produces semi-Dirac nodes at two critical points","Nodal structure evolution in bismuth cuprates via supermodulation","Supermodulation turns d-wave nodes semi-Dirac in cuprates","Two critical couplings turn cuprate nodes into semi-Dirac"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the supermodulation acts as a constant, momentum-independent potential $V$ coupling the main band to the $\\pm Q$ translated bands, with those translated $d$-wave gaps equal in magnitude and opposite in sign to the main gap at the crossing points; the analytic dispersion $E=\\pm V\\pm\\sqrt{\\epsilon^2+\\Delta^2}$ and the critical values $V_{c1}\\approx\\Delta_c$ and $V_{c2}\\approx2\\Delta_c$ follow from that exact relation.","fun_headline_variants_meta":{"raw":{"variants":["Supermodulation creates semi-Dirac nodes in Bi2212 and Bi2201","BiO buckling produces semi-Dirac nodes at two critical points","Nodal structure evolution in bismuth cuprates via supermodulation","Supermodulation turns d-wave nodes semi-Dirac in cuprates","Two critical couplings turn cuprate nodes into semi-Dirac"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000509,"raw_usage":{"total_tokens":2469,"prompt_tokens":928,"completion_tokens":1541,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":544,"completion_tokens_details":{"reasoning_tokens":1456}},"tokens_in":544,"tokens_out":1541,"duration_ms":12204,"temperature":1.0,"reasoning_tokens":1456,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T05:23:06.245227+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Angle-resolved photoemission on optimally doped Bi2212, focused at the crossing point where the main Fermi surface meets its superstructure replica, should reveal a Bogoliubov minimum touching zero near $V\\approx12.5$ meV with a $\\sqrt{E}$ density of states; if the minimum stays gapped or the detected $V$ is far from $0.3\\Delta_0$, the central claim is contradicted.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"These photoemission data show the d-wave gap changes sign at the crossing points, the condition that makes the nodal evolution possible."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"This work supplies the tight-binding dispersion used for the Fermi surfaces and crossing-point locations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"This work sets the optimal-doping gap $\\Delta_0\\approx40$ meV for Bi2212, fixing the scale for $V_{c1}$."},{"cited_title":"Valla, I","cited_arxiv_id":null,"evidence_quote":"This work determines the supermodulation potential $V\\approx12.5$ meV from avoided crossings, placing Bi2212 at $V_{c1}$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"This work confirms the selective hybridization between main and superstructure bands that supports the constant-$V$ model."}],"review_version":1}