{"id":"e1011974-9854-4f1d-a83a-f5fc6bbd54ed","arxiv_id":"1908.04820","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Gap openings at Dirac points away from the Fermi energy are shown to be a direct signature of non-unitary multiorbital superconductivity, characterized by the non-unitarity parameter Υ = Tr[τz ΔΔ†].","lead":"Superconductors with pairing that is not the same on every orbital should open a gap at high-energy Dirac points in their band structure, while ordinary unitary pairing leaves those points untouched. This theoretical analysis gives ARPES a new observable to identify orbital-selective, non-unitary superconducting order in materials like iron chalcogenides and twisted bilayer graphene.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unitary pairing involving a third band at the Dirac energy can gap the Dirac point, so the claimed non-unitary signature requires an isolated two-band crossing that the paper does not justify.","rationale":"The reader's weakest_assumption lists 'other bands cross near the same energy' among several failure modes; my concern isolates that mechanism as the most load-bearing and makes it concrete: even without any competing order, a third band at the same energy coupled by unitary pairing can break the two-band closure that Eq. (4) assumes. The paper's own examples and lattice model are strictly two-band, so they provide no evidence against this failure. This does not overturn the paper's derivation for genuinely isolated Dirac crossings, but it does mean the abstract's general ARPES diagnostic is conditional on an unstated and untested isolation requirement. Since the reader already graded the paper CONDITIONAL and identified related concerns, my read does not move the verdict; it sharpens the condition that must be stated for the claim to hold.","tokens_in":7697,"tokens_out":17447,"duration_ms":180250,"concrete_test":"Construct a three-orbital tight-binding model on a honeycomb lattice: sublattices A and B form the Dirac point at energy -Λ, and a third orbital C is added at each site with on-site energy -Λ and a small normal-state hybridization to A (so C crosses the Dirac energy). Choose a pairing matrix Δ in the full three-orbital basis that is unitary (ΔΔ† ∝ I) but not block-diagonal, e.g., Δ = Δ0 exp(iθ(|A⟩⟨C| + |C⟩⟨A|)) in the A/B/C basis. Compute the BdG electron spectral function A(ω,k) around the original Dirac momentum. If the Dirac point gaps for this globally unitary pairing, the paper's criterion fails; if it remains gapless, the isolated two-band assumption is sufficient for the central claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on the Schrieffer-Wolff reduction in Eq. (4), which projects onto a strictly two-orbital Dirac subspace and treats the Dirac point as an isolated crossing. In real multiorbital systems, the Dirac point need not be isolated: a third band at the same energy can be coupled to the Dirac orbitals by the superconducting pairing. Even if the pairing restricted to the two crossing bands is unitary (Δ ∝ τ0), interorbital unitary coupling to the third band generates a second-order self-energy for the Dirac subspace that is not proportional to identity; after downfolding the third band, the effective pairing in the A/B subspace becomes non-unitary and acquires a τz mass, opening a gap at the Dirac point. The lattice validation in Sec. III uses only a two-band honeycomb model, so it cannot detect this failure. Without an argument that the Dirac point is isolated in the target materials (iron chalcogenides, twisted bilayer graphene), the categorical conclusion that a Dirac gap observed by ARPES implies non-unitary pairing is not established. This isolated-crossing assumption is the load-bearing weak spot: the one-to-one mapping between Υ and the mass γz in Eq. (6) is derived and demonstrated only within a closed two-band subspace.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a spectroscopic diagnostic for non-unitary multiorbital superconductivity. Starting from a two-orbital k·p model with a Dirac point at energy -Λ away from the Fermi level, the authors include superconductivity in a BdG Hamiltonian and perform a Schrieffer-Wolff / second-order perturbative reduction. The effective correction is Δ(K)Δ(K)†/(2Λ); its τ_z component is controlled by Υ = Tr[τ_z ΔΔ†], Eq. (6), so that for Υ ≠ 0 the Dirac point acquires a mass γ_z = Υ/(4Λ) and opens a gap, while for unitary pairing (ΔΔ† ∝ I) it remains gapless. The mechanism is validated with a two-band honeycomb lattice model and spectral-function calculations. The paper further shows that a domain wall between two non-unitary regions hosts a chiral state inside the gap, and concludes that ARPES observations of remote Dirac-point gaps can identify non-unitary multiorbital pairing, with iron chalcogenides and twisted bilayer graphene mentioned as targets.","tokens_in":7960,"tokens_out":11691,"duration_ms":120199,"significance":"If the correspondence is taken in the two-band subspace for which it is derived, the paper gives a clean and potentially useful result: a finite-energy spectral feature whose gap is controlled by an orbital non-unitarity parameter rather than by the Fermi-surface gap. The derivation is explicit, the lattice calculation supports the analytic formula in the stated limit, and the ARPES spectral-function prediction is falsifiable. The topological interface statement is also a concrete additional consequence. The main limitation is that the diagnostic is established only for an isolated two-band crossing; whether it survives in the multiband environment of the cited materials is not addressed, which tempers the significance of the broader claim.","major_comments":[{"comment":"The signature is derived under an isolated two-band assumption that the paper does not establish for the materials it names. The Schrieffer-Wolff reduction in Eq. (4) projects onto the two crossing orbitals and treats all other bands as irrelevant; the lattice check in Sec. III is itself a two-band honeycomb model, so it cannot detect a third band at the same energy. If a third band at the Dirac energy is coupled to the crossing orbitals by an otherwise unitary pairing, downfolding that band produces a second-order self-energy in the A/B subspace whose ΔΔ† is not proportional to the identity; the τ_z part of that self-energy would open the Dirac gap even though the full pairing is unitary. The abstract and conclusion claim that a gap opening is a signature of non-unitary multiorbital order is therefore too strong unless the paper adds an explicit argument that such third-band processes are absent for the iron-chalcogenide or twisted-bilayer cases, or restricts the claim to a strictly isolated two-band crossing.","section":"Sec. II, Eq. (4); Sec. III"},{"comment":"The text states that the numerical Dirac gap agrees with Υ/(2Λ) only when |Λ| ≫ |Δ0|, but the parameters listed for Fig. 2(c,d) use Δ0 = 0.71Λ, which is not in that regime. The comparison panels in Fig. 2(e,f) should state the parameter range and show the condition explicitly, so the reader can see where the analytic result is expected to hold; as written, the agreement for the strong-pairing panels is outside the stated validity domain and does not by itself validate the asymptotic formula.","section":"Sec. III, Fig. 2"}],"minor_comments":[{"comment":"There are several typographical errors: \"gneral\" in the introduction, \"onsite\" where \"onset\" is meant, \"thorugh\" and \"frecuency\" in Sec. IV, and \"spectra function\" for \"spectral function\" in Sec. III.","section":"Throughout"},{"comment":"The phrase \"γx = γy = γz = 0 defines a unitary pairing state\" should be phrased as \"ΔΔ† proportional to the identity\" to avoid ambiguity about whether Δ itself is unitary as a matrix.","section":"Sec. II, after Eq. (5)"},{"comment":"The caption should clarify that Δ0 = 0.71Λ applies only to panels (c,d), and that panels (e,f) scan parameters with the condition |Λ| ≫ |Δ0| enforced where the analytic formula is claimed.","section":"Sec. III, Fig. 2 caption"},{"comment":"The sentence in the conclusions about an incipient charge density wave should be clarified, because a CDW-induced gap at the same Dirac point would be a false positive for the non-unitarity diagnostic unless the spectral function is compared above and below the superconducting transition.","section":"Sec. V"}],"recommendation":"major_revision","confidential_remarks":"The central two-band derivation is sound and the paper is clearly written, but the abstract and conclusions overstate the applicability. The third-band confounding scenario is not exotic: any unitary pairing that couples a degenerate third band to the Dirac orbitals will alter the low-energy two-band subspace in a way that mimics the non-unitary gap. I would ask the authors to add an explicit statement of the isolated-crossing assumption and either provide a band-structure argument for the named materials or limit the conclusions to the minimal model."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThis paper gives a clean, minimal argument: a Dirac point away from the Fermi level in a multiorbital superconductor acquires a mass term gamma_z = Upsilon/(4Lambda), where Upsilon = Tr[tau_z Delta Delta^dagger] is an orbital non-unitarity measure. Unitary pairing leaves the Dirac crossing intact; pairing that breaks orbital symmetry (e.g., pairing on one sublattice only) opens a gap. The Schrieffer-Wolff step is straightforward, and the lattice model in Sec. III matches the analytic prediction well in the stated |Lambda| >> |Delta| regime. The ARPES spectral-function panels are convincing, and the interface-mode section is a nice bonus. The paper is worth serious engagement.\n\nThe soft spots are real but not fatal. First, the abstract overgeneralizes: only non-unitary states with Upsilon != 0 open the gap, while other non-unitary matrices merely shift the Dirac point. The wording should be tightened. Second, and more importantly, the derivation assumes an isolated two-band crossing. The stress-test note is right: if a third band sits near the Dirac energy and is coupled to the crossing orbitals by pairing, downfolding can generate an effective tau_z mass even from nominally unitary pairing. The paper does not address this, and the honeycomb-lattice validation uses only two bands, so it cannot detect the failure. For real systems like iron chalcogenides or twisted bilayer graphene, the one-to-one mapping between Upsilon and the observed gap is not guaranteed without checking the isolation of the crossing. Competing orders (CDW, magnetic) that also gap the same point are mentioned only in passing.\n\nStill, the central two-band result is correct, and the Upsilon parameter is a useful diagnostic in the appropriate regime. The paper would be strengthened by a discussion of how to verify the isolated-crossing condition and how to distinguish pairing-induced gaps from other symmetry-breaking orders. I would cite this paper for the Upsilon formula and would bring it to a reading group focused on unconventional superconductivity. A serious referee should see it; with revisions that clarify scope and address the isolation caveat, it is a solid contribution.","headline":"A clean, useful two-band criterion for detecting orbital-selective pairing via remote Dirac gaps, with an overbroad abstract and a real isolated-crossing caveat.","tokens_in":8499,"tokens_out":1668,"would_cite":true,"duration_ms":19462,"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":"This paper shows that a Dirac crossing away from the Fermi energy that develops a gap upon entering the superconducting state is a direct signature of non-unitary multiorbital pairing, with the gap size controlled by the non-unitarity…","keywords":["non-unitary superconductivity","multiorbital pairing","Dirac points","ARPES","orbital-selective pairing","Bogoliubov–de Gennes","pair density wave","chiral interface states"],"falsifier":"Measure the gap at a remote Dirac crossing in a superconductor as a function of the energy distance $\\Lambda$ (for example by doping or by tuning spin-orbit splitting): the paper predicts $\\gamma_z=\\Upsilon/(4\\Lambda)$, so the gap should shrink as $1/\\Lambda$ and vanish when $\\Upsilon=0$. A remote Dirac gap that does not scale inversely with $\\Lambda$, or that appears in a material whose gap matrix is independently known to be unitary, would falsify the diagnosis.","tokens_in":7483,"feed_emoji":"🔬","tokens_out":8956,"duration_ms":83553,"temperature":0.7,"pith_summary":"This paper argues that the fate of Dirac crossings located well away from the Fermi energy can diagnose the orbital structure of a superconductor. In a normal-state band structure with a Dirac point below the chemical potential, a unitary pairing state leaves the crossing intact, whereas a non-unitary multiorbital pairing state opens a gap there. The magnitude of that gap is set by the non-unitarity parameter $\\Upsilon = \\mathrm{Tr}[\\tau_z \\Delta\\Delta^\\dagger]$ divided by the energy distance $\\Lambda$ to the crossing, so an angle-resolved photoemission measurement of a remote Dirac gap would directly reveal orbital-selective pairing. The same mechanism also predicts chiral interface states between domains with opposite non-unitary order. A sympathetic reader would care because ARPES is a standard tool and this gives a high-energy spectroscopic fingerprint of an otherwise hard-to-determine pairing symmetry.","feed_headline":"A remote Dirac gap reveals non-unitary superconductivity","feed_subtitle":"Angle-resolved photoemission reads the orbital structure of a superconductor from high-energy Dirac gaps.","key_machinery":"The object that carries the argument is the non-unitarity parameter $\\Upsilon = \\mathrm{Tr}[\\tau_z \\Delta\\Delta^\\dagger]$ for the $2\\times2$ superconducting gap matrix $\\Delta$ in orbital space. It appears as the coefficient of the $\\tau_z$ Pauli matrix in the effective Dirac Hamiltonian obtained by a Schrieffer-Wolff / second-order perturbation treatment of the pairing at the Dirac momentum, giving the mass term $\\gamma_z = \\Upsilon/(4\\Lambda)$. This identity maps the orbital structure of the order parameter directly onto the existence and size of a spectral gap at a remote Dirac crossing: $\\Upsilon\\neq0$ opens the gap, $\\Upsilon=0$ does not. The same parameter also fixes the topological response of the gapped Dirac cone through the Chern number $C=\\frac12\\,\\mathrm{sign}(\\Upsilon)\\mathrm{sign}(\\Lambda)$, which predicts chiral states at interfaces between domains with opposite non-unitary order.","core_discovery":"The paper's central claim, stated on its own terms, is that non-unitarity in orbital space is the controlling quantity for Dirac crossings away from the Fermi level in a superconductor. Starting from a minimal two-orbital $k\\cdot p$ Dirac Hamiltonian $H_0^{\\mathrm{DP}}(k)=-\\Lambda\\tau_0+\\tau_x k_x+\\tau_y k_y$ and coupling it to a $2\\times2$ gap matrix $\\Delta$, a Schrieffer-Wolff reduction gives the effective correction $\\Delta\\Delta^\\dagger/2\\Lambda$. Decomposing this correction into Pauli components, the coefficient $\\gamma_z = \\Upsilon/(4\\Lambda)$ of $\\tau_z$ acts as a mass term, opening a gap at the Dirac point precisely when the non-unitarity parameter $\\Upsilon=\\mathrm{Tr}[\\tau_z\\Delta\\Delta^\\dagger]$ is nonzero. Unitary pairing ($\\Upsilon=0$) leaves the Dirac point gapless; non-unitary pairing such as pairing on only one orbital ($\\Delta=\\Delta_0(\\tau_0+\\tau_z)/2$) opens it. The claim is verified in a honeycomb-lattice model, including the predicted $1/\\Lambda$ scaling, and extended to domain walls where the gapped Dirac points carry half-integer Chern numbers and produce chiral interface excitations.","pith_inferences":["Beyond the paper, the $\\Upsilon/(4\\Lambda)$ scaling suggests a quantitative route: once the normal-state Dirac position is known from ARPES, the measured remote gap size could be converted into an estimate of the orbital-selectivity strength $\\Upsilon$, something the paper does not attempt.","Beyond the paper, the diagnosis is safest in compounds where charge-density-wave, magnetic, or lattice-symmetry-breaking orders are absent, since any competing order that gaps the same crossing would mimic the non-unitary signature; the paper notes the pair-density-wave example itself carries a sublattice charge imbalance.","Beyond the paper, the half-integer Chern number per gapped Dirac cone implies that a network of non-unitary domains could act as a chiral conductor at remote energies, decoupled from the low-energy superconducting condensate; the paper discusses a single interface but not such networks."],"forward_implications":["An ARPES scan of the superconducting state that shows a gap at a Dirac point well below the Fermi energy is direct evidence of non-unitary multiorbital pairing, even if the low-energy gap looks conventional.","The remote Dirac gap should scale as $\\Upsilon/(4\\Lambda)$, so sweeping the chemical potential or spin-orbit splitting changes the gap in a predictable way.","Domain walls between two degenerate non-unitary orders produce a chiral state inside the remote Dirac gap, while leaving the Fermi-level superconducting gap unchanged.","The criterion is independent of the momentum structure of the pairing, so it applies to $s$-wave, $d$-wave, or other pairing channels on equal footing.","In three-dimensional Dirac semimetals the mechanism changes character: a $2\\times2$ Dirac cone cannot open a gap, only shift, so the diagnostic is dimension-specific."],"supporting_citations":[{"why":"Report ARPES gap openings at remote Dirac points in iron chalcogenides, the experimental motivation the mechanism is built to explain.","marker":"[8,9]"},{"why":"Reports orbital-selective pairing in iron-based superconductors, the multiorbital phenomenon the non-unitarity criterion targets.","marker":"[7]"},{"why":"Establishes that the iron chalcogenide band structure hosts remote Dirac crossings, making the ARPES test feasible.","marker":"[19]"},{"why":"Supplies the recursive Green's function algorithm used to compute the interface spectral function exactly.","marker":"[16]"},{"why":"Index theorem connecting the half-integer Chern numbers to the presence of a single interface state in the Dirac gap.","marker":"[17]"},{"why":"Provide the spin-orbit coupling form used in the honeycomb-lattice model to create remote Dirac points at half filling.","marker":"[12,13]"},{"why":"Prior discussion of remote Dirac points in superconductors that the paper adopts and generalizes.","marker":"[10]"},{"why":"Show that twisted bilayer graphene hosts Dirac points about 10 meV below the Fermi level, cited as a candidate application.","marker":"[23–26]"}],"fun_headline_variants":["Dirac gap opening exposes non-unitary superconductivity","Finite-energy Dirac gaps flag non-unitary pairing","ARPES reads orbital structure from Dirac gaps","Non-unitary order leaves a Dirac-gap fingerprint","High-energy Dirac points reveal orbital pairing"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the Dirac point is an isolated two-band crossing whose energy separation from the Fermi level is much larger than the pairing energy, so the pairing acts as a weak perturbation; if other bands or a competing order such as a charge-density wave or magnetism also live at that crossing, the measured gap can no longer be attributed uniquely to non-unitary pairing.","fun_headline_variants_meta":{"raw":{"variants":["Dirac gap opening exposes non-unitary superconductivity","Finite-energy Dirac gaps flag non-unitary pairing","ARPES reads orbital structure from Dirac gaps","Non-unitary order leaves a Dirac-gap fingerprint","High-energy Dirac points reveal orbital pairing"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000646,"raw_usage":{"total_tokens":2962,"prompt_tokens":934,"completion_tokens":2028,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":550,"completion_tokens_details":{"reasoning_tokens":1957}},"tokens_in":550,"tokens_out":2028,"duration_ms":16119,"temperature":1.0,"reasoning_tokens":1957,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:32:59.944157+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the gap at a remote Dirac crossing in a superconductor as a function of the energy distance $\\Lambda$ (for example by doping or by tuning spin-orbit splitting): the paper predicts $\\gamma_z=\\Upsilon/(4\\Lambda)$, so the gap should shrink as $1/\\Lambda$ and vanish when $\\Upsilon=0$. A remote Dirac gap that does not scale inversely with $\\Lambda$, or that appears in a material whose gap matrix is independently known to be unitary, would falsify the diagnosis.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports orbital-selective pairing in iron-based superconductors, the multiorbital phenomenon the non-unitarity criterion targets."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes that the iron chalcogenide band structure hosts remote Dirac crossings, making the ARPES test feasible."},{"cited_title":"Komendov\\'a , author A","cited_arxiv_id":null,"evidence_quote":"Prior discussion of remote Dirac points in superconductors that the paper adopts and generalizes."}],"review_version":1}