{"id":"62567d3a-ea75-443b-b271-8b29f84e070d","arxiv_id":"2608.10902","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Flexural eigenfrequencies of a skyrmion string in an N-layer ferromagnetic film are |ω_m| = (4J'S/ħ) sin²(πm/(2N)), derived analytically and matched by spin-lattice simulations.","lead":"A skyrmion string in a ferromagnetic film of N atomic layers can bend and oscillate at discrete frequencies given by a simple formula involving the interlayer exchange and spin length. The paper derives this spectrum and confirms it in numerical simulations, giving experiments a concrete target for detecting and exciting these modes.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (22) is proven only for Belavin–Polyakov skyrmions with D=4πQ; the abstract and §IV present it as universal, and the numerical check does not quantify D or report peak positions.","rationale":"The reader's weakest-assumption analysis correctly identifies that Eq. (22) is the BP-skyrmion special case of Eq. (21). The derivation in §II is otherwise transparent: the Thiele-equation treatment with interlayer tension T=J′a²Dσ² and gyrovector G0=4πℏQσ is internally consistent, and the mode-counting with the discrete Laplacian is standard. However, the central claim as stated in the abstract—a parameter-free, DMI-independent resonance spectrum for a ferromagnetic film—depends on D/(4πQ) being exactly one. For a DMI-stabilized skyrmion this ratio can deviate from unity, and the paper provides no quantitative evidence that the deviation is negligible. The numerical section says only that peaks are 'in accordance with the theory' and that frequencies are weakly dependent on anisotropy and field; without the computed D factor, the measured peak positions, or error bars, the reader cannot verify whether the simulation is testing Eq. (22) or merely Eq. (21) with a nearly-BP profile. This is a fixable presentation and verification gap, not a fatal flaw: the underlying method and the general formula Eq. (21) appear sound, and the qualitative mode structure is reproduced. The paper should be conditionally accepted with the requirement that the D-dependence be stated explicitly and the numerical comparison be made quantitative. I therefore agree with the reader's conditional verdict and do not see grounds to move to rejection.","tokens_in":12143,"tokens_out":11162,"duration_ms":106857,"concrete_test":"Repeat the numerical protocol of §III with N=10, J′/J=0.2, A/J=0.2, H/J=-0.03, and the same anisotropy parameter D/J as used in the paper (which must be reported). For the relaxed single-skyrmion string, evaluate D = (1/2)∫[(∂m/∂x)²+(∂m/∂y)²]d²r from one layer's spin configuration using Eq. (5), and confirm Q=1. Then extract the first six resonance peaks from FS(ω) in Fig. 2 and compare them against Eq. (21) (with the measured D) and Eq. (22) (with D=4πQ). If the peaks match Eq. (21) but not Eq. (22), the universality claim in the abstract and §IV is refuted and must be revised; if D/(4πQ) is within the peak-width tolerance, reporting the numerical values settles the issue in the paper's favor.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central result |ω_m|=(4J′S/ℏ)sin²(πm/2N) is obtained by inserting D=4πQ into the general expression for the mode frequencies. Combining Eq. (21) with T=J′a²Dσ² and G0=4πℏQσ gives |Ω_m|=(J′S D/(πℏ|Q|))sin²(πm/2N), so Eq. (22) follows only when D=4π|Q|. For DMI-stabilized skyrmions in a conventional ferromagnet, D is generally larger than 4π|Q|, so the frequencies are not universal and carry an additional dependence on the skyrmion profile. The abstract states the formula without this restriction, and §IV calls it a 'rigorous analytical result' applicable to 'conventional ferromagnets possessing skyrmions,' which overstates the derivation. The simulations use a DMI-stabilized skyrmion, but the text reports only qualitative agreement and 'weak dependence' on anisotropy and field; it does not list the numerical peak frequencies, the value of D/J in the Hamiltonian (Eq. 55), or the actual D of the relaxed skyrmion. Without these numbers, the claim that Eq. (22) holds for the non-BP skyrmion is not quantitatively established; Eq. (21) is the defensible general statement, and Eq. (22) should be presented as its BP-skyrmion special case.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives flexural oscillation modes of a skyrmion string spanning N ferromagnetic atomic layers. Starting from Thiele equations for skyrmion centers and a nearest-layer coupling energy, the authors obtain Eq. (21), a discrete spectrum Ω_m = -(4T/G0) sin²(πm/(2N)), and then specialize to Belavin-Polyakov (BP) skyrmions, for which the gradient energy D equals 4π|Q|, yielding the headline formula |ω_m| = (4J'S/ℏ) sin²(πm/(2N)) (Eq. (22)). They also derive driven oscillations, resonant growth, damping saturation, and spin-orbit-torque excitation. Numerical spin-lattice simulations for a DMI-stabilized skyrmion string and skyrmion-string lattices (N=10) are presented as support, with fluctuation spectra of excited displacements shown in Figs. 2-4, 6-7.","tokens_in":12381,"tokens_out":7430,"duration_ms":70278,"significance":"The analytical derivation is transparent and uses standard tools: Thiele dynamics, a harmonic interlayer tension, and diagonalization of the discrete Laplacian. No parameters are fitted to the predicted frequencies, and the spin-lattice simulations are an independent numerical test. The general result, Eq. (21), is a useful and testable dispersion relation for flexural modes of skyrmion strings. The BP special case, Eq. (22), is clean, parameter-free, and potentially important for experiments on synthetic multilayers. However, the paper's headline claim applies Eq. (22) to 'conventional ferromagnets possessing skyrmions' in the abstract and Section IV, which is not justified by the derivation unless the skyrmion profile is of BP type with D=4π|Q|. The numerical confirmation is also reported only qualitatively. With appropriate restriction of the claim and quantitative comparisons, the paper would be a solid contribution.","major_comments":[{"comment":"Eq. (21) reads Ω_m = -4T/G0 sin²(πm/(2N)) = -(J'S D/(πℏ Q)) sin²(πm/(2N)) after inserting T=J'a²Dσ², G0=4πℏQσ, and σ=S/a². Eq. (22) follows only when D=4π|Q|, which, as stated after Eq. (5), holds for a pure-exchange Belavin-Polyakov skyrmion. DMI-stabilized skyrmions in a conventional ferromagnet have a different profile, so D generally differs from 4π|Q| and the frequencies are not universal. The abstract and Section IV therefore overstate the result: Eq. (22) should be presented as the BP-skyrmion special case, with Eq. (21) as the general expression, and the phrase 'rigorous analytical result ... applies to conventional ferromagnets possessing skyrmions' should be amended accordingly.","section":"Section IIB, Eqs. (21)-(22); Section IV; Abstract"},{"comment":"The numerical verification is described only as 'good agreement'; the manuscript does not report the numerical peak frequencies, does not give the value of the uniaxial anisotropy D/J appearing in Eq. (55) (the main parameter set lists only J'/J, A/J, and H/J), and does not measure D of the relaxed skyrmion profile. Since the simulated skyrmion is DMI-stabilized and therefore not a BP skyrmion, agreement with Eq. (22) could instead be agreement with Eq. (21) for a D close to 4πQ. The authors should provide a table comparing computed peak positions with Eqs. (21) and (22) and report a numerical estimate of D for the relaxed skyrmion; without these data the confirmation claim is not quantitatively established.","section":"Section III, Eq. (55) and Figures 2-4, 6-7"},{"comment":"The abstract states that the result is 'confirmed numerically for ... skyrmion-string lattices at finite temperature,' but Section IV states that at elevated temperature 'only the lowest flexural resonances remain clearly identifiable, while higher modes are obscured by hybridization with lattice and spin-wave excitations,' and Fig. 7 indeed shows additional peaks to the right of the theoretical positions. This is an inconsistency between the abstract's confirmation claim and the actual numerical evidence; the abstract should be qualified to state that at finite temperature only the low-m modes are confirmed for the lattice case.","section":"Abstract and Section IV; Figures 6-7"}],"minor_comments":[{"comment":"The sentence 'Bending of skyrmion strings under a thermal gradient has been observed by Ran et al. [30]' is duplicated.","section":"Introduction"},{"comment":"The Introduction includes m=0 in the range m=0,1,...,N-1, while the abstract and Eq. (22) state m=1,...,N-1; since the m=0 mode is the zero-frequency rigid translation, the range should be stated consistently.","section":"Introduction vs. Abstract and Eq. (22)"},{"comment":"The text 'k_m = mπ/ℏ' appears to be a typo; the continuous wave number for a film of thickness b should be k_m = mπ/b, and with a=b=1 this is mπ/N.","section":"After Eq. (28)"},{"comment":"The skyrmion-locator formula is typeset in a way that obscures the intended ratio; it should read R = Σ_{s_{i,z}>0} r_i s_{i,z}^2 / Σ_{s_{i,z}>0} s_{i,z}^2.","section":"Eq. (58)"},{"comment":"The notation in Eq. (43) uses j both as the current-density vector and as a free index; the expression should be written with an explicit summation over the Cartesian index, e.g., I_{ij}[e_n × j(t)]_j.","section":"Eq. (43)"}],"recommendation":"major_revision","confidential_remarks":"The central derivation is sound but the presentation overgeneralizes the BP-skyrmion result, and the numerical validation is not quantitative. Both issues are fixable within the manuscript's scope. The self-citation pattern is not problematic for this field."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper gives a clean analytical spectrum for flexural modes of a skyrmion string in an N-layer ferromagnetic film: |omega_m| = (4J'S/hbar) sin^2(pi m/2N). That formula is new for finite N, as far as I can tell from the cited literature. The derivation is standard Thiele dynamics plus the discrete Laplacian, and it is done carefully. The generalization to driven oscillations, including the spin-orbit-torque excitation and the power absorption estimate, is a useful extra that makes the result concrete for experiments.\n\nWhat I like: the two-layer case is worked out first, the discrete spectrum follows from a textbook eigenvalue problem, and the continuous limit is checked against the helical-wave picture. The numerics use an independent spin-lattice model with DMI and finite temperature, which is the right way to test the theory. The agreement shown in the figures looks genuine.\n\nThe soft spots are real but not fatal. The abstract and Section IV present Eq. (22) as universal, but it only follows from Eq. (21) when the skyrmion profile is Belavin-Polyakov, D = 4pi|Q|. For DMI-stabilized skyrmions, which is what the simulations actually use, the frequencies carry the factor D/(4piQ). Eq. (21) is the defensible general statement; Eq. (22) should be labeled as the BP special case. That is a presentation problem, not a mathematical one.\n\nThe numerical section is also under-quantified. The text says there is good agreement but does not list the measured peak frequencies, give error bars, or state the anisotropy constant D/J used in the Hamiltonian. The relaxed skyrmion's D is never computed. Without those numbers, the claim that Eq. (22) holds for the non-BP DMI skyrmion is not established. Given the theory is parameter-free, a simple comparison table would settle it.\n\nMinor issues: the anisotropy value is missing but the model is reproducible from the other parameters; the self-citations to the authors' earlier methods are appropriate and do not carry the derivation.\n\nBottom line: this is a solid subfield contribution. The derivation is transparent, the simulations are real, and the overstatement in the abstract is fixable. A serious referee can sort out the D dependence and ask for quantitative comparisons. I would send it to review.","headline":"Useful analytical result for skyrmion-string flexural modes; the clean BP-skyrmion formula in the abstract needs a caveat, and the numerics need quantitative reporting, but the core derivation is solid.","tokens_in":12984,"tokens_out":1041,"would_cite":true,"duration_ms":12116,"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 skyrmion string in an N-layer ferromagnetic film bends at eigenfrequencies $(4J'S/\\hbar)\\sin^2(\\pi m/(2N))$, as derived and confirmed by simulation.","keywords":["skyrmion string","flexural oscillations","ferromagnetic film","interlayer exchange","Thiele dynamics","Belavin-Polyakov skyrmion","microwave absorption","spin-lattice simulation"],"falsifier":"Measure the flexural resonance peaks of a single skyrmion string while tuning the DMI strength (or any parameter that changes $D/Q$) with $J'$ and $N$ fixed: Eq. (22) predicts no shift, whereas Eq. (21) predicts every peak moving proportionally to $D/(4\\pi|Q|)$. A numerical run with a skyrmion profile whose $D/Q$ is, say, half of $4\\pi$ should show all peaks at half the frequencies of Eq. (22).","tokens_in":11849,"feed_emoji":"🌀","tokens_out":7149,"duration_ms":67471,"temperature":0.7,"pith_summary":"The paper asks how a skyrmion string, meaning a vertical stack of one two-dimensional skyrmion per atomic layer, bends when its layers shift sideways. It derives a closed spectrum for the flexural eigenmodes of an $N$-layer string: $|\\omega_m|=(4J'S/\\hbar)\\sin^2[\\pi m/(2N)]$ for $m=1,\\dots,N-1$, where $J'$ is the interlayer exchange and $S$ the spin length. The frequencies depend on the layer count and the interlayer coupling but not on the skyrmion's size, anisotropy, or magnetic field, so the spectrum is a clean fingerprint of string bending stiffness. The paper reports spin-lattice simulations for individual strings and for strings in a lattice that reproduce the predicted peaks, including at finite temperature. If the result holds, it gives an analytically exact and nearly parameter-free description of a degree of freedom previously treated mostly numerically.","feed_headline":"Flexural modes of skyrmion strings follow a sin-squared ladder","feed_subtitle":"In an N-layer ferromagnet, bending resonances depend only on interlayer exchange and N, and spin-lattice runs agree.","key_machinery":"The machinery is a discrete chain of skyrmion-center displacements $d_n=R_n-R_{n+1}$ governed by linearized Thiele equations. Interlayer exchange contributes an energy $U=(T/2)\\sum_n|d_n|^2$, with tension $T=J'a^2D\\sigma^2$, and the gyrovector $G_0=4\\pi\\hbar Q\\sigma$ turns the restoring force into circular motion rather than simple relaxation. Diagonalizing the discrete Laplacian with free boundaries gives eigenvectors $\\sin(\\pi mn/N)$ and eigenvalues $\\lambda_m=4\\sin^2[\\pi m/(2N)]$; substituting those into the Thiele equation yields $\\Omega_m=-T\\lambda_m/G_0$, and using $D=4\\pi|Q|$ produces the final formula. This maps the many-spin problem onto the normal modes of a one-dimensional lattice of coupled gyroscopes.","core_discovery":"The central claim is that flexural oscillations of a skyrmion string in a ferromagnetic film are discrete sinusoidal modes of the relative displacement between neighboring layers, with eigenfrequencies given by Eq. (22): $|\\omega_m|=(4J'S/\\hbar)\\sin^2[\\pi m/(2N)]$. The derivation starts from the general expression $|\\Omega_m|=(J'S D/(\\pi \\hbar Q))\\sin^2[\\pi m/(2N)]$, where $D$ is the gradient energy of the in-layer skyrmion and $Q$ its topological charge. For a pure-exchange Belavin-Polyakov skyrmion, $D=4\\pi|Q|$, and the general expression reduces to the parameter-free form. The paper verifies this analytic ladder numerically for up to ten atomic layers, both for a single skyrmion string and for a lattice of skyrmion strings, using the skyrmion-center displacement as the dynamical variable.","pith_inferences":["An implication the paper leaves implicit is that any measured deviation from the $\\sin^2$ ladder in a DMI-stabilized skyrmion directly measures the ratio $D/(4\\pi|Q|)$, effectively turning the flexural spectrum into a probe of the skyrmion's internal profile.","The same coupled-gyroscope chain should describe bending of other columnar magnetic textures, such as bubble strings, with the only input changing being the ratio $D/Q$ that sets the effective tension.","One could extend the analysis to asymmetric multilayers with different spin densities in different layers; the derivation shows the effective circling frequency then depends on $1/G_1+1/G_2$, so the symmetric formula Eq. (22) would acquire a layer-dependent correction."],"forward_implications":["Every flexural mode $m=1,\\dots,N-1$ can be excited by a surface spin-orbit torque applied to the bottom layer, with resonant growth of the displacement and resonant microwave absorption at the predicted frequencies.","In synthetic multilayers the resonance ladder can be tuned by choosing the interlayer coupling $J'$, and the spacing between successive modes is set by the layer number $N$ through $\\sin^2[\\pi m/(2N)]$.","The flexural modes provide an interface-addressable channel of microwave absorption distinct from perpendicular standing spin waves, with absorbed power per unit area depending only weakly on skyrmion size.","At low temperature the analytic spectrum survives in a skyrmion-string lattice; at elevated temperature the low-lying modes remain identifiable while higher modes hybridize with lattice and spin-wave excitations.","Undamped resonant driving makes the oscillation amplitude grow linearly in time and the pumped energy grow as $t^2$, which can destabilize the string unless dissipation limits the amplitude."],"supporting_citations":[{"why":"Supplies the Belavin-Polyakov skyrmion profile for which $D=4\\pi|Q|$, the identity that converts the general frequency formula into the parameter-free Eq. (22).","marker":"[34]"},{"why":"Supplies the Thiele equations of skyrmion motion used to write the layer dynamics and the flexural mode equation.","marker":"[35]"},{"why":"Provides the sine eigenvectors of the discrete Laplacian that are the normal modes of the $N$-layer chain.","marker":"[36]"},{"why":"The displacement-of-skyrmion-centers method whose generalization yields the string flexural equations.","marker":"[24]"},{"why":"Companion treatment of skyrmion-center displacements in ferrimagnetic films, part of the same derivation chain.","marker":"[25]"},{"why":"Supplies the skyrmion-locator formula used to track the layer-by-layer centers in the numerical simulations.","marker":"[40]"},{"why":"Supplies the formula for the effective skyrmion size $\\lambda_{\\rm eff}$ used to characterize the numerically studied strings.","marker":"[42]"}],"fun_headline_variants":["Skyrmion string flexing locks to sin^2(mπ/2N)","Layer number N controls bending frequencies of skyrmion strings","Exact sin^2 ladder for skyrmion string vibrations in films","Ferromagnetic film skyrmion strings: eigenfrequencies from exchange only","Skyrmion string bending modes depend only on N and interlayer exchange"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The clean spectrum relies on every layer's skyrmion being a pure-exchange Belavin-Polyakov profile whose gradient energy is exactly $4\\pi$ times its topological charge; for conventional DMI-stabilized skyrmions that relation fails and the resonance frequencies shift.","fun_headline_variants_meta":{"raw":{"variants":["Skyrmion string flexing locks to sin^2(mπ/2N)","Layer number N controls bending frequencies of skyrmion strings","Exact sin^2 ladder for skyrmion string vibrations in films","Ferromagnetic film skyrmion strings: eigenfrequencies from exchange only","Skyrmion string bending modes depend only on N and interlayer exchange"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000341,"raw_usage":{"total_tokens":1829,"prompt_tokens":847,"completion_tokens":982,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":463,"completion_tokens_details":{"reasoning_tokens":886}},"tokens_in":463,"tokens_out":982,"duration_ms":8823,"temperature":1.0,"reasoning_tokens":886,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T14:40:47.019018+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the flexural resonance peaks of a single skyrmion string while tuning the DMI strength (or any parameter that changes $D/Q$) with $J'$ and $N$ fixed: Eq. (22) predicts no shift, whereas Eq. (21) predicts every peak moving proportionally to $D/(4\\pi|Q|)$. A numerical run with a skyrmion profile whose $D/Q$ is, say, half of $4\\pi$ should show all peaks at half the frequencies of Eq. (22).","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Belavin-Polyakov skyrmion profile for which $D=4\\pi|Q|$, the identity that converts the general frequency formula into the parameter-free Eq. (22)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Thiele equations of skyrmion motion used to write the layer dynamics and the flexural mode equation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the sine eigenvectors of the discrete Laplacian that are the normal modes of the $N$-layer chain."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The displacement-of-skyrmion-centers method whose generalization yields the string flexural equations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Companion treatment of skyrmion-center displacements in ferrimagnetic films, part of the same derivation chain."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the skyrmion-locator formula used to track the layer-by-layer centers in the numerical simulations."},{"cited_title":"Kachkachi and D","cited_arxiv_id":null,"evidence_quote":"Supplies the formula for the effective skyrmion size $\\lambda_{\\rm eff}$ used to characterize the numerically studied strings."}],"review_version":1}