{"id":"69ae0c8c-37c9-4289-8d8f-0fc8a35c734c","arxiv_id":"2607.09881","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":5.5,"correctness_risk":"low","formal_verification":"none","parameter_count":3,"one_line_summary":"Piecewise elementary-function approximations to Stoner–Wohlfarth hysteresis loops (fixed-angle and randomly oriented) match numerical solutions to within ~0.05 except near switching.","lead":"Approximate closed-form formulas are derived for Stoner–Wohlfarth hysteresis loops of single particles and random ensembles. They replace numerical integration with elementary functions usable for nanoparticle and magnet modeling.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"No significant objection identified","rationale":"The reader correctly isolates the low-field expansion near the astroid and the hand-tuned matching points as the weakest assumptions, yet also correctly notes that both are quantified and do not invalidate the claim of “reasonably small” error. The paper’s own error plots (Figs. 6–7) already supply the concrete evidence that |δ| stays below 0.05 almost everywhere and that the chosen joins are near-optimal. Because the comparison is strictly internal (analytic approximation versus numerical solution of the identical SW equations) there is no external-data circularity. Consequently no load-bearing concern arises that would move the verdict away from ACCEPT; the recommended verification is merely a reproducibility check that the published error bounds survive independent re-implementation.","tokens_in":19112,"tokens_out":471,"duration_ms":5179,"concrete_test":"Independently recompute the absolute deviation |δ| of Eqs. 50–55 versus a high-resolution numerical quadrature of Eq. 24 on a dense h0 grid (Δh0≤0.01) over [−2,2]; confirm that max |δ| remains ≲0.05 except in a narrow neighborhood of h0=0.5 and that the hand-chosen join at |h0|=0.75 still yields the global minimum of the L∞ error among candidate joins in [0.6,0.9].","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the piecewise elementary formulas (fixed-θ_H: Eqs. 16–20; random ensemble: Eqs. 50–55) reproduce the exact numerical SW magnetization to absolute error ≲0.05 over most of the field range, with the largest (still “reasonably small”) deviation near h0≈0.5 for the averaged loop. The low-field small-angle expansions and the hand-chosen matching abscissae (h0=±1.5 fixed-angle, ±0.75 averaged) are the only soft points, but both are openly quantified in Figs. 6–7 and do not undermine the stated accuracy. No circular fitting to external data occurs; the comparison is internal to the same model. The argument therefore holds under scrutiny.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The paper derives approximate explicit piecewise formulas for Stoner–Wohlfarth hysteresis loops, both for a single particle at fixed easy-axis angle θ_H (Eqs. 16–20) and for an ensemble with randomly oriented anisotropy axes (Eqs. 50–55). The derivations rest on low-field small-angle expansions of the equilibrium condition (moment near the easy axis, Eqs. 8–9) and a high-field expansion (moment near the field direction, Eq. 12), with matching abscissae chosen to minimize global error against numerical solutions of the same model. Absolute deviations remain ≲0.05 over most of the field range, with the largest (still described as reasonably small) deviation near h_0≈0.5 for the averaged loop.","tokens_in":19329,"tokens_out":715,"duration_ms":11263,"significance":"If the stated accuracy holds, the formulas supply a practical, elementary-function alternative to numerical root-finding and angular integration for the classic SW model. They rest on a transparent physical picture rather than pure least-squares parametrization, recover known limits (initial susceptibilities, approach-to-saturation law), and cover both fixed-orientation and random-ensemble cases. The internal validation against exact numerical solutions of the SW equations, together with the explicit inverse-astroid formulas in Appendix A, makes the results immediately usable for modeling nanoparticles, thin films, and polycrystalline magnets without statistical sampling.","major_comments":[],"minor_comments":[{"comment":"In §3.1–3.2 the small-angle approximations are introduced without an a-priori error bound; a short remark quantifying the maximum |θ′| or |θ*| on the astroid (already plotted in Fig. 5) would help the reader judge the domain of validity before the numerical checks of §6.","section":null},{"comment":"The matching fields |h_0|=1.5 (fixed-θ_H) and |h_0|=0.75 (ensemble) are stated to be chosen for best fit, yet no sensitivity analysis is given. A single sentence or inset showing how the max |δ| changes if the join is moved by ±0.1 would strengthen the claim that the choice is robust.","section":null},{"comment":"Figure 2 caption and the surrounding text refer to “dotted curves” for the analytic formulas, but the line styles are hard to distinguish in grayscale; adding a legend or using dashed/solid consistently would improve readability.","section":null},{"comment":"Eq. (23) for the initial curve at fixed θ_H is left in a rather cumbersome product-of-sines form; a brief expansion for small h_0 would make the connection to the known χ_max = M_0²/(2K_a) more transparent.","section":null},{"comment":"Typographical: in the abstract and highlights the arXiv identifier appears as 2607.09881 while the header uses the same; consistency with the final journal citation style should be checked at proof stage.","section":null}],"recommendation":"accept","confidential_remarks":"The manuscript is a solid, self-contained contribution that fills a genuine practical gap. No novelty or citation concerns; it is a natural fit for a condensed-matter or magnetism journal that values usable analytic approximations."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"This is a practical paper that does exactly what it says. Savin and Koksharov give piecewise elementary formulas for Stoner–Wohlfarth hysteresis—both fixed easy-axis angle and the random-ensemble average—that stay within absolute error ~0.05 of the numerical solution over most of the field range. The largest deviation sits at h0≈0.5 on the averaged loop, right where the slope jumps, and they show it clearly in Fig. 7.\n\nWhat is new is the construction itself. They start from the low-field (moment near easy axis) and high-field (moment near field) limits of the energy stationarity condition, get closed expressions for θ*, then integrate analytically for the random case using the inverse-astroid roots they derive in the appendix. No least-squares fitting of the magnetization curves, no Monte-Carlo sampling. The matching points (|h0|=1.5 for fixed-angle, 0.75/0.5 for the ensemble) are chosen by hand for best visual fit, and the ln expansion is truncated at third order; those are free parameters, but they are few, openly stated, and the error plots quantify the residual. Prior analytic work either stayed at low field, used pure numerical parametrization, or stayed more general and heavier. This is lighter and covers the whole range.\n\nThe soft spots are real but proportional. The small-angle expansions are stretched near the astroid, and the matching abscissae are post-hoc. Both are discussed in Sections 3–6 and Figs. 5–7; the authors do not hide them. The comparison is internal to the same model, so there is no circularity with external data. Math checks out: the integrals are elementary and correctly evaluated, the high-field piece recovers the classical approach-to-saturation law, and the initial susceptibilities match known limits.\n\nThis is for people who write nanoparticle or permanent-magnet codes and want a fast, elementary SW kernel without calling a root finder every time. It does not open a new physical question, but it is a solid, reproducible convenience inside a mature subfield. I would send it to peer review without hesitation; a referee can ask them to label the matching points as optimized rather than uniquely derived, and that is about it. Worth citing if you need the formulas.","headline":"Clean, usable elementary approximations to SW loops that actually work; soft spots are openly quantified and do not sink the claim.","tokens_in":19896,"tokens_out":553,"would_cite":true,"duration_ms":5455,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Approximate elementary formulas reproduce Stoner–Wohlfarth hysteresis loops for single particles and random ensembles across the full field range.","keywords":["Stoner–Wohlfarth model","hysteresis loop","magnetic nanoparticles","single domain","analytical formulas","astroid","easy axis","random anisotropy"],"falsifier":"Direct numerical comparison of the piecewise formulas against the exact energy-minimization solution of the Stoner–Wohlfarth equations over a dense grid of field values and easy-axis angles; if the absolute magnetization error systematically exceeds ~0.05 outside the immediate vicinity of switching, or if the averaged-loop coercivity error grows beyond a few percent, the claim fails.","tokens_in":19999,"feed_emoji":"🧲","tokens_out":767,"duration_ms":7244,"temperature":0.7,"pith_summary":"The paper derives closed-form approximate expressions for the magnetization curves of Stoner–Wohlfarth particles. It covers both a particle whose easy axis is fixed relative to the applied field and an ensemble whose axes are randomly oriented. The authors start from a simple physical picture: at low field the moment stays near the easy axis, while at high field it stays near the field direction. Those small-angle expansions, joined at carefully chosen intermediate field values, produce piecewise formulas built only from elementary functions. The resulting expressions match numerical solutions of the classic model to within a few percent over most of the field range, the largest (still modest) discrepancy occurring near the sharp slope change of the averaged loop. The formulas therefore give a practical, integration-free alternative for calculating SW hysteresis without numerical root-finding or angular averaging.","feed_headline":"Elementary formulas give full SW hysteresis loops","feed_subtitle":"Low- and high-field expansions match numerical Stoner–Wohlfarth magnetization to a few percent","key_machinery":"The low-field expansions for the equilibrium angle θ* (Eqs. 8–9) and the high-field expansion (Eq. 12), joined at hand-chosen matching fields (h0 = ±1.5 for fixed-angle loops, h0 = ±0.75 for the averaged loop) that minimize global deviation from the numerical solution.","core_discovery":"Piecewise elementary formulas, obtained by matching low-field (moment near easy axis) and high-field (moment near applied field) expansions at optimized intermediate points, reproduce the full Stoner–Wohlfarth hysteresis loops for both fixed easy-axis orientation and randomly oriented ensembles, with absolute magnetization error remaining below roughly 0.05 except near the switching field.","pith_inferences":["The matching-field choice could be automated by a simple least-squares or minimax criterion, turning the present hand-tuned intervals into a fully algorithmic construction.","Because the formulas are elementary, they can be embedded directly in larger micromagnetic or finite-element codes that currently call numerical SW solvers at every mesh point.","The same low-/high-field matching strategy may extend to other uniaxial anisotropy models whose energy landscapes share a similar two-regime structure."],"forward_implications":["Hysteresis loops for fixed-angle and randomly oriented SW particles can be evaluated with only elementary functions, without root-finding or angular integration.","Initial susceptibilities of both single particles and random ensembles recover the known analytic limits of the model.","The same piecewise construction supplies an explicit initial-magnetization curve for the random ensemble.","The formulas remain usable near the switching field, where previous low-field approximations break down."],"fun_headline_variants":["Piecewise expansions match full SW hysteresis for fixed and random axes","Low-high field matching yields explicit Stoner-Wohlfarth loop formulas","Elementary piecewise formulas track SW magnetization with error below 0.05","Approximations capture SW loops near easy axis and applied field","Matched expansions reproduce Stoner-Wohlfarth hysteresis for ensembles"],"cache_read_input_tokens":16512,"weakest_assumption_plain":"The small-angle expansions that keep the magnetic moment near the easy axis remain accurate enough even near the Stoner–Wohlfarth switching curve, where the reduced field is no longer small.","fun_headline_variants_meta":{"raw":{"variants":["Piecewise expansions match full SW hysteresis for fixed and random axes","Low-high field matching yields explicit Stoner-Wohlfarth loop formulas","Elementary piecewise formulas track SW magnetization with error below 0.05","Approximations capture SW loops near easy axis and applied field","Matched expansions reproduce Stoner-Wohlfarth hysteresis for ensembles"]},"model":"grok-4.5","effort":"low","cost_usd":0.003984,"raw_usage":{"total_tokens":1166,"prompt_tokens":705,"num_sources_used":0,"completion_tokens":92,"cost_in_usd_ticks":39840000,"prompt_tokens_details":{"text_tokens":705,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":369,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":705,"tokens_out":92,"duration_ms":4161,"temperature":1.0,"reasoning_tokens":369,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-14T14:51:00.942464+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Direct numerical comparison of the piecewise formulas against the exact energy-minimization solution of the Stoner–Wohlfarth equations over a dense grid of field values and easy-axis angles; if the absolute magnetization error systematically exceeds ~0.05 outside the immediate vicinity of switching, or if the averaged-loop coercivity error grows beyond a few percent, the claim fails.","supporting_citations":[],"review_version":1}