{"id":"89f874ef-8e2c-4bac-9d79-47a965f58156","arxiv_id":"2412.12517","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A voltage bias can drive a spin spiral through rigid rotation, quasi-periodic motion, and chaos, with the rotation slowing down as the wire grows.","lead":"Spin spirals in a one-dimensional conductor, when a voltage is applied, can rotate rigidly, wobble quasi-periodically, or tumble chaotically. This voltage control could let electrical currents manipulate spin polarization in nanoscale magnetic devices without external fields.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Thermodynamic-limit claim rests on a frozen nonequilibrium occupation: time-dependent b± could alter the conical-state fate.","rationale":"The reader's verdict of CONDITIONAL with moderate confidence is appropriate, and the identified weakest assumption (adiabatic separation) is real. My stress-test agrees that the electron-distribution treatment is the most delicate part of the argument, but I single out a more specific and self-admitted gap: the infinite-system calculation freezes the nonequilibrium occupation factors b± at a single time, even though the text concedes they should be time-dependent. This is the most load-bearing concern because the 'static conical state persists to arbitrarily long times' claim depends directly on this infinite-system model, while the finite-size t1 ∝ L^1.75 scaling only establishes metastability in finite chains. The finite-size RR/QP/CP phase diagram (Fig. 2 and Fig. 3) is supported by direct integration of Eq. (4) and is less affected by this issue. The proposed self-consistent b± update is a concrete, feasible computation that would either validate or invalidate the thermodynamic-limit conclusion. I therefore recommend keeping the reader's CONDITIONAL verdict unchanged, with the condition explicitly including this self-consistency check. The paper also has independent support in the form of direct time integration and explicit simulations at multiple L values, which I credit; the concern here is about the extrapolation to the infinite system, not about the internal consistency of the finite-size numerics.","tokens_in":20369,"tokens_out":13201,"duration_ms":128494,"concrete_test":"Recompute the infinite-system dynamics self-consistently: at each integration step of Eq. (S24), re-extract b± (or the full n+(k)) from the instantaneous uniform spin orientation (m∥, mz, φ) using the scattering calculation of Fig. S7, and compare the trajectory against the frozen-b± result. If after the transient mz remains constant for times far exceeding t1 (e.g., 10^4 t0 at the same parameters), the static-conical-state claim holds; if mz drifts, oscillates, or shows slow secular growth, the claim fails. A complementary finite-size check is to simulate L ≈ 2000 and monitor bulk mz for several t1; a deviation from the frozen-b± prediction would corroborate the concern.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim that a static conical state persists to arbitrarily long times in the thermodynamic limit is established in the Supplemental via the 'infinite system' calculation. There, the electronic occupation n+(k) is assumed to have the double-step form of Eq. (S25) with amplitudes b± extracted once, at t = 500 w^-1, from a finite-size simulation (Fig. S7), then held fixed while integrating Eq. (S24). The Supplemental explicitly acknowledges that 'in principle, since the spin texture changes in time, the values b± should also be taken as time-dependent.' This is not a minor technicality: b± are determined by contact scattering, which depends on the instantaneous spin texture. As mz grows from 0 to a finite value during the transient, the scattering changes; if b± continue to drift on timescales comparable to t1, the torque in Eq. (S24) changes and the conical state could slowly evolve (e.g., mz keeps increasing or oscillates) rather than remain stationary. The finite-size result t1 ∝ L^1.75 (Fig. S7b) only shows that in finite chains the metastable conical state survives until the RR instability sets in; it does not prove that the exact infinite-system fixed point (with time-dependent b±) is a stable attractor. Moreover, the red dashed curve in Fig. 6 is generated with the same frozen b± used to model the finite system, so the agreement in Fig. 6 is partly circular and cannot by itself validate the infinite-time persistence of the conical state.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a one-dimensional chain of classical Heisenberg spins exchange-coupled to itinerant electrons, driven out of equilibrium by a voltage bias between two reservoirs. Using the non-equilibrium Green's function formalism for the electron steady state and a Landau-Lifshitz-type equation with a phenomenological damping term for the spins, the authors report three dynamical regimes as the bias increases: a rigidly rotating (RR) spiral, a quasi-periodic (QP) phase, and a chaotic phase (CP). They further show that the RR phase produces a finite average spin current, analyze the linear-response precession vector Omega via Eq. (7), and derive the asymptotic scalings |Omega|/deltaV ~ L^-1 in region I and |Omega|/deltaV ~ eta^{-1} L^-3 and eta L^-3 in regions II and III. Based on this scaling and on a finite-size/infinite-system comparison, the paper claims that in the thermodynamic limit the spin texture freezes into a static conical state that persists to arbitrarily long times.","tokens_in":20613,"tokens_out":7577,"duration_ms":79762,"significance":"If the central claims hold, the paper demonstrates a new mechanism by which a voltage bias alone can destabilize spiral magnetic order and generate collective spin dynamics, including a time-crystalline-like rigid rotation and eventual chaos, with experimentally accessible spin-current signatures. The qualitative three-regime picture is supported by direct numerical integration of Eq. (4), and the phase classification is made explicit through criteria such as Eq. (S19) and transport fingerprints. The linear-response analysis leading to Eq. (8) involves no fitted constants. However, the two strongest claims - the L^-3 scaling and the infinite-lifetime conical state - rest on a linearized equation that is not derived in the manuscript or the Supplemental Material, and on an infinite-system calculation with frozen occupation steps. These issues are load-bearing and need to be addressed before the quantitative conclusions can be accepted.","major_comments":[{"comment":"The central scaling law (8) is obtained from Eq. (7), but Eq. (7) is stated without a derivation. The text says that the partial derivatives in Eq. (6) can be computed explicitly, but the actual linearization of Eq. (4) around the equilibrium spiral, the explicit form of A, the definition of the vector or scalar omega, and the treatment of the constraint |m_i|=1 are not given. The Supplemental Material does not contain this derivation either; it jumps to the infinite-system calculation. Since Eq. (8) is one of the main quantitative results, the full derivation of Eq. (7) must be provided, including how the numerical derivatives with respect to m_j and deltaV are evaluated and how the resulting linear system is solved without introducing additional approximations.","section":"Linear response regime, Eq. (7)"},{"comment":"The conclusion that a static conical state persists to arbitrarily long times in the thermodynamic limit relies on the infinite-system calculation in which the occupation steps b+ and b- are extracted once at t = 500 w^{-1} from a finite-size simulation and then held fixed while integrating Eq. (S24). The Supplemental Material explicitly acknowledges that 'in principle, since the spin texture changes in time, the values b± should also be taken as time-dependent.' This is not a minor caveat: b± are determined by contact scattering, which depends on the instantaneous spin texture and on mz. If b± drift on timescales comparable to t1, the torque in Eq. (S24) changes and the conical state could slowly evolve rather than remain stationary. The finite-size result t1 ~ L^1.75 shows that the metastable state survives until the RR instability in finite chains, but it does not prove that the exact infinite-system fixed point is stable in the fully self-consistent dynamics. Moreover, the red dashed curve in Fig. 6 is generated with b± extracted from the same L=1000 simulation, so the agreement in Fig. 6 is partly a consistency check of the tilted-spiral ansatz, not an independent validation of infinite-time persistence. The authors should either compute b± self-consistently as a function of the instantaneous spin configuration, or provide a quantitative bound on their time dependence and a stability argument for the frozen-occupation fixed point.","section":"Infinite system and Supplemental Material, Eq. (S25)"},{"comment":"The finite-size scaling underlying Eq. (8) uses system sizes selected at the centers of the beating-period intervals (the red dots in Fig. S2(b)). The Supplemental Material states that this selection is made to obtain smoother dependence because physical properties jump when two additional beating periods fit into the chain. Since Fig. 5(b) uses the same selected values of L, the apparent L^-3 collapse could be partly an artifact of the commensuration selection. The authors should present the scaling for generic values of L, or at least quantify the spread around the selected centers, in order to establish that the L^-3 law is not a consequence of choosing special system sizes.","section":"Supplemental Material, Fig. S2(b) and Fig. 5(b)"},{"comment":"The entire numerical scheme assumes the adiabatic separation tau_e << tau_m for the electron and spin timescales. The paper argues a posteriori that tau_m grows as L^3 in the RR phase, which justifies the approximation for large L. However, the QP and CP phases are studied at finite L and finite bias, and no estimate of tau_m in those regimes is given. Since the three-regime phase diagram and the corresponding transport predictions are central claims, the authors should provide a quantitative check of the adiabatic condition (for example, a comparison of the largest spin-precession frequency with the electron relaxation or dwell time) or state clearly where the approximation is expected to break down.","section":"Methodology and phase diagram"}],"minor_comments":[{"comment":"The expression for Omega_i in Eq. (S18) appears to be a scalar, while the text refers to it as the common rotating vector Omega_i. Please provide the correct vector form; this is important because the criterion Eq. (S19) compares these quantities.","section":"Supplemental Material, Eq. (S18)"},{"comment":"The notation for the precession vector is inconsistent: Eq. (5) uses Omega, while Eq. (7) uses omega. It should be made explicit whether omega is a scalar precession frequency or the vector Omega, and the algebra in Eq. (7) should be written with unambiguous vector notation.","section":"Linear response regime, Eqs. (5) and (7)"},{"comment":"The claim that Fig. 3 is representative of the thermodynamic limit is based on the convergence of deltaVc1 for three selected values of eta. It would be useful to comment on the convergence of the other phase boundaries, especially the RR-QP boundary in the small-eta region where the authors state that the asymptotic behavior is difficult to characterize.","section":"Supplemental Material, Fig. S4(b)"},{"comment":"The re-entrant behavior of the RR phase near eta/w ~ 0.15 is mentioned but not discussed. A brief explanation, or at least a statement of whether this feature persists in the thermodynamic limit, would help the reader interpret the phase diagram.","section":"Phase diagram, Fig. 3"},{"comment":"The statement that 'even an infinitesimal bias can induce non-trivial dynamics' should be reconciled with the thermodynamic-limit result that the conical state is static. The order of limits (L -> infinity before t -> infinity) is clear in the main text, but the conclusion would benefit from repeating this order-of-limits caveat explicitly.","section":"Conclusion"}],"recommendation":"major_revision","confidential_remarks":"The manuscript reports interesting physics and the qualitative phase diagram appears supported by direct numerical integration. However, the two headline quantitative claims - the L^-3 scaling and the infinite-lifetime conical state - are not fully supported as written: Eq. (7) lacks a derivation, and the infinite-system calculation assumes time-independent occupation steps in a regime where the authors themselves note they should be time-dependent. These are fixable with additional analysis, so I do not recommend rejection, but the manuscript should not be accepted without a complete derivation of Eq. (7), a robustness check of the finite-size scaling at generic L, and a self-consistent treatment or quantitative justification of the frozen b± approximation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know: this paper has a genuinely new result. Voltage bias on a 1D conductor with spiral magnetic order drives the spiral into a rigidly rotating state at small bias, then quasi-periodic, then chaotic. The phase diagram and the L^-3 precession-frequency scaling are not in the prior literature, and the spin-current pumping in the RR state is a concrete observable. It is a solid subfield contribution, not a paradigm shift.\n\nWhat is done well: the three-regime picture comes from direct integration of Eq. (4), not from fitting. The linear response analysis leading to Eq. (7) is a sensible self-consistent expansion, and the scaling collapses in Fig. 5(b) are convincing within the chosen L values. The authors are also unusually honest: they flag the weak residual 2k_F modulation of the equilibrium spiral, they acknowledge region I has large finite-size effects, and they explicitly note that the b± occupation steps should in principle be time-dependent in the infinite-system calculation.\n\nThe soft spots are real but mostly minor. Eq. (7) is stated without a full derivation; the finite-size scaling deliberately picks L at the centers of beating periods, which is reasonable but needs justification; phase boundaries have no error bars and the RR/CP classification uses threshold heuristics. The bigger question is the stress-test concern about frozen b± in the infinite-system calculation. The paper's claim that a static conical state persists arbitrarily long in the thermodynamic limit rests on that calculation. The authors' defense—mz stays small, so b± are nearly constant—is plausible but not rigorous. If b± drift on a timescale comparable to t1, the conical state could slowly evolve. That said, this does not sink the qualitative picture; the L^-3 scaling is independent evidence that precession freezes, and the infinite-system calculation is clearly labeled as an approximation.\n\nWho gets value: anyone in mesoscopic spintronics or nonequilibrium magnetism. It deserves a serious referee. I would send it to review, and ask for derivations, error estimates, and ideally code/data. The central claims are interesting enough that the quantitative gaps should be patched rather than desk-rejected.","headline":"Real new phase diagram and L^-3 spin precession freeze in voltage-biased spirals, but the 'arbitrarily long times' conical state rests on frozen occupation steps that the authors concede should drift.","tokens_in":21197,"tokens_out":3510,"would_cite":true,"duration_ms":31511,"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 voltage bias across a one-dimensional conductor can set an equilibrium spin spiral into rigid rotation, then quasi-periodic wobbling, then chaos, with the rotation freezing in the thermodynamic limit.","keywords":["spin spiral","voltage bias","adiabatic approximation","nonequilibrium Green's functions","spin torque dynamics","spin pumping","chaotic spin dynamics","one-dimensional conductor"],"falsifier":"A direct numerical integration of the full time-dependent Schrödinger equation for the electrons (no adiabatic freeze) on a chain of $L\\approx 50$–$100$ at $\\delta V/w\\approx 0.03$–$0.08$ would settle it: if the RR→QP→CP sequence and the $|\\boldsymbol\\Omega|\\propto L^{-3}$ scaling do not survive, the adiabatic separation is the load-bearing approximation. On the experimental side, the spin-current noise spectrum of an engineered atomic chain under bias should show a single sharp precession peak in the RR regime whose frequency drops as the chain length grows; its absence would falsify the prediction.","tokens_in":20116,"feed_emoji":"🧲","tokens_out":8353,"duration_ms":70368,"temperature":0.7,"pith_summary":"This paper studies a one-dimensional chain of classical localized spins exchange-coupled to itinerant electrons and connected to two reservoirs held at different chemical potentials. It claims that a finite voltage bias alone, without any external oscillating field, destabilizes the equilibrium spiral order and produces collective spin dynamics: at small bias the whole spiral rotates rigidly, at larger bias it wobbles quasi-periodically, and at still larger bias it evolves chaotically. The three regimes are distinguishable in transport, and the rigidly rotating state pumps spin polarization out of the reservoirs. The authors further argue that the rotation frequency vanishes as the chain length grows, so in the thermodynamic limit a static conical spin state persists for arbitrarily long times. If correct, this is a minimal setting in which a DC voltage alone drives spontaneous time-dependent order from an equilibrium symmetry-broken state.","feed_headline":"Bias alone sets a spin spiral rotating, then chaotic","feed_subtitle":"The rotating state pumps spin polarization, and longer chains freeze the motion—DC voltage as a control knob for magnetic order.","key_machinery":"The machinery is a Landau-Lifshitz-type torque equation (Eq. 4) in which the instantaneous electronic spin polarization $\\langle\\mathbf s_i\\rangle$ is computed by the nonequilibrium Green's function formula (Eq. 3) for the frozen spin configuration, so the spins evolve under a self-consistent, highly nonlinear torque. Two auxiliary devices carry the argument: a rigid-rotation ansatz $\\partial\\mathbf m_i/\\partial t = \\boldsymbol\\Omega\\times \\mathbf m_i$ that reduces the small-bias dynamics to a linear-response eigenvalue problem (Eq. 7), and a gauge transformation to a uniform spin frame that lets the authors treat the infinite system as a single-site dynamics with a double-step occupation function. The damping parameter $\\eta$ is what selects the three orientation regimes of $\\boldsymbol\\Omega$ and controls the crossover scales $\\eta\\sim J$ and $\\eta\\sim L^{-2}$.","core_discovery":"Within an adiabatic approximation—electrons reach their steady state for each instantaneous spin configuration via the nonequilibrium Green's function formalism—and with a phenomenological damping term $\\eta$, the classical spins evolve by $\\partial \\mathbf m_i/\\partial t = J\\langle \\mathbf s_i\\rangle\\times \\mathbf m_i + \\eta(\\langle \\mathbf s_i\\rangle\\times \\mathbf m_i)\\times \\mathbf m_i$. The paper shows that the planar spiral is not stable under bias: electron polarization develops a perpendicular component that exerts a torque, and feedback from the deformed spin texture yields three dynamical phases—rigid rotation (RR), quasi-periodic (QP), and chaotic (CP)—with re-entrant RR regions in the damping-vs-bias phase diagram. A rotating-frame ansatz $\\partial \\mathbf m_i/\\partial t = \\boldsymbol\\Omega\\times\\mathbf m_i$ captures the small-bias states and yields the angular velocity $\\boldsymbol\\Omega$ as a function of damping and system size; in the linear-response regime the precession vector lies perpendicular to the spiral plane for weak damping, parallel to it for strong damping, and tilted in between. The central quantitative claim is the finite-size scaling $|\\boldsymbol\\Omega|/\\delta V \\propto L^{-3}$ in regions II and III, which implies $\\boldsymbol\\Omega\\to 0$ in the thermodynamic limit and hence a metastable conical state with diverging lifetime.","pith_inferences":["An extension the paper leaves implicit: the same bias-induced torque mechanism should operate on any collective mode with a chiral or spiral character—helical magnetic chains, and possibly spin-density waves or superconducting phase textures—so voltage alone might drive nonequilibrium dynamics in those orders too.","The $L^{-3}$ freezing suggests a practical probe: in finite chains the crossover from precession to frozen conical order could be observed as a sharp length dependence of spin-current noise, which would also test the adiabatic assumption.","One could engineer a two-terminal device in which the RR state acts as a DC-voltage-driven spin battery; coupling two such chains with opposite chirality would produce a pure spin current with no moving parts.","The reported re-entrant RR region at $\\eta/w \\sim 0.15$ implies that damping does not simply destroy the rotating phase; tuning dissipation—for instance by coupling to a substrate—could switch the chain between chaos and rigid rotation at fixed bias."],"forward_implications":["At small bias the RR state is a voltage-controlled spin-current source: the rotating spiral transfers spin polarization along the precession axis, measurable as a finite time-averaged $\\bar J^z_{sL}\\propto\\Omega$.","The RR, QP, and CP phases are distinguishable in transport: charge current is constant in RR, oscillates with incommensurate frequencies in QP, and fluctuates irregularly in CP.","Because $|\\boldsymbol\\Omega|\\propto L^{-3}$ in regions II and III, longer chains precess more slowly; for $L\\to\\infty$ before $t\\to\\infty$, the system freezes into a static conical state whose lifetime diverges with system size.","The qualitative phase sequence and the $\\boldsymbol\\Omega$ orientation regimes persist for generic spiral parameters, so the phenomenon is a property of spiral order under bias rather than a fine-tuned point.","The results imply that a DC bias can act directly on the Goldstone mode of a symmetry-broken magnetic state, without oscillating drives, which the authors propose as a general route to nonequilibrium dynamical phases."],"supporting_citations":[{"why":"Supplies the nonequilibrium Green's function steady-state formalism (Eq. 3) used to compute the electronic spin polarization for a frozen spin configuration.","marker":"[16–18]"},{"why":"Provides the magnetic Archimedes screw analogy and shows that a rigidly precessing spiral pumps spin polarization.","marker":"[57, 58]"},{"why":"Predicts spiral nuclear order in interacting nanowires, the equilibrium order whose bias-induced fate this paper analyzes.","marker":"[23, 29]"},{"why":"Reports experimental evidence of the predicted spiral order in nanowires, connecting the model to realized systems.","marker":"[30]"},{"why":"Gives the equilibrium spiral configurations and phase diagram used as the starting state for the dynamics.","marker":"[26, 28]"},{"why":"Motivates the phenomenological damping term in Eq. (4) as arising from the electronic system beyond the adiabatic approximation.","marker":"[59, 60]"},{"why":"Supplies the formula used to compute the spin current $J_{sL}$ that distinguishes the dynamical phases.","marker":"[62]"},{"why":"Gives the double-step occupation function structure used in the infinite-system conical-state calculation.","marker":"[63]"}],"fun_headline_variants":["Bias spins spiral into chaos","Voltage bias turns spiral into rotor, then chaos","DC bias drives spin spiral from rotation to chaos","Spin spiral under bias: rigid, quasi, chaotic"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes that the electrons reach their steady state for each instantaneous spin configuration much faster than the spins move, so the torques are computed from a frozen spin texture; if precession is fast enough to drag the electron distribution out of equilibrium, both the torques and the phase boundaries could change.","fun_headline_variants_meta":{"raw":{"variants":["Bias spins spiral into chaos","Voltage bias turns spiral into rotor, then chaos","DC bias drives spin spiral from rotation to chaos","Spin spiral under bias: rigid, quasi, chaotic"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000493,"raw_usage":{"total_tokens":2441,"prompt_tokens":983,"completion_tokens":1458,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":599,"completion_tokens_details":{"reasoning_tokens":1400}},"tokens_in":599,"tokens_out":1458,"duration_ms":13190,"temperature":1.0,"reasoning_tokens":1400,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T14:00:30.597925+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct numerical integration of the full time-dependent Schrödinger equation for the electrons (no adiabatic freeze) on a chain of $L\\approx 50$–$100$ at $\\delta V/w\\approx 0.03$–$0.08$ would settle it: if the RR→QP→CP sequence and the $|\\boldsymbol\\Omega|\\propto L^{-3}$ scaling do not survive, the adiabatic separation is the load-bearing approximation. On the experimental side, the spin-current noise spectrum of an engineered atomic chain under bias should show a single sharp precession peak in the RR regime whose frequency drops as the chain length grows; its absence would falsify the prediction.","supporting_citations":[{"cited_title":"Braunecker, P","cited_arxiv_id":null,"evidence_quote":"Reports experimental evidence of the predicted spiral order in nanowires, connecting the model to realized systems."},{"cited_title":"Dynamics of spin spirals in a voltage biased 1D conductor","cited_arxiv_id":null,"evidence_quote":"Gives the double-step occupation function structure used in the infinite-system conical-state calculation."}],"review_version":1}