{"id":"1e97daef-ead7-485b-8bff-0249dcf3f0b0","arxiv_id":"2506.08714","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"SU(2) non-Abelian hopping in a 2D Hatano-Nelson cylinder yields Hopf-link band braiding and a GBZ-based polarization parameter that classifies left, right, and bipolar skin modes.","lead":"This paper studies a two-dimensional lattice model with spin-dependent, non-reciprocal hopping and reports that its complex energy bands form braided loops with a measurable topological winding using only nearest-neighbor couplings. It proposes a polarization measure built from the generalized Brillouin zone that classifies left-, right-, and bipolar skin localization, and analyzes unusual degenerate boundary states.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The I_p-based NHSE classifier is under-specified: Eq. (11) counts roots of Eq. (10) without stating which energy/root set is used, its phase maps are admitted to be numerically unconverged, and no direct GBZ-to-exact-OBC validation is shown.","rationale":"I read the paper in good faith and checked the parts within my competence. The braiding-degree derivation is internally consistent: after subtracting (1/2)Tr H, the σ0 terms cancel and Eq. (8) is the integral of d/dk_x ln[J_R^2 sin^2β e^{-2ik_x}+J_L^2 sin^2α e^{2ik_x}], whose winding is indeed ±2 away from the EP surface Eq. (9); the cylinder reduction to the 1D model of Ref. [52] is openly credited. The genuinely new and load-bearing claim is the I_p polarization classifier for NHSE. The reader's weakest assumption correctly identifies the auxiliary-GBZ and root-counting construction as the shaky link. My stress-test sharpens this: I_p is not merely under-validated, it is under-specified in the manuscript, because Eq. (11) does not say which β_n set is used when the characteristic equation depends on E and the GBZ is a continuum contour. The paper's own admission of serrated, 15-day-converged phase maps and the absence of code/data mean a referee cannot reproduce or falsify the central classification. This does not prove the physics wrong: direct exact-diagonalization density plots at selected points (Figs. 4, 7, 8) are real evidence and no internal contradiction is apparent. The appropriate response is therefore conditional acceptance requiring a clear definition of I_p, a convergence study, and a direct GBZ-to-exact-OBC comparison, not rejection. I therefore leave the reader's CONDITIONAL verdict unchanged.","tokens_in":16039,"tokens_out":25523,"duration_ms":319221,"concrete_test":"Recompute I_p and the C/Q eigenstate maps on the same J_L-J_R grid as Fig. 4(a) (α=π/3, β=π/2, k_y=1), using exact diagonalization of Hamiltonian (3) for N=50 and N=200. Then evaluate Eq. (11) with the four roots of Eq. (10) at each exact OBC eigenvalue, using high-precision arithmetic, and overlay the I_p=±1 and I_p=0 region boundaries on the Q boundaries from the exact eigenstates. Also plot individual eigenstate densities (not just the average P(j)) for every state in the nominal bipolar region. If the I_p boundaries move by more than one grid step under the per-eigenvalue versus averaged convention, or if no individual eigenstate in the bipolar region is localized at both ends, the classifier and the bipolar-skin claim would need revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central new result is the I_p-based classification of left/right/bipolar skin modes. That result is load-bearing and under-specified. Eq. (11) defines I_p as a count of all 2M GBZ roots β_n relative to |β|=1, but for the quartic characteristic equation (10) the four roots depend on the energy E; the GBZ condition |β_M|=|β_{M+1}| selects a contour of β values, not a finite set {β_1,...,β_4}. The manuscript never states whether I_p is evaluated at a single OBC eigenvalue, averaged over the OBC spectrum, or computed from a θ-grid of the auxiliary resultant; different choices can change the root count and hence the I_p sign or zero. The phase maps in Fig. 4(a) are explicitly conceded to be numerically unconverged ('serrated patterns... stems from insufficient computational precision... requires a computational duration of 15 days'), and no code or data is released. The only validation against real-space eigenstates is at a handful of sampling points (Figs. 4, 9, 10). Thus the central claim that I_p 'quantitatively discerns' left-, right-, and bipolar skin modes is not established as stated. By contrast, the Hopf-link braiding part reduces to the 1D result of Ref. [52] after subtracting (1/2)Tr H, and the braiding-degree computation modulo convergence checks out.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript studies a two-dimensional Hatano-Nelson model with SU(2) non-Abelian gauge fields on a cylinder geometry. The authors claim two main results: (i) the complex-energy bands form a Hopf-link braiding with braiding degree nu = +-2, characterized by an exceptional-point boundary J_R^2 sin^2(beta) = J_L^2 sin^2(alpha), achievable with only nearest-neighbor x-direction couplings; and (ii) a generalized-Brillouin-zone polarization parameter I_p, defined in Eq. (11), quantitatively classifies left-, right-, and bipolar skin modes, including a bipolar skin effect at J_L = J_R with alpha = beta. The Hopf-link computation is presented as exact, while the I_p classification is supported by numerical phase diagrams and a limited set of real-space eigenstate comparisons. The paper also discusses zero-imaginary-energy eigenstates at the topological boundary, with claims of degeneracy and size-independent bipolar localization.","tokens_in":16173,"tokens_out":5805,"duration_ms":70574,"significance":"If the central claims hold, the paper would be a useful contribution to non-Hermitian topological physics: it would show that non-Abelian gauge fields can generate spectral braiding and bipolar skin localization without long-range couplings, which is a notable conceptual step beyond previous Abelian-gauge and long-range-hopping constructions. The exact derivation of the braiding degree is a genuine strength: Eq. (8) follows from the model Hamiltonian with no fitted parameters, and the cancellation of k_y in the traceless part correctly reduces the cylinder problem to a one-dimensional model, matching the known construction of Ref. [52]. The cross-validation of I_p against directly encoded real-space eigenstates (Section IV C) is good practice and is convincing at the selected sampling points. However, the central quantitative claim about I_p is currently under-specified and is not supported by fully converged numerical data, so the paper requires major revision before the classification result can be accepted.","major_comments":[{"comment":"The definition of I_p is not operationally well specified. The characteristic equation (10) is quartic in beta for each complex energy E, and the GBZ condition |beta_M| = |beta_{M+1}| stated in Section IV A selects a continuum of (beta, E) pairs, not a fixed finite set {beta_1, ..., beta_{2M}}. The manuscript does not state whether I_p is evaluated at a single OBC eigenvalue, averaged over the OBC spectrum, or computed from a theta-grid of the auxiliary resultant, and these different choices can change the root count and therefore the sign or zero of I_p. As a result, the phase maps in Figures 4, 9, and 10 are not reproducible from the text as it stands.","section":"Section IV B, Eq. (11)"},{"comment":"The manuscript itself concedes that the serrated patterns in the central I_p phase diagram 'stems from insufficient computational precision' and that reaching the shown precision 'requires a computational duration of 15 days.' Since Figure 4(a) is the primary evidence for the claim that I_p 'quantitatively discerns' left-, right-, and bipolar skin modes, the admitted numerical artifacts undermine the quantitative claim until convergence is demonstrated, the computations are repeated with a more robust method, or the claim is restricted to parameter regions that are demonstrably converged.","section":"Section IV B, Figure 4(a)"},{"comment":"The verification of I_p against real-space eigenstates is limited to a small number of marker points (pentagram, square, circle) in each parameter plane, and the paper provides no direct comparison between GBZ-predicted OBC spectra and exact OBC spectra. The auxiliary-GBZ construction of Section IV A is load-bearing for the entire I_p classification, yet its fidelity for this spinful quartic characteristic equation is nowhere checked by an exact spectrum overlay. A systematic agreement measure over the full parameter plane, rather than a handful of selected points, is needed to support the claimed correspondence between I_p and the real-space skin modes.","section":"Section IV C, Figures 5, 9, 10"},{"comment":"The abstract claims that the bipolar localization at the topological boundary is 'unaffected by size effects,' but the supporting evidence in Figure 7(c1) tracks only the most edge-localized eigenstate for system sizes up to N = 300. A size-scaling analysis of the full eigenstate distribution, or at least of the fraction of eigenstates exhibiting the bipolar profile, is needed to substantiate this claim.","section":"Section V B, Figure 7(c1)"}],"minor_comments":[{"comment":"The spinor notation for c^dagger_{x,y} and c_{x,y} is typeset in a garbled way; the components of the spinor and the Pauli-matrix structure should be written out explicitly.","section":"Section II, Eq. (1)"},{"comment":"The auxiliary-GBZ construction via the resultant of f(beta, E) and f(beta e^{i theta}, E) is introduced without derivation or a self-contained explanation; readers unfamiliar with Ref. [22] will not be able to follow the subsequent I_p computation.","section":"Section IV A"},{"comment":"The color scale for I_p is not defined in the figure or caption, so statements such as 'yellow area' and 'yellow-blue-middle area' cannot be interpreted quantitatively.","section":"Section IV B, Figure 4(a)"},{"comment":"The placeholder 'citevideo' appears in the sentence describing Figure 3(d)-(f); this should be replaced with a proper citation or reference to the supplementary video.","section":"Section III, text near Eq. (7)"},{"comment":"There are multiple typographical issues, including 'systematic-ally' in the abstract, 'y-direction cylinder-types' in Section V, and inconsistent notation for J_L and J_R in the text and figures.","section":"Various"},{"comment":"The phase diagrams in Figures 9 and 10 are described as reflecting 'rotation' of skin effects, but no quantitative measure of the degree of rotation is provided; the qualitative language should be aligned with the actual I_p values.","section":"Appendix C"}],"recommendation":"major_revision","confidential_remarks":"The Hopf-link braiding derivation is sound and could stand on its own, but the paper's headline claim of an I_p-based quantitative classification of skin modes is not yet established: the definition of I_p is underspecified, the central phase diagram is admitted to be numerically non-converged, and the validation is sparse. I recommend major revision rather than rejection, because these issues are fixable with a clearer specification, converged numerics, and a direct GBZ-versus-exact-OBC comparison. The paper might benefit from releasing code or data to support reproducibility, especially given the disclosed 15-day computation time."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the braiding computation is solid, and the paper is honest about crediting Ref. [52] for the nearest-neighbor Hopf-link result. What is genuinely new is the I_p polarization parameter, the real-space encoding cross-checks (C and Q), and the boundary-state degeneracy/IPR analysis. That new content is worth engaging, but I_p as currently defined does not support the claim that it \"quantitatively discerns\" left-, right-, and bipolar skin modes.\n\nI rechecked the exact part. Eq. (8) winds ±2 away from the EP surface Eq. (9), and the sigma_0 term (including 2 cos k_y) drops out of the traceless Hamiltonian, so the cylinder reduces to the 1D model of Ref. [52]. That is correct, and the paper openly says so. R(alpha,beta) is a clean ratio measure, and the cross-validation at selected sampling points (Figs. 4, 9, 10) is good practice. The boundary analysis at |J_R| = J_L, with half the spectrum purely real and size-robust bipolar localization, is internally consistent.\n\nThe stress-test note lands. Eq. (11) counts roots of Eq. (10), but Eq. (10) depends on E; the GBZ condition |beta_M| = |beta_{M+1}| gives a contour, not a finite 2M-element set. The manuscript never states whether I_p is evaluated at one OBC eigenvalue, averaged over the spectrum, or computed from a theta grid. Different choices can change the root count and flip the sign of I_p. That is a load-bearing ambiguity, because I_p is the headline new tool. The phase maps in Fig. 4(a) are explicitly admitted to be numerically unconverged (\"15 days\"), and the convergence problems sit exactly in the regions claimed to distinguish bipolar from one-way skin. There is no code or data, no direct GBZ-versus-exact-OBC spectrum comparison, and only a handful of real-space validation points. Also, k_y is fixed to 1 throughout, and the text has production artifacts (dangling \"citevideo,\" a URL-style reference) that should be cleaned.\n\nThese are serious but fixable issues. The exact topology part holds up; I_p needs a precise definition and better validation. The paper is not incoherent, and the authors engage prior literature, especially Ref. [52]. If I were the editor, I would send this to peer review rather than desk reject: a good referee can pin down the I_p definition and ask for code/data or a direct GBZ check. As is, I would not cite the I_p claim without redoing the numerics myself.","headline":"The exact topology part is correct but largely inherited from Ref. [52]; the new I_p skin classifier is under-specified and needs revision before I would trust its quantitative claim.","tokens_in":16939,"tokens_out":2272,"would_cite":false,"duration_ms":29270,"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":"A non-Abelian SU(2) gauge field can make a 2D non-Hermitian Hatano-Nelson cylinder show Hopf-link band braiding and bipolar skin localization using only nearest-neighbor x-direction hoppings.","keywords":["non-Hermitian skin effect","non-Abelian gauge field","Hatano-Nelson model","generalized Brillouin zone","Hopf-link braiding","exceptional point","bipolar skin effect","polarization parameter"],"falsifier":"Exact-diagonalize Hamiltonian (3) at $N=50$ on a fine grid in $(J_L,\\alpha)$ with $J_R=0.6$, $\\beta=\\pi/2$, compute the real-space center of mass $C(J_L,\\alpha)$ of all eigenstates, and overlay it on the $I_p$ map of Figure 4(a); any region where the two disagree on the skin direction (left, right, or bipolar) would refute the claim that $I_p$ quantitatively classifies the non-Hermitian skin effect.","tokens_in":15652,"feed_emoji":"🔗","tokens_out":8260,"duration_ms":87094,"temperature":0.7,"pith_summary":"The paper studies a two-dimensional Hatano-Nelson model whose leftward and rightward hoppings carry different SU(2) gauge phases, and it argues that these non-Abelian phases by themselves create two effects usually thought to require extra ingredients. First, the two complex energy bands of the periodic cylinder form Hopf-link braiding, with braiding degree $\\nu=\\pm2$, even though the model contains only nearest-neighbor couplings in the $x$-direction. Second, the open-boundary eigenstates can localize on both ends of the cylinder — a bipolar skin effect — even when the hopping amplitudes and phases are reciprocal, $J_L=J_R$ and $\\alpha=\\beta$. To capture that localization, the paper introduces a polarization parameter $I_p$ computed from the generalized Brillouin zone and checks it against directly encoded real-space eigenstates. If correct, the work shows that non-Abelian gauge fields are a standalone control knob for non-Hermitian topology and skin localization.","feed_headline":"Non-Abelian gauge creates Hopf-link bands and two-sided skin modes","feed_subtitle":"SU(2) phases turn nearest-neighbor hoppings into spectral braiding and controllable skin localization.","key_machinery":"The load-bearing object is the generalized Brillouin zone (GBZ) of the spinful lattice. For a finite chain with open boundaries, the GBZ is the locus of complex momenta $\\beta=e^{ik}$ selected by the condition $|\\beta_M|=|\\beta_{M+1}|$ on the roots of the characteristic equation, computed here through the auxiliary-GBZ resultant method; the directions in which the GBZ contour lies outside or inside the unit circle determine where eigenstates accumulate. The paper's new diagnostic is the polarization parameter $I_p = (N_>(R_d)-N_<(R_d))/(N_>(R_d)+N_<(R_d))$, the normalized imbalance between GBZ roots with modulus larger and smaller than the Brillouin-zone radius $R_d=1$. The non-Abelian gauge enters through the two Pauli-matrix phases $e^{i\\alpha\\sigma_y}$ and $e^{i\\beta\\sigma_x}$, whose non-commutation couples the two spin sectors and is the reason the characteristic equation acquires the structure that supports both Hopf-link braiding and the bipolar skin effect. The braiding degree $\\nu=\\int_0^{2\\pi}\\frac{dk_x}{2\\pi i}\\frac{d}{dk_x}\\ln(J_R^2\\sin^2\\beta\\,e^{-2ik_x}+J_L^2\\sin^2\\alpha\\,e^{2ik_x})$ supplies the topological bookkeeping.","core_discovery":"On the paper's own terms, the central discovery is that replacing the usual scalar or Abelian phases in a 2D Hatano-Nelson cylinder with non-commuting SU(2) phases ($\\sigma_y$ on leftward hops, $\\sigma_x$ on rightward hops) makes the complex energy bands braid into a Hopf link in $(\\mathrm{Re}E,\\mathrm{Im}E,k_x)$ space without long-range hoppings. The braiding degree is $\\nu=\\pm2$, and the transition between the two braiding types occurs at the exceptional-point surface $J_R^2\\sin^2\\beta = J_L^2\\sin^2\\alpha$, which is also the line on which the skin direction reverses. The paper further claims that the generalized Brillouin zone contains enough information to classify the skin effect quantitatively: the polarization parameter $I_p$, defined by the imbalance of GBZ roots outside and inside the unit circle, separates left-skin, right-skin, and bipolar skin regions, and it does so even at the reciprocal point $J_L=J_R$, $\\alpha=\\beta$, where ordinary non-reciprocity is absent. At the topological phase boundary half the eigenstates have purely real energies, and those zero-imaginary-energy states show pronounced degeneracy and a bipolar localization that persists as the system size grows.","pith_inferences":["One extension the paper does not pursue is using $I_p$ as a general-purpose diagnostic for any multiband non-Hermitian lattice; if the GBZ root counting is faithful there, it would replace expensive exact-diagonalization scans with a cheap scalar phase map.","The bipolar skin effect at $J_L=J_R$, $\\alpha=\\beta$ suggests a testable statement the authors leave implicit: the non-Abelian phases themselves act as an effective non-reciprocity, so tuning $\\alpha$ or $\\beta$ alone should rotate the skin direction continuously through the bipolar regime.","Because the paper's own phase maps show serrated, numerically unconverged regions, a natural next check is to recompute $I_p$ with higher-precision or root-polishing methods on a coarse grid; if the serration collapses into smooth boundaries, the bipolar regimes become sharper experimental predictions.","In a synthetic platform such as a topolectrical circuit, the predicted size-independent bipolar localization could be detected by measuring impedance or voltage profiles at the topological boundary, without needing time-resolved dynamics."],"forward_implications":["Hopf-link spectral braiding in this model requires no couplings beyond nearest-neighbor $x$-direction hoppings, so the effect should be realizable in the same kinds of short-range synthetic lattices used for Hatano-Nelson systems.","The single equation $J_R^2\\sin^2\\beta=J_L^2\\sin^2\\alpha$ gives a two-parameter phase diagram in which the braiding degree jumps between $\\nu=+2$ and $\\nu=-2$, providing a direct target for experimental phase-boundary searches.","The polarization parameter $I_p$ gives a scalar, momentum-space route to telling left-, right-, and bipolar skin regimes apart, which means the skin classification no longer depends only on inspecting real-space density profiles.","At the critical line $|J_R|=J_L$, half of the eigenstates become purely real and remain bipolarly localized as the system size grows, so the boundary itself carries robust, size-independent spatial structure that can be probed by dynamics and inverse participation ratio."],"supporting_citations":[{"why":"showed that non-Abelian gauge fields can produce Hopf-link spectral topology with nearest-neighbor hopping, the effect this paper extends to the 2D cylinder Hatano-Nelson model","marker":"[52]"},{"why":"supplies the auxiliary generalized Brillouin zone and resultant method used here to compute the GBZ from the characteristic equation","marker":"[22]"},{"why":"established that the GBZ contour and its deviation from the unit circle encode the non-Hermitian skin effect, the basis for the polarization parameter","marker":"[18]"},{"why":"ties the local radius of the GBZ to the direction of skin localization, which motivates counting roots inside versus outside the unit circle","marker":"[57]"},{"why":"introduced the GBZ framework for open-boundary spectra of non-Hermitian systems, on which the paper's boundary-condition analysis depends","marker":"[6]"},{"why":"defines the Hopf-link braiding degree used for the topological phase diagram, even though that work required longer-range couplings","marker":"[11]"}],"fun_headline_variants":["SU(2) gauge twists 2D Hatano-Nelson into Hopf-link spectrum","Non-Abelian phases drive Hopf-link braiding and bipolar skin effect","Hopf-link bands emerge from non-Abelian gauge in non-Hermitian lattice","Non-Abelian gauge yields Hopf-link bands and two-sided skin localization"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole skin classification rests on the assumption that the auxiliary-GBZ construction — the ordered-root condition $|\\beta_M|=|\\beta_{M+1}|$ applied to the quartic characteristic equation — faithfully predicts where open-boundary eigenstates sit, and that counting roots inside versus outside the unit circle tracks real-space localization even in parameter regions where the paper's own numerical maps are unconverged.","fun_headline_variants_meta":{"raw":{"variants":["SU(2) gauge twists 2D Hatano-Nelson into Hopf-link spectrum","Non-Abelian phases drive Hopf-link braiding and bipolar skin effect","Hopf-link bands emerge from non-Abelian gauge in non-Hermitian lattice","Non-Abelian gauge yields Hopf-link bands and two-sided skin localization"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00065,"raw_usage":{"total_tokens":3042,"prompt_tokens":1063,"completion_tokens":1979,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":679,"completion_tokens_details":{"reasoning_tokens":1891}},"tokens_in":679,"tokens_out":1979,"duration_ms":17478,"temperature":1.0,"reasoning_tokens":1891,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T05:06:39.656559+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Exact-diagonalize Hamiltonian (3) at $N=50$ on a fine grid in $(J_L,\\alpha)$ with $J_R=0.6$, $\\beta=\\pi/2$, compute the real-space center of mass $C(J_L,\\alpha)$ of all eigenstates, and overlay it on the $I_p$ map of Figure 4(a); any region where the two disagree on the skin direction (left, right, or bipolar) would refute the claim that $I_p$ quantitatively classifies the non-Hermitian skin effect.","supporting_citations":[{"cited_title":"Liang, D","cited_arxiv_id":null,"evidence_quote":"showed that non-Abelian gauge fields can produce Hopf-link spectral topology with nearest-neighbor hopping, the effect this paper extends to the 2D cylinder Hatano-Nelson model"},{"cited_title":"Zhang, Z","cited_arxiv_id":null,"evidence_quote":"supplies the auxiliary generalized Brillouin zone and resultant method used here to compute the GBZ from the characteristic equation"},{"cited_title":"Kawabata, K","cited_arxiv_id":null,"evidence_quote":"established that the GBZ contour and its deviation from the unit circle encode the non-Hermitian skin effect, the basis for the polarization parameter"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"ties the local radius of the GBZ to the direction of skin localization, which motivates counting roots inside versus outside the unit circle"},{"cited_title":"Ashiba, Z","cited_arxiv_id":null,"evidence_quote":"introduced the GBZ framework for open-boundary spectra of non-Hermitian systems, on which the paper's boundary-condition analysis depends"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"defines the Hopf-link braiding degree used for the topological phase diagram, even though that work required longer-range couplings"}],"review_version":1}