{"id":"3d23ac14-5605-4978-b064-f2cf94209566","arxiv_id":"2508.17084","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A Floquet-engineered spin model with flip-flip and flop-flop interactions produces polarization-dependent chiral transport and, in the large-spin limit, chiral solitons whose direction is set by polarization.","lead":"This paper proposes a new class of spin models, called spin Kitaev models, where spin excitations with different polarizations travel in opposite directions along a chain. The work also predicts chiral solitons and solitonic molecules in these models, with a concrete Floquet protocol for experimental realization.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The chiral-soliton and solitonic-molecule claims are derived only in the semiclassical large-S limit (Eqs.","rationale":"The paper's strongest-claim decomposition is fair: the polarization-dependent chiral transport is supported by exact numerics for S=1/2 to 5 and by a short-time pathway argument that does not require the semiclassical limit. The chiral solitons and molecules, however, are the genuinely new advertised phenomena, and they live entirely in the semiclassical approximation. I checked the mean-field derivation and the construction of the traveling profiles; internally, Eqs. (7)-(9) and the boundary-condition analysis in the Supplement are consistent with each other and with the stated domain-wall boundary conditions. The missing piece is external validation against the original quantum model at finite S. This matters because the soliton states are domain walls rather than small spin excitations, and because the dynamical equations contain directional amplification, so stability is not automatic. The Floquet protocol and the experimental-feasibility statement are secondary; even if the protocol has 1/T corrections, the model and the semiclassical mapping are what the claims rest on. Since no finite-S quantum test of the solitons is reported, and the abstract states the soliton results without this caveat, the reader's conditional verdict is appropriate. My concern does not move the verdict.","tokens_in":18191,"tokens_out":21361,"duration_ms":209844,"concrete_test":"Simulate the full quantum spin-S dynamics of H_sK at theta=phi=pi/2 with Delta/J=0.1, h=0 and h=0.0001J, using e.g. time-evolving block decimation (TEBD) on a chain of about 80-120 sites for S=2 and S=3, initialized as the product of spin-coherent states corresponding to the semiclassical soliton profile (e.g., vs=3J, where the width is tens of sites). Compare the co-moving density Jx_n(t) and Jy_n(t) with Eq. (8) integrated over Jt about 50-100. If the peak position shifts or the width doubles by more than the 1/S estimate, the semiclassical soliton claim is not robust; if the profile propagates shape-preserving, the concern is resolved. As a cheaper analytical check, linearize the 1/S quantum corrections (truncated Wigner or Bogoliubov) around the soliton and look for unstable modes with growth rates comparable to J.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central new results beyond the linear Bosonic Kitaev model are the chiral solitons and solitonic molecules advertised in the abstract and shown in Figs. 2 and 3. These follow exclusively from the semiclassical equations of motion (Eqs. 7-9), in which spin operators are replaced by expectation values and fluctuations of order 1/S are dropped (Supplement, 'Mean field approach', Eq. S14). The main text states only that 'we consider a large but finite S, where the semiclassical approach applies'; no estimate of the 1/S corrections, no linear stability analysis of the soliton profile against quantum fluctuations, and no finite-S quantum simulation of the soliton dynamics are provided. The exact diagonalization benchmarks in Fig. 1 and Fig. S2 cover only small-α chiral transport (S up to 5), not the soliton states. Two further features make the gap consequential. First, the soliton profiles are not small excitations of the fully polarized state used for transport: they are domain walls with Jz(n→-∞)=-1 and Jz(n→+∞)=+1, so the relevant quantum state is strongly inhomogeneous, and entanglement growth near the domain wall could dephase the semiclassical trajectory on the time scales shown. Second, the demonstration that Eq. (8) translates without dispersion shows that the traveling-wave ansatz is a solution; it does not establish stability under perturbations, which the inherent directional amplification of the underlying Hatano-Nelson dynamics makes a nontrivial requirement. Thus the advertised soliton phenomenology in the quantum spin model is the least secure part of the central claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"Lv et al. propose Floquet protocols to realize a 'spin Kitaev model' whose Hamiltonian contains both exchange (flip-flop) and pair-creation-like (flip-flip/flop-flop) interactions. They study the resulting dynamics in one dimension. The paper's main results are: (i) a Floquet pulse sequence, Eq. (2), that generates the model with arbitrary phases and spin S; (ii) a demonstration, via a short-time interference calculation and exact diagonalization for S=1/2 through 5, that spin excitations polarized along x and y propagate chirally in opposite directions; (iii) a semiclassical large-S mapping of the model to a nonlinear Hatano-Nelson equation, Eqs. (7)-(8); and (iv) the construction of chiral solitons and solitonic molecules in that semiclassical limit, with dynamics shown in Figs. 2 and 3. The abstract advertises polarization-dependent chiral transport and chiral solitons as the central physics.","tokens_in":18498,"tokens_out":4483,"duration_ms":47255,"significance":"If the results hold, the paper introduces a new class of spin models with directional, polarization-dependent transport that is accessible in cold-molecule, magnetic-atom, and Rydberg-atom platforms. The chiral transport claim is well supported: the interference argument leading to Eq. (6) is simple and transparent, and the exact diagonalization data for S=1/2 up to 5 show clear asymmetric propagation consistent with the short-time prediction. The paper is also self-contained: the Floquet sequence is explicit, the effective Hamiltonian is derived, and the mean-field equations are obtained from the Heisenberg equations of motion. The soliton and solitonic-molecule results are intriguing and represent the main novel physics, but they currently rest entirely on the semiclassical large-S approximation, with no finite-S quantum test and no stability analysis. The paper is clearly written and the supplementary material is thorough, which strengthens its usefulness to the community even though the soliton predictions need additional support.","major_comments":[{"comment":"The chiral-soliton and solitonic-molecule results are derived and verified only in the semiclassical large-S limit, where spin operators are replaced by expectation values and 1/S quantum fluctuations are dropped (Supplement, Eq. S14). The manuscript does not provide a finite-S quantum check (e.g., exact diagonalization or truncated Wigner for moderate S) nor a linear stability analysis of the soliton profiles against perturbations. Since the abstract presents chiral solitons as a central result, this missing support is load-bearing. The traveling-wave ansatz shows that the profile is a solution of the mean-field equation, but it does not establish that it is robust in the quantum model; the directional amplification of the underlying Hatano-Nelson dynamics makes stability a nontrivial requirement.","section":"Chiral Solitons, Eqs. (7)-(8) and Supplemental 'Mean field approach'"},{"comment":"Equation (8) is derived from Eq. (7) by replacing J^z_n with sqrt(1-(J^x_n)^2), which is only valid where J^z_n >= 0. The soliton profiles in Fig. 2(a) satisfy J^z_n -> -1 as n -> -infinity and cross zero, so the positive square root with a single sign cannot describe the full profile. The Supplement (Eq. S14 onward) correctly uses a '±' sign, but the main text omits it. As written, Eq. (8) is inconsistent with the boundary conditions used to construct the solitons; the sign of the square root must be made piecewise explicit, or the derivation must be restricted to regions where J^z > 0 and the boundary condition adjusted accordingly.","section":"Chiral Solitons, Eq. (8) and Fig. 2"},{"comment":"The solitonic molecule is a bound state of x- and y-solitons in a magnetic field h, and its travel direction is shown to depend on the molecule's orientation. The demonstration is numerical solution of the semiclassical equations only. In addition, the value h = -0.0001J used in Fig. 3 is extremely small; while the Supplement shows larger h for one orientation (Fig. S7), the dependence of the binding on h for both orientations and the molecule's stability under quantum fluctuations are not addressed. These points should be discussed or supplemented for the molecule claim to be at the same level of support as the transport claim.","section":"Chiral Solitonic Molecules, Eq. (9) and Fig. 3"}],"minor_comments":[{"comment":"The same symbol U is used for the Floquet propagator and for the unitary rotation operator U; this is confusing and should be disambiguated (e.g., U_F for the Floquet operator).","section":"Floquet Implementation, Eq. (2)"},{"comment":"The quadrature operator is written as 'ˆP = i(ˆa† - ˆa)/√2', but the text later uses P_n for its expectation value while also using p for the molecule's polarization; consider renaming one of these to avoid confusion.","section":"After Eq. (4)"},{"comment":"The replacement J^z_n = sqrt(1-(J^x_n)^2) should carry an explicit ± sign, as the Supplement does, to indicate that J^z can be negative.","section":"Chiral Solitons, Eq. (8)"},{"comment":"The soliton velocities vs = ±50J in Fig. 2 and vs = ±200J in Fig. 3 are much larger than the microscopic coupling scale; a brief remark on the experimental regime in which such fast, sharply localized solitons could be observed would be helpful.","section":"Fig. 2 and Fig. 3"},{"comment":"In the text near Fig. S3, the sentence 'In Fig. S3(a-d)' refers to the dependence of soliton profiles on Δ, but Fig. S3 has only two panels (a,b); the cross-reference should be corrected.","section":"Supplement, 'Boundary conditions of the solitons'"}],"recommendation":"major_revision","confidential_remarks":"The chiral transport result is convincing and well-supported, and the paper is well organized. The main risk is the soliton and solitonic-molecule claims, which currently rest on the semiclassical large-S limit without finite-S verification or stability analysis. I would ask the authors to either provide a finite-S quantum test (even on small chains) or clearly construe the solitons as predictions of the semiclassical limit. The sign issue in Eq. (8) should also be fixed; it is a genuine inconsistency that affects the soliton construction. If these points are addressed, the paper would be a solid contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You'll want to know this one for the Floquet protocol and the transport result. The soliton half is real but thinner than the abstract lets on.\n\nWhat is actually new: the paper proposes a concrete pulse sequence, Eq. (2), that engineers a spin Kitaev model with flip-flip/flop-flop terms and arbitrary phases, starting from exchange interactions available in molecules, magnetic atoms, and Rydberg arrays. That construction is explicit and credible. The main physics claim, polarization-dependent chiral transport, is well supported: a simple first-order interference argument (Eq. 6) shows x- and y-polarized excitations amplifying in opposite directions, and exact diagonalization for S=1/2 through 5 confirms it. The infinite-S limit reducing to the known bosonic Kitaev model is acknowledged, not hidden. The genuinely new content is the finite-S transport and the nonlinear density-dependent mass in Eq. (8).\n\nThe soft spots are concentrated in the soliton section. The chiral solitons and solitonic molecules in Figs. 2 and 3 come exclusively from the semiclassical large-S equations of motion, with operators replaced by expectation values and 1/S fluctuations dropped. The paper gives no estimate of 1/S corrections, no linear stability analysis of the soliton profile, and no finite-S quantum dynamics for the solitons. This matters because the soliton profiles are domain walls, Jz going from -1 to +1, so they are not small fluctuations of the polarized state used for the transport benchmarks. Entanglement near the domain wall could dephase the semiclassical trajectory on the times shown. The traveling-wave ansatz shows these are solutions of the semiclassical equations, but does not establish stability, and stability is nontrivial in a Hatano-Nelson-type setting with directional amplification.\n\nThat said, the semiclassical limit is a legitimate place to look for solitons, and for large S one might expect the quantum corrections to be small. The issue is that the paper does not make that case. The abstract's statement that the results are \"directly accessible in current experiments\" is also stronger than what is demonstrated; the protocol is plausible but the feasibility discussion is brief.\n\nThe citation pattern is fine. The curved-space interpretation leans on the authors' own earlier papers, but that is interpretation, not a load-bearing derivation.\n\nWho is this for? The Floquet engineering and quantum simulation community, and people working on non-Hermitian transport. The transport part should survive scrutiny; the soliton part needs qualification or additional evidence.\n\nMy recommendation: send it to a serious referee. Conditional acceptance after the authors quantify or soften the soliton claims would be appropriate.","headline":"A clean Floquet construction and solid ED evidence for polarization-dependent chiral transport, with the chiral-soliton story resting on unverified semiclassical dynamics.","tokens_in":19043,"tokens_out":1839,"would_cite":true,"duration_ms":21656,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper proposes a Floquet-realizable spin Kitaev model in which x- and y-polarized spin excitations travel in opposite chiral directions, and whose large-spin limit hosts polarization-dependent chiral solitons.","keywords":["spin Kitaev model","chiral transport","Floquet engineering","Hatano-Nelson model","non-Hermitian skin effect","spin solitons","quantum simulation","polar molecules"],"falsifier":"Track the full quantum dynamics for a spin chain with $S=3/2$ and $N=7$, set $\\theta=\\phi=\\pi/2$ and $\\Delta=0.8J$, and start from $\\hat{U}_x$ applied to the fully polarized state. If $S^x_n(t)$ does not grow preferentially at $n_0+1$ while $S^y_n(t)$ grows at $n_0-1$ in the early-time window, or if a localized profile initialized from the nonlinear Hatano-Nelson equation disperses instead of translating at constant velocity greater than $2J$, the central transport and soliton claims fail.","tokens_in":1983,"feed_emoji":"🧲","tokens_out":3098,"duration_ms":83766,"temperature":0.7,"pith_summary":"The paper proposes a one-dimensional spin model, the spin Kitaev model, whose interactions include the usual flip-flop exchange together with flip-flip and flop-flop terms that do not conserve total magnetization. It shows how a sequence of Floquet pulses can realize this model for arbitrary spin S with controllable phases, using current ultracold molecule, magnetic atom, or Rydberg platforms. In a fully polarized chain, a local x-polarized excitation amplifies and travels to the right while a y-polarized excitation travels to the left, a polarization-dependent chiral skin effect. In the large-S limit the model becomes a nonlinear Hatano-Nelson equation whose density-dependent mass supports chiral solitons, and a magnetic field binds these into chiral solitonic molecules whose travel direction is controlled by their internal orientation. The authors argue these effects are directly accessible in experiments and open routes to spintronics, information processing, and quantum sensing.","feed_headline":"Spin Kitaev model sends x and y excitations opposite ways","feed_subtitle":"A Floquet protocol engineers flip-flip spin terms, producing chiral transport and one-way solitons in large-spin chains.","key_machinery":"The central object is the spin Kitaev Hamiltonian $\\hat{H}_{sK}=\\hat{H}_{ex}+\\hat{H}_{ff}$, where $\\hat{H}_{ex}$ contains exchange terms $J e^{-i\\theta}\\hat{S}^-_n\\hat{S}^+_{n+1}$ and $\\hat{H}_{ff}$ contains pair-flip terms $\\Delta e^{-i\\phi}\\hat{S}^+_n\\hat{S}^+_{n+1}$. The carrying identity is the quadrature decomposition: at $S\\to\\infty$, the $X$ and $P$ quadratures obey $\\dot X_n=(\\Delta+J)X_{n-1}+(\\Delta-J)X_{n+1}$ and $\\dot P_n=(J-\\Delta)P_{n-1}-(\\Delta+J)P_{n+1}$, two Hatano-Nelson chains with opposite nonreciprocity. At large finite $S$, the semiclassical equation $\\dot J_n^x=[(\\Delta+J)J_{n-1}^x+(\\Delta-J)J_{n+1}^x]J_n^z$ with $J_n^z=\\sqrt{1-(J_n^x)^2}$ turns the nonlinear term into a density-dependent mass that halts exponential amplification and stabilizes chiral solitons.","core_discovery":"The central claim is that adding flip-flip and flop-flop terms to a spin-exchange chain creates a new class of spin models with transport properties beyond conventional spin dynamics. For $\\theta=\\phi=\\pi/2$, exact diagonalization shows that x-polarized excitations are amplified toward the right and y-polarized excitations toward the left for spin S from 1/2 to 5. In the $S\\to\\infty$ limit, Holstein-Primakoff transformation reduces the model to a bosonic Kitaev model whose X and P quadratures each satisfy a Hatano-Nelson chain with opposite nonreciprocity. For large finite S, retaining the nonlinear factor in the semiclassical equations yields a nonlinear Hatano-Nelson model whose density-dependent mass stabilizes localized one-way solitons moving faster than 2J; with a magnetic field, x- and y-solitons bind into molecules whose direction of travel is set by the relative ordering of the two components.","pith_inferences":["If the solitonic molecule's direction is governed by its internal polarization $p$, then reversing the internal order of the two solitons should flip the propagation direction without reversing the lattice or the field, suggesting an internal-state-controlled spin diode.","The same density-dependent mass mechanism could be tested in platforms other than large-spin chains, including momentum-state or cavity implementations, where the nonlinear Hatano-Nelson equation may be realized without engineering collective spins.","The early-time interference prediction $S^x_{n_0\\pm 1}=S\\alpha t(\\Delta\\pm J)$ offers a quantitative observable for isolating the contribution of each pathway; testing it at larger $\\alpha$ would probe how the interference evolves beyond the perturbative regime.","The boundary-condition dependence of the soliton direction means that domain walls in $J_z$ could act as tunable launching points: flipping the global $J_z$ boundary condition reverses both soliton and molecule velocities."],"forward_implications":["With $\\theta=\\phi=\\pi/2$, a local x-polarized excitation in the fully polarized chain is amplified to the right, while a y-polarized excitation is amplified to the left, for spin S from 1/2 to 5 in exact diagonalization.","The semiclassical large-S equations reduce to a nonlinear Hatano-Nelson equation, while the $S\\to\\infty$ limit recovers linear Hatano-Nelson chains for the X and P quadratures with couplings $\\Delta\\pm J$.","Chiral soliton solutions exist only for $|v_s|>2J$; as $v_s$ approaches $2J$ the soliton narrows to zero width, and for $v_s\\le 2J$ the localized solutions disappear.","A magnetic field binds one x-soliton and one y-soliton into a solitonic molecule whose travel direction is set by the internal order of the two solitons, i.e. by the sign of $p=n_x-n_y$.","The proposed Floquet pulse sequence realizes the spin Kitaev model with arbitrary phases $\\theta$ and $\\phi$, and arbitrary effective spin S built from ensembles of spin-1/2 particles."],"supporting_citations":[{"why":"Supplies the bosonic Kitaev-Majorana chain whose quadrature decomposition into two Hatano-Nelson chains is the $S\\to\\infty$ limit the paper builds on.","marker":"[41]"},{"why":"Provides the non-Hermitian Hatano-Nelson model and the skin-effect amplification language used to characterize the chiral transport.","marker":"[37–40]"},{"why":"Gives the hyperbolic-space-curvature interpretation of asymmetric non-Hermitian couplings that the paper uses to explain the origin of chiral solitons.","marker":"[54]"},{"why":"Shows how paired spin flips in dipolar multilayers generate bosonic Kitaev dynamics, supporting the collective-spin route to large-S implementations.","marker":"[52]"},{"why":"Demonstrates that flip-flip terms can arise from interband coupling in a double-well lattice, evidence that such non-number-conserving spin terms are physically accessible.","marker":"[35]"},{"why":"Provides the experimentally demonstrated Floquet-engineered XYZ spin model that the proposed pulse sequence extends to the spin Kitaev model.","marker":"[32]"},{"why":"Establishes the Dicke-manifold protection used to make many spin-1/2 particles behave as a collective large spin.","marker":"[50]"}],"fun_headline_variants":["Spin Kitaev model: x goes right, y goes left","Opposite polarization, opposite spin transport","Floquet-engineered spin chain yields one-way solitons","Chiral solitons on demand via spin polarization"],"cache_read_input_tokens":21120,"weakest_assumption_plain":"The chiral solitons and solitonic molecules are derived and simulated in the large-S semiclassical limit, where spin operators are replaced by their expectation values and $1/S$ quantum fluctuations are neglected; if those fluctuations destabilize the localized traveling waves, the advertised soliton phenomenology would not occur in the quantum spin model.","fun_headline_variants_meta":{"raw":{"variants":["Spin Kitaev model: x goes right, y goes left","Opposite polarization, opposite spin transport","Floquet-engineered spin chain yields one-way solitons","Chiral solitons on demand via spin polarization"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000474,"raw_usage":{"total_tokens":2356,"prompt_tokens":949,"completion_tokens":1407,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":565,"completion_tokens_details":{"reasoning_tokens":1342}},"tokens_in":565,"tokens_out":1407,"duration_ms":12646,"temperature":1.0,"reasoning_tokens":1342,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T17:07:49.775702+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Track the full quantum dynamics for a spin chain with $S=3/2$ and $N=7$, set $\\theta=\\phi=\\pi/2$ and $\\Delta=0.8J$, and start from $\\hat{U}_x$ applied to the fully polarized state. If $S^x_n(t)$ does not grow preferentially at $n_0+1$ while $S^y_n(t)$ grows at $n_0-1$ in the early-time window, or if a localized profile initialized from the nonlinear Hatano-Nelson equation disperses instead of translating at constant velocity greater than $2J$, the central transport and soliton claims fail.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows how paired spin flips in dipolar multilayers generate bosonic Kitaev dynamics, supporting the collective-spin route to large-S implementations."},{"cited_title":"Christakis, J","cited_arxiv_id":null,"evidence_quote":"Provides the experimentally demonstrated Floquet-engineered XYZ spin model that the proposed pulse sequence extends to the spin Kitaev model."}],"review_version":2}