{"id":"b15b7e3d-0ea9-4a3b-8347-045c683decab","arxiv_id":"2506.08502","paper_version":1,"verdict":"REJECT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":1,"one_line_summary":"The displacement of a soliton after one nonlinear pumping cycle is claimed to equal the non-Abelian Chern number of the braiding nonlinear bands divided by the number of braiding bands.","lead":"This paper proposes a unified topological invariant, a non-Abelian Chern number of braiding nonlinear bands, to predict soliton displacement in nonlinear Thouless pumping. If correct, the formula D = C_NL divided by the number of braiding bands would unify integer and fractional nonlinear transport and guide the design of nonlinear topological devices.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (3) defines C_NL with a non-orthogonal frame; a U(N) rotation does not preserve normalization or the connection, so C_NL is not a gauge-invariant topological integer and Eq. (1) lacks a well-defined numerator.","rationale":"The reader's weakest assumption identifies the same load-bearing defect: the non-Abelian Berry connection in Eq. (3) is applied to non-orthogonal nonlinear eigenvectors, and the stated normalization argument does not supply a well-defined U(N) gauge structure. I make the failure more explicit: with S=<Phi_a|Phi_b> not equal to I, the connection transforms with i U^dagger S dU instead of i U^dagger dU, so the trace of the curvature is not gauge invariant. Since C_NL is the numerator of the central claim D=C_NL/N, this flaw is decisive for the paper's main equation as stated. The missing derivation of D from an adiabatic transport calculation is a second, independent gap, but the gauge problem alone is sufficient justification for the reader's rejection. The proposed test using the spectral projector would settle the issue computationally because the first Chern number of the rank-N subbundle is unambiguous and does not depend on choosing a non-orthogonal eigenvector frame.","tokens_in":11305,"tokens_out":9178,"duration_ms":112749,"concrete_test":"Reuse the g=4.0 three-site AAH data of Fig. 2(g): at each (k,t) build the spectral projector P(k,t) onto the three lowest nonlinear eigenstates. Compute the gauge-invariant first Chern number C_proj = (1/(2*pi*i)) integral_{BZ} Tr(P dP ^ dP) over the 2D torus, and separately compute C_NL from Eq. (3) with the paper's non-orthogonal frame. If C_proj differs from C_NL, or if C_NL changes under a smooth U(3) rotation U(k,t) applied to the three vectors with periodic boundary conditions, then Eq. (3) is not the topological invariant of the degenerate subspace and the numerator of Eq. (1) is not a well-defined quantized quantity.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central formula D=C_NL/N depends on C_NL being a genuine non-Abelian Chern number. For the N lowest nonlinear bands the paper takes eigenvectors Phi_a(k,t) that are normalized individually, ||Phi_a||=1, but not mutually orthogonal, so the overlap matrix S_ab=<Phi_a|Phi_b> differs from delta_ab. The connection in Eq. (3), [A_k]_ab=i<Phi_a|d_k Phi_b>, is not a U(N) connection: under Phi_a' = sum_b U_ba Phi_b one gets A_k' = U^dagger A_k U + i U^dagger S d_k U rather than U^dagger A_k U + i U^dagger d_k U, with an analogous t-component. The overlap term is not a pure gauge term unless S=I, so Tr F and the integrated (1/(2*pi*i)) integral Tr F are not invariant under smooth U(N) transformations; also the rotated vectors no longer have unit norm. Thus C_NL is not shown to be quantized or even a property of the degenerate subspace. No adiabatic derivation links the soliton center-of-mass displacement to Tr F/N; Eq. (1) is asserted rather than derived. The numerical matches in Figs. 2-4 and S2 therefore cannot certify Eq. (1) while the invariant entering it is frame-dependent.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a unified topological invariant for nonlinear Thouless pumping of solitons in one-dimensional lattices with Kerr nonlinearity. The central claim, Eq. (1), is that the center-of-mass displacement D of a soliton after one pump cycle equals C_NL/N, where C_NL is a non-Abelian Chern number of the lowest braiding (degenerate) nonlinear bands and N is the number of such bands. In the weak-nonlinearity limit the formula is claimed to reduce to the Abelian Chern number and give integer pumping; in the strong-nonlinearity limit, for braiding bands, it is claimed to give fractional pumping such as D = -1/2, -1/3, or 0. The paper supports this with nonlinear band-structure calculations and dynamical simulations for three-site, five-site, and seven-site Aubry-André-Harper models and for a nonlinear Rice-Mele model.","tokens_in":11507,"tokens_out":9824,"duration_ms":134117,"significance":"If valid, the proposed formula would be a valuable bridge between nonlinear band topology and soliton transport, unifying integer, zero, and fractional pumping in a single expression. The numerical work covers several lattice sizes and model families, and the predictions are concrete and falsifiable. However, the central definition of C_NL in Eq. (3) is not gauge invariant for the non-orthogonal nonlinear eigenvectors explicitly used in the paper, so the claimed topological meaning of the numerator of Eq. (1) is not established. Moreover, Eq. (1) is asserted rather than derived, and the numerical evaluation of C_NL is not specified in a way that can be reproduced or checked. These are load-bearing issues, not presentation concerns.","major_comments":[{"comment":"The definition of the non-Abelian Berry connection uses non-orthogonal eigenvectors: [A_k]_ab = i <Φ_a|∂_k Φ_b>. Under a change of frame Φ_a' = Σ_b U_ba Φ_b, one obtains A_k' = U† A_k U + i U† S ∂_k U, where S_ab = <Φ_a|Φ_b>. Since the paper explicitly states that the nonlinear eigenvectors are not orthogonal, S is not the identity and the extra term is not a pure gauge term. Consequently Tr F and the integral in Eq. (3) are not invariant under smooth U(N) transformations, and the rotated vectors are not in general normalized. The paper's statement that 'they still preserve U(N) gauge due to the normalization of eigenvectors' is therefore incorrect. C_NL as defined is not a gauge-invariant topological invariant, and its quantization is not established. This invalidates the topological meaning of the numerator in Eq. (1).","section":"Section 'Chern number of nonlinear bands', Eq. (3)"},{"comment":"Equation (1) is asserted without derivation. The text states that the pumping charge is the average nonlinear Chern number over the dimension of the degenerate subspace, but no adiabatic or projection argument connects the center-of-mass displacement of a single soliton initial condition to Tr F / N. In the linear non-Abelian case, one would need to specify carefully how a single filled band or a single initial state in an N-fold degenerate subspace relates to the trace of the Berry curvature; here the nonlinear self-consistency makes that relation even less immediate. The numerical agreement in Figs. 2-4 cannot certify Eq. (1) because the quantity C_NL entering it is not well-defined.","section":"Section 'Nonlinear AAH model', Eq. (1)"},{"comment":"The selection of the integer N and of which bands form the 'braiding' subspace is not governed by a clear, a priori criterion. For the five-site model, N=2 is used at g=1.6 and N=5 at g=3.0; for the seven-site model, N=3 at g=1.5 and N=5 at g=3.0. The paper does not provide an algorithmic rule for choosing N in terms of a spectral gap or a well-defined degeneracy structure, particularly because in the 'ill-defined' regions the lowest bands also touch loop bands. If N is chosen after inspecting the dynamical displacement D, the agreement between D and C_NL/N is partly circular.","section":"Section 'Nonlinear AAH model', Figs. 3 and S2"},{"comment":"The numerical computation of C_NL is not described in a reproducible way. The main text and the Supplement describe the modified Newton method for finding eigenvectors, but do not give the lattice discretization or the smooth-gauge procedure used to evaluate the integral in Eq. (3). Since Eq. (3) is gauge dependent for non-orthogonal frames, the reported values C_NL = -1, 0 cannot be independently checked. Reference [50] is a standard lattice-gauge formula for orthonormal tight-binding Bloch states; the paper does not explain how it is applied to non-orthogonal nonlinear eigenvectors.","section":"Section 'Chern number of nonlinear bands' and Supplemental Material"}],"minor_comments":[{"comment":"In the concluding paragraph, 'fictionalization' should be 'fractionalization'.","section":"Summary"},{"comment":"The phrase 'a system with linear interaction' appears to be a typo; in context it should likely read 'a system with nonlinear interaction' or 'with a staggered potential'.","section":"Section 'Nonlinear Rice-Mele model'"},{"comment":"The sign convention in H = H_lin - diag(g|Ψ|^2) and the self-consistency condition H(Ψ)Ψ = EΨ are not stated explicitly; adding one sentence would improve clarity.","section":"Eq. (2)"},{"comment":"The panels of Fig. 2(b) for g=1.5, 2.3, and 3.0 are not labeled individually, making it hard to connect the dynamical evolution to the phase diagram in Fig. 2(c).","section":"Fig. 2(b)"}],"recommendation":"reject","confidential_remarks":"The central mathematical definition of C_NL appears to be incorrect for the non-orthogonal frame used in the paper, and Eq. (1) is not derived. This is not a local presentation issue; rescuing the manuscript would require redefining the invariant (for example with an explicit orthogonalization of the degenerate subspace or a metric-corrected connection), redoing all numerical evaluations of C_NL, and providing a dynamical derivation of Eq. (1). I would not recommend revision in the present form, though a substantially rewritten manuscript with a well-defined invariant and a non-circular comparison might be worth reconsidering."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: the paper's central invariant, C_NL, is not a well-defined topological number. They define a non-Abelian Berry connection using nonlinear eigenvectors that are normalized but not orthogonal. Under a U(N) rotation, the connection picks up an overlap-matrix term that is not a pure gauge, so the trace of the curvature and its integral are not invariant. That breaks the claimed quantization and leaves Eq. (1) with an ill-defined numerator. This isn't a minor gap; it's the load-bearing math.\n\nThat said, the paper does something genuinely new. They propose a unified formula D=C_NL/N connecting soliton displacement to a non-Abelian Chern number of braiding nonlinear bands, and they test it across 3-, 5-, and 7-site AAH models and a Rice-Mele model. The numerical phase diagrams and fractional charges (-1/2, -1/3) are interesting, and the band-structure calculations using Newton downhill are thorough. The honest exclusion of ill-defined regions is a plus.\n\nBut the problems beyond the gauge issue: there's no adiabatic derivation linking the COM displacement to Tr F/N. The formula is asserted, then checked numerically at a few points. The choice of which bands form the degenerate subspace and the division by N appears post hoc. And there are no error bars or convergence checks on the Chern-number integrals. The ill-defined regions are shaded and set aside, which is fine, but it means the theory only applies where the numerics already work.\n\nIn short, this is a promising conjecture with attractive numerics, not a proof. The intended audience is people working on nonlinear topological pumping who want a design rule for fractional charges; they'll find the numerics suggestive but should not cite Eq. (1) as an established result. I'd send it to a careful referee—the claim is important and the flaw is subtle enough that a good report could redirect the authors to a proper formulation, e.g., via orthogonalization or a different invariant. But I wouldn't accept it in current form.","headline":"Central non-Abelian Chern number is not gauge invariant for non-orthogonal nonlinear bands, so the paper's main formula D=C_NL/N is a conjecture, not a proven result.","tokens_in":12131,"tokens_out":3498,"would_cite":false,"duration_ms":38052,"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":"This paper proposes a single topological invariant, the ratio of a non-Abelian Chern number to the number of braiding nonlinear bands, that it claims governs soliton displacement in nonlinear Thouless pumps across weak and strong…","keywords":["nonlinear Thouless pumping","soliton transport","non-Abelian Chern number","nonlinear bands","fractional pumping","Kerr nonlinearity","Aubry-André-Harper model","Rice-Mele model"],"falsifier":"Compute $C_{\\mathrm{NL}}$ for the five-site AAH model at $g=1.6$ after applying a smooth U(2) gauge rotation to the two degenerate eigenvectors at every $(k,t)$ point; if the value changes under that rotation, the invariant is gauge-dependent and Eq. (1) fails. Alternatively, run the full nonlinear time evolution in a well-defined braiding phase and check whether the center-of-mass displacement equals $C_{\\mathrm{NL}}/N$; a mismatch would falsify the formula.","tokens_in":11018,"feed_emoji":"⚛️","tokens_out":11819,"duration_ms":114652,"temperature":0.7,"pith_summary":"This paper claims that the quantized transport of solitons in a nonlinear Thouless pump is governed by one formula across all nonlinearity strengths: the displacement per cycle equals a non-Abelian Chern number divided by the number of braiding nonlinear bands, $D = C_{\\mathrm{NL}}/N$. In the weak-nonlinearity limit, where the bands stay separated, this reduces to the Abelian Chern number of the occupied nonlinear band, recovering integer pumping. In the strong-nonlinearity regime, where nonlinear Kerr loops expand and the lowest bands intertwine into a degenerate braiding subspace, the same formula yields fractional displacements such as $-1/2$ and $-1/3$. The authors test this invariant numerically in several Aubry-André-Harper and Rice-Mele lattices and report agreement with dynamical simulations. A reader who wants a unified topological picture of nonlinear soliton transport would find that here a single integer-like quantity organizes both the weak and strong regimes.","feed_headline":"One formula predicts integer and fractional soliton pumping","feed_subtitle":"Weak nonlinearity gives integer transport; strong nonlinearity gives fractional steps like -1/2 and -1/3.","key_machinery":"The non-Abelian Wilczek-Zee Berry connection on the degenerate nonlinear band subspace, with matrix elements $[A_k]_{ab} = i \\langle \\Phi_a | \\partial_k \\Phi_b \\rangle$, and the associated non-Abelian Chern number $C_{\\mathrm{NL}}$ obtained by integrating the trace of the Berry curvature over the $1+1$D Brillouin zone. This object replaces the Abelian Chern number of linear band theory: it is defined even when nonlinear eigenvectors are non-orthogonal, and the dimension $N$ of the degenerate subspace enters as the denominator in the pumping formula, encoding the braiding of bands.","core_discovery":"In the paper's own formulation, the discovery is a generalized TKNN expression for nonlinear soliton pumping: $D = C_{\\mathrm{NL}}/N$, where $C_{\\mathrm{NL}}$ is the non-Abelian Chern number of the lowest braiding nonlinear bands and $N$ is the number of those bands. The nonlinear bands are obtained by solving the instantaneous nonlinear Schrödinger equation via a modified Newton method; they are generically non-orthogonal and, at strong nonlinearity, form degenerate subspaces in which the standard Abelian Chern number is undefined. The paper defines a non-Abelian Berry connection on this subspace and computes $C_{\\mathrm{NL}}$ as the integral of the trace of the non-Abelian Berry curvature over the two-dimensional Brillouin zone. It shows that in well-separated bands this reduces to the Abelian Chern number, while for $N$ braiding bands it gives a fractional pumped charge; numerical simulations of the three-, five-, and seven-site Aubry-André-Harper models and the Rice-Mele model match the formula in the well-defined phases.","pith_inferences":["Beyond the paper, if $D = C_{\\mathrm{NL}}/N$ survives closer scrutiny, it implies a nonlinearity-tuned sequence of rational pumping values set purely by the size of the braiding multiplet, which would allow engineering fractional transport without interparticle interactions.","The gauge-invariance of the non-orthogonal eigenvector connection is the point most likely to need tightening; a biorthogonal or metric-corrected connection could preserve quantization even if the raw Wilczek-Zee trace does not.","The ill-defined region, where loops touch non-filled bands, may correspond to a non-adiabatic Landau-Zener regime; measuring how sharply the integer-to-fractional transition occurs could probe the band-touching structure the paper identifies."],"forward_implications":["In the weak nonlinear regime, the pumped displacement of a soliton is quantized to the integer Abelian Chern number of the occupied nonlinear band.","In the strong nonlinear regime, when $N$ lowest bands braid, the displacement is fractional, $C_{\\mathrm{NL}}/N$; the paper reports $-1/2$ for two braiding bands in the five-site model and $-1/3$ for three braiding bands in the seven-site model.","There is an intermediate ill-defined region where loop bands touch non-filled bands, the Chern number is undefined, and soliton pumping is non-adiabatic and unquantized.","The invariant applies to both the nonlinear Aubry-André-Harper model and the nonlinear Rice-Mele model, yielding a phase diagram with well-defined topological regions separated by ill-defined ones.","The paper states that the correspondence between the topological invariant and soliton displacement can be checked in ultracold-atom and photonic-waveguide experiments by tuning the nonlinearity and measuring the displacement."],"supporting_citations":[{"why":"Reports fractional soliton pumping in an AAH lattice at moderate nonlinearity, the phenomenon the paper's unified invariant aims to explain.","marker":"[33]"},{"why":"Introduces the Kerr-loop structures in nonlinear energy bands that the paper uses to define the braiding-band regime.","marker":"[47]"},{"why":"Supplies the Wilczek-Zee non-Abelian Berry connection used to build the invariant for degenerate nonlinear bands.","marker":"[49]"},{"why":"The generalized TKNN formula connecting pumped charge to Chern number, which Eq. (1) extends to the nonlinear braiding case.","marker":"[51]"},{"why":"Establishes the weak-nonlinearity nonlinear-band Chern number picture that the new invariant reduces to in the weak regime.","marker":"[46]"},{"why":"Demonstrates quantized pumping of a single soliton in an AAH model with Kerr nonlinearity, the physical setting the paper generalizes.","marker":"[32]"},{"why":"Provides the numerical lattice-gauge method the paper uses to evaluate the non-Abelian Chern number in the strong-braiding cases.","marker":"[50]"}],"fun_headline_variants":["Unified invariant unifies integer and fractional soliton pumping","Nonlinear soliton pumping: one invariant for all regimes","From Abelian to non-Abelian: a single topological formula","One invariant links weak and strong nonlinear pumping","Non-Abelian topology tunes soliton pump fractions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The formula stands on the claim that the non-Abelian Berry connection defined from the non-orthogonal nonlinear eigenvectors is a well-defined, gauge-invariant object with a quantized trace; if that trace is not a true topological integer, the numerator $C_{\\mathrm{NL}}$ in Eq. (1) loses its topological meaning.","fun_headline_variants_meta":{"raw":{"variants":["Unified invariant unifies integer and fractional soliton pumping","Nonlinear soliton pumping: one invariant for all regimes","From Abelian to non-Abelian: a single topological formula","One invariant links weak and strong nonlinear pumping","Non-Abelian topology tunes soliton pump fractions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00015,"raw_usage":{"total_tokens":1178,"prompt_tokens":905,"completion_tokens":273,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":521,"completion_tokens_details":{"reasoning_tokens":194}},"tokens_in":521,"tokens_out":273,"duration_ms":4187,"temperature":1.0,"reasoning_tokens":194,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T05:09:20.610476+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute $C_{\\mathrm{NL}}$ for the five-site AAH model at $g=1.6$ after applying a smooth U(2) gauge rotation to the two degenerate eigenvectors at every $(k,t)$ point; if the value changes under that rotation, the invariant is gauge-dependent and Eq. (1) fails. Alternatively, run the full nonlinear time evolution in a well-defined braiding phase and check whether the center-of-mass displacement equals $C_{\\mathrm{NL}}/N$; a mismatch would falsify the formula.","supporting_citations":[{"cited_title":"Jürgensen, S","cited_arxiv_id":null,"evidence_quote":"Reports fractional soliton pumping in an AAH lattice at moderate nonlinearity, the phenomenon the paper's unified invariant aims to explain."},{"cited_title":"Wu and Q","cited_arxiv_id":null,"evidence_quote":"Introduces the Kerr-loop structures in nonlinear energy bands that the paper uses to define the braiding-band regime."},{"cited_title":"Wilczek and A","cited_arxiv_id":null,"evidence_quote":"Supplies the Wilczek-Zee non-Abelian Berry connection used to build the invariant for degenerate nonlinear bands."},{"cited_title":"Thouless, M","cited_arxiv_id":null,"evidence_quote":"The generalized TKNN formula connecting pumped charge to Chern number, which Eq. (1) extends to the nonlinear braiding case."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the weak-nonlinearity nonlinear-band Chern number picture that the new invariant reduces to in the weak regime."},{"cited_title":"Jürgensen, S","cited_arxiv_id":null,"evidence_quote":"Demonstrates quantized pumping of a single soliton in an AAH model with Kerr nonlinearity, the physical setting the paper generalizes."},{"cited_title":"Fukui, Y","cited_arxiv_id":null,"evidence_quote":"Provides the numerical lattice-gauge method the paper uses to evaluate the non-Abelian Chern number in the strong-braiding cases."}],"review_version":1}