{"id":"7b15904d-534e-4ab9-9b8f-1e5911e8e975","arxiv_id":"2501.14169","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Non-Hermitian skin and topological modes can drive all oscillators in a coupled chain into a single-frequency synchronized state, shown in simulations and nonlinear topoelectrical circuits.","lead":"This paper shows that non-reciprocal coupling between nonlinear oscillators, known as the non-Hermitian skin effect, can make every oscillator in a chain lock to a single frequency no matter how the system starts. The authors built small electric circuits that reproduce this behavior and showed it survives weak disorder.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The min-IPR selection rule is the load-bearing link between non-Hermitian mode coupling and the predicted final state, but Section 2 only asserts it: the effective Hamiltonian is complex, not a Lyapunov function, and the stated gradient form mismatches Eq. (1) by a factor 1/2.","rationale":"The reader's weakest assumption and my concern coincide: the effective-potential/min-IPR argument is the unproved load-bearing link. I agree with CONDITIONAL because the empirical simulations and two circuit demonstrations are independent evidence that some robust single-frequency synchronized state exists; what is at stake is the paper's proposed mechanism and its predictive specificity, not the existence of the phenomenon. If the fixed-point enumeration finds stable non-minimal-IPR synchronized states, then the stated selection rule is false and the central claim as a mechanism would need major revision. If it confirms the min-IPR states are the only stable single-frequency attractors, the concern is resolved. The derivative-factor mismatch in the effective Hamiltonian is easily corrected, but it illustrates that the variational argument was not checked; it should be fixed or removed. I therefore keep the reader's CONDITIONAL verdict, with the additional explicit check as the condition.","tokens_in":16494,"tokens_out":9416,"duration_ms":88432,"concrete_test":"Enumerate all single-frequency synchronized solutions of Eq. (1) for the N=15 parameters of Fig. 1 by writing Z_l(t)=a_l e^{iωt} and solving the resulting 2N real algebraic equations for (ω, Re a_l, Im a_l) with multi-start Newton from every linear eigenmode and random phase initial guesses. For each converged solution, compute the normalized IPR and linear stability. If any stable synchronized solution has IPR above the value of the two reported in-phase/anti-phase modes, the minimal-IPR selection rule is refuted; if the only stable solutions are exactly those two modes, the heuristic is supported. As a control, integrate Eq.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The decisive step is the claim in Section 2 that, because eigenmodes of the Hatano-Nelson chain are non-orthogonal, 'the system is expected to evolve into the linear eigenstate with the minimal effective potential — namely, the eigenstate with the minimal IPR.' This is the only argument connecting non-Hermitian mode coupling to the predicted final state, and it is not a derivation. The effective Hamiltonian H = Σ 0.5(−α|Z_l|^2 + 0.5β|Z_l|^4) − i Σ J_kl Z_l^*Z_k is complex; the dynamical flow Ż = −∂H/∂Z^* is therefore not a gradient descent on a real potential, and the coupling term is neither positive-definite nor sign-definite. In addition, the stated H does not reproduce Eq. (1): with the standard Wirtinger derivative, −∂H/∂Z_l^* yields (iω0 + α/2 − β|Z_l|^2/2)Z_l − i(J_+Z_{l+1}+J_-Z_{l−1}), so the gain and nonlinearity are too small by a factor 2. Even ignoring this mismatch, the quartic term 0.25βΣ|Z_l|^4 favors small IPR only at fixed Σ|Z_l|^2, but the gain term changes that norm during evolution. Thus the minimal-IPR selection rule is a heuristic, and the paper's prediction of which mode wins, and hence which spatial profile and frequency appear, rests entirely on this rule.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes and experimentally demonstrates that in chains of Stuart-Landau oscillators with non-reciprocal couplings, the dynamics converge to a single-frequency synchronized state whose spatial profile matches the minimum-IPR eigenstate of the linear Hatano-Nelson or non-Hermitian SSH model, for essentially all initial conditions over a range of parameters and weak disorder. The authors call this non-Hermitian global synchronization, identify linear skin-state, nonlinear skin-state, topological, and skin-topological synchronized regimes, and support the claims with ODE simulations, disorder statistics, phase diagrams, and circuit experiments for N=5 versus N=9 (skin) and N=9 versus N=15 (topological).","tokens_in":16905,"tokens_out":5539,"duration_ms":50357,"significance":"If the proposed mechanism is correct, the work offers a novel and appealing route to robust global synchronization that is insensitive to initial conditions and scalable in size, and the experimental realization in nonlinear topoelectrical circuits is a substantial contribution. The paper contains nontrivial, falsifiable predictions: the final oscillation frequency and spatial profile are compared with independently computed minimum-IPR eigenstates, and the size-dependent onset of synchronization (N=5 vs N=9 and N=9 vs N=15) was predicted before the circuit measurements. These are genuine strengths. The central limitation is that the theoretical selection rule identifying the min-IPR eigenstate as the global attractor is asserted rather than derived, and the experimental evidence, while qualitatively consistent, covers only a small number of samples.","major_comments":[{"comment":"The statement that the system evolves into the linear eigenstate with the minimal effective potential, namely the minimal-IPR eigenstate, is not supported by the equations as written. The effective Hamiltonian H = sum 0.5(-alpha|Z_l|^2 + 0.5 beta|Z_l|^4) - i sum J_kl Z_l^* Z_k is complex, so it is not a real potential or Lyapunov function for the flow Z_dot = -dH/dZ*. Moreover, with standard Wirtinger derivatives this H does not reproduce Eq. (1): the gain and nonlinearity terms come out with coefficients alpha/2 and beta/2 rather than alpha and beta. Thus the 'minimal effective potential' argument, which is the only link between non-Hermitian mode non-orthogonality and the predicted final state, is a heuristic at this point. This matters because the paper's central predictions, namely which eigenstate wins and therefore which frequency and spatial profile appear, rest entirely on this rule.","section":"Section 2 (effective Hamiltonian paragraph)"},{"comment":"Even if the effective Hamiltonian were corrected, the claim that the quartic term sum 0.25 beta |Z_l|^4 favors the minimal-IPR state requires an implicit fixed-norm assumption. During the evolution, the Stuart-Landau gain term alpha|Z_l|^2 changes the total norm, so the comparison of IPR at fixed norm is not valid. The paper does not provide a derivation of gradient-descent-like dynamics on the mode coefficients C_n(t), nor a numerical test showing that these coefficients follow such a landscape. I would like to see either an analytical argument (for example, a Lyapunov function or an adiabatic elimination of the mode amplitudes) or direct numerical evidence that the expansion coefficients evolve toward the minimal-IPR mode in a way that depends on an effective potential. Without this, the mechanism remains a conjecture even though the simulations and experiments are consistent with it.","section":"Section 2 (convergence to min-IPR)"},{"comment":"The order parameter R_o is defined as max[R(t)] - min[R(t)] for t>t0, which was chosen after observing that the two attractors are in-phase and anti-phase. Because R_o is calibrated to the two observed synchronized states, the phase diagrams in Figs. 1f, 1g, 3j, and 3k partly codify the classification rather than independently detecting synchronization. This does not invalidate the reported profile and frequency matches, but it weakens the claim that the diagrams demonstrate a synchronization transition without prior knowledge of the attractors. Please either justify R_o on independent grounds (for example, by showing it is equivalent to a standard measure after a sublattice rotation) or present an additional order parameter that detects anti-phase synchronization in a principled way.","section":"Section 2 (order parameter R_o)"}],"minor_comments":[{"comment":"There is a typo in the sentence 'the system is excepted to evolve' which should read 'expected to evolve'.","section":"Section 2"},{"comment":"The definition IPR(epsilon) = sum |phi_l(epsilon)|^4 should specify the normalization convention for the eigenstates, for example sum |phi_l|^2 = 1, since the IPR value depends on the normalization.","section":"Section 2 (IPR definition)"},{"comment":"The caption states 'at the time marked by blue lines in (b1)-b(2)' but the blue lines are dashed vertical lines in the main panels; the notation is inconsistent and should be cleaned up.","section":"Figure 1 caption"},{"comment":"The imaginary unit is denoted j in the circuit equations but i in Eq. (1); while this is a common convention, it should be stated explicitly to avoid confusion.","section":"Equations (3)-(4)"},{"comment":"The factor-of-two mismatch between the effective Hamiltonian and Eq. (1) may be a typographical error, but it should be corrected and the derivation rechecked so that the gradient-descent statement is at least dimensionally consistent.","section":"Section 2 (effective Hamiltonian)"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely to be of interest to the journal's readership, and the experimental demonstrations are valuable. The main concern is that the min-IPR selection rule, which is the core theoretical mechanism, is not derived and the effective potential argument is internally inconsistent as presented. This is fixable within the manuscript's scope by either supplying a rigorous derivation or substantially reframing the theoretical claims as an empirical rule supported by numerics and experiments. I would also encourage the authors to provide more extensive experimental repetitions with varied initial conditions, given the strength of the initial-state-immunity claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this paper claims something genuinely new: a non-reciprocal Stuart-Landau chain whose dynamics converge to a single global oscillation regardless of initial conditions, selected by the non-Hermitian skin effect. Prior work did boundary or phase synchronization; this targets the bulk and shows size-dependent onset, with a topological zero-mode version. The predictions are not back-fitted: frequencies and spatial profiles are compared to independently computed linear eigenstates, and the N=5 vs N=9 vs N=15 behavior was predicted before the circuit measurements. The topological circuit with pseudospins to get real-valued couplings is a solid piece of experimental work, even if the samples are small.\n\nThe soft spots are real but not fatal. The central theoretical claim—that the system evolves to the minimum-IPR eigenstate—rests on a single paragraph invoking an effective potential. That argument is not a derivation. The effective Hamiltonian is complex, so it is not a Lyapunov function; the quartic term selects delocalized states only at fixed norm, and the gain term changes the norm. Worse, as written, Ż = −∂H/∂Z* gives 0.5α and 0.5β, not α and β from Eq. (1), so the Hamiltonian does not reproduce the equations of motion. That is a concrete technical error in the supporting argument, not just a missing proof. The simulations and experiments may still be right, but the \"why\" is only a heuristic. The order parameter R_o is defined after seeing the two attractors, so it is a classification tool rather than a predictive measure; that is acceptable but should be acknowledged. Experiments use one initial condition per circuit and no repeated samples or error bars; the qualitative match carries the weight.\n\nFor a reader in non-Hermitian photonics or laser arrays, this is worth engaging with: the phenomenon is plausible, the circuit demonstration is the first of its kind, and the size-dependent transition from unsynchronized to synchronized is an interesting effect. The paper deserves a serious referee: the flaws are in the theoretical framing, not in the core observation. I would send it to review, asking for the selection rule to be either proven or explicitly reframed as a numerical observation, and for the effective-Hamiltonian mismatch to be fixed.","headline":"Plausible new mechanism for global oscillator synchronization via skin-effect mode convergence, with real circuit evidence, but the min-IPR selection rule that carries the theory is a heuristic, not a derivation.","tokens_in":17385,"tokens_out":3452,"would_cite":true,"duration_ms":32909,"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":"Non-reciprocally coupled Stuart-Landau oscillators globally synchronize through the non-Hermitian skin effect, independent of initial conditions and weak disorder.","keywords":["global synchronization","non-Hermitian skin effect","Stuart-Landau oscillators","non-reciprocal coupling","topological zero modes","topoelectrical circuits","inverse participation ratio","Kuramoto order parameter"],"falsifier":"Start a chain from an initial condition whose overlap with the two most concentrated patterns is exactly zero, for example by exciting only a weakly concentrated pattern. The paper's initial-state-independence claim predicts that this trajectory must still lock to the frequency and spatial profile of the most concentrated pattern; observing a multi-frequency steady state, or locking to a different pattern, would refute the central mechanism.","tokens_in":16276,"feed_emoji":"🔄","tokens_out":15973,"duration_ms":126167,"temperature":0.7,"pith_summary":"The paper claims that the non-Hermitian skin effect—the boundary pile-up of eigenstates in chains with unequal left/right couplings—can be converted into a synchronization mechanism. In a Hatano-Nelson chain (a one-dimensional lattice with unequal left/right couplings) of Stuart-Landau oscillators (self-sustained nonlinear oscillators), the dynamics is reported to converge, for one thousand random initial states, to one of two single-frequency oscillations whose frequencies match the eigenenergies of the two most concentrated linear eigenstates. The proposed selection rule is that the nonlinear effective potential is minimized by the eigenstate with the smallest inverse participation ratio, $\\mathrm{IPR} = \\sum_l |\\varphi_l(\\varepsilon)|^4$, and the strong non-orthogonality of skin-localized modes lets that state absorb all the others. The same rule is shown to extend to the zero-energy midgap state of a non-Hermitian SSH chain (a two-band lattice with alternating couplings), producing topological global synchronization, and to nonlinear eigenstates in longer chains. The authors validate the mechanism in fabricated nonlinear topoelectrical circuits, with a 9-site circuit locking at 23.5 kHz to a skin mode and a 15-site circuit locking at 92.6 kHz to a topological zero mode.","feed_headline":"All oscillators lock to one frequency with non-reciprocal couplings","feed_subtitle":"Non-reciprocal couplings make every site oscillate at one shared frequency, confirmed in 9- and 15-site circuits.","key_machinery":"The load-bearing object is the minimum-IPR selection rule: the eigenstate with the smallest inverse participation ratio, $\\mathrm{IPR} = \\sum_l |\\varphi_l|^4$, is the state that minimizes the nonlinear part of the effective Hamiltonian and therefore becomes the global synchronized state. The machinery has three parts. First, non-reciprocal couplings ($J_+\\neq J_-$) in a Hatano-Nelson chain produce the non-Hermitian skin effect, localizing every eigenstate at one boundary and breaking the orthogonality of the eigenbasis; the biorthogonal right/left eigenstate expansion then shows that initially excited modes couple strongly. Second, the effective Hamiltonian $H = \\sum_l 0.5(-\\alpha|Z_l|^2 + 0.5\\beta|Z_l|^4) - i\\sum_{kl}J_{kl}Z_l^*Z_k$ assigns a lower potential to more concentrated states, so the dynamics should transfer population to the minimum-IPR mode. Third, in the topological extension, the non-Hermitian SSH chain's zero-energy midgap state can become the minimum-IPR state for suitable intercell coupling, making the topological mode the global attractor; when the linear expansion breaks down, nonlinear eigenstates play the same role.","core_discovery":"The central discovery is that non-reciprocal coupling makes collective synchronization a property of the lattice's most concentrated eigenstate rather than of carefully prepared initial conditions. For a Hatano-Nelson chain (a one-dimensional tight-binding chain with unequal hopping amplitudes $J_+\\neq J_-$), all eigenstates are localized at one boundary, and the authors find that the long-time dynamics always settles into the linear eigenstate with the minimal inverse participation ratio, $\\mathrm{IPR} = \\sum_l |\\varphi_l(\\varepsilon)|^4$, with the common oscillation frequency given by $|\\varepsilon+\\omega_0|$. The explanation is that the effective Hamiltonian $H = \\sum_l 0.5(-\\alpha|Z_l|^2 + 0.5\\beta|Z_l|^4) - i\\sum_{kl} J_{kl}Z_l^*Z_k$ has a nonlinear potential term $0.25\\beta\\sum_l|Z_l|^4$ that is smallest for the most concentrated state, while the skin effect's breakdown of eigenstate orthogonality couples all initially excited modes so the system drains into that minimum. The same minimum-IPR logic is then applied to a non-Hermitian SSH chain, where tuning the intercell coupling can make the midgap topological zero mode the minimum-IPR state and thereby synchronize the whole lattice; in longer chains the role passes to nonlinear eigenstates obtained numerically.","pith_inferences":["Beyond the paper's examples, the minimum-IPR rule suggests a general design recipe: to make an array synchronize at a target frequency and spatial profile, engineer the non-reciprocal lattice so that the target eigenstate has the smallest IPR; mode competition should then select it automatically.","The same logic should transfer to other non-Hermitian lattices with strongly non-orthogonal eigenmodes, such as disordered or quasiperiodic chains and two-dimensional skin-effect geometries; direct simulations of those systems would test the generality of the mechanism.","The transition regions where synchronization fails, and the observation that disorder can partially rescue synchronization there, imply that the basin of attraction is controlled by effective mode-coupling strength rather than by nonlinearity alone, which could be tested by tuning coupling asymmetry at fixed gain.","For laser arrays, the paper's picture points to a concrete engineering target: non-reciprocal coupling alone could enforce single-mode, single-frequency operation in a large array without external injection locking or careful pump shaping."],"forward_implications":["With the paper's parameters, any of one thousand random initial states of a 15-site Hatano-Nelson chain converges to one of two single-frequency oscillations, so global synchronization is achieved without preparing a particular initial condition.","Weak disorder in the onsite frequencies (up to $W=0.1\\omega_0$) leaves the synchronization probability at 100 percent in the linear-skin region and above 90 percent in the nonlinear-skin region, so the effect tolerates structural perturbations.","System size acts as a control knob: longer chains first turn on linear skin-state synchronization, then pass through an unsynchronized transition region, and then recover synchronization through nonlinear skin states.","In the non-Hermitian SSH chain, tuning the intercell coupling selects among linear skin-state, nonlinear skin-state, topological, and skin-topological cluster synchronization, with the midgap zero mode synchronizing the full lattice when it has the minimum IPR.","Fabricated circuits confirm the two headline cases: a 9-site circuit shows single-frequency voltage oscillations at 23.5 kHz matching the minimum-IPR skin mode, and a 15-site circuit synchronizes at 92.6 kHz matching the topological zero mode."],"supporting_citations":[{"why":"Defines the non-Hermitian skin effect as boundary localization whose number of states scales with volume; the paper's mode-selection mechanism depends on this effect.","marker":"[30]"},{"why":"Supplies the original non-reciprocal chain (Hatano-Nelson) that the paper uses as its base lattice.","marker":"[35]"},{"why":"Provides the biorthogonal right/left eigenstate formalism used to expand the dynamics and quantify non-orthogonality.","marker":"[38]"},{"why":"Earlier demonstration that non-reciprocal phase transitions can drive synchronization; the paper positions its result as a distinct, global form of synchronization.","marker":"[27]"},{"why":"Shows how the skin effect can shape topological in-gap states, the mechanism behind the paper's topological global synchronization.","marker":"[41]"},{"why":"Defines the Kuramoto order parameter that the paper extends (as $R_o$) to quantify both in-phase and anti-phase synchronization.","marker":"[9]"},{"why":"A Hatano-Nelson laser array exhibiting non-Hermiticity and nonlinearity, cited as the main application platform for the proposed mechanism.","marker":"[74]"},{"why":"Establishes the circuit-lattice mapping that justifies using nonlinear topoelectrical circuits to realize the lattice models experimentally.","marker":"[55]"}],"fun_headline_variants":["Non-Hermitian physics makes all oscillators lock together","Skin effect synchronizes all oscillators in a chain","Non-reciprocal couplings force one shared frequency","Topological zero mode synchronizes entire lattice","Robust global sync from non-Hermitian skin effect"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes, without proof, that the nonlinear terms create a genuine energy landscape whose lowest point is always the most concentrated vibration pattern, and that the interactions between patterns reliably push every starting condition into that lowest point.","fun_headline_variants_meta":{"raw":{"variants":["Non-Hermitian physics makes all oscillators lock together","Skin effect synchronizes all oscillators in a chain","Non-reciprocal couplings force one shared frequency","Topological zero mode synchronizes entire lattice","Robust global sync from non-Hermitian skin effect"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000407,"raw_usage":{"total_tokens":2134,"prompt_tokens":984,"completion_tokens":1150,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":600,"completion_tokens_details":{"reasoning_tokens":1089}},"tokens_in":600,"tokens_out":1150,"duration_ms":7345,"temperature":1.0,"reasoning_tokens":1089,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T15:19:33.305692+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Start a chain from an initial condition whose overlap with the two most concentrated patterns is exactly zero, for example by exciting only a weakly concentrated pattern. The paper's initial-state-independence claim predicts that this trajectory must still lock to the frequency and spatial profile of the most concentrated pattern; observing a multi-frequency steady state, or locking to a different pattern, would refute the central mechanism.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the non-Hermitian skin effect as boundary localization whose number of states scales with volume; the paper's mode-selection mechanism depends on this effect."},{"cited_title":"Hatano and D","cited_arxiv_id":null,"evidence_quote":"Supplies the original non-reciprocal chain (Hatano-Nelson) that the paper uses as its base lattice."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the biorthogonal right/left eigenstate formalism used to expand the dynamics and quantify non-orthogonality."},{"cited_title":"Fruchart et al., Nature 2021, 592, 363-369","cited_arxiv_id":null,"evidence_quote":"Earlier demonstration that non-reciprocal phase transitions can drive synchronization; the paper positions its result as a distinct, global form of synchronization."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows how the skin effect can shape topological in-gap states, the mechanism behind the paper's topological global synchronization."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the Kuramoto order parameter that the paper extends (as $R_o$) to quantify both in-phase and anti-phase synchronization."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"A Hatano-Nelson laser array exhibiting non-Hermiticity and nonlinearity, cited as the main application platform for the proposed mechanism."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the circuit-lattice mapping that justifies using nonlinear topoelectrical circuits to realize the lattice models experimentally."}],"review_version":1}