{"id":"9f9af1d7-cbe6-446b-9b37-9c57b9bb1484","arxiv_id":"2411.18282","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"For a magnetar with a toroidal magnetic field confined to its crust, the paper derives the Zeeman-split magneto-elastic oscillation spectrum and gives simple formulas and fit constants for frequencies across stellar masses.","lead":"Magnetars are neutron stars that sometimes flare and show quasi-periodic oscillations, and this paper computes how a magnetic field reshapes their vibration spectrum. The authors show that a purely toroidal crustal field produces a simple Zeeman splitting that depends on the azimuthal quantum number squared, which may change how observed oscillations are identified.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"First-order perturbation theory is extrapolated to B0 ~ 10^15–3×10^15 G where the magnetic term is not small, and the one-parameter variational check in Sec. III E provides only an upper bound, not a certified error, so Eq. (31) is not quantitatively established in the QPO-fitting regime.","rationale":"I read the paper as making a well-defined claim within the stated toy model: for a purely toroidal crustal field vanishing at both boundaries, FOPT yields Eq. (31). The algebra behind Eq. (23) is internally consistent, and the angular integration leading to the m^2 factor checks out. The untouched m=0 elastic branch is a natural consequence of u×B=0 for axisymmetric torsional displacement. The paper is candid about the toy-model status of the field geometry and about the Alfvén-leakage limitation. The load-bearing weakness is not the physical realism of the field configuration but the quantitative regime of the perturbative expansion: the QPO fits use B0 values at which the magnetic term in Eq. (31) is of order the elastic term, so the first-order treatment is not controlled. The restricted variational estimate is a legitimate Rayleigh-quotient bound but cannot certify the FOPT error because both estimates are upper bounds, so the gap does not bound the error from below. This is an internal-consistency concern that a concrete non-perturbative eigenvalue calculation would settle. It does not change the CONDITIONAL verdict but sharpens the condition: Eq. (31) should be validated numerically before the quoted frequencies and B0 intervals are used for QPO identification.","tokens_in":21209,"tokens_out":10853,"duration_ms":99406,"concrete_test":"Perform a non-perturbative calculation of the fundamental (n=0) torsional magneto-elastic modes for the same 1.4 M_sun BSk21 crust model with ψ0(x)=16x^2(1−x)^2, imposing Y,r=0 at R1 and R2 and no external field, by directly solving the coupled radial eigenvalue problem obtained from ρH ω^2 u = T_µ[u] + T_B[u] (Eq. 3) with the full B0 sinθ ψ0(x) toroidal field. Compare the resulting ν_{0ℓm} at B0 = 1, 2, and 3×10^15 G with Eq. (31) for, e.g., (ℓ,m)=(2,2), (4,3), and (7,7). If any frequency differs by more than 5–10%, the FOPT extrapolation is not quantitatively valid, and the QPO/B0 claims in Figs. 8–9 should be revised to remain within a validated regime.","verdict_should_be":"UNCHANGED","load_bearing_attack":"At the field strengths used for the QPO fits (B0 ~ 2×10^15–3×10^15 G in Sec. IV, Figs. 8 and 9), the squared magnetic correction m^2 ν_B0*^2 B_*^2 is comparable to or larger than the elastic term (1/4)(ℓ+2)(ℓ−1)ν_µ0^2 in Eq. (31). For the ℓ=2, m=2 mode at B*=3, the magnetic term is 632 Hz^2 versus the elastic 532 Hz^2; for ℓ=7, m=7 at B*=2.7, the magnetic term is roughly 6300 Hz^2 versus 7200 Hz^2. Thus the usual smallness condition for first-order perturbation theory fails in exactly the regime where the paper claims quantitative frequencies and B0 constraints. The check in Sec. III E does not repair this: the trial family Y=1+w sin^2(πx/2) cannot change the angular dependence of the mode, couples to no other ℓ, and—since it is used in the Rayleigh quotient—produces only an upper bound on the true eigenvalue. The FOPT value Y=const is another upper bound, so the gap between them is not a lower bound on the error, and 'reasonable qualitative agreement' does not certify Eq. (31) to any stated accuracy. The paper itself acknowledges that mode crossings become quasi-crossings not captured by FOPT (Sec. III D), a separate manifestation of the same invalid-perturbation regime. The m^2 Zeeman structure is likely robust, but the specific frequencies, and therefore the illustrative identifications in Figs. 8–9, are not.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a first-order perturbation theory (FOPT) treatment of non-axisymmetric magneto-elastic oscillations of magnetar crusts containing a purely toroidal magnetic field confined to the crust. For fundamental modes (n=0) the authors derive the simple Zeeman-type formula Eq. (31), ν^2_{0ℓm} = (1/4)(ℓ+2)(ℓ−1)ν^2_{µ0} + m^2 ν^2_{B0*} B_*^2, with ν_{µ0}=23.06 Hz and ν_{B0*}=4.192 Hz for a 1.4 M_sun BSk21 model, and an analogous expression Eq. (35) for ordinary modes with radial nodes. They provide self-similarity fits over stellar masses, variational checks, and illustrative comparisons with QPOs from SGR 1900+14 and SGR 1806–20. The central qualitative claim is that the m^2 Zeeman term greatly enriches the oscillation spectrum and affects QPO interpretation.","tokens_in":21626,"tokens_out":5140,"duration_ms":67262,"significance":"The derivation of Eqs. (22)–(23) is algebraically consistent and transparent: the angular integral eliminates the Legendre dependence and leaves a simple radial integral. The m^2 dependence of the magnetic frequency is a clean, falsifiable prediction, and the paper is commendably explicit about the toy-model nature of the assumed field geometry. The self-similarity fits and the small set of auxiliary frequencies in Tables I and II are useful practical tools. If the FOPT result were quantitatively valid at the fields used for QPO fitting, the paper would give a very economical description of a large portion of the magneto-elastic spectrum. However, the quantitative accuracy of Eq. (31) in the regime B_0 ∼ 2–3×10^15 G, where the magnetic term is comparable to or larger than the elastic term, is not established by the checks presented in the manuscript.","major_comments":[{"comment":"The paper itself acknowledges in §III D that mode crossings at sufficiently high B_0 become quasi-crossings that are 'not described by our FOPT approach.' Since the illustrative fits in §IV use B_0 = (1.7–1.9)×10^15 G and B_0 = (2.7±0.12)×10^15 G, where many such crossings occur in Figs. 1 and 2, the specific ℓ,m identifications and B_0 constraints in Figs. 8–9 are not reliable within the stated approximations. In addition, the claim in §II E that 'under the formulated assumptions the sum rule (12) is exact' is misleading: Eq. (12) is derived by evaluating a Rayleigh quotient with unperturbed (zero-order) eigenfunctions, which gives only the first-order correction in B^2, not an exact result for the full problem. Clarifying the status of Eq. (12) is important because the paper explicitly extends it to fields where ω_B is comparable to ω_µ.","section":"§III A, Eq. (21), and §IV"}],"minor_comments":[{"comment":"Typo: 'nimber of radial nodes' should read 'number of radial nodes.'","section":"§II D"},{"comment":"Typo: 'with minis sign' should read 'with minus sign.'","section":"§II C"},{"comment":"Typo: 'These can can be numbered' should read 'These can be numbered.'","section":"§III F"},{"comment":"The frequency '8.38.6' should read '838.6'.","section":"Fig. 6 caption"},{"comment":"The name 'SGR 1860–20' appears twice in §IV; this should be 'SGR 1806–20' as used elsewhere.","section":"§IV"},{"comment":"The denominator notation '|u2|' is unclear; it should be the modulus squared of the vector displacement, e.g., |u|².","section":"§II E, Eq. (13)"},{"comment":"The QPO frequencies are shown as dotted lines on the figures, but the individual frequencies are not labeled on the plots; the text is clear, but labels would improve readability.","section":"§III D and Figs. 8–9"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid and clearly written extension of the authors' previous work on Zeeman splitting of magnetar oscillations. The main technical gap is the quantitative validity of FOPT at the high field strengths used in the QPO fits; this is fixable either by restricting the quantitative claims to the regime ω_B ≪ ω_µ, or by supplementing the variational check with a genuine error estimate (e.g., a full coupled-mode calculation or a two-parameter trial family that can mix ℓ). The paper's own caveats are honest, but the central quantitative formula, Eq. (31), is used at fields where its derivation does not apply, so a major revision is appropriate rather than acceptance in the present form."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take on Yakovlev & Fedorov, arXiv:2411.18282. The genuinely new thing is the m^2 Zeeman-splitting formula for torsional crustal modes when the field is purely toroidal and confined to the crust. They reduce the spectrum to two constants per stellar model, and the result — m=0 unsplit, magnetic frequency independent of ell for fundamentals — is not in the cited dipole-field papers. The algebra in Sec. III is consistent; I checked the theta integration and the reduction to the radial integral (Eqs. 22-23) is right. They are also honest about the toy model: toroidal field vanishing at both boundaries, no Alfven leakage, and they say it plainly. The self-similarity fits for 1-2.2 solar masses are a practical add-on.\n\nThe soft spots are exactly where the reader and stress-test put them. First, Eq. (31) is used at B0 ~ 2-3 x 10^15 G, where m^2 nu_B0*^2 B*^2 is comparable to or larger than the elastic term. First-order perturbation theory is not certified there. The variational check in Sec. III E is a restricted family (one parameter, no angular coupling) and gives only an upper bound; the gap between FOPT and variational values is not a rigorous error bar. The paper acknowledges quasi-crossings, which is the same regime problem showing up in another place. So the m^2 structure is probably robust, but the specific predicted frequencies in the QPO-fitting regime are not quantitatively established.\n\nSecond, the confinement of the field to the crust is an assumption designed to avoid Alfven leakage. Real magnetars will have core fields and non-zero boundary currents; that will change frequencies and introduce damping/leakage. The authors say this too, in Sec. V, but it means the QPO identifications in Figs. 8-9 are illustrative, not constraints. They pick B0 ranges by hand and leave some QPOs unexplained. That's fine as a demonstration, but it shouldn't be read as a measurement.\n\nWho's this for? Anyone working on magnetar QPOs or crustal oscillations. The paper gives a compact parametrization that observers can test against better data. It deserves a serious referee: the core derivation is clean, the limitations are disclosed, and the m^2 Zeeman effect is a real addition to the dipole-field results. I'd send it to peer review with a request to either soften the strong-field claims or add a better validity test, and to keep the toy-model framing in the abstract.","headline":"Clean analytic result for Zeeman splitting in toroidal-field magnetars, honest about its toy-model status; the QPO fits lean on perturbation theory beyond its validity, but the core formula deserves serious referee time.","tokens_in":22142,"tokens_out":1912,"would_cite":true,"duration_ms":17651,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["97.60.Jd","04.40.Dg","98.70.Rz","32.60.+i"],"model":"deepseek-v4-flash","headline":"A purely toroidal crustal magnetic field splits each magnetar oscillation mode into m-dependent Zeeman frequencies.","keywords":["magnetars","quasi-periodic oscillations","magneto-elastic oscillations","Zeeman splitting","toroidal magnetic field","neutron star crust","soft gamma repeaters","torsional oscillations"],"falsifier":"Compute the full non-axisymmetric magneto-elastic spectrum of a crust-confined toroidal field at $B_0 \\approx 2\\times10^{15}$ G in full general relativity without first-order perturbation theory; if any low-$\\ell$ mode frequency deviates from Eq. (31) by more than the variational spread shown in the paper's Fig. 5, the sum rule fails. Observationally, a clean Zeeman multiplet from one flare must obey the $m^2$ progression $\\nu^2_{0,2,2}-\\nu^2_{0,2,0}=4\\nu^2_{B0*}\\,B_*^2$ and $\\nu^2_{0,2,1}-\\nu^2_{0,2,0}=\\nu^2_{B0*}\\,B_*^2$; an observed spacing pattern inconsistent with these ratios would rule out the model.","tokens_in":20973,"feed_emoji":"🌌","tokens_out":14703,"duration_ms":102848,"temperature":0.7,"pith_summary":"The paper establishes that a purely toroidal magnetic field confined to the crust of a magnetar splits each torsional oscillation mode into a Zeeman family of $\\ell+1$ frequencies labeled by azimuthal number $m$, with the magnetic correction growing as $m^2$. For fundamental (nodeless) modes the whole spectrum is compressed into a two-number formula, $\\nu^2_{0\\ell m} = \\frac14(\\ell+2)(\\ell-1)\\nu^2_{\\mu0} + m^2\\nu^2_{B0*}\\,B_*^2$, with $\\nu_{\\mu0} = 23.06$ Hz and $\\nu_{B0*} = 4.192$ Hz for a $1.4\\,M_\\odot$ BSk21 star. The $m$-dependent splitting, the authors argue, makes the theoretical QPO spectrum much richer than the axisymmetric $m=0$ case and changes how observed quasi-periodic oscillations from magnetar flares should be matched to stellar models. The result matters because low-frequency QPOs are the main observable window into magnetar crusts and magnetic fields.","feed_headline":"Toroidal fields split magnetar oscillation modes into Zeeman families","feed_subtitle":"Each crustal mode gains ℓ+1 m-components, so observed QPOs can match theory at new frequencies.","key_machinery":"The central object is the first-order perturbation matrix $T_{m'm}$ built from the zero-order torsional eigenfunctions of Eqs. (6)$-$ (7); for axially symmetric fields it is diagonal in $m$, which converts the problem into the sum rule $\\omega^2 = \\omega_\\mu^2 + \\omega_B^2$ and reduces the magnetic correction to a ratio of radial integrals. For a toroidal field $B_\\phi = B_0\\sin\\theta\\,\\psi(x)$ that vanishes at both crust boundaries, the angular integrals eliminate the Legendre structure and leave $\\omega_B^2 = (1-x_g)\\, m^2B_0^2 \\int \\psi^2 |Y|^2 r^2\\,dr \\,/\\, (4\\pi \\int \\rho_H |Y|^2 r^4\\,dr)$, so the Zeeman term is exactly proportional to $m^2$ and to the crustal magnetic energy. The formula's power is that all angular dependence disappears: for fundamental modes the whole spectrum is set by two auxiliary frequencies, $\\nu_{\\mu0}$ and $\\nu_{B0*}$, and for $n>0$ by three constants per radial-node family, with self-similar mass scalings from a companion study carrying the results across $M=1$ $-$ $2.2\\,M_\\odot$.","core_discovery":"According to the paper, when the magnetic field is purely toroidal and vanishes at both crust boundaries, the first-order perturbation result is exact in the restricted problem: the squared frequency of any fundamental magneto-elastic mode is the sum of a pure shear term and a magnetic term, $\\omega^2 = \\omega_\\mu^2 + m^2\\omega_{B0*}^2\\,B_*^2$. For $B_\\phi = B_0\\sin\\theta\\,\\psi(x)$ with $\\psi(x)=16x^2(1-x)^2$, this yields the closed form Eq. (31), and the same structure extends to ordinary modes with radial nodes, where the zero-field fine structure $\\nu^2_{\\mu n\\ell} = \\nu^2_{\\mu n} + (\\ell+2)(\\ell-1)\\delta\\nu^2_{\\mu n}$ acquires an added $m^2\\nu^2_{Bn*}\\,B_*^2$ term. Because the perturbation matrix is diagonal in $m$ for axially symmetric fields, the quantum numbers $(n,\\ell,m)$ remain good labels, and modes with opposite signs of $m$ stay degenerate. The authors show that at $B_0$ of a few times $10^{15}$ G the Zeeman families of different $\\ell$ overlap and cross, producing dense frequency regions where QPO identification becomes ambiguous, and they use this to interpret observed QPOs from SGR 1900+14 and SGR 1806$-$20.","pith_inferences":["The $m^2$ scaling offers an observational lever to separate toroidal from poloidal fields: for poloidal crustal fields the magnetic frequency depends on $\\ell$ and does not vanish at $m=0$, whereas a purely toroidal field forces $\\omega_B=0$ at $m=0$, so a detected $m=0$ mode that shifts with field strength would point to a poloidal component.","If a future flare reveals a full Zeeman multiplet from a single $\\ell$, the derived product $\\nu_{B0*}\\,B_*$ could be combined with an independent crust mass and radius estimate to constrain the crustal magnetic energy and the radial profile $\\psi(x)$ through the coefficient $\\kappa$.","The analogy with the Paschen$-$Back effect suggests that at $B_0 \\gtrsim 10^{15}$ G the high-frequency fine structure reorganizes from ordering by $\\ell$ to ordering by $m$; searching observed QPO catalogues for this reordering would be a direct test, though the paper leaves this step to future work."],"forward_implications":["For a crust-confined toroidal field, every fundamental mode frequency is set by just two constants, $\\nu_{\\mu0}$ and $\\nu_{B0*}$, through Eq. (31), so measuring one member of a multiplet predicts the others.","The $m=0$ frequencies are independent of $B_0$, while $m>0$ frequencies grow as $m^2B_0^2$; at $B_0 \\gtrsim 10^{15}$ G the bunches for different $\\ell$ overlap and cross, creating dense allowed frequency regions.","High-frequency QPOs ($\\nu \\gtrsim 150$ Hz) fall into densely covered Zeeman domains even at lower fields, so identifying them with specific $(n,\\ell,m)$ modes becomes highly ambiguous.","Low-frequency QPOs from SGR 1900+14 can be matched with a $1.4\\,M_\\odot$ model at $B_0 \\approx (1.7$ $-$ $1.9)\\times10^{15}$ G, while the SGR 1806$-$20 hyperflare requires a more massive ($\\approx 2.2\\,M_\\odot$) star with $B_0 \\approx (2.7\\pm0.12)\\times10^{15}$ G."],"supporting_citations":[{"why":"First quantitative estimate of Zeeman splitting of magnetar QPO frequencies; the effect this paper extends to toroidal fields.","marker":"[48]"},{"why":"Provides the first-order perturbation-matrix formalism and the quadrature sum rule for squared frequencies used throughout.","marker":"[49]"},{"why":"Companion first-order study for dipole crustal fields; its QPO interpretation is compared with the toroidal-field results.","marker":"[50]"},{"why":"Supplies the self-similarity relations and fitted auxiliary frequencies for torsional oscillations used to extend results across stellar masses.","marker":"[26]"},{"why":"Supports the constant radial wave function for fundamental modes and the crustal microphysics used for the base torsion frequency.","marker":"[25]"},{"why":"Source of the general-relativistic elastic-mode equation used for zero-order torsional frequencies.","marker":"[13]"},{"why":"Provides the BSk21 unified equation of state and analytic stellar models used for the mass series.","marker":"[54]"},{"why":"Supplies the Coulomb-solids shear modulus expression used to compute the crustal elastic energy.","marker":"[55]"},{"why":"Summarizes the observed SGR QPO frequencies used for illustrative comparison in the paper.","marker":"[46]"},{"why":"Standard formalism for magneto-elastic oscillations that motivates the baseline wave equations.","marker":"[38]"}],"fun_headline_variants":["Zeeman families of crustal modes explain magnetar QPOs","Exact mode splitting in magnetars with toroidal fields","Toroidal field yields exact m-splitting in magnetar oscillations","Magnetar QPOs: Zeeman effect makes spectra richer","Toroidal fields: magnetar modes split into Zeeman multiplets"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the toroidal field drops to zero at both the outer surface and the crust-core interface, locking all magnetic energy in the solid crust; if real magnetars have fields that thread the core or carry boundary currents, the calculated frequencies, the absence of Alfvén leakage, and the QPO identifications would all change.","fun_headline_variants_meta":{"raw":{"variants":["Zeeman families of crustal modes explain magnetar QPOs","Exact mode splitting in magnetars with toroidal fields","Toroidal field yields exact m-splitting in magnetar oscillations","Magnetar QPOs: Zeeman effect makes spectra richer","Toroidal fields: magnetar modes split into Zeeman multiplets"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000606,"raw_usage":{"total_tokens":2839,"prompt_tokens":971,"completion_tokens":1868,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":587,"completion_tokens_details":{"reasoning_tokens":1779}},"tokens_in":587,"tokens_out":1868,"duration_ms":12914,"temperature":1.0,"reasoning_tokens":1779,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:21:39.460164+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the full non-axisymmetric magneto-elastic spectrum of a crust-confined toroidal field at $B_0 \\approx 2\\times10^{15}$ G in full general relativity without first-order perturbation theory; if any low-$\\ell$ mode frequency deviates from Eq. (31) by more than the variational spread shown in the paper's Fig. 5, the sum rule fails. Observationally, a clean Zeeman multiplet from one flare must obey the $m^2$ progression $\\nu^2_{0,2,2}-\\nu^2_{0,2,0}=4\\nu^2_{B0*}\\,B_*^2$ and $\\nu^2_{0,2,1}-\\nu^2_{0,2,0}=\\nu^2_{B0*}\\,B_*^2$; an observed spacing pattern inconsistent with these ratios would rule out the model.","supporting_citations":[{"cited_title":"What Magnetar Seismology can Teach us about the Magnetic Fields","cited_arxiv_id":"0903.3319","evidence_quote":"First quantitative estimate of Zeeman splitting of magnetar QPO frequencies; the effect this paper extends to toroidal fields."},{"cited_title":"Zeeman splitting of torsional oscillation frequencies of magnetars","cited_arxiv_id":"2312.10022","evidence_quote":"Provides the first-order perturbation-matrix formalism and the quadrature sum rule for squared frequencies used throughout."},{"cited_title":"Powerful flares and magneto-elastic oscillations of magnetars","cited_arxiv_id":"2409.11178","evidence_quote":"Companion first-order study for dipole crustal fields; its QPO interpretation is compared with the toroidal-field results."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Source of the general-relativistic elastic-mode equation used for zero-order torsional frequencies."},{"cited_title":"Magnetar Oscillations II: spectral method","cited_arxiv_id":"1110.2107","evidence_quote":"Standard formalism for magneto-elastic oscillations that motivates the baseline wave equations."}],"review_version":1}