{"id":"29e34b53-786e-4476-9887-1783d21da5f9","arxiv_id":"2507.02154","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Micromagnetic simulations show that disordered iron-oxide nanoflowers exhibit a secondary coercivity maximum in the vortex state, at the transition between core-dominated and flux-closure-dominated reversal.","lead":"Iron-oxide nanoflowers used in magnetic hyperthermia were studied in large-scale micromagnetic simulations across 10 to 400 nm sizes. The simulations reveal a size 'sweet spot' around 110 nm where magnetic switching changes mode and coercivity peaks, suggesting a design rule for more efficient nanoheaters.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'tailorable hyperthermia sweet spot' rests on equating zero-temperature DC coercivity with AC loss; no hysteresis-loop area or SAR is computed, so the practical peak may be an artifact of the proxy.","rationale":"The coercivity-size map and the two-mode reversal description are internally consistent and supported by baseline comparisons (k = 1, ideal nanospheres). My concern is not with the micromagnetic result itself but with its translation into hyperthermia guidance. The paper's title, abstract, introduction, and conclusion all present the coercivity peak as the 'hyperthermia sweet spot,' and the Introduction explicitly equates heating performance with coercivity. That identification is load-bearing: if it fails, the practical claims fail even though the fundamental coercivity map may be correct. The reader's weakest_assumption identifies the same issue, so I agree with the conditional verdict. I also considered the statistical overlap of the two reversal modes noted in footnote 68 — both processes occur in an overlapping size range, with the peak identified with whichever process dominates — which could complicate the sharp V_core/V_NF = 1/3 transition language, but the hyperthermia proxy is more consequential because it is the stated purpose of the paper. The proposed test is feasible: the authors already have full hysteresis loops, so integrating loop area versus field amplitude is a post-processing step, and a full AC/thermal check is a natural extension of their Mumax3 framework. No change to the reader's conditional verdict is needed.","tokens_in":14633,"tokens_out":6923,"duration_ms":81909,"concrete_test":"Compute loop areas from the already-simulated M(H) curves (or run stochastic-LLG AC simulations at T = 300 K) for d = 70, 90, 110, 130, 150, and 200 nm with k = 0.25 and 15 nm grains, integrating ∮M·dH over sinusoidal fields with amplitudes μ0H0 = 5, 10, 15, 20, and 30 mT at f = 100–300 kHz. If the loop-area peak diameter differs from the d ≈ 110 nm coercivity peak by more than ~20% for any clinically relevant (H0, f), or if no area peak exists, the 'hyperthermia sweet spot' claim must be revised. A minimal version of this test is to integrate minor-loop areas from the existing zero-T full hysteresis loops up to ±H0, which directly checks whether area tracks coercivity.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central practical claim — that the coercivity peak tailors the hyperthermia optimum — depends on identifying heating performance with coercivity. The text states this identification explicitly (Introduction: 'heating performance, represented by the coercivity') and the abstract/conclusion generalize it to a 'hyperthermia sweet spot.' But MHT heating is proportional to the hysteresis-loop area ∮M·dH at the applied frequency and field amplitude, not to the zero-temperature quasi-static coercivity. The paper reports only μ0H_C from DC, 0 K simulations (Fig. 1) and never computes loop area, SAR, or finite-temperature AC response. This matters because the two vortex reversal modes have very different loop shapes: in the flux-closure regime (Fig. 2B) the magnetization reverses gradually over a wide field interval, so a lower-coercivity loop could still enclose substantial area; conversely, at the claimed sweet spot d ≈ 110 nm, μ0H_C ≈ 14 mT, and clinical hyperthermia fields are typically μ0H0 ≈ 10–30 mT at 100–500 kHz. If H0 < H_C, only minor loops are traversed and the enclosed area can be small, while a larger diameter with lower H_C but a squarer or broader minor loop could dissipate more. Thermal fluctuations at 300 K will further suppress the zero-temperature coercive barrier, especially near 100 nm, potentially shifting or erasing the peak. Thus the DC coercivity maximum and the SAR maximum need not coincide, and the practical nanoheater guidance is currently unsupported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents large-scale micromagnetic simulations of iron-oxide nanoflowers from d = 10 nm to 400 nm, using a Voronoi grain structure with random anisotropy axes and a reduced inter-grain exchange factor k to model spin disorder. It reports a coercivity-versus-size map with three regimes: a Stoner-Wohlfarth coherent-rotation regime below roughly 50 nm, a domain-wall-nucleation regime between about 50 and 75 nm, and a vortex regime above about 70 nm. In the vortex regime the authors find a non-monotonic coercivity with a secondary maximum near d = 110 nm for k = 0.25, which is absent in defect-free nanospheres. They attribute the peak to the interplay of grain anisotropy and grain-boundary pinning, identify a Vcore/VNF = 1/3 threshold separating core-dominated from flux-closure-dominated reversal, and conclude that the grain size can tailor this coercivity 'sweet spot' for magnetic hyperthermia.","tokens_in":14935,"tokens_out":3860,"duration_ms":47444,"significance":"If the central mechanism is correct, the work is a valuable step beyond macrospin descriptions of nanoflowers: it provides a coercivity-size map over a wide diameter range, uses 15-25 random realizations per size, includes a clear baseline comparison with fully coupled nanoflowers and defect-free nanospheres, and proposes a concrete geometric threshold (Vcore/VNF = 1/3) for the reversal-mode crossover. The separation of reversal modes and the demonstration that grain disorder rather than mere size creates the secondary coercivity peak are physically interesting and testable. However, the practical claim that the coercivity peak marks a 'hyperthermia sweet spot' is not supported by the simulations as presented, because heating performance in magnetic hyperthermia depends on the hysteresis-loop area under the applied AC field and frequency, not on the zero-temperature DC coercivity.","major_comments":[{"comment":"The paper's central practical claim — that the coercivity peak located at the reversal-mode transition tailors the 'hyperthermia sweet spot' — is unsupported because heating performance is identified with zero-temperature DC coercivity without computing the AC hysteresis-loop area or specific absorption rate. The Introduction states that heating performance is 'represented by the coercivity,' but MHT losses are proportional to ∮M·dH at the applied field amplitude and frequency, and clinical fields (typically μ0H0 ≈ 10–30 mT at 100–500 kHz) may be below the reported peak coercivity of about 14 mT, so minor-loop and loop-shape effects can decouple coercivity from dissipated power. Please compute loop areas or SAR at relevant AC conditions (or at least at finite temperature for representative sizes), or restrict the conclusions to coercivity rather than hyperthermia efficiency.","section":"Section II.A and Conclusion"},{"comment":"The Vcore/VNF = 1/3 threshold is asserted for nanoflowers, but the quantitative core-volume data shown in Fig. 2(E) are for the ideal nanosphere only. The sentence 'In the NF, we find the same Vcore/VNF = 1/3 threshold' is not backed by a corresponding plot or tabulated NF core-volume analysis in the main text or the available supplement. Because this threshold is load-bearing for the explanation of the coercivity peak location, the authors should present NF core-volume fractions versus diameter, or explicitly state where the supporting data appear.","section":"Section II.B and Fig. 2(E)"},{"comment":"The absolute coercivity scale and the location of the sweet spot depend on the inter-grain coupling k = 0.25, which is selected to match experimental coercivity values [19,44], but no quantitative comparison between simulated and measured coercivities is shown, and no sensitivity analysis for k is reported. Since the abstract and Fig. 1 emphasize 'excellent agreement' and the practical field-amplitude comparisons use the simulated peak height, please provide the experimental comparison plot, the statistical spread over the 15-25 random realizations, and at least a brief statement of how the peak height and location vary over the plausible range of k and material parameters.","section":"Section II.A and Methods"}],"minor_comments":[{"comment":"The text contains a typo: 'Stoner-Wolhfarth' should be 'Stoner-Wohlfarth'.","section":"Section II.A"},{"comment":"The caption contains 'nmnm' and should read 'nm'.","section":"Fig. 1(A) caption"},{"comment":"The word 'offeres' should be 'offers'.","section":"Conclusion"},{"comment":"The d^-1.85 power law is fit to a very small number of points in the 50-75 nm range; please report the fit range, R^2, and standard error, and avoid overinterpreting the difference from the d^-1.5 expectation.","section":"Section II.A"},{"comment":"The text describes Fig. 2(C) for d = 120 nm, whereas the caption states d = 110 nm; this inconsistency should be corrected.","section":"Fig. 2 caption and Section II.B"},{"comment":"The statement that 'all data supporting the findings are included' would be more useful if simulation scripts or a repository link were provided, given that the results depend on specific random Voronoi tessellations and Mumax3 settings.","section":"Data Availability"}],"recommendation":"major_revision","confidential_remarks":"The core micromagnetic study is sound and the reversal-mode classification is a useful contribution, but the title and conclusions overreach by making an uncomputed claim about hyperthermia performance. Adding AC loop-area simulations or consistently reframing the results as a coercivity map without heating-efficiency claims would address the main concern. The paper is otherwise within scope for cond-mat.mes-hall."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThis one is worth reading: it maps coercivity versus size for disordered iron-oxide nanoflowers from 10 to 400 nm, and it finds something genuinely new. Vortex states do not simply give monotonically decreasing coercivity. The coercivity rises to a secondary peak around 110 nm (for the k=0.25 case) and then falls again. The paper ties that peak to a change in reversal mechanism: core-immediate reversal when the vortex core volume exceeds about one-third of the particle volume, and flux-closure rotation when the core is smaller. The Vcore/VNF = 1/3 threshold is a clean, physical criterion, and the comparison with defect-free nanospheres is exactly the right control to show that the peak comes from grain disorder.\n\nThe simulations look careful: 15-25 random Voronoi realizations per size, reasonable material parameters, and a check that the grain structure does not affect the small-size single-domain regime. The k-sweep in Figure 3 gives a nice mechanistic story for why intermediate inter-grain coupling maximizes pinning without destroying vortex coherence. I believe the coercivity map itself.\n\nWhere the paper goes soft is in two places. First, the abstract attributes the peak to grain-boundary pinning, saying weak inter-grain exchange allows the core to be pinned by random anisotropy axes. The body text says the opposite: the peak height is the same for k=1 and k=0.25, so random grain anisotropy drives the peak, and the inter-grain exchange only shifts its location. That inconsistency should be fixed.\n\nSecond, the hyperthermia sweet-spot language overreaches. The paper equates heating performance with zero-temperature DC coercivity and never computes a hysteresis loop area, SAR, or minor-loop response at clinical field amplitudes and frequencies. The two reversal modes have different loop shapes, and if the applied AC field is below the coercive field, the enclosed area can be small or vanish. The DC coercivity maximum is a plausible heuristic for where maximum losses might occur, but it is not established. I would soften the practical claims and present this as a coercivity predictor, not a heating predictor.\n\nThe d^-1.85 power law between 50 and 75 nm is fit to very few points; treat it as illustrative rather than quantitative.\n\nOverall, the central reversal-mechanism result holds up and moves beyond macrospin treatments. A good referee can fix the abstract and the proxy claims. I would send it to peer review rather than desk reject, and I would tell the authors to either compute AC loop areas or explicitly limit their claims to coercivity. The intended audience is the magnetic hyperthermia simulation community, and they will cite the size map regardless.","headline":"A genuinely useful simulation study of vortex reversal in disordered nanoflowers, with a clear secondary coercivity peak and a clean reversal-mode criterion; the hyperthermia sweet-spot claims overreach because no AC loop area or SAR is computed.","tokens_in":15478,"tokens_out":2437,"would_cite":true,"duration_ms":27946,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["75.60.Jk","75.75.-m","75.50.Tt"],"model":"deepseek-v4-flash","headline":"Disordered iron-oxide nanoflowers reach a maximum coercivity at a diameter of about 110 nm, because grain anisotropy pins the vortex core until the flux-closure region takes over.","keywords":["magnetic nanoflowers","coercivity","vortex state","magnetic hyperthermia","micromagnetic simulation","spin disorder","grain-boundary pinning","iron-oxide nanoparticles"],"falsifier":"Measure the specific absorption rate (heating power per gram) of monodisperse iron-oxide nanoflowers with diameters from about 70 nm to 150 nm under a clinical alternating field of roughly 20–30 mT at 100–300 kHz: if the SAR maximum is not at or near the diameter where the zero-temperature coercivity peaks (about 110 nm for the calibrated microstructure), the coercivity-proxy assumption fails.","tokens_in":14421,"feed_emoji":"🧲","tokens_out":11018,"duration_ms":116269,"temperature":0.7,"pith_summary":"Iron-oxide nanoflowers are promising hyperthermia heaters because their interior is magnetically disordered, and this paper argues that the disorder is what creates a favorable size window. Simulating particles from 10 to 400 nm in diameter, it maps the coercive field—the reverse field needed to switch the magnetization—and finds that above the single-domain limit the magnetization folds into a vortex. In the vortex regime the coercivity is non-monotonic: it rises to a secondary maximum at about 110 nm for a grain-boundary coupling of $k = 0.25$, then falls as the particle grows. The peak marks the size at which the reversal mechanism changes from an immediate flip of the vortex core to a gradual rotation led by the surrounding flux-closure moments. If this is right, nanoflower synthesis can target a concrete diameter to maximize magnetic losses, and the paper explains why disordered multicore particles outperform smooth spheres.","feed_headline":"Spin disorder creates a ~110-nm nanoflower heating sweet spot","feed_subtitle":"Vortex reversal switches from core flip to flux rotation at the size where disordered grains heat best.","key_machinery":"The load-bearing object is the disordered vortex state in a simulated nanoflower: a near-spherical cluster of grains with randomly oriented uniaxial easy axes, coupled across grain boundaries by an exchange stiffness reduced by a factor $k$. The reversal analysis tracks the polar angle of the net magnetization and the relative vortex-core volume $V_{\\mathrm{core}}/V_{\\mathrm{NF}}$. The key structural threshold is $V_{\\mathrm{core}}/V_{\\mathrm{NF}} \\approx 1/3$, the point at which flux-closure moments outnumber core moments and the reversal mode changes from immediate core reversal to a gradual perpendicular rotation. The competing energy terms are demagnetizing energy, which favors the vortex; grain anisotropy, which creates pinning potentials for the core; and inter-grain exchange, which keeps the vortex profile coherent—too little coupling lets grains switch independently and kills the peak.","core_discovery":"The central discovery is a mechanism-based coercivity–size map for realistic iron-oxide nanoflowers. Above the single-domain threshold the flower hosts a magnetic vortex: an inhomogeneous spin texture in which moments circulate around a small, strongly magnetized core. In the disordered flower the coercivity rises with size just after vortex nucleation because the randomly oriented easy axes of the grains pin the vortex core, then reaches a maximum of about 14 mT at $d \\approx 110$ nm for inter-grain exchange factor $k = 0.25$, and decays as $d^{-1.5}$ as the core shrinks. The reversal switches at a vortex-core volume fraction $V_{\\mathrm{core}}/V_{\\mathrm{NF}} \\approx 1/3$: below that threshold the core reverses immediately along the field, above it the surrounding flux-closure moments rotate the core perpendicular to the field before switching. The same threshold in a perfect sphere leads to vanishing coercivity, so the peak is a signature of grain anisotropy and grain-boundary pinning, not of geometry alone. The location of the peak moves to larger sizes when grains are fully exchange-coupled (about 150 nm for $k = 1$), showing that disorder tunes the sweet spot while the peak height comes from the grains themselves.","pith_inferences":["The paper's coercivity proxy could be tested directly by a finite-temperature alternating-field simulation that computes hysteresis-loop area, and by calorimetry on monodisperse nanoflowers near 70–150 nm; nothing in the text guarantees the loop area peaks at exactly the coercivity peak.","The $V_{\\mathrm{core}}/V_{\\mathrm{NF}} \\approx 1/3$ threshold may serve as a cheap screening rule for other soft-magnetic particle geometries, since it is derived from the vortex profile rather than from the specific grain model; this universality is not claimed by the paper.","Because the peak height is nearly independent of the exchange-reduction factor while the peak position is not, grain-size control offers an experimentally accessible dial for the sweet spot; this suggests a similar tuning strategy for other multicore particle architectures such as chains or clusters, a transfer the paper leaves implicit."],"forward_implications":["At the calibrated microstructure (15 nm grains, $k = 0.25$), nanoflowers near $d \\approx 110$ nm should display a coercivity maximum of roughly 14 mT, bracketed by lower values on both sides.","Tuning the grain size shifts the sweet spot: fully exchange-coupled grains place the peak at about 150 nm, so a larger number of weaker grain-boundary links behaves like a smaller number of stronger ones.","Defect-free nanospheres of the same material should show a monotonic coercivity decline beyond the single-domain limit, with no secondary peak, because they lack the random easy axes that pin the core.","Beyond the peak, the coercivity decays as a power law: approximately $d^{-1.5}$ at $k = 0.25$ and $d^{-1.3}$ at $k = 1$, consistent with earlier vortex-state scalings.","The switch between reversal modes is controlled by the vortex-core volume fraction $V_{\\mathrm{core}}/V_{\\mathrm{NF}} \\approx 1/3$, a geometric threshold that applies to both disordered and ideal spheres, while only the pinning landscape decides whether that threshold yields a coercivity peak."],"supporting_citations":[{"why":"the experimental nanoflower measurements that the $k = 0.25$ simulation parameters are calibrated against.","marker":"[19]"},{"why":"reports vortex states in nanoflowers, the experimental phenomenon the simulations are built to explain.","marker":"[44]"},{"why":"gives the ideal-sphere vortex nucleation threshold of $7.2\\ell_{\\mathrm{ex}}$ that is compared with the disordered threshold.","marker":"[29]"},{"why":"provides the theoretical context that intra-particle disorder can raise heating efficiency by an order of magnitude, motivating the disorder-centered analysis.","marker":"[22]"},{"why":"supplies the random-anisotropy factor 0.48 used to interpret the single-domain coercivity.","marker":"[61]"},{"why":"the micromagnetic solver used to compute all hysteresis loops and reversal trajectories.","marker":"[72]"},{"why":"the geometric grain-generation procedure used to build the irregular polycrystalline nanoflower shape.","marker":"[73]"},{"why":"the vortex-regime coercivity scaling law against which the $d^{-1.5}$ decay is compared.","marker":"[30]"}],"fun_headline_variants":["Heating sweet spot at 110-nm nanoflowers pinned by grain disorder","Vortex reversal switch at 110 nm tunes nanoflower heating","Disordered nanoflowers peak heating at 110 nm","Coercivity peak at 110 nm from vortex core pinning","Grain disorder sets nanoflower heating sweet spot at 110 nm"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that coercivity from static hysteresis is a valid proxy for heating under clinical alternating fields, so that the size maximizing coercivity also maximizes heat release, together with the modelling premise that a random-grained particle with one reduced inter-grain exchange factor represents a real nanoflower.","fun_headline_variants_meta":{"raw":{"variants":["Heating sweet spot at 110-nm nanoflowers pinned by grain disorder","Vortex reversal switch at 110 nm tunes nanoflower heating","Disordered nanoflowers peak heating at 110 nm","Coercivity peak at 110 nm from vortex core pinning","Grain disorder sets nanoflower heating sweet spot at 110 nm"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000646,"raw_usage":{"total_tokens":3070,"prompt_tokens":1150,"completion_tokens":1920,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":766,"completion_tokens_details":{"reasoning_tokens":1841}},"tokens_in":766,"tokens_out":1920,"duration_ms":15501,"temperature":1.0,"reasoning_tokens":1841,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T20:37:29.436356+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the specific absorption rate (heating power per gram) of monodisperse iron-oxide nanoflowers with diameters from about 70 nm to 150 nm under a clinical alternating field of roughly 20–30 mT at 100–300 kHz: if the SAR maximum is not at or near the diameter where the zero-temperature coercivity peaks (about 110 nm for the calibrated microstructure), the coercivity-proxy assumption fails.","supporting_citations":[{"cited_title":"Berns, U","cited_arxiv_id":null,"evidence_quote":"the experimental nanoflower measurements that the $k = 0.25$ simulation parameters are calibrated against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"reports vortex states in nanoflowers, the experimental phenomenon the simulations are built to explain."},{"cited_title":"Goiriena-Goikoetxea, A","cited_arxiv_id":null,"evidence_quote":"gives the ideal-sphere vortex nucleation threshold of $7.2\\ell_{\\mathrm{ex}}$ that is compared with the disordered threshold."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the theoretical context that intra-particle disorder can raise heating efficiency by an order of magnitude, motivating the disorder-centered analysis."},{"cited_title":"Gavilán, K","cited_arxiv_id":null,"evidence_quote":"supplies the random-anisotropy factor 0.48 used to interpret the single-domain coercivity."},{"cited_title":"Superparamagnetic Superparticles for Magnetic Hyperthermia Therapy: Overcoming the Particle Size Limit","cited_arxiv_id":"2411.17172","evidence_quote":"the micromagnetic solver used to compute all hysteresis loops and reversal trajectories."},{"cited_title":"Erokhin and D","cited_arxiv_id":null,"evidence_quote":"the geometric grain-generation procedure used to build the irregular polycrystalline nanoflower shape."},{"cited_title":"Goiriena-Goikoetxea, D","cited_arxiv_id":null,"evidence_quote":"the vortex-regime coercivity scaling law against which the $d^{-1.5}$ decay is compared."}],"review_version":1}