{"id":"5f61dc32-2a19-451e-8393-d5ec5d18120d","arxiv_id":"2411.14160","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Simulations of a magnetized wire reveal new localized states, including asymmetric drifting double breathers and an intermittent 3-5-7 soliton complex, with multistability in the drive parameter plane.","lead":"This paper maps the localized magnetic states supported by a driven magnetic wire, including stationary solitons, breathing solitons, and drifting double breathers. The interest is in how balanced energy loss and external magnetic driving create stable multi-soliton structures.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Reported drifting-breather regime conflicts with the paper's own fast-mode cutoff; boundary/radiation contamination may explain the drift.","rationale":"The reader's weakest assumption was that the numerical setup distinguishes true attractors from long-lived transients and resolves the weak radiation controlling the drift. I agree with that concern, and the most concrete, internally checkable version of it is the inconsistency between the stated fast-mode cutoff (v>0.005) and the Zone-4 drifting breathers with speeds up to 0.032. The central claim would be true if the simulations were converged and boundary effects negligible, but the paper provides no dedicated resolution or domain-size study for the new states, and the boundary condition is reflective. The transient of 5x10^3 plus a 10^3 steady window is short relative to the 1/LLE timescales for the weakly chaotic breathers (LLE ~1e-4 implies a 10^4 timescale). The intermittent state is shown for a single parameter point with no quantitative period or LLE, so its 'very robust and regular' characterization rests on one 10^3-unit window. These are not grounds for rejection, because the reported states are plausible and the numerical framework is standard, but they are grounds for requiring additional verification before the stability diagram is treated as definitive. The proposed test directly checks whether drift velocity, direction, and soliton-count cycling survive higher resolution, larger domains, and longer runs; if they do, the concern is resolved, and if they do not, the paper's claim would need revision.","tokens_in":19964,"tokens_out":7728,"duration_ms":74232,"concrete_test":"Take a representative Zone-4 point, e.g. (nu,h0)=(-0.35,0.65), and the intermittent point (-0.20,0.45). Rerun each with dz=1/12 or 1/24 (or a pseudo-spectral spatial discretization), a domain extended to 2L=500 (or absorbing/periodic boundaries), and total integration time 5x10^4, measuring zCM(t), drift velocity, and soliton-count cycle in successive windows of length 10^3. Accept the intrinsic-drift claim only if velocity and direction are stable to <10% across windows and independent of domain size and resolution, and if the same states are not excluded by the v>0.005 cutoff. Also recompute the largest Lyapunov exponent over 2x10^4-10^5 with bootstrap uncertainty estimates to confirm the sign.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 4 sets a cutoff: 'Modes with the drift velocity exceeding 0.005 are considered as fast-moving ones and dismissed ... due to artifacts expected from their collisions with the boundaries.' Section 5.3 then classifies Zone 4, the only non-negligible-drift zone, as asymmetric breathing double solitons with average speeds in 10^-4 to 10^-1 and maximum 0.032. Most of these states violate the paper's own boundary-contamination cutoff, so the headline phenomenon of symmetry-broken drift is at risk of being a boundary/radiation-recycling artifact rather than an intrinsic attractor property. The steady-state diagnostic lasts only 10^3 time units after a 5x10^3 transient; with Neumann boundaries at z=±125, radiation emitted during the transient reflects and returns on timescales short compared with the run, and the radiation-recoil asymmetry that is supposed to produce the drift is exactly the quantity most sensitive to such reflections. The intermittent 3-5-7 soliton cycle is likewise documented at one parameter point over 10^3 units and could be a long-lived transient. The map is therefore a one-initial-condition basin snapshot, not a demonstrated existence/stability diagram.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript presents a numerical parameter-space study of a one-dimensional magnetic wire described by the Landau-Lifshitz-Gilbert equation with a DC plus AC perpendicular magnetic field. The authors scan the detuning-amplitude plane (ν, h0) with one symmetric two-soliton initial condition and classify the resulting localized states into standard solitons, breathers, double breathers, drifting asymmetric double breathers, fast-moving modes, and an intermittent multi-soliton complex exhibiting 3-5-7 soliton cycles. They characterize these states by power spectral densities, center-of-mass drift, micromagnetic energy, and the largest Lyapunov exponent, and report multistability among localized and delocalized responses.","tokens_in":20213,"tokens_out":4359,"duration_ms":43645,"significance":"If validated, the results would extend the known phenomenology of parametrically driven dissipative solitons from the PDNLS approximation to the full LLG model, and would document a spontaneous symmetry-breaking drift mechanism for double breathers as well as a regular intermittent multi-soliton cycle. The paper's strengths are that it uses the standard LLG equation without fitted parameters, performs a systematic scan with multiple diagnostics, and cross-checks the solver with an independent Python/Numba implementation. The main significance is contingent, however, on the numerical convergence and boundary-independence of the drifting and intermittent states, which are currently not established at the level needed for the paper's central claims.","major_comments":[{"comment":"The fast-mode exclusion criterion in Sec. 4 states that modes with drift velocity exceeding 0.005 are dismissed as artifacts of boundary collisions, yet Sec. 5.3 reports that Zone 4, the only zone with non-negligible drift, has average speeds in the range 10^-4 to 10^-1 with a maximum of 0.032. Since Zone 4 is the sole evidence for the spontaneously drifting asymmetric double breathers, this internal contradiction is load-bearing. The authors should either exclude these states by their own criterion, which would eliminate the headline phenomenon, or provide boundary-converged evidence that the drift is intrinsic, for example using longer runs, absorbing boundary conditions, or the ring geometry mentioned in the paper.","section":"Secs. 4 and 5.3"},{"comment":"The largest Lyapunov exponents in the range of about 10^-4 are reported without error bars or convergence tests; the value is averaged over the short window t = 1.8e4 to 2.0e4 of a single trajectory, with a sampling step of about 0.98. For a dissipative system with slow internal modulation and emitted radiation, finite-time exponents of this magnitude can be numerical artifacts or transient contamination. The claim that breathing solitons are weakly chaotic needs convergence checks with respect to the averaging window, the initial perturbation, and multiple trajectories.","section":"Sec. 5.5, Eq. (12)"},{"comment":"The phase diagram is generated from a single symmetric initial condition, Eq. (6), at each parameter point, so the colored regions are response maps for that particular initial condition rather than existence or stability regions of the states. The paper nevertheless labels Fig. 4 as 'Existence regions' and the abstract claims 'existence and stability areas.' This overstates the result; the authors should either rephrase the maps as basin-of-attraction maps for the specified initial condition or demonstrate that the states persist for other initial conditions.","section":"Sec. 5.1, Fig. 4(a)"},{"comment":"The intermittent 3-5-7 soliton complex is documented at one parameter point, (ν, h0) = (-0.20, 0.45), over a steady-state window of only 10^3 time units. Because the cycle period is not reported relative to this window, a long-lived transient cannot be excluded. The 'very robust and regular' characterization requires a run covering many cycles, preferably with additional initial conditions, and a comparison with the duration of the transient stage.","section":"Sec. 5.1, Fig. 5"},{"comment":"The steady-state diagnostics use spatial discretization dz = 1/6 and a steady-state window of 10^3 time units, but no convergence test at finer mesh or longer time is provided for the new states. This is especially relevant for the slow drift (speeds down to 10^-4) and for the weak radiation that is claimed to control the drift, because reflected radiation from the Neumann boundaries at z = ±125 can return on timescales comparable to the reported windows. Representative convergence checks for each dynamical zone would substantiate the claim that these are true attractors rather than long-lived numerical transients.","section":"Sec. 3 and Sec. 5"}],"minor_comments":[{"comment":"The caption labels panels as '3.c' and '3.d'; these should be '(c)' and '(d)' for consistency with the other panels.","section":"Fig. 3 caption"},{"comment":"The construction line is written as h0(ν) = 2.857ν + 1.729 in Fig. 4 but as h0(νh0) in Sec. 5.4; please define the parameterization consistently.","section":"Sec. 5.4 and Fig. 4 caption"},{"comment":"The peak-count criterion restricts attention to peaks at ω < 3Ω, but the number of peaks will depend on the noise floor and the peak-detection algorithm; these implementation details should be stated for reproducibility.","section":"Sec. 4, PSD counting"},{"comment":"There is a duplicated word in 'University of of Tarapacá'; it should be 'University of Tarapacá.'","section":"Acknowledgments"},{"comment":"The critical amplitude equation is introduced with hcrit0 and then used as h_0 in the following sentence; please make the notation uniform.","section":"Sec. 2, Eq. (5)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of nlin.PS and addresses a topic of current interest in parametrically driven magnetic solitons. The main risk is that the central claims about drifting asymmetric double breathers and the intermittent 3-5-7 soliton cycle are not yet supported by convergent, boundary-controlled numerics. The internal inconsistency between the fast-mode cutoff in Sec. 4 and the drift speeds in Zone 4 of Sec. 5.3 should be resolved before acceptance. I do not see evidence of any inappropriate practice, but the one-initial-condition nature of the maps and the short diagnostic windows should be clarified to the reader."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThis is a systematic numerical survey of localized states in a 1D Landau-Lifshitz-Gilbert wire with parametric drive. The genuinely new pieces are drifting, symmetry-broken double breathers and an intermittent 3-5-7 soliton complex, which I don't think have been reported before for the LLG setting. The authors map existence regions in the (ν, h0) plane on a fine grid and characterize the states with PSD, center-of-mass velocity, micromagnetic energy, and largest Lyapunov exponents. They also cross-check the algorithm with two independent codes. That is solid work, and the paper is a legitimate incremental extension of the earlier PDNLS studies and the authors' own single-soliton results.\n\nThe soft spot is where the stress-test note points. Section 4 dismisses modes with drift velocity above 0.005 as fast-moving and boundary-contaminated. Section 5.3 then presents a drift zone with average speeds up to 0.032, including the drifting double breathers that are the headline new result. The paper never explains why the cutoff doesn't apply to Zone 4; it only says the breathers stay far from the edges. But the drift mechanism is radiation recoil, and with Neumann boundaries at ±125 and a 5000-unit transient, emitted radiation has plenty of time to reflect off the boundaries and contaminate exactly the left-right recoil imbalance that supposedly produces the drift. The steady-state window is just 1000 units, short for separating true attractor drift from a slow transient. The intermittent 3-5-7 complex is documented at one parameter point over the same 1000-unit window, so a long-lived transient is not excluded.\n\nSmaller issues: the map is a basin snapshot for one initial condition rather than a general existence/stability diagram; the LLE values near 10^-4 come from a 2000-unit averaging window with no error bar; the classification thresholds are hand-set. None of these is fatal, but together with the cutoff inconsistency they make the drift and intermittency claims provisional.\n\nThis paper is for people working on dissipative magnetic solitons or the PDNLS-to-LLG correspondence. The parameter map and the diagnostics are useful, and the paper deserves a serious referee. I would send it out, but ask the authors to reconcile the cutoff inconsistency, run longer simulations (ideally with periodic boundary conditions) to rule out radiation-reflection artifacts, and test a few extra initial conditions in the disputed regions. That is a conditional acceptance, not a rejection.","headline":"Useful numerical map of LLG wire states, but the drifting-breather headline is undercut by the paper's own fast-mode cutoff and radiation-reflection concerns.","tokens_in":20768,"tokens_out":6007,"would_cite":false,"duration_ms":50742,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q51","35B36","37D45","82D40"],"pacs":["05.45.-a","75.78.-n"],"model":"deepseek-v4-flash","headline":"A one-dimensional magnetic wire driven by an AC field can hold soliton pairs in non-stationary bound states: symmetric and asymmetric breathers, drifting bound pairs, and an intermittent state that cycles through three, five, and seven…","keywords":["Landau-Lifshitz-Gilbert equation","breathers","soliton bound states","magnetic wire","parametric resonance","Lyapunov exponents","multistability","spontaneous symmetry breaking"],"falsifier":"Repeat the same parameter scan with a finer spatial mesh (for instance dz≤1/12), steady-state windows of at least $10^{4}$ time units, and a ring geometry with periodic boundary conditions; if the drifting asymmetric double breathers or the cyclic three-five-seven soliton state disappear or change character, they were numerical transients or boundary artifacts rather than intrinsic attractors.","tokens_in":19802,"feed_emoji":"🧲","tokens_out":9268,"duration_ms":76549,"temperature":0.7,"pith_summary":"The paper aims to establish that a one-dimensional magnetic wire described by the Landau-Lifshitz-Gilbert equation, with damping balanced by a perpendicular AC drive near the 2:1 parametric resonance, supports a richer family of localized bound states than previously catalogued. It reports that two solitons can bind into stationary pairs, breathing pairs, and asymmetric breathing pairs that drift steadily, and that a small parameter window hosts an intermittent multi-soliton complex whose soliton count cycles periodically through three, five, and seven. The authors also map out where each state exists in the plane of drive amplitude and frequency detuning, quantify the drift and the largest Lyapunov exponent of each regime, and attribute the drift to radiation recoil after spontaneous symmetry breaking. If true, this extends soliton-bound-state phenomenology from the abstract parametrically driven nonlinear Schrödinger equation to a concrete magnetic-wire model, and it shows that multiple stable attractors can coexist in the same parameter window.","feed_headline":"Driven magnetic wire cycles through 3, 5, and 7 solitons","feed_subtitle":"The same wire hosts drifting asymmetric breathers and a periodic state that switches between three, five, and seven solitons.","key_machinery":"The machine carrying the argument is the one-dimensional Landau-Lifshitz-Gilbert equation with a damping term, exchange and anisotropy energies, and a spatially uniform field containing DC and AC parts; near the 2:1 resonance this LLG model reduces to the parametrically driven damped nonlinear Schrödinger equation, which is the standard universal model for such forced dissipative systems. The numerical apparatus scans the (ν,h0) plane at fine steps, classifies states by their envelope dynamics and power-spectral density, measures center-of-mass drift by linear fits, and computes the largest Lyapunov exponent. The explanatory mechanism for the newly reported drift is radiation recoil: breathers emit low-amplitude dispersive waves, and once the bound pair's internal symmetry is broken, the left- and right-moving radiation fluxes differ, producing a net thrust.","core_discovery":"The central claim is that the driven damped LLG wire does not merely reproduce the known single- and double-soliton states; it also hosts non-stationary bound states of its own. Specifically, the paper reports symmetric double breathers, spontaneously asymmetric double breathers that emit unequal radiation and therefore drift, fast-moving double breathers that collide with the domain edges, and a localized intermittent complex, found near (ν,h0)=(−0.20,0.45), whose soliton number repeatedly cycles through seven, five, and three as solitons collide, merge, and split. The existence and stability regions for all of these states are laid out in the (ν,h0) plane under the first Arnold tongue, and multistability is documented by the coexistence of single-soliton, double-soliton, uniform, and pattern states for the same parameters. Lyapunov analysis assigns negative exponents to the regular double-soliton states and small positive exponents to breathing states, marking the latter as weakly chaotic.","pith_inferences":["Editorial inference: because the LLG wire reduces to the parametrically driven damped nonlinear Schrödinger equation near resonance, the same family of breathing bound states and the three-five-seven soliton cycle should be sought in other parametrically driven systems, such as Faraday-wave experiments, microcavity soliton combs, or coupled pendulum arrays.","Editorial inference: the drifting asymmetric double breather behaves like a nano-scale 'soliton motor' whose direction and speed encode the symmetry-breaking state; it could serve as a sensitive experimental probe of radiation emission in magnetic nanowires.","Editorial inference: the fuzzy boundaries between breathing-state regions hint at fractal basin boundaries; a basin-entropy or initial-condition-scan study would test whether the observed multistability is organized by a riddled basin structure."],"forward_implications":["The (ν,h0) existence map gives concrete drive-amplitude and detuning windows where each bound-state species is the attractor reached from a symmetric two-soliton initial condition.","The radiation-recoil mechanism makes drift a diagnostic: an asymmetric double breather's average velocity reports the degree of internal symmetry breaking.","The intermittent complex is periodic on the long timescale and cycles through three-, five-, and seven-soliton configurations without global drift, so the soliton number itself becomes a dynamical variable.","Multistability means the same wire parameters can sustain single-soliton, double-soliton, uniform, and pattern states; the initial condition selects which one is observed.","Along the studied parameter line, standard double solitons have negative largest Lyapunov exponents while breathing solitons and subharmonic patterns have positive ones, so the breathing regimes are weakly chaotic rather than quasiperiodic."],"supporting_citations":[{"why":"Supplies the magnetic-wire LLG setup, parameter values, and the two-soliton precession state from which the present bound states grow.","marker":"Urzagasti et al. [2012]"},{"why":"Defines the single-breather soliton regime in the same wire and the power-spectral-density criterion used to distinguish standard solitons from breathers.","marker":"Urzagasti et al. [2013]"},{"why":"Establishes the numerical accuracy of the spatial discretization and Runge-Kutta integration, and the soliton-antisoliton interaction baseline.","marker":"Urzagasti et al. [2014b, 2013]"},{"why":"Provides the parametrically driven damped nonlinear Schrödinger equation, the universal model whose phase-locked solitons underlie the 2:1 resonance physics.","marker":"Barashenkov et al. [1991]"},{"why":"Demonstrates stable two-soliton complexes in that model, the reference bound-state analog for the double-soliton states reported here.","marker":"Barashenkov and Zemlyanaya [1999]"}],"fun_headline_variants":["Magnetic wire breathers cycle through 3,5,7 solitons","Parametric drive makes magnetic wire breathers drift","Driven magnetic wire hosts chaotic drifting breathers","Soliton complexes in magnetic wire switch counts","Breather bound states show 3-5-7 soliton cycle"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The results rest on the assumption that the numerical scan, run on a mesh with dz≈1/6 and steady-state windows of about $10^{3}$ to 2×$10^{4}$ time units, distinguishes true long-lived attractors from transients and resolves the weak radiation that controls the breathers' drift.","fun_headline_variants_meta":{"raw":{"variants":["Magnetic wire breathers cycle through 3,5,7 solitons","Parametric drive makes magnetic wire breathers drift","Driven magnetic wire hosts chaotic drifting breathers","Soliton complexes in magnetic wire switch counts","Breather bound states show 3-5-7 soliton cycle"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000572,"raw_usage":{"total_tokens":2747,"prompt_tokens":1029,"completion_tokens":1718,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":645,"completion_tokens_details":{"reasoning_tokens":1634}},"tokens_in":645,"tokens_out":1718,"duration_ms":12746,"temperature":1.0,"reasoning_tokens":1634,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:27:46.714853+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the same parameter scan with a finer spatial mesh (for instance dz≤1/12), steady-state windows of at least $10^{4}$ time units, and a ring geometry with periodic boundary conditions; if the drifting asymmetric double breathers or the cyclic three-five-seven soliton state disappear or change character, they were numerical transients or boundary artifacts rather than intrinsic attractors.","supporting_citations":[],"review_version":1}