{"id":"65b7093d-8656-43c8-8c04-ad70d23077a6","arxiv_id":"2509.05424","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Fe10Dy10 is predicted to host a maximally toroidal ground doublet with a large toroidal moment and a finite-temperature response measurable via a shaped-laser protocol.","lead":"This paper predicts that the Fe10Dy10 molecular ring hosts a record-large toroidal magnetic moment in its ground state and proposes a laser-based protocol to detect it. A generalist might read it because it offers a concrete route to the first direct observation of a long-predicted molecular toroidal moment, which could enable new magnetoelectric control schemes.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Toroidal relaxation time from the paper's own microSQUID data (U=11.5 K, τ0=8×10^-13 s) is ~10^-11 s at 4 K, contradicting the 'slow toroidal relaxation' assumption central to the proposed detection protocol.","rationale":"The reader's conditional verdict identified the slow-relaxation assumption and the missing magnetoelectric tensor as weaknesses. I agree with those but focus on the relaxation timescale as the single most load-bearing issue because it is directly contradicted by the paper's own experimental data. The reported microSQUID barrier U=11.5 K and attempt time τ0=8×10^-13 s yield a relaxation time of ~10^-11 s at 4 K, while the proposed protocol requires the toroidal polarization to survive for at least the measurement time of the μSQUID (milliseconds). The AC susceptibility data similarly indicate fast relaxation. This invalidates the 'accumulation' step of the protocol. The missing magnetoelectric tensor further undermines the readout step, but the relaxation contradiction alone suffices to make the direct-observation claim unsupported. The theoretical model of the toroidal ground state and the finite-temperature toroidal susceptibility may still be valid, so the paper deserves conditional acceptance rather than outright rejection. My recommendation is therefore unchanged from the reader's conditional verdict, with a request for additional relaxation measurements or calculations and the missing magnetoelectric tensor derivation.","tokens_in":27419,"tokens_out":6639,"duration_ms":71039,"concrete_test":"Using the reported microSQUID parameters (U=11.5 K, τ0=8×10^-13 s from SI Note 5), evaluate τ(T)=τ0 exp(U/T) at the proposed operating temperature (e.g., 4.2 K) and compare it with the μSQUID integration time and the required laser repetition period. If τ is shorter than ~1 μs, the accumulation step fails. Additionally, require the authors to provide the full derivation of the magnetoelectric tensor and the predicted electric-field-induced magnetic moment for the proposed laser geometry; without this, the readout claim is untestable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that toroidal polarisation can be 'prepared, accumulated and read out under realistic experimental conditions' relies on the Section 5 assumption that toroidal relaxation is 'sufficiently slow.' The paper's own data contradict this. Supplementary Note 5 reports an Arrhenius fit U=11.5 K, τ0=8×10^-13 s from microSQUID magnetisation decay. At liquid-He temperature (T≈4 K), τ≈τ0 exp(U/T)≈1.4×10^-11 s—orders of magnitude shorter than any realistic μSQUID integration time (~100 μs or longer) and shorter than the interval between successive laser pulses needed for cumulative accumulation. AC susceptibility (SI Note 4) places the out-of-phase maximum above 1500 Hz even below 1.8 K, implying τ<10^-4 s at those temperatures and even faster at 4 K. Thus the proposed population imbalance would decay before it can be measured. Furthermore, the abstract promises an ab initio-informed magnetoelectric tensor for electric-field readout, but no such tensor or calculation appears in the main text or Supplementary Information. Without a quantitative readout model and a measured or computed toroidal relaxation time, the detection protocol is not established. The static splitting estimate (ΔE≈1.5 cm^-1) may be correct, but it alone does not demonstrate that polarisation can be accumulated and observed.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the icosanuclear ring Fe10Dy10 and constructs an ab initio-informed transfer-matrix model that projects the 62-billion-dimensional low-energy Hilbert space onto 2^10 Dy(III) Ising configurations dressed by Fe(III) spin excitations, with Fe–Fe exchange treated perturbatively. The model is compared with powder magnetization, χT, and specific-heat measurements and shows good agreement. On this basis the authors predict a maximally toroidal ground doublet with a toroidal moment of order 1500 μB Å, introduce a toroidal susceptibility ξ as a linear response to ∇×B, and propose a laser-pulse protocol in which a shaped optical waveform induces a curl of the magnetic field that removes the degeneracy of the counter-rotating toroidal states. The central claim is that Fe10Dy10 allows toroidal polarization to be prepared, accumulated, and read out under realistic experimental conditions.","tokens_in":27749,"tokens_out":4794,"duration_ms":54241,"significance":"If the result holds, the paper provides a significant theoretical framework for finite-temperature toroidal response in a large molecular ring, going beyond the archetypal Dy3 triangle. The transfer-matrix treatment of a 62-billion-dimensional space with ab initio parameters is an impressive technical achievement, and the reproduction of magnetization, susceptibility, and specific heat gives nontrivial support to the underlying spin Hamiltonian. The definition of ξ as a thermodynamic response to a magnetic-field curl is a useful conceptual step. However, the paper's headline claim—that toroidal polarization can be prepared, accumulated, and read out under realistic conditions—is not quantitatively established. The detection protocol rests on an unsupported assumption about toroidal relaxation times, and a promised ab initio-informed magnetoelectric tensor for electric-field readout is absent from the manuscript and supplementary information.","major_comments":[{"comment":"The proposed preparation-and-detection protocol requires that toroidal relaxation be 'sufficiently slow' (Section 5), but no toroidal relaxation time is computed or measured. The available dynamics data indicate the opposite: Supplementary Note 5 reports an Arrhenius fit U = 11.5 K and τ0 = 8×10^-13 s from microSQUID magnetization decay, giving τ ≈ 1.4×10^-11 s at 4 K, the temperature invoked for the protocol. Supplementary Note 4 places the AC out-of-phase maximum above 1500 Hz even below 1.8 K, implying τ < 10^-4 s at those temperatures and even shorter at 4 K. The population imbalance generated by a laser pulse would thus decay orders of magnitude faster than any realistic measurement or inter-pulse accumulation interval. The authors may argue that magnetization relaxation differs from toroidal relaxation, but then they must provide an explicit estimate or calculation of the toroidal","section":"Section 5, Eq. (10); Supplementary Notes 4 and 5"},{"comment":"The abstract promises that 'an ab initio-informed magnetoelectric tensor predicts an electric-field-induced magnetic moment within μSQUID detectability,' but no such tensor is derived, tabulated, or used anywhere in the main text or the Supplementary Information. The detection section (Section 5) discusses only the τ·(∇×B) coupling and gives a static splitting estimate; it does not model the electric-field readout, the induced magnetic moment, or the μSQUID response. Without a quantitative readout model, the claim that the toroidal polarization is readable within μSQUID detectability is not established.","section":"Abstract; Section 5"},{"comment":"The finite-temperature toroidal response is characterized by the equilibrium linear-response function ξ, but the proposed laser protocol is a strongly time-dependent, far-from-equilibrium process involving repeated 10 fs pulses. The paper does not connect the equilibrium ξ or the static splitting ΔE = 2τ·(∇×B) to the accumulation dynamics under the full optical waveform. In particular, no account is given of the optical electric field acting directly on the molecular charges, sample heating, or the repetition-rate constraint imposed by the relaxation time. The static splitting estimate may be correct, but it does not by itself demonstrate that a measurable toroidal polarization can be accumulated.","section":"Section 4, Eq. (6)-(7); Section 5"}],"minor_comments":[{"comment":"The manuscript title in the main text is 'Giant Molecular Toroidal Moment Amenable to Direct Observation in a Fe10Dy10 Ring,' whereas the submitted arXiv title is 'Finite-Temperature Toroidal Moment Amenable to Direct Observation in an Fe10Dy10 Molecular Ring.' Please harmonize the title.","section":"Title"},{"comment":"The text says 'Using Eq. (17, 18)' before those equations are introduced in Methods. Please renumber equations or adjust the cross-reference.","section":"Section 2.4"},{"comment":"The model is repeatedly described as 'parameter free,' but Eq. (5) contains two fitted parameters a and b for the lattice specific heat. This is not a problem for the magnetic model, but the wording should be qualified to avoid confusion.","section":"Section 3, Eq. (5)"},{"comment":"Typo: 'cohercive field' should be 'coercive field.'","section":"Supplementary Note 5"}],"recommendation":"major_revision","confidential_remarks":"The transfer-matrix modeling and the comparison with thermodynamic data are solid and could justify a strong paper about the predicted toroidal ground state and finite-temperature toroidal susceptibility. The main risk is overclaiming: the abstract and Section 5 promise a complete preparation-and-readout protocol, but the dynamics and readout parts are not worked out. If the authors can either supply a quantitative toroidal relaxation estimate and a magnetoelectric-tensor calculation, or substantially soften the 'direct observation under realistic conditions' claim, the paper would be defensible. I would not recommend rejection, because the underlying static/thermodynamic results appear sound and valuable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The Fe10Dy10 paper is worth reading for the model, not for the detection claim. Soncini et al. build a parameter-free transfer-matrix description of a 62-billion-dimensional space using ab initio crystal-field and broken-symmetry DFT exchange couplings, and they reproduce powder magnetisation, chiT, and heat capacity across a wide temperature range. The introduction of the toroidal susceptibility xi as a finite-temperature linear response to curl B is a genuine new formalism, and the prediction of a ground doublet with tau ~ 1500 muB angstrom is striking if the model holds. That part is careful and honest; the SI even notes that the microSQUID data show no direct signature of the toroidal ground state, and explains why.\n\nThe soft spots are all on the direct-observation protocol. First, their own relaxation data contradict the assumption of sufficiently slow toroidal relaxation. AC susceptibility shows out-of-phase maxima above 1500 Hz below 1.8 K, and microSQUID relaxation gives U = 11.5 K, tau0 = 8e-13 s, implying tau ~ 1e-11 s at 4 K. That is orders of magnitude faster than any cumulative laser-pulse scheme or microSQUID readout. Unless toroidal relaxation is somehow decoupled from magnetic relaxation (not shown), the proposed population imbalance will decay before it can be measured. Second, the abstract promises an ab initio-informed magnetoelectric tensor for electric-field readout, but no such tensor appears in the main text or SI. Third, the estimated laser-induced splitting has a numerical slip: 2 x 1500 muB A x 5e-4 T/A = 1.5 muB T ~ 0.7 cm-1, not 1.5 cm-1. That does not break the qualitative point but suggests the quantitative case is not polished.\n\nNet: the modelling is a real contribution and the toroidal susceptibility concept should survive peer review, but the central claim of amenable to direct observation is currently a hope, not a demonstration. This deserves a serious referee—the formalism and the prediction are important enough—but the authors should be pushed to compute or measure toroidal relaxation, provide the missing tensor, fix the arithmetic, and soften the title and abstract.\n\nRecommendation: send to peer review, but expect major revision.","headline":"Solid ab initio model and a new toroidal susceptibility, but the paper's own relaxation data undermine the central direct-observation claim.","tokens_in":28250,"tokens_out":4061,"would_cite":true,"duration_ms":39991,"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":"The paper predicts that Fe10Dy10, a 20-ion 3d–4f molecular ring, hosts a maximally toroidal ground doublet with a toroidal moment of order 1500 μB Å that a shaped near-infrared laser pulse can split and make directly detectable.","keywords":["toroidal moment","single-molecule toroics","molecular magnetism","Fe10Dy10","toroidal susceptibility","magnetoelectric coupling","transfer-matrix model","lanthanide wheel"],"falsifier":"A single-crystal micro-SQUID magnetization measurement below 1 K should show the predicted nonmagnetic ground state and sharp level-crossing steps at fields of roughly 6–40 mT, depending on orientation, where the toroidal ground doublet crosses the magnetic excited state. If instead the zero-field ground state is magnetic at these temperatures, the claimed maximally toroidal ground doublet is falsified; similarly, the proposed 10-fs pulse experiment would falsify the detection scheme if no electric-field-induced magnetization appears at the estimated pulse parameters.","tokens_in":27337,"feed_emoji":"🧲","tokens_out":14105,"duration_ms":123336,"temperature":0.7,"pith_summary":"Single-molecule toroics are magnetic wheels whose spin arrangement forms a closed vortex, carrying a toroidal moment. Because opposite vortex chiralities are degenerate under ordinary magnetic fields, toroidal order has so far been inferred indirectly rather than seen directly. This paper predicts that the icosanuclear ring Fe10Dy10 changes that: an ab initio-parameterised transfer-matrix model, which reproduces measured magnetization, susceptibility, and specific heat, places a maximally toroidal ground doublet with |τ| ≈ 1500 μB Å—about 25 times larger than the Dy3 archetype. To make the prediction testable, the paper introduces the toroidal susceptibility ξ, a linear-response function giving the toroidal polarization induced by a magnetic-field curl, and shows the response remains substantial up to about 10 K. It then estimates that a focused 10-fs near-infrared pulse creates a curl splitting of ≈1.5 cm−1 between the two chiralities and predicts a magnetoelectric signal readable by micro-SQUID magnetometry, proposing a concrete path to the first direct observation of a molecular toroidal moment.","feed_headline":"Molecular ring's magnetic vortex now laser-detectable","feed_subtitle":"A 10-fs laser pulse should split the two vortex states by 1.5 cm−1, making molecular toroidal order measurable.","key_machinery":"The central object is the toroidal moment operator τ = Σ_i r_i × M_i, evaluated in the low-energy manifold of the ring. The argument is carried by three tools. First, an ab initio-parameterised transfer-matrix model: the 2^10 × 6^10 ≈ 62-billion-dimensional Hilbert space of ten DyIII Ising doublets and ten FeIII S = 5/2 spins is compressed into a product of ten 24×24 transfer matrices, with FeIII–FeIII exchange treated in first-order perturbation theory. This model reproduces powder magnetization, χT, and field-dependent specific heat without fitting parameters (aside from two lattice terms). Second, the newly introduced toroidal susceptibility ξαβ = −∂²F/∂(∇×B)α∂(∇×B)β, the curl-field analo","core_discovery":"The paper's central claim is that Fe10Dy10 realizes a maximally toroidal ground state. In the model, the ten DyIII Ising spins form a zero-noded s-wave vortex around the elliptical wheel; the ground Kramers doublet consists of the two time-reversed chiralities of this vortex and carries essentially no magnetic moment but a toroidal moment τ = Σ_i r_i × M_i of order 1500 μB Å. A magnetic state of ~83 μB lies only 0.2 cm−1 higher, and most of the dense low-energy spectrum carries both toroidal and magnetic character. The authors define the molar toroidal susceptibility ξ as the second derivative of the free energy with respect to ∇×B, so that T→0 gives the square of the ground-state toroidal m","pith_inferences":["Editorial inference: because the toroidal moment grows linearly with ring radius and per-ion moments, applying the same transfer-matrix-plus-ξ machinery to even larger 3d–4f wheels should push the laser-induced splitting well beyond 1.5 cm−1 and make direct detection easier.","Editorial inference: the paper estimates the response to a single pulse but does not simulate the full optical waveform; a natural two-pulse pump-probe variant would use the first pulse to create the curl-induced imbalance and a delayed second pulse to read it out via the magnetoelectric tensor, with pulse shape and timing matched to the slow toroidal relaxation suggested by the AC-susceptibility ","Editorial inference: the proximity of the magnetic state at 0.2 cm−1 suggests chemical control—ligand substitution or strain that tunes the Fe–Dy exchange constants—could push the magnetic excitation higher, purifying the toroidal ground state and extending its finite-temperature window.","Editorial inference: if the predicted magnetoelectric readout works, the electric-dipole symmetry of τ implies a route to electrically write toroidal chirality, since the two chiralities are time-reversed partners; the paper motivates but does not demonstrate such switching."],"forward_implications":["Fe10Dy10 becomes a candidate platform for the first direct observation of a molecular toroidal moment; a shaped near-infrared pulse should split the ground doublet by ~1.5 cm−1, and the resulting population imbalance should be readable in micro-SQUID magnetometry.","The toroidal susceptibility ξ establishes a standard thermodynamic response function for single-molecule toroics, letting future work compute and compare finite-temperature toroidal polarization across molecules and field configurations.","The parameter-free transfer-matrix strategy makes the full 62-billion-state low-energy spectrum of large 3d–4f wheels computationally accessible, so the same machinery can be reused for related heterometallic rings.","The prediction that a uniform field stabilizes a constant toroidal polarization up to ~10 K implies toroidal states can be addressed at experimentally accessible temperatures, not only in the millikelvin limit.","A magnetic excited state only 0.2 cm−1 above the toroidal ground doublet means that even a ~5 mT field can make the ground state magnetic; this competition is captured by the model and is relevant for interpreting low-field magnetization."],"supporting_citations":[{"why":"Defines toroidal moments in condensed matter and their coupling to magnetic-field curls, the basis for the H_Tor interaction and for ξ.","marker":"[3]"},{"why":"Derives molecular response to non-uniform magnetic fields, grounding the τ·(∇×B) splitting formula used in the detection estimate.","marker":"[8]"},{"why":"Establishes the nonmagnetic Kramers doublet of the Dy3 triangle as a toroidal state, the benchmark the paper compares against.","marker":"[16]"},{"why":"Reports a net toroidal moment in a {Dy6} ring, the previous molecular reference that Fe10Dy10 is argued to surpass.","marker":"[17]"},{"why":"Demonstrates a ferrotoroidic ground state in a heterometallic CrIII–DyIII6 wheel, motivating the 3d–4f design.","marker":"[24]"},{"why":"Reports the synthesis and family of Fe10Ln10 nano-toruses, the source of the Fe10Dy10 compound studied here.","marker":"[29]"},{"why":"Supplies the ab initio multiconfigurational code used to compute the Dy(III) crystal-field doublets and magnetic axes.","marker":"[30]"},{"why":"Provides the broken-symmetry formula used to extract the Fe–Dy and Fe–Fe exchange couplings from DFT.","marker":"[32]"},{"why":"Shows how to express response properties as analytical derivatives of the Helmholtz free energy, the basis for the toroidal susceptibility ξ.","marker":"[34]"},{"why":"Introduces the transfer-matrix representation of decorated Ising rings used to compute the 62-billion-state partition function.","marker":"[43]"}],"fun_headline_variants":["Direct laser readout of molecular toroidal moment","Fe10Dy10 ring: toroidal polarization seen at finite T","Light pulse imprints and reads molecular vortex","Toroidal order in molecular ring made measurable","Laser protocol detects toroidal polarization in ring"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The detection estimate hinges on the assumption that the ring's toroidal moment couples to the magnetic-field curl of the shaped laser pulse with the linear form H = τ·(∇×B) and the estimated |∇×B| ≈ 5×10^-4 T/Å, and that toroidal relaxation is slow enough for the induced population imbalance to accumulate and be measured; the paper does not model the full optical waveform or heating.","fun_headline_variants_meta":{"raw":{"variants":["Direct laser readout of molecular toroidal moment","Fe10Dy10 ring: toroidal polarization seen at finite T","Light pulse imprints and reads molecular vortex","Toroidal order in molecular ring made measurable","Laser protocol detects toroidal polarization in ring"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000224,"raw_usage":{"total_tokens":1359,"prompt_tokens":866,"completion_tokens":493,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":610,"completion_tokens_details":{"reasoning_tokens":420}},"tokens_in":610,"tokens_out":493,"duration_ms":5335,"temperature":1.0,"reasoning_tokens":420,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T05:23:41.688118+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A single-crystal micro-SQUID magnetization measurement below 1 K should show the predicted nonmagnetic ground state and sharp level-crossing steps at fields of roughly 6–40 mT, depending on orientation, where the toroidal ground doublet crosses the magnetic excited state. If instead the zero-field ground state is magnetic at these temperatures, the claimed maximally toroidal ground doublet is falsified; similarly, the proposed 10-fs pulse experiment would falsify the detection scheme if no electric-field-induced magnetization appears at the estimated pulse parameters.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines toroidal moments in condensed matter and their coupling to magnetic-field curls, the basis for the H_Tor interaction and for ξ."},{"cited_title":", author Ligabue, A","cited_arxiv_id":null,"evidence_quote":"Derives molecular response to non-uniform magnetic fields, grounding the τ·(∇×B) splitting formula used in the detection estimate."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the nonmagnetic Kramers doublet of the Dy3 triangle as a toroidal state, the benchmark the paper compares against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports a net toroidal moment in a {Dy6} ring, the previous molecular reference that Fe10Dy10 is argued to surpass."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Demonstrates a ferrotoroidic ground state in a heterometallic CrIII–DyIII6 wheel, motivating the 3d–4f design."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the synthesis and family of Fe10Ln10 nano-toruses, the source of the Fe10Dy10 compound studied here."},{"cited_title":", author Piccardo, M","cited_arxiv_id":null,"evidence_quote":"Supplies the ab initio multiconfigurational code used to compute the Dy(III) crystal-field doublets and magnetic axes."},{"cited_title":", author Fukui, H","cited_arxiv_id":null,"evidence_quote":"Provides the broken-symmetry formula used to extract the Fe–Dy and Fe–Fe exchange couplings from DFT."},{"cited_title":"& author Soncini, A","cited_arxiv_id":null,"evidence_quote":"Shows how to express response properties as analytical derivatives of the Helmholtz free energy, the basis for the toroidal susceptibility ξ."},{"cited_title":"& author Chibotaru, L","cited_arxiv_id":null,"evidence_quote":"Introduces the transfer-matrix representation of decorated Ising rings used to compute the 62-billion-state partition function."}],"review_version":1}