Pith. sign in

REVIEW 4 major objections 4 minor 7 references

Realizing Blume-Capel Degrees of Freedom with Toroidal Moments in a Ruby Artificial Spin Ice

T0 review · 4 major / 4 minor · reviewed 2026-08-05 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Magnetic loops in a Ruby lattice become Blume-Capel spins

desk verdict A genuinely new ASI geometry that maps toroidal plaquettes to Blume-Capel states, with solid experiments but a too-weak test of the mapping itself. read the letter →

arxiv 2509.01522 v1 pith:TRAN33FZ submitted 2025-09-01 cond-mat.mes-hall

Luca Berchialla , Gavin M. Macauley , Flavien Museur , Tianyue Wang , Armin Kleibert , Peter M. Derlet , Laura J. Heyderman This is my paper · ORCID
classification cond-mat.mes-hall
keywords artificialspiniceBlume-CapelmodeltoroidalmomentferrotoroidicorderRubylatticephasetransitionsnanomagnetarraythree-state
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper claims that the three possible magnetic states of a closely packed plaquette in a Ruby artificial spin ice—a clockwise head-to-tail loop, an anticlockwise loop, or no closed loop—are exactly the three states of the Blume-Capel spin variable t = +1, −1, 0. In this mapping, toroidal moments of hexagons or triangles become the spins, with the two lattice parameters a and b setting the effective interaction strength J and the anisotropy Δ. Monte Carlo simulations of the full dipolar model, of the Blume-Capel model, and real-space magnetic imaging together show the same sequence: a high-temperature paramagnetic regime, a paratoroidic crossover where plaquette loops form without long-range order, and a second-order transition into a ferrotoroidic ground state. If correct, this is the first real-space realization of Blume-Capel degrees of freedom and suggests that artificial spin ice superstructures can implement other multi-state spin Hamiltonians by design.

What carries the argument

The load-bearing object is the plaquette toroidal moment: for each hexagon or triangle, the sense of circulation of its macrospins (clockwise, anticlockwise, or absent). These toroidal moments sit on a triangular lattice (hexagonal plaquettes) or a hexagonal lattice (triangular plaquettes) and act as the Blume-Capel variables. The mapping carries the argument by replacing the many-spin dipolar Hamiltonian with the two-parameter Blume-Capel Hamiltonian; J is set by a two-configuration energy difference (ground state versus one reversed toroidal moment), and Δ by the single-plaquette energy gap. The model then predicts the two-step ordering sequence, including the paratoroidic crossover and th

What would settle it

Measure the heat capacity or order-parameter distribution of a large Ruby ASI through the paratoroidic-to-ferrotoroidic transition: a discontinuous jump in the toroidal order parameter, or a heat-capacity peak that grows faster than log L with system size, would falsify the second-order Blume-Capel branch.

Watch

Extended reading notes

Core claim

For a plaquette of nanomagnets, the macroscopically distinct configurations are: all moments head-to-tail clockwise, all head-to-tail anticlockwise, or not all head-to-tail. These map onto t = +1, t = −1, and t = 0 of the Blume-Capel Hamiltonian H = −J Σ t_i t_j − Δ Σ t_i². J is fixed by the energy difference between two adjacent plaquettes with same-sign toroidal moments and one with a reversed moment; Δ is fixed by the gap between the average excited and ground plaquette energies. With J(a,b) and Δ(a,b) extracted from micromagnetic calculations, the Blume-Capel model reproduces the crossover and transition temperatures of the full dipolar Ruby ASI, and the experiments access a previously u

Load-bearing premise

The load-bearing premise is that each plaquette's micromagnetic configuration is reducible to one three-state variable with two fixed parameters J and Δ taken from local energy differences, and that as-grown arrays all sit at one effective annealing temperature; if multiplet correlations renormalize J and Δ with temperature, the mapping fails.

Editorial extensions

If this is right

  • A Ruby ASI with tunable a/b is a real-space laboratory for Blume-Capel physics: imaging the signs of plaquette toroidal moments directly measures the Blume-Capel order parameter Ψ.
  • For large Δ/J the ferrotoroidic transition temperature is exactly known from the Ising model on triangular and hexagonal lattices (4J/ln 3 and 2J/ln(2+√3)), giving a quantitative check of the mapping.
  • The same superstructure principle—using the collective state of a group of nanomagnets as a single spin variable—can be applied to other multi-state Hamiltonians.
  • The positive-Δ region explored here excludes the tricritical point; reaching it would require superstructures with antiferromagnetic toroidal interactions.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A direct test would be to measure the heat capacity of large Ruby ASI arrays under controlled thermal cycling: if the paratoroidic-to-ferrotoroidic peak sharpens with system size as log L, the second-order classification is confirmed; if it develops a jump, the Blume-Capel mapping needs a first-order branch or additional couplings.
  • The coarse-graining to J and Δ implicitly assumes intra-plaquette correlations can be integrated out once and for all; a natural extension is to compute temperature-dependent renormalized J(Δ) from the full dipolar model and compare with the fixed-parameter Blume-Capel predictions.
  • The manuscript cites a reference 28 for the deposition-annealing effective-temperature argument, but reference 28 does not appear in the bibliography; the effective-temperature comparison should be checked once the source is restored.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 4 minor

Summary. The manuscript reports the fabrication and characterization of a Ruby artificial spin ice (ASI) in which the lattice parameters a and b are independently tuned. It shows experimentally (MFM and PEEM) and via Monte Carlo simulations that the system can reach a ferrotoroidic ground state either directly or through an intermediate paratoroidic regime, depending on the ratio of interaction strengths within triangular and hexagonal plaquettes. The central claim is that fully-formed toroidal moments on closely-packed plaquettes realize the three states of the Blume-Capel (BC) spin variable t∈{-1,0,+1}, with effective parameters J and Δ extracted from micromagnetic energy differences, and that the observed two-step ordering and phase diagram are captured by the BC model. The paper reports a comparison of the order parameter Ψ between dipolar simulations of the Ruby ASI and BC simulations, including a 'modified' BC model with adjusted high-temperature state probabilities.

Significance. If the proposed mapping is valid, this would be a notable advance: a real-space, directly imageable realization of Blume-Capel degrees of freedom and a tunable toroidal-moment system with a Blume-Capel phase diagram. The experimental dataset is substantial, combining MFM and PEEM imaging over a wide range of lattice parameters, and the finite-size scaling analysis in Supplementary Figure 3 is a useful check of the critical behavior. The micromagnetic extraction of J and Δ from Mumax3 is clearly described and is a non-trivial step toward quantitative comparison. However, the validation of the central mapping is the weakest link, and the paper’s own modified BC construction partially fits the high-temperature behavior rather than predicting it.

major comments (4)
  1. [§'Blume-Capel degrees of freedom', Fig. 5b, Eq. (3)] The validation of the BC mapping rests on comparison of Ψ, which depends only on the single-site populations n_{+1}, n_{-1} and their imbalance. It does not test the interaction term -JΣt_i t_j or the Boltzmann structure of the coarse-grained variables. Agreement between the upper and lower panels of Fig. 5b therefore does not establish that neighboring toroidal moments interact through the two-parameter BC Hamiltonian. A direct test is needed, e.g. conditional probabilities P(t_i | t_j) or spin-spin correlators from the full dipolar Monte Carlo, compared with BC simulations using the extracted J and Δ without adjusting state probabilities.
  2. [Methods, 'Monte Carlo simulations' (modified Blume-Capel model)] The modified BC model sets the high-temperature state probabilities to 62/64, 1/64 and 1/64 for hexagonal plaquettes (and 6/8, 1/8, 1/8 for triangular plaquettes). Therefore the high-temperature plateau Ψ≈2/64 or Ψ≈2/8 is reproduced by construction rather than derived from the Hamiltonian. This makes the agreement in the crossover region partly circular. The authors should either use the unmodified BC model and quantify the mismatch, or subtract the fitted baseline and show that deviations from the fit are captured by the model. As written, the best-agreement claim is weakened by this fitting.
  3. [Methods, 'Monte Carlo simulations' (proposal distribution)] The modified BC simulation uses a non-symmetric proposal distribution: from state 0 the proposed states are drawn with probabilities 62/64, 1/64, 1/64, while from ±1 they are drawn with 1/3 each. No Metropolis-Hastings proposal-ratio correction is given. If the acceptance probability is the standard min(1, exp(-ΔE/T)), this non-symmetric proposal changes the equilibrium measure unless compensated. The authors need to specify the detailed-balance condition or state explicitly that the proposal asymmetry is corrected. Without this, the reported BC results may not correspond to the intended Hamiltonian.
  4. [Methods, 'Effective temperature'] The effective temperature is obtained by matching the experimental vertex and flux-closed plaquette populations to the Monte Carlo simulation of the same model that is being validated. This fitting can mask systematic errors in the Hamiltonian or in the mapping. To support the claim that the experiments traverse a constant-T_eff line, an independent estimate of T_eff (e.g., from blocking temperature physics) or a consistency check across many lattice-parameter ratios should be provided. As it stands, the quantitative comparison along the dashed line in Fig. 3a is not fully predictive.
minor comments (4)
  1. [Methods, 'Monte Carlo simulations'] Typo: 'having a positive fully-formed toroidal moment (state +1) is 1/64 and of having a positive fully-formed toroidal moment (state -1)' should read 'negative' for the -1 state.
  2. [References] Reference 28 is cited in the text ('effective thermal annealing ... 28,29') but is missing from the reference list; the numbering jumps from 27 to 29. Please add the missing reference or renumber.
  3. [Methods, 'Micromagnetic simulations'] Minor typos: 'and and the cell sizes' should be 'and the cell sizes'; also 'the probabilities of the proposed state are 1/3 for each of 0, +1 and -1' is clear but could be rephrased for readability.
  4. [Fig. 5b and Supplementary Figs. S7/S8] The color maps of Ψ are useful but the color scale and contour levels are not defined in the caption. Adding a color bar and explicit contour lines for Ψ=2/3, 1, 2 would make the comparison quantitative.

Circularity Check

2 steps flagged · score 5.0 of 10

High-temperature plateau in the modified Blume-Capel model is set by fitted state probabilities, and the effective temperature is calibrated to the same populations used for validation.

  1. fitted input called prediction [Methods, 'Monte Carlo simulations'; Section 'Blume-Capel degrees of freedom for toroidal moments and phase transitions']
    "In the modified Blume-Capel model, we adjust the probabilities of each state (+1, -1 and 0) so that they reflect the probability of obtaining the corresponding state in the Ruby ASI. In the Ruby ASI, considering hexagonal plaquettes on the triangular lattice, the probability of having a non-fully-formed toroidal moment (state 0) is 62/64, of having a positive fully-formed toroidal moment (state +1) is 1/64 and of having a positive fully-formed toroidal moment (state - 1) is 1/64."

    The modified Blume-Capel model fixes the single-site probabilities to the experimentally measured high-temperature toroidal populations. The order parameter Ψ = (n_+1 + n_-1) + |n_+1 - n_-1| is a function only of those single-site populations, so with p_0=62/64, p_+1=p_-1=1/64, the high-temperature value is Ψ = 2/64 by arithmetic, not by the Blume-Capel Hamiltonian. The paper then presents the modified model as providing the 'best agreement' with the Ruby ASI at high temperature; that agreement is imposed by construction.

  2. fitted input called prediction [Methods, 'Effective temperature']
    "The effective temperature is determined by identifying the Monte Carlo simulation temperature where the vertex and flux-closed plaquette populations observed in the experimental data best match those in the simulation. Specifically, for each temperature, an error function is computed as the weighted sum of the differences between the experimental and Monte Carlo populations."

    The effective temperature is fitted to the experimental vertex and flux-closed plaquette populations. Figure 4 then compares the same type of populations (t_H±, t_T±, Φ) against Monte Carlo curves at the fitted temperature k_BT/D≈0.469, and uses that comparison to support the traversal of the phase diagram. The agreement between the dashed Monte Carlo curve and the experimental points is therefore partly a restatement of the fitting criterion rather than an independent prediction. This is a calibration step, so it is only partially circular.

full rationale

Most of the central derivation is not circular. J and Δ are extracted from micromagnetic two-configuration energy differences (Figure 5a, Methods), not from the target Blume-Capel phase diagram, and the unmodified Blume-Capel simulations with those parameters independently produce the paratoroidic-to-ferrotoroidic crossover/transition structure; the exact Ising critical temperatures are cited from Baxter. No load-bearing self-citation or imported-uniqueness argument is used. The circularity is confined to two fitted inputs presented as validation. First, the 'modified Blume-Capel model' fixes single-site state probabilities to the measured high-temperature toroidal populations; because Ψ is a function of exactly those populations, the high-temperature plateau Ψ ≈ 2/64 is imposed by construction. Second, the effective temperature used to overlay experiment and simulation is calibrated to the same vertex/plaquette populations that Figure 4 then plots, so that agreement is partially a fit-quality check. These are partial circularities; they do not make the entire derivation equivalent to its inputs, because the phase-transition temperatures and the Δ/J dependence of the paratoroidic and ferrotoroidic regions remain emergent model predictions.

Assumptions & free parameters 2 free parameters · 4 assumptions · 0 invented entities

No new physical particles, forces, or conserved quantities are introduced. The Blume-Capel variable is an effective mapping of existing toroidal moments, and the toroidal moments themselves were established in prior ASI work (refs 6-8).

free parameters (2)
  • Effective temperature T_eff = Approximately kBT/D = 0.469 for as-grown samples
    Determined per lattice parameter by minimizing the weighted difference between experimental and Monte Carlo vertex and flux-closed plaquette populations (Methods 'Effective temperature'). Used to place as-grown states on the phase diagram.
  • Modified Blume-Capel state probabilities = Hexagons: 62/64 (0), 1/64 (+1), 1/64 (-1); Triangles: 6/8 (0), 1/8 (+1), 1/8 (-1)
    Set to the high-temperature populations of toroidal moments in the Ruby ASI Monte Carlo simulation (Methods 'Monte Carlo simulations'). This modification is what makes the modified Blume-Capel model match the Ruby ASI high-temperature behavior, i.e., fitting to the data being compared.
assumptions (4)
  • domain assumption Each nanomagnet is a single-domain Ising macrospin with dipole-dominated interactions
    Standard artificial spin ice assumption, used throughout and in the Monte Carlo simulations (Methods).
  • domain assumption As-grown configurations approximate thermal equilibrium at an effective temperature set by deposition annealing
    Invoked to compare as-grown MFM images with Monte Carlo snapshots (Intro and Methods 'Effective temperature'), with citations to refs 28-30.
  • ad hoc to paper A plaquette with a fully formed toroidal moment is a valid three-state degree of freedom interacting through a nearest-neighbor Blume-Capel Hamiltonian
    Central mapping in Section 'Blume-Capel degrees of freedom for toroidal moments and phase transitions'; J and Delta are defined via local two-configuration energy differences (Figure 5a). This is the load-bearing coarse-graining assumption.
  • domain assumption The paratoroidic to ferrotoroidic transition is in the 2D Ising universality class with nu = 1
    Supported by finite-size scaling analysis in Supplementary Figure 3 for three parameter sets; used to assign critical temperatures via exact triangular and hexagonal lattice Ising results (Baxter, ref 31).

how reviews work

0 comments
Cite this review

Pith. "Pith review of Realizing Blume-Capel Degrees of Freedom with Toroidal Moments in a Ruby Artificial Spin Ice." pith.science (2026). https://pith.science/paper/TRAN33FZ

@misc{pith2026250901522,
  author       = {Pith},
  title        = {Pith review of: Realizing Blume-Capel Degrees of Freedom with Toroidal Moments in a Ruby Artificial Spin Ice},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/TRAN33FZ}},
  note         = {Machine review of arXiv:2509.01522}
}
read the original abstract

Realizing exotic Hamiltonians beyond the Ising model is a key pursuit in experimental statistical physics. One such example is the Blume-Capel model, a three-state spin model, whose phase diagram features a tricritical point where second-order and first-order transition lines converge, leading to a coexistence of paramagnetic, ferromagnetic, and disordered phases. Here, we realize an artificial crystal of single-domain nanomagnets, placed on the links of the Ruby lattice, enabling real-space observation of the Blume-Capel degrees of freedom. These Blume-Capel degrees of freedom are represented by the presence, sign and interactions of the toroidal moments that emerge naturally in plaquettes of nanomagnets in the Ruby artificial spin ice. By precisely tuning the lattice parameters of the Ruby artificial spin ice, we demonstrate control over the two-step ordering process of the toroidal moments, whereby there is a high-temperature crossover from a paramagnetic phase to an intermediate paratoroidic regime, followed by a second-order phase transition to a ferrotoroidic ground state. This sequence of toroidal phases and transitions is accurately captured by the Blume-Capel framework and provides a direct realization of a substantial portion of the phase diagram associated with the model. This establishes a new platform for exploring exotic Hamiltonians in terms of artificial spin ice superstructures, here with groups of nanomagnets forming toroidal moments. The success of this mapping paves the way for an entirely new frontier in artificial spin ice: intentionally engineering lattice designs whose effective Hamiltonians mediate unconventional forms of magnetic order, with new behaviors and functionalities.

Figures

Figures reproduced from arXiv: 2509.01522 by the authors.

Figure 4
Figure 4. Signatures of the two-step ordering process and toroidic phases in as-grown configurations determined with MFM. (a) Average effective temperature of as-grown states as a function of the lattice parameter ratio. (b)-(c) Average fraction of formed toroidal moments associated with (b) hexagons 𝑡𝑡𝐻𝐻± = 𝑡𝑡𝐻𝐻+ + 𝑡𝑡𝐻𝐻− and (c) triangles 𝑡𝑡𝑇𝑇± = 𝑡𝑡𝑇𝑇+ + 𝑡𝑡𝑇𝑇− in as-grown states as a function of the lattice parameter ratio. … view at source ↗
Figure 5
Figure 5. Assigning Blume-Capel degrees of freedom to the Ruby ASI. (a) The Hamiltonian of the Blume-Capel model is composed of two terms: a nearest-neighbour interaction of strength J, and an anisotropy term ∆, which is the energy difference between the 0 and the ±1 states. In our effective model for the Ruby ASI, the nearest-neighbour interaction J is given by the energy difference between a ferrotoroidic ground state and a… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

7 extracted references · 3 canonical work pages

  1. [1]

    H.; Marrows, C

    (1) Skjærvø, S. H.; Marrows, C. H.; Stamps, R. L.; Heyderman, L. J. Advances in Artificial Spin Ice. Nat. Rev. Phys. 2020, 2 (1), 13–28. https://doi.org/10.1038/s42254-019-0118-3. (2) Berchialla, L.; Macauley, G. M.; Heyderman, L. J. Focus on Three-Dimensional Artificial Spin Ice. Appl. Phys. Lett. 2024, 125 (22), 220501. https://doi.org/10.1063/5.0229120...

    doi:10.1038/s42254-019-0118-3 2020

  2. [5]

    Snakes in the Plane: Controllable Gliders in a Nanomagnetic Metamaterial

    https://doi.org/10.48550/arXiv.2505.01116. (35) Imre, A.; Csaba, G.; Ji, L.; Orlov, A.; Bernstein, G. H.; Porod, W. Majority Logic Gate for Magnetic Quantum -Dot Cellular Automata. Science 2006, 311 (5758), 205–208. https://doi.org/10.1126/science.1120506. (36) Arava, H.; Derlet, P. M.; Vijayakumar, J.; Cui, J.; Bingham, N. S.; Kleibert, A.; Heyderman, L....

    work page Pith review arXiv doi:10.48550/arxiv.2505.01116 2006

  3. [964]

    (39) Gartside, J

    https://doi.org/10.1038/s41467-024-45319-7. (39) Gartside, J. C.; Stenning, K. D.; Vanstone, A.; Holder, H. H.; Arroo, D. M.; Dion, T.; Caravelli, F.; Kurebayashi, H.; Branford, W. R. Reconfigurable Training and Reservoir Computing in an Artificial Spin -Vortex Ice via Spin- Wave Fingerprinting. Nat. Nanotechnol. 2022, 17 (5), 460–469. https://doi.org/10....

    doi:10.1038/s41467-024-45319-7 2022

  4. [1982]

    (32) Graner, F.; Glazier, J. A. Simulation of Biological Cell Sorting Using a Two-Dimensional Extended Potts Model. Phys. Rev. Lett. 1992, 69 (13), 2013–2016. https://doi.org/10.1103/PhysRevLett.69.2013. (33) Reichardt, J.; Bornholdt, S. Detecting Fuzzy Community Structures in Complex Networks with a Potts Model. Phys. Rev. Lett. 2004, 93 (21), 218701. ht...

    doi:10.1103/physrevlett.69.2013 1992

  5. [1987]

    P.; Akerman, J.; Stein, A.; Phatak, C.; Evans, R

    (29) Morgan, J. P.; Akerman, J.; Stein, A.; Phatak, C.; Evans, R. M. L.; Langridge, S.; Marrows, C. H. Real and Effective Thermal Equilibrium in Artificial Square Spin Ices. Phys. Rev. B 2013, 87 (2), 024405. https://doi.org/10.1103/PhysRevB.87.024405. (30) Hofhuis, K.; Hrabec, A.; Arava, H.; Leo, N.; Huang, Y.-L.; Chopdekar, R. V.; Parchenko, S.; Kleiber...

    work page doi:10.1103/physrevb.87.024405 2013

  6. [2025]

    Magnetic ordering in out-of-plane artificial spin systems based on the Archimedean lattices

    https://doi.org/10.48550/arXiv.2503.08462. (27) Grünbaum, B.; Shephard, G. C. (Geoffrey C. Tilings and Patterns; New York : W.H. Freeman,

    work page Pith review arXiv doi:10.48550/arxiv.2503.08462

  7. [2562]

    (41) Kitaev, A

    https://doi.org/10.1038/s41467-023-38286-y. (41) Kitaev, A. Anyons in an Exactly Solved Model and Beyond. Annals of Physics 2006, 321 (1), 2–111. https://doi.org/10.1016/j.aop.2005.10.005. (42) Vansteenkiste, A.; Leliaert, J.; Dvornik, M.; Helsen, M.; Garcia -Sanchez, F.; Van Waeyenberge, B. The Design and Verification of MuMax3. AIP Advances 2014, 4 (10)...

    doi:10.1038/s41467-023-38286-y 2006

Pith tools

Reviewed August 5, 2026 · model on record in the stance chip above.