{"id":"584b54d6-5731-4acc-aefc-21cad51b6a93","arxiv_id":"2412.16760","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Chiral double-layer edge modes map to massive Dirac edge states, with a formal bulk-edge correspondence and boundary-phase control of edge currents.","lead":"This paper derives edge and bulk quasiparticle modes for a chiral electron double layer superconductor using a mapping to a massive Dirac equation. It shows a formal 'bulk-edge correspondence' and proposes boundary-condition phases as a way to control edge currents.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The lower Dirac block in Eq. (B6) is not an eigenfunction of h_D(-|Δ|): the sign of ρ^{-1} is wrong, and this swaps the two middle components of Ψ in Eq. (10), reversing which layer's edge current ψ controls.","rationale":"I read the paper as a formal construction of edge and bulk eigenfunctions for a chiral double layer via the unitary map to two massive Dirac Hamiltonians. The Dirac part itself, including the K_n(cr) solutions and the analytic continuation c→±ik to J_n/Y_n bulk waves, is a legitimate and clean observation. The load-bearing weak point is not the analytic continuation as such, but the double-layer eigenfunction used to make physical predictions. Direct substitution shows that the lower Dirac block in Eq. (B6) is not an eigenfunction of h_D(-|Δ|) for the in-gap energies: the required coefficient is -ρ^{-1}, not +ρ^{-1}. This error propagates through W† into Eq. (10) and reverses the layer assignments in the ψ-controlled currents, so the central claim that ψ separately controls the two edge currents is unsupported as printed. The reader's boundary-condition concern is also valid: no self-adjoint boundary condition at r0 is specified, so the continuum of in-gap modes is not fully determined. I do not see the two concerns as identical; the sign error is independent and would survive even after a boundary condition is chosen. Because the error is a correctable algebraic mistake rather than a flaw in the overall strategy, the reader's conditional verdict remains appropriate, but the requested revisions should explicitly require re-deriving Eq. (10) and the current formulas with the corrected lower-block sign.","tokens_in":7535,"tokens_out":28602,"duration_ms":246410,"concrete_test":"Substitute the printed lower spinor (K_n(cr), iρ^{-1}K_{n+1}(cr)e^{iα})e^{inα} into h_D(-|Δ|) and use Eq. (B5): the first radial equation demands cρ^{-1}=E+|Δ|, while Eq. (B5) gives cρ^{-1}=-(E+|Δ|). Repeat with the corrected coefficient -iρ^{-1} and recompute Ψ_E,n = W†Φ_E,n; verify that Eq. (10)'s middle components become i(ρ+e^{iψ}ρ^{-1}) and i(ρ-e^{iψ}ρ^{-1}) in the opposite order. If this check confirms the swap, the ψ=0 layer suppression and C↑,C↓ in Eqs. (14)-(15) switch layers, falsifying the printed control claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"For the lower block h_D(-|Δ|), write the second spinor component as i b K_{n+1}(cr)e^{iα}. The first radial equation for mass -|Δ| is c b = E+|Δ|, so b=(E+|Δ|)/c. Equation (B5) gives ρ=(E-|Δ|)/c, hence the paper's choice b=ρ^{-1}=c/(E-|Δ|)=-(E+|Δ|)/c. This satisfies c b = E+|Δ| only at E=-|Δ| with c=0, so the printed lower-block spinor is not an eigenfunction in the in-gap regime. The correct coefficient is b=-ρ^{-1}. Repeating the W† construction with g=-iρ^{-1}K_{n+1}e^{iα} interchanges the two middle components of Eq. (10): the terms i(ρ-e^{iψ}ρ^{-1}) and i(ρ+e^{iψ}ρ^{-1}) appear in the opposite order (for A'_n=e^{iψ}A_n). At E=0, ψ=0, this reverses which layer's zero mode is suppressed, and the layer current coefficients C↑,C↓ in Eqs. (14)-(15) trade their ψ-dependence. The claimed independent control of the two edge currents is therefore not supported by the printed calculation. The boundary-condition issue raised by the reader is real, but this sign error is more immediately load-bearing because it changes the physical predictions even within the family of K_n solutions.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies quasiparticle excitations of a chiral electron double layer with interlayer pairing, using a unitary mapping (App. A) from the BdG Hamiltonian H_EEDL to a block-diagonal massive Dirac Hamiltonian H_Dirac. In a circular geometry with a hole, the authors write edge-mode eigenfunctions in terms of modified Bessel functions K_n(cr) with in-gap energies E = ±sqrt(|Δ|^2 - c^2), then pass to bulk modes by the formal replacement c → ±ik and identify E = ±|Δ| as a localization-delocalization transition. They also compute layer-resolved edge currents and claim that the relative boundary-condition phase ψ controls the two layers' edge currents separately, allowing individual switching off of zero modes. The paper's central claims are the edge/bulk dichotomy, the analytic-continuation bulk-edge correspondence, and the ψ-controlled layer currents.","tokens_in":7892,"tokens_out":8685,"duration_ms":73024,"significance":"If the claims were established, the paper would offer a useful duality-based route to edge-state wavefunctions and currents in a paired chiral double layer, with a concrete experimental handle (boundary-condition phase ψ) for controlling layer-resolved currents. The strengths include an explicit unitary transformation in App. A, a clean separation-of-variables calculation, and a clear formal connection between the edge-mode functions K_n and the bulk functions J_n,Y_n through analytic continuation. However, two load-bearing issues prevent acceptance in the present form: a sign error in the lower Dirac block that changes the physical predictions for layer currents, and the absence of an imposed boundary condition at the hole, which leaves the in-gap spectrum underdetermined. The claimed transition at E = ±|Δ| is also not established as a genuine transition. These issues are local in the sense that the calculational framework is promising, but the printed results and the headline physical conclusions require substantial correction.","major_comments":[{"comment":"The lower-block spinor in Eq. (B6) is not an eigenfunction of h_D(-|Δ|). Substituting (A'K_n, A'iρ^{-1}K_{n+1}e^{iα}) into the first radial Dirac equation for mass -|Δ| gives cρ^{-1} = E + |Δ|, whereas the upper-block relation cρ = E - |Δ| combined with c^2 = |Δ|^2 - E^2 implies cρ^{-1} = -(E + |Δ|); the printed expression is therefore only consistent at E = -|Δ| in the c → 0 limit. The correct lower-block coefficient is -iρ^{-1}. Repeating the W† construction with the correct sign interchanges the middle two components of Eq. (10): the terms i(ρ + e^{iψ}ρ^{-1}) and i(ρ - e^{iψ}ρ^{-1}) appear in the opposite order. This changes the physical predictions already within the family of K_n solutions, e.g., at E = 0, ψ = 0 the suppressed zero mode switches from one layer to the other, and the ψ-dependence of C↑ and C↓ in Eqs. (14)-(15) is exchanged. Since the central claim of independent layer control rests on these formulas, this sign error is load-bearing and must be corrected.","section":"Appendix B, Eq. (B6); Eq. (10)"},{"comment":"The edge-mode problem is underdetermined because no boundary condition at the circular hole r = r0 is imposed. Equation (5) is a decaying solution of the radial Dirac equation for any c in (0, |Δ|); without specifying the self-adjoint boundary condition at r0 (or the matching condition for a finite disk), the in-gap spectrum is not discrete and every c is admissible. The text acknowledges this indirectly when it says that 'the actual values of the boundary conditions are determined by the coupling of the double layer to its environment', but that statement does not replace the missing derivation. As a result, the paper's 'edge-mode spectrum' and the precise edge/bulk distinction are not yet established.","section":"§2, Eqs. (5)-(7) and text after Eq. (11)"},{"comment":"The claims of a 'bulk-edge correspondence' and of a 'transition' at E = ±|Δ| are stronger than what is demonstrated. The replacement c → ±ik is a formal analytic continuation that turns the modified Bessel functions K_n into combinations of J_n and Y_n; it does not, by itself, identify which bulk modes satisfy the same boundary condition as a given edge mode, and no scattering problem is solved. The divergence of the decay length 1/c as E^2 → |Δ|^2 shows that the edge modes become extended at the spectral boundary, but this is a continuum limit rather than evidence of a thermodynamic or spectral phase transition. The authors should either prove a precise statement of the correspondence (e.g., a bijection between edge eigenstates and bulk scattering states with the same boundary condition) or substantially soften the transition language.","section":"§3 and Conclusion, text after Eqs. (8)-(9)"}],"minor_comments":[{"comment":"The four-component wavefunction is written as Ψ(r) = (ψ↑1(r), ψ↑2(r), ψ↓1(r)ψ↓2(r))^T; a comma is missing between ψ↓1(r) and ψ↓2(r).","section":"§1, after Eq. (1)"},{"comment":"For the full disk the text says to replace K_n by I_n, but no boundary condition is given that would select the allowed values of c (or the analog for I_n). Please clarify which geometry (infinite plane with a hole versus finite disk) is used in each formula and where the boundary condition enters.","section":"§2, paragraph after Eq. (7)"},{"comment":"The asymptotic expansion of K_n(-iz) is used to motivate the identification of bulk modes; it would be helpful to cite a uniform identity such as the DLMF relation between modified and ordinary Bessel functions, since the asymptotic form shown is not valid uniformly for all n and z.","section":"Eq. (8)"}],"recommendation":"major_revision","confidential_remarks":"The sign error in Eq. (B6) is easily fixable but it directly changes the headline claim about ψ-controlled layer currents. The missing boundary-condition analysis is a more fundamental gap, but it could be addressed in a revision by imposing a concrete boundary condition at r0 and rederiving the quantization condition for c. The duality references [7,8] are published and independent, so there is no circularity concern. The paper is not ready for acceptance in its present form."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: the paper's advertised new result—phase ψ control of layer-selective edge currents—is not supported by the calculation as printed. The reader's boundary-condition concern is valid, but the stress-test note found a more immediate problem: Eq. (B6) writes the lower-block Dirac spinor with ρ^{-1}, and that spinor is not an eigenfunction of h_D(−|Δ|) except at c=0. The correct coefficient is −ρ^{-1}. This is not cosmetic. Repeating the W† transformation with the correct sign swaps the two middle components of Eq. (10), and for ψ=0, E=0 it is the bottom layer's zero mode that is suppressed, not the top layer's. The claimed independent control of the two edge currents therefore rests on a sign error that reverses the physical prediction.\n\nWhat the paper does well: the unitary mapping from the EEDL to two massive Dirac Hamiltonians is a legitimate route, and the circular ansatz with modified Bessel functions is a clean way to write edge-mode profiles. The relative phase ψ between layer amplitudes is a real, modest extension of the authors' earlier duality work, and the current formulas are explicit enough to check.\n\nThe soft spots, in order of severity. First, the sign error above. It is load-bearing because Eqs. (14)–(15) and the ψ=0 example change when the middle components are swapped. Second, no boundary condition is ever imposed at r0. The decay constant c in (0,|Δ|) is never fixed, so the in-gap spectrum is underdetermined; any c is allowed. The \"transition\" at E=±|Δ| is only the limit c→0 and is not established as a physical transition. Third, the analytic continuation c→±ik is formal: replacing K_n(cr) by K_n(−ikr) gives J_n+iY_n, but no boundary matching at the hole is performed, so the \"bulk-edge correspondence\" is asserted rather than demonstrated.\n\nWho should read this: someone working on chiral double-layer superconductors might find the setup useful once corrected, and the paper is a good teaching example of how a sign error in a transformed eigenfunction can flip the physics. But as printed, the central predictions are not reliable, and I would not cite it yet.\n\nRecommendation: send it to peer review; it deserves a serious referee, but expect major revision. The boundary-condition gap can be fixed by imposing explicit boundary conditions at r0, and the sign error is a one-line correction—but that correction changes the paper's main claim.","headline":"A sign error in the lower Dirac block flips which layer's edge current is controlled; the boundary-condition omission is real but secondary.","tokens_in":8344,"tokens_out":5599,"would_cite":false,"duration_ms":44197,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"In chiral electron double layers, quasiparticle modes split into edge and bulk spectra linked by an analytic continuation.","keywords":["chiral double layer","Bogoliubov–de Gennes","edge modes","bulk-edge correspondence","analytic continuation","quasiparticle currents","Dirac Hamiltonian","superconductivity"],"falsifier":"Solving the full Bogoliubov–de Gennes equations in a disk with a hole and an explicit boundary condition (e.g., $\\Psi(r_0)=0$) would show whether a continuum of in-gap edge modes exists; if the spectrum is discrete or empty, the claimed transition at $E=\\pm|\\Delta|$ and the analytic-continuation bulk-edge correspondence would not hold in that geometry.","tokens_in":7351,"feed_emoji":"🌀","tokens_out":5998,"duration_ms":47449,"temperature":0.7,"pith_summary":"The paper studies quasiparticles in a chiral electron double layer with interlayer pairing, described by a Bogoliubov–de Gennes Hamiltonian. It argues that in a circular geometry with an edge, the quasiparticle modes split into edge modes localized at the edge and bulk modes, and the two are connected by an analytic continuation $c\\to \\pm ik$. It further claims that a phase $\\psi$ in the boundary conditions between the two layers controls the edge current in each layer separately, allowing one layer's edge mode to be switched off independently. If these claims hold, the edge currents can be manipulated by external fields and are sustained by the supercurrents.","feed_headline":"Chiral double-layer quasiparticles split into edge and bulk modes","feed_subtitle":"A boundary phase between the two layers switches each edge current independently, offering experimental control.","key_machinery":"The central object is the Bogoliubov–de Gennes Hamiltonian $H_{\\rm EEDL}(\\Delta)$ of the chiral electron double layer, which a unitary transformation $W$ maps to a block-diagonal pair of massive Dirac Hamiltonians $h_D(\\pm|\\Delta|)$. The eigenfunctions of each Dirac layer in a circular geometry are built from modified Bessel functions $K_n(cr)$ for edge modes; the analytic continuation $c\\to \\pm ik$ converts them via the identity $K_n(-iz)=\\frac{\\pi i}{2}e^{i\\pi(n+1)/2}[J_n(z)+iY_n(z)]$ into Bessel combinations describing bulk modes. The boundary phase $\\psi$ between the layer amplitudes $A_n$ and $A'_n=e^{i\\psi}A_n$ is the control parameter that fixes the separate edge currents in the two layers.","core_discovery":"The central claim is that the quasiparticle spectrum of the chiral double layer separates into in-gap edge modes with energies $|E|<|\\Delta|$, described by modified Bessel functions $K_n(cr)$ with $c=\\sqrt{|\\Delta|^2-E^2}$, and bulk modes with $|E|>|\\Delta|$ obtained by the analytic continuation $c\\to \\pm ik$. The paper presents this continuation as a bulk-edge correspondence: the same eigenfunctions interpolate between exponentially localized and extended behavior as $E$ crosses $\\pm|\\Delta|$, which it likens to a localization-delocalization transition. The phase $\\psi$ entering through $A'_n=e^{i\\psi}A_n$ determines the relative weight of the two Dirac sectors in the double-layer wavefunction, and through it the magnitude, direction, and phase of the edge currents in the two layers. In particular, $\\psi=0$ suppresses the zero mode in the top layer and $\\psi=\\pi$ suppresses it in the bottom layer. The edge quasiparticle currents are balanced by supercurrents, so they can be maintained without external sources.","pith_inferences":["If one imposes a concrete boundary condition at the hole radius (for instance, vanishing of the wavefunction at $r_0$), the in-gap energies would become discrete; whether the analytic-continuation transition survives this is a direct test of the paper's edge-mode picture.","The mapping to two Dirac Hamiltonians suggests the same edge-mode structure should appear in any paired double-layer system dual to the electron–hole bilayer, so the predictions could be checked in dipolar exciton bilayers.","Because the boundary phase $\\psi$ controls the currents, a gated edge could be used as a switch for the edge-channel conductance in a single layer; this is an experimentally testable consequence the paper does not develop.","The analytic continuation $c\\to \\pm ik$ is formal; a lattice simulation of the tight-binding BdG equations in a disk would clarify whether the bulk-edge correspondence persists beyond the Dirac continuum approximation."],"forward_implications":["If the bulk-edge correspondence holds, edge modes and bulk bands are connected through $c\\to \\pm ik$, with $E=\\pm|\\Delta|$ as the critical energies where the localization length diverges.","The boundary phase $\\psi$ can independently suppress the zero mode in either layer, implying that edge channels can be switched layer-selectively.","The edge currents in the two layers flow in opposite directions, and their magnitude depends on $\\psi$, so the phase can tune the current.","Charge conservation requires the quasiparticle edge currents to be balanced by supercurrents, so the edge modes are sustained rather than radiated away.","The transition at $E=\\pm|\\Delta|$ is reminiscent of a localization-delocalization transition, suggesting that varying the energy drives the quasiparticles from edge-localized to bulk-extended states."],"supporting_citations":[{"why":"Supplies the pairing mechanism in the electron double layer via duality to the electron-hole bilayer, establishing the model's physical basis.","marker":"[7]"},{"why":"Provides the BdG Hamiltonian of the chiral double layer used here.","marker":"[8]"},{"why":"Defines the bulk-edge correspondence concept that the paper claims to realize in this system.","marker":"[9]"},{"why":"Another formulation of bulk-edge correspondence that the paper's construction extends.","marker":"[10]"},{"why":"Supplies the Bessel-function identity used in the analytic continuation from edge to bulk modes.","marker":"[16]"}],"fun_headline_variants":["Phase controls edge currents in chiral double layers","Chiral double-layer edge modes show localization transition","Independent edge current control in chiral double layers via phase","Bulk-edge correspondence from analytic continuation in double layers","Switch either layer's edge current by tuning phase in chiral double layer"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper assumes that the in-gap edge-mode eigenfunction $K_n(cr)$ is valid for every $c$ in $(0,|\\Delta|)$ without imposing a boundary condition at the hole radius, so the edge spectrum is taken as given rather than derived from boundary data.","fun_headline_variants_meta":{"raw":{"variants":["Phase controls edge currents in chiral double layers","Chiral double-layer edge modes show localization transition","Independent edge current control in chiral double layers via phase","Bulk-edge correspondence from analytic continuation in double layers","Switch either layer's edge current by tuning phase in chiral double layer"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000745,"raw_usage":{"total_tokens":3282,"prompt_tokens":869,"completion_tokens":2413,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":485,"completion_tokens_details":{"reasoning_tokens":2337}},"tokens_in":485,"tokens_out":2413,"duration_ms":16586,"temperature":1.0,"reasoning_tokens":2337,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T10:14:56.690041+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solving the full Bogoliubov–de Gennes equations in a disk with a hole and an explicit boundary condition (e.g., $\\Psi(r_0)=0$) would show whether a continuum of in-gap edge modes exists; if the spectrum is discrete or empty, the claimed transition at $E=\\pm|\\Delta|$ and the analytic-continuation bulk-edge correspondence would not hold in that geometry.","supporting_citations":[{"cited_title":"Sinner, Yu","cited_arxiv_id":null,"evidence_quote":"Supplies the pairing mechanism in the electron double layer via duality to the electron-hole bilayer, establishing the model's physical basis."},{"cited_title":"Ziegler, A","cited_arxiv_id":null,"evidence_quote":"Provides the BdG Hamiltonian of the chiral double layer used here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the bulk-edge correspondence concept that the paper claims to realize in this system."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Another formulation of bulk-edge correspondence that the paper's construction extends."},{"cited_title":"Abramowitz and I","cited_arxiv_id":null,"evidence_quote":"Supplies the Bessel-function identity used in the analytic continuation from edge to bulk modes."}],"review_version":1}