{"id":"9eed4027-5a05-4432-976c-1c4d8c47f233","arxiv_id":"2502.09764","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Atomistic simulations show that current-driven mixtures of skyrmions and skyrmioniums form tilted lanes and eventually all skyrmioniums collapse into skyrmions, giving a three-phase dynamic diagram.","lead":"This paper uses computer simulations of two kinds of magnetic whirlpools, skyrmions and skyrmioniums, to show that electric currents make them separate into tilted traffic lanes, and strong currents crush the larger whirlpools into the smaller ones. It maps the current and density conditions for each phase, connecting the behavior to lane formation seen in pedestrian crowds and colloids.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Single-realization, fixed-parameter simulations leave the JS–LS–MSk phase boundaries unverified; a finite-size and ensemble reproducibility check would settle whether the claimed phases are robust.","rationale":"The Reader's weakest assumption concerns the representativeness of the fixed material parameters, particularly α = 0.4 and D/J. That is a legitimate concern about generalizability, but the more load-bearing issue is internal to the simulations: because each state point is a single 200 ns run at one system size with no order parameter, error bars, or initial-condition averaging, the paper does not yet demonstrate that the observed JS/LS/MSk labels are reproducible phase behavior rather than finite-size or initial-condition effects. These two concerns are related: a sensitivity sweep would also implicitly test robustness, but it would not establish reproducibility. The proposed concrete test addresses the internal validity first: if the same phases appear at larger system size and across random initial configurations, then the central qualitative claim is secure and parameter sensitivity becomes a quantitative calibration issue; if the labels change, the phase diagram needs revision regardless of material parameters. The Discussion section explicitly acknowledges that thermal effects and alternative current-drive schemes are not treated, which further supports a conditional rather than unconditional reading. I do not see an internal inconsistency or a fatal flaw; the work is a plausible computational study that needs reproducibility checks, so the Reader's CONDITIONAL verdict should stand unchanged.","tokens_in":15533,"tokens_out":6941,"duration_ms":81669,"concrete_test":"Use the same Hamiltonian and SOT dynamics (Sec. II) to recompute the Nsk=14, Nskium=1 and Nsk=25, Nskium=1 state points at L = 272 nm (four times the area, matched densities) and at the original L = 136 nm for at least 10 independently randomized initial configurations. For each run, record the time-averaged θskium and a quantitative lane order parameter, e.g., the anisotropy of the two-species density correlation along the lane direction. If the JS-to-LS sign reversal persists, the phase boundary shifts by less than one current bin, and the phase labels agree across realizations and system sizes, the concern is resolved; if labels vary or the boundary moves with L, the claimed phase diagram is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that a driven skyrmion/skyrmionium mixture exhibits three dynamic phases (partially jammed, laned, and skyrmion-only moving crystal) with a sign reversal of the skyrmionium deflection at the JS–LS boundary. The evidence for each phase boundary is a single 200 ns time average for a single initial configuration at L = 136 nm (Sec. II), mostly with Nskium = 1, with phase labels assigned by visual inspection of trajectories and by the sign of the time-averaged velocity component ⟨Vy⟩_skium (Figs. 2, 5, 8, 10, and 13; phase diagram in Fig. 7). No quantitative order parameter is defined that would let an independent researcher identify JS versus LS versus MSk states; no error bars or initial-condition ensemble are reported; no system-size variation is shown; and α = 0.4 is fixed, although the Hall-angle and mobility differences that produce the tilted lanes scale with damping in the Thiele-like dynamics. The load-bearing assumption is therefore that the observed transitions are stationary, reproducible, thermodynamic-like phase behavior rather than finite-size artifacts or run-specific metastable trajectories. This assumption is currently untested, so the phase diagram in Fig. 7 may not be robust.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents atomistic Landau-Lifshitz-Gilbert simulations of a thin ferromagnetic film containing coexisting skyrmions and skyrmioniums driven by spin-orbit torque. It claims that such mixtures exhibit three dynamic phases as a function of applied current and skyrmion density: a partially jammed state (JS) in which the skyrmionium is dragged by the surrounding skyrmions and moves at a finite angle, a laning state (LS) in which the skyrmionium opens a tilted lane and its angle of motion reverses or approaches zero, and a skyrmion-only moving lattice (MSk) in which all skyrmioniums collapse into skyrmions. A phase diagram is constructed for a single skyrmionium among various skyrmion densities, and additional simulations with multiple skyrmioniums show that lane formation is enhanced or suppressed depending on the density ratio. The results are compared to tilted laning in chiral pedestrian flows.","tokens_in":1591,"tokens_out":4509,"duration_ms":80159,"significance":"If the claimed phases are robust, the paper introduces a new and potentially important class of driven binary systems: mixtures of magnetic textures with different Hall angles and mobilities. The work uses a standard atomistic LLG+SOT model with literature-based material parameters and does not fit any parameter to force the laning result, which is a strength. The connection to tilted laning in pedestrian dynamics provides a fresh interdisciplinary link. However, the significance is bounded by the fact that the phase diagram rests on single-realization, fixed-parameter simulations; without reproducibility and sensitivity checks, the central claim of three distinct phases is not yet fully established.","major_comments":[{"comment":"The JS-LS-MSk phase boundaries are assigned from a single 200 ns time average for one initial configuration per (j, nsk) point, with phase labels based on visual inspection of trajectories and on the sign or zero-crossing of <Vy>_skium. No quantitative order parameter (for example a lane order parameter, a collision frequency, or a transverse displacement measure) is defined, and no initial-condition ensemble or error bar is reported. Since the central claim is the existence of three distinct dynamic phases and a phase diagram, the phases should be reproducible over multiple random initializations and characterized by a well-defined observable; otherwise the boundaries in Fig. 7 are not falsifiable.","section":"Section III, Figs. 2, 5, 7"},{"comment":"All simulations use a single system size L = 136 nm with periodic boundary conditions and small numbers of textures (Nsk from 8 to 26, Nskium from 1 to 12). No system-size variation is shown. In driven binary mixtures, laning can be strongly affected by finite-size effects and periodic boundary conditions, so the existence and width of the LS region in Fig. 7 should be checked at least at one other L (or at the same density with a doubled cell) to establish that the phases are thermodynamic-like rather than artifacts of the box size or of the specific realization.","section":"Section II"},{"comment":"The phase diagram is computed for a single material parameter set (J = 1 meV, D = 0.2J, K = 0.01J, alpha = 0.4, T = 0). The skyrmion Hall angle, the mobility difference between skyrmions and skyrmioniums, and the current at which skyrmioniums collapse all depend on alpha and D/J; the paper itself notes that low damping reduces skyrmionium stability (ref. 39). No sensitivity analysis is given, so it is unclear whether the JS-LS-MSk boundaries would shift or disappear for parameters typical of other skyrmion-hosting films. A small parameter sweep or at least a quantitative discussion of the expected dependence is needed to support the generality implied by the phase diagram.","section":"Section II (parameters) and Section VI"},{"comment":"The definitions of the ID (independent dynamics) and LS phases appear to be mutually inconsistent. The text states that the LS appears in the range 0.00027 < nsk < 0.00189 and that the ID phase is stable for nsk <= 0.00038, and then says that lane formation is still possible in the ID state. If ID and LS can overlap in a non-trivial way, the phase diagram in Fig. 7 needs a clearer separation criterion; if they are meant to be exclusive, the reported density ranges need to be reconciled. This ambiguity affects the interpretation of the low-density part of the phase diagram.","section":"Section III, Fig. 7 and text around it"}],"minor_comments":[{"comment":"Reference 44 contains a typo in the title: 'skyrmionum' should be 'skyrmionium'.","section":"References"},{"comment":"The notation in the SOT torque term is confusing: m x (j x z) x m uses the same symbol for the current unit vector and the current density j; please distinguish the unit vector (for example j-hat) from the scalar density.","section":"Eq. (2)"},{"comment":"The dashed horizontal line marks the sign change of <Vy>_skium, but the figures would benefit from error bars or at least a statement that the plotted values are time averages whose statistical uncertainty is small compared to the symbol size.","section":"Figs. 2 and 5"},{"comment":"The text says that at nsk = 0.00065 lane formation is optimized and the LS extends up to j = 4.50 x 10^9 A/m^2, but no simulation snapshot or velocity data are shown for this particular density; a reference to a figure or a short discussion would help the reader connect this point to the other results.","section":"Section III"},{"comment":"The abstract states that 'the skyrmionium angle of motion is reversed' in the laned state, but for the Nsk = 25 sample (Fig. 5) the angle only approaches zero without becoming positive. The abstract should be qualified to indicate that the sign reversal occurs for moderate skyrmion densities but not for all laned states.","section":"Abstract and Section VII"}],"recommendation":"major_revision","confidential_remarks":"The central physical idea is attractive and the simulations are standard, but the evidence for the phase diagram needs strengthening. I would ask the authors to add either (a) an ensemble of random initializations with error bars on the phase boundaries, or (b) a quantitative order parameter and a finite-size check, before publication. The current manuscript is a solid starting point, but the 'three phases' claim is too central to rest on single trajectories. The overlap between ID and LS definitions should also be clarified."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The core physics here is new and the simulations look honest. Showing that a driven skyrmion/skyrmionium mixture separates into jammed, laned, and skyrmion-only phases, with tilted lanes and a sign reversal of the skyrmionium deflection, is a real step beyond the prior Magnus-dispersion particle work by the same group. The build is clean: standard LLG+SOT, material parameters taken from the literature, no fitting to force the laning. The velocity and angle curves in Figs. 2, 5, 8, and 10 show clear qualitative transitions, and the trajectory snapshots back up the phase labels.\n\nThe soft spots are real but addressable. The phase diagram in Fig. 7 rests on single simulations per parameter point, with no error bars, no ensemble averaging over initial configurations, and no system-size variation. The phase labels themselves are assigned by visual inspection and by the sign of the time-averaged <Vy>, not by a quantitative order parameter. That would be fine for identifying interesting behavior, but it is not yet a thermodynamic-like phase boundary map. There is also no sensitivity analysis for alpha or D/J, even though those control the Hall angle and mobility contrast that produce the lane tilt and the collapse current. The stress-test note is fair: a few runs at different damping, a larger cell, and a repeat with a different random seed would go a long way.\n\nThat said, the limitations are stated plainly in the text (T=0, no thermal effects, the Discussion flags the STT alternative), and the analogy to the Bacik pedestrian lanes is appropriately qualitative. The paper does not oversell; it presents a phase diagram but acknowledges the regimes in prose rather than overclaiming precision.\n\nMy verdict is conditional, not skeptical. I would send this to peer review, but I would ask the authors to add an ensemble or size check and at least one alpha sensitivity point. The central idea is likely correct, but the boundary locations in Fig. 7 currently have unknown error bars, and the paper would be stronger if the authors showed the transitions do not depend on one lucky trajectory.","headline":"A believable atomistic demonstration of laning in skyrmion/skyrmionium mixtures, but the phase diagram is built from single runs at fixed parameters and needs reproducibility checks before it can stand.","tokens_in":16330,"tokens_out":2022,"would_cite":true,"duration_ms":23238,"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":"Under a spin-orbit-torque current, a mixture of skyrmions and skyrmioniums forms three dynamical phases—partially jammed, tilted laning, and skyrmion-only moving crystal—with the laning transition signaled by a reversal of the…","keywords":["skyrmion","skyrmionium","laning transition","phase separation","Hall angle","spin-orbit torque","atomistic simulation","magnetic textures"],"falsifier":"At skyrmion density $n_{sk} = 0.00081\\ \\mathrm{nm}^{-2}$ with one skyrmionium, the paper predicts the skyrmionium transverse velocity reverses sign near $j = 1.25 \\times 10^{9}\\ \\mathrm{A/m^{2}}$ when laning begins. A trajectory-resolved measurement or simulation that finds no such sign reversal, or no lane for densities up to $0.00189\\ \\mathrm{nm}^{-2}$, would contradict the core phase diagram.","tokens_in":15350,"feed_emoji":"🌀","tokens_out":6509,"duration_ms":59301,"temperature":0.7,"pith_summary":"Skyrmions and skyrmioniums are two kinds of swirling magnetic textures that can coexist in a thin magnetic film. This paper uses atomistic simulations to show that when an electric current drives such a mixture, the two species do not just move past each other: they form three distinct dynamical phases. At low current the skyrmionium is dragged along by the slower, Hall-deflected skyrmions; at higher current it opens a tilted lane through them and its sideways motion reverses direction; at even higher current it collapses into an ordinary skyrmion. The authors map these phases as a function of current and skyrmion density. The point is that binary mixtures of magnetic textures with different topology are a new setting for laning transitions, with tilted lanes arising naturally from the skyrmion Hall angle.","feed_headline":"Skyrmions and skyrmioniums separate into tilted lanes under current","feed_subtitle":"Atomistic simulations map three flow phases: jamming, laning, and collapse into a skyrmion-only lattice.","key_machinery":"The load-bearing objects are the two magnetic textures themselves: a skyrmion, a swirling spin texture with topological charge $Q = \\pm 1$, and a skyrmionium, a concentric double-twist texture with $Q = 0$. Their dynamics are integrated with the atomistic Landau-Lifshitz-Gilbert equation augmented by a spin-orbit-torque term. The mechanism that produces the phases is the combination of the mobility difference (skyrmioniums move roughly twice as fast) with the skyrmions' intrinsic Hall angle of about $67^\\circ$ in these parameters, while skyrmioniums have no Hall angle. Lanes appear when the current is large enough for the skyrmionium to separate from the skyrmion bath; their tilt direction is controlled by the skyrmion Hall angle, and their stability window is bounded by skyrmion density and by current-induced collapse of the skyrmionium.","core_discovery":"The central claim is that a current-driven skyrmion/skyrmionium mixture is a genuine binary active system with a controllable phase diagram. Because skyrmions carry a nonzero topological charge they move slowly and at a finite skyrmion Hall angle, while skyrmioniums have zero net charge, move about twice as fast, and travel straight along the drive. The simulations show that at low drives the fast, straight-moving skyrmionium is repeatedly scattered by the slow, tilted-moving skyrmions, producing a partially jammed state in which the skyrmionium itself moves at a finite positive angle. Above a threshold current the skyrmionium pushes open a lane through the skyrmions; because the lane is set by the skyrmions' Hall deflection it is tilted opposite to the skyrmion Hall angle, and once the skyrmionium locks to this lane its transverse velocity reverses sign. If the current is high enough, or the skyrmion density too high, the pressure on the skyrmionium collapses its inner skyrmion and the system becomes a lattice of ordinary skyrmions only.","pith_inferences":["Beyond the paper: any other pair of magnetic textures with different Hall angles or mobilities, such as skyrmion-antiskyrmion mixtures, should show analogous lane formation with the lane tilt set by the sign of the Hall angle.","Beyond the paper: the sharp collapse boundary suggests a current-pulse scheme that converts a skyrmionium into a skyrmion and back, acting as a binary switch; the paper does not demonstrate the reverse conversion.","Beyond the paper: thermal fluctuations might widen the laning window by helping the skyrmionium penetrate the skyrmion lattice, but they could also lower the collapse current; the paper notes this trade-off but does not simulate temperature.","Beyond the paper: a direct experimental test would track coexisting skyrmions and skyrmioniums with time-resolved X-ray microscopy and look for the predicted sign reversal of the skyrmionium transverse velocity at the jamming-laning boundary."],"forward_implications":["A single skyrmionium in a skyrmion bath undergoes a dynamic transition from partially jammed to laned flow at a well-defined current, signaled by a reversal of the skyrmionium's transverse velocity.","The laning state exists only in a finite window of skyrmion densities: below it the textures move independently, and above it the skyrmionium is dragged until it collapses.","Multiple skyrmioniums can reinforce one another's lanes, widening the current window for laning when the skyrmion density is moderate.","At sufficiently high currents all skyrmioniums collapse into skyrmions, producing a skyrmion-only moving lattice.","The tilted-lane geometry is opposite to the skyrmion Hall angle, connecting this magnetic system to chiral pedestrian flows."],"supporting_citations":[{"why":"Supplies the chiral pedestrian laning result that the tilted lanes are compared to.","marker":"[14]"},{"why":"Provides the binary active-matter laning framework and the Magnus-dispersion treatment this work extends to magnetic textures.","marker":"[13]"},{"why":"Establishes the intrinsic skyrmion Hall angle that sets the lane tilt.","marker":"[28]"},{"why":"Supplies the skyrmionium property of high velocity with no Hall effect, which creates the mobility contrast.","marker":"[36]"},{"why":"Provides the comparison of skyrmionium and skyrmion dynamics and interactions underpinning the density-dependent behavior.","marker":"[46]"},{"why":"Documents current-driven skyrmionium collapse into skyrmions, the process that terminates the laning state.","marker":"[45]"},{"why":"Supplies the atomistic spin-dynamics method used for the simulations.","marker":"[55]"},{"why":"Provides the Pt/Co/MgO material parameters used to stabilize the simulated Neel skyrmions.","marker":"[60]"}],"fun_headline_variants":["Skyrmion mixtures form tilted lanes under current","Chiral lanes emerge in skyrmion-skyrmionium flows","Spin-orbit torque drives skyrmion-skyrmionium laning","Skyrmion phase separation: jamming to laning to collapse"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The results rest on the assumption that the specific set of magnetic-film parameters used in the simulations, including the damping value 0.4 and the ratio of Dzyaloshinskii-Moriya to exchange strength 0.2, is representative of real films where skyrmions and skyrmioniums coexist.","fun_headline_variants_meta":{"raw":{"variants":["Skyrmion mixtures form tilted lanes under current","Chiral lanes emerge in skyrmion-skyrmionium flows","Spin-orbit torque drives skyrmion-skyrmionium laning","Skyrmion phase separation: jamming to laning to collapse"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000265,"raw_usage":{"total_tokens":1685,"prompt_tokens":1103,"completion_tokens":582,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":719,"completion_tokens_details":{"reasoning_tokens":508}},"tokens_in":719,"tokens_out":582,"duration_ms":5416,"temperature":1.0,"reasoning_tokens":508,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T20:32:54.986117+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"At skyrmion density $n_{sk} = 0.00081\\ \\mathrm{nm}^{-2}$ with one skyrmionium, the paper predicts the skyrmionium transverse velocity reverses sign near $j = 1.25 \\times 10^{9}\\ \\mathrm{A/m^{2}}$ when laning begins. A trajectory-resolved measurement or simulation that finds no such sign reversal, or no lane for densities up to $0.00189\\ \\mathrm{nm}^{-2}$, would contradict the core phase diagram.","supporting_citations":[{"cited_title":"Reichhardt , author J","cited_arxiv_id":null,"evidence_quote":"Provides the binary active-matter laning framework and the Magnus-dispersion treatment this work extends to magnetic textures."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the skyrmionium property of high velocity with no Hall effect, which creates the mobility contrast."},{"cited_title":"Comparing Dynamics, Pinning and Ratchet Effects for Skyrmionium, Skyrmions, and Antiskyrmions","cited_arxiv_id":"2412.02001","evidence_quote":"Provides the comparison of skyrmionium and skyrmion dynamics and interactions underpinning the density-dependent behavior."},{"cited_title":"Xia , author X","cited_arxiv_id":null,"evidence_quote":"Documents current-driven skyrmionium collapse into skyrmions, the process that terminates the laning state."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the atomistic spin-dynamics method used for the simulations."},{"cited_title":"Boulle , author J","cited_arxiv_id":null,"evidence_quote":"Provides the Pt/Co/MgO material parameters used to stabilize the simulated Neel skyrmions."}],"review_version":1}