Single-Molecule Toroics: Magnetic Vortices

• Single-molecule toroics (SMTs) are molecular clusters with toroidal magnetic states formed by head-to-tail spin arrangements that cancel net dipoles.
• They leverage strong single-ion anisotropy and competing inter-ion couplings, particularly in lanthanide-based and mixed 3d–4f systems, to stabilize counter-rotating vortex configurations.
• Recent advances combine polarized neutron diffraction, ab initio modeling, and microwave-control protocols to achieve and probe coherent toroidal-state manipulation.

Single-molecule toroics (SMTs) are molecular clusters whose low-energy magnetic states carry a toroidal, or anapole, moment generated by head-to-tail spin and orbital-current textures arranged around a closed loop, so that the net magnetic dipole vanishes while a magnetic vortex remains finite (Soncini et al., 5 Sep 2025). In lanthanide-based archetypes such as Dy$_3$, coupled double triangles such as MDy$_6$ ($M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$), and larger mixed $3d$–$4f$ wheels such as Fe$_{10}$Dy$_{10}$, SMT behavior arises from strong single-ion anisotropy combined with weak, often competing, inter-ion couplings that stabilize counter-rotating ground doublets and suppress uniform-field magnetization (Hymas et al., 11 Apr 2025). Recent work has extended the field from indirect assignment by magnetometry and ab initio modeling to atomic-scale visualization by polarized neutron diffraction and to explicit protocols for coherent toroidal-state preparation and field-curl manipulation (Jenkins et al., 8 Jun 2026).

1. Definition of the toroidal degree of freedom

For a discrete molecular cluster, the toroidal moment is defined from the spatial distribution of local magnetic moments. In the Dy$_3$ neutron-diffraction study it is written as

$$\mathbf{T}=\sum_i \mathbf{r}_i \times \boldsymbol{\mu}_i,$$

whereas the Fe$_{10}$Dy$_6$0 work uses the molecular toroidal-moment operator

$_6$1

The microwave-control study notes that standard SMT conventions sometimes include a factor $_6$2,

$_6$3

and attributes these differences to multipolar convention; in the Fe$_6$4Dy$_6$5 analysis, $_6$6 is the operative definition (Soncini et al., 5 Sep 2025).

The toroidal moment is odd under time reversal and odd under spatial inversion. Under time reversal, $_6$7, hence $_6$8; under inversion, $_6$9, again implying $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$0. The microwave study describes $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$1 as a polar vector of magnetic origin, with time-odd and space-odd character, while the neutron-diffraction study emphasizes that toroidal moments can stack ferro- or antiferrotoroidically in crystals and thereby enable magnetoelectric couplings and fourth-order ferroic phenomena (Hymas et al., 11 Apr 2025).

Physically, a toroidal state is a head-to-tail loop of magnetic dipoles, or equivalently a vortex of spin and orbital currents, producing no net magnetic poles and typically a vanishing dipole moment. Its natural conjugate field is not a uniform magnetic field but a magnetic field with nonzero curl. In the Fe$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$2Dy$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$3 formulation, the coupling is

$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$4

so the degeneracy of counter-rotating toroidal partners is lifted by

$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$5

Because $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$6 is simultaneously odd under $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$7 and $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$8, toroidal polarization can mediate linear magnetoelectric coupling; in achiral systems this requires breaking $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$9 with magnetic-field curls, whereas in chiral molecules broken $3d$0 can allow mixing so that uniform $3d$1 contributes as well (Soncini et al., 5 Sep 2025).

2. Microscopic origin in lanthanide molecular clusters

The archetypal SMT is the Dy$3d$2 triangle. Its three Dy$3d$3 ions carry strongly axial local moments whose easy axes are nearly tangential to the triangular rim, so the in-plane vector sum of dipoles vanishes while the vortex pattern generates a toroidal moment normal to the molecular plane. In the non-collinear Ising description used for Dy-based SMTs, pseudo-spins $3d$4 define a toroidal ground doublet $3d$5, consisting of time-reversal-related counter-rotating vortex textures with vanishing net magnetic moment (Hymas et al., 11 Apr 2025).

The trimeric complex $3d$6 provides a particularly explicit realization. Single-crystal neutron diffraction established trigonal $3d$7 symmetry with a $3d$8-OH bridge, intramolecular Dy–Dy distance $3d$9, and four molecules per unit cell. The Dy$4f$0 core has $4f$1 symmetry through its center, and the three local easy axes lie approximately tangential to the triangle with a slight out-of-plane tilt, producing a head-to-tail arrangement in the $4f$2 plane and a toroidal moment along the crystallographic $4f$3 axis. Within the $4f$4 lattice, two molecular toroidal domains rotate clockwise and two counterclockwise per unit cell, yielding antiferrotoroidic stacking along $4f$5 and zero net toroidization in the single crystal (Jenkins et al., 8 Jun 2026).

Ab initio calculations on this Dy$4f$6 system show strongly axial ground Kramers doublets at each Dy site with $4f$7–$4f$8, $4f$9–$_{10}$0, $_{10}$1–$_{10}$2, and near-pure $_{10}$3 character $_{10}$4. The principal magnetic axes have out-of-plane tilt $_{10}$5 and in-plane canting $_{10}$6–$_{10}$7, consistent with the experimentally inferred $_{10}$8 canting. The lowest eight Kramers doublets span $_{10}$9–$_{10}$0, with the first excited doublet at $_{10}$1, yielding a computed Orbach barrier $_{10}$2 (Jenkins et al., 8 Jun 2026).

In coupled double triangles MDy$_{10}$3, two Dy$_{10}$4 toroidal moieties are stacked and linked by a central ion. The Dy sites retain approximately tangential anisotropy axes with weak out-of-plane canting angles $_{10}$5 and $_{10}$6, and the low-energy configurations can be ferrotoroidic (FT, con-rotating) or antiferrotoroidic (AFT, counter-rotating) depending on geometry, exchange, and dipolar interactions. The microwave-control study emphasizes that weak inversion-symmetry breaking, expressed as $_{10}$7, becomes a crucial spectroscopic handle because perfect inversion symmetry prevents a static field from distinguishing inversion-related toroidal configurations (Hymas et al., 11 Apr 2025).

3. Scaling from Dy$_{10}$8 to Fe$_{10}$9Dy$_3$0

A central development in SMT research is the move from small triangles to larger rings in order to amplify the toroidal response. Since $_3$1, both larger molecular radii and larger local moments increase $_3$2. The Fe$_3$3Dy$_3$4 icosanuclear wheel alternates Fe and Dy ions around an elliptical ring of average radius $_3$5, compared with $_3$6 for Dy$_3$7. Using $_3$8 and $_3$9, the estimated toroidal moment scale rises from

$$\mathbf{T}=\sum_i \mathbf{r}_i \times \boldsymbol{\mu}_i,$$0

to

$$\mathbf{T}=\sum_i \mathbf{r}_i \times \boldsymbol{\mu}_i,$$1

which the Fe$$\mathbf{T}=\sum_i \mathbf{r}_i \times \boldsymbol{\mu}_i,$$2Dy$$\mathbf{T}=\sum_i \mathbf{r}_i \times \boldsymbol{\mu}_i,$$3 study presents as an approximate $$\mathbf{T}=\sum_i \mathbf{r}_i \times \boldsymbol{\mu}_i,$$4 enhancement (Soncini et al., 5 Sep 2025).

The compound $$\mathbf{T}=\sum_i \mathbf{r}_i \times \boldsymbol{\mu}_i,$$5 contains alternating $$\mathbf{T}=\sum_i \mathbf{r}_i \times \boldsymbol{\mu}_i,$$6–$$\mathbf{T}=\sum_i \mathbf{r}_i \times \boldsymbol{\mu}_i,$$7 centers around an elliptical wheel. Ab initio CAHF/CASCI-SO calculations show that the Dy$$\mathbf{T}=\sum_i \mathbf{r}_i \times \boldsymbol{\mu}_i,$$8 sites possess well-isolated ground Kramers doublets composed of nearly pure $$\mathbf{T}=\sum_i \mathbf{r}_i \times \boldsymbol{\mu}_i,$$9 states with axial $_{10}$0-tensors,

$_{10}$1

and first excited Kramers doublets at $_{10}$2–$_{10}$3. Their principal axes are arranged tangentially around the ring, with planar projections following a vortex-like pattern. Fe$_{10}$4 ions contribute orbitally nondegenerate $_{10}$5 ground multiplets with negligible zero-field splitting and $_{10}$6. Broken-symmetry DFT yields asymmetric nearest-neighbor ferromagnetic Fe–Dy couplings $_{10}$7–$_{10}$8 and $_{10}$9–$_6$00, together with weak antiferromagnetic Fe–Fe couplings $_6$01 to $_6$02 (Soncini et al., 5 Sep 2025).

These interactions generate a dense manifold of toroidal excitations. Even when the description is restricted to Dy ground Kramers doublets and Fe $_6$03 manifolds, the product space contains

$_6$04

states, described in the paper as a “$_6$05 billion dimensional toroidal space.” The spectrum is organized by an Ising-band structure with $_6$06-, $_6$07-, $_6$08-, $_6$09-, $_6$10-, and $_6$11-wave distortions of the vortex pattern. A magnetic state with moment $_6$12 lies only $_6$13 above the nonmagnetic toroidal ground doublet, so a uniform field of $_6$14 along an easy axis is sufficient to interchange toroidal and magnetic ground states; the same model reproduces powder magnetization at $_6$15 and specific-heat data up to $_6$16 (Soncini et al., 5 Sep 2025).

4. Theoretical descriptions and toroidal response functions

SMT modeling typically starts from a non-collinear Ising or anisotropic spin Hamiltonian in which strong crystal-field anisotropy projects the lanthanide sites onto local Kramers-doublet pseudospins. For the Dy$_6$17 neutron study the spin-Hamiltonian context is written

$_6$18

with local response expressed through

$_6$19

so that polarized neutron diffraction can refine the site-resolved susceptibility tensors $_6$20 directly (Jenkins et al., 8 Jun 2026).

For Fe$_6$21Dy$_6$22, the effective Hamiltonian is decomposed as

$_6$23

where the Dy sites are treated as Ising spins $_6$24 with $_6$25, the Fe ions remain explicit quantum $_6$26 spins, $_6$27 contains intrafragment dipolar interactions, $_6$28, and $_6$29. Because the weak Fe–Fe exchange is treated perturbatively to first order, the zeroth-order Hamiltonian factorizes into local Fe-site Hamiltonians $_6$30, and the partition function can be compressed into traces of products of $_6$31 transfer matrices,

$_6$32

This reduction turns the $_6$33-state Hilbert space into a computationally manageable thermodynamic problem while retaining excellent agreement with experiment (Soncini et al., 5 Sep 2025).

The same formulation yields standard magnetic observables,

$_6$34

and also the toroidal expectation and response,

$_6$35

The Fe$_6$36Dy$_6$37 paper introduces $_6$38 as the toroidal susceptibility tensor, and its molar form $_6$39, as a thermodynamic linear-response function that measures finite-temperature toroidal polarization induced by a weak magnetic-field curl. In analogy with Van Vleck susceptibility, $_6$40 contains a diagonal fluctuation term proportional to $_6$41 and an off-diagonal virtual-transition term. In the $_6$42 limit,

$_6$43

so $_6$44 directly encodes the ground-state toroidal moment (Soncini et al., 5 Sep 2025).

Direct evaluation of $_6$45 for Fe$_6$46Dy$_6$47 predicts a sizable ground-state toroidal moment and an enhanced finite-temperature toroidal polarization up to $_6$48 under modest uniform fields, attributed to field-induced spreading of toroidal levels across the spectrum. This suggests that toroidal thermodynamics in larger wheels is not limited to a single isolated doublet but reflects the structure of a broad, densely packed vortex manifold (Soncini et al., 5 Sep 2025).

5. Experimental identification of toroidal order

Historically, SMT assignments relied largely on magnetometry and ab initio calculations. The Dy$_6$49 neutron-diffraction study identifies this as a limitation, because bulk magnetization suppression and simulated local axes cannot directly resolve site-specific non-collinearity inside a crystal. Its central methodological advance is the combined use of polarized neutron diffraction (PND), variable-field single-crystal neutron diffraction, ab initio calculations, and magnetometry as a quantitative framework for probing molecular toroidal order (Jenkins et al., 8 Jun 2026).

In PND, flipping ratios

$_6$50

are measured at Bragg reflections and fitted to obtain local susceptibility tensors $_6$51 through $_6$52. For Dy$_6$53-2, the measurements were performed on DEMAND at $_6$54, neutron polarization $_6$55, $_6$56, and $_6$57, using $_6$58 good-quality flipping ratios and yielding a goodness of fit $_6$59. The refined Dy-site susceptibility ellipsoids are strongly anisotropic, nearly tangential to the triangle, and contain nonzero off-diagonal terms that directly evidence non-collinearity. This constitutes an atomic-scale visualization of the head-to-tail arrangement responsible for toroidicity (Jenkins et al., 8 Jun 2026).

Variable-field neutron diffraction complements PND by resolving field-induced order. With $_6$60 at $_6$61, the $_6$62 Bragg reflection shows the onset of field-induced magnetic order at $_6$63, peaking near $_6$64. At $_6$65, magnetic scattering appears below $_6$66. Magnetic symmetry refinement gives the best description in $_6$67, with an in-plane ordered moment $_6$68 per Dy and a fixed $_6$69-axis component $_6$70 chosen to match bulk data. The refined pattern combines antiferromagnetic alignment of toroidal domains along $_6$71 with head-to-tail in-plane alignment, thereby establishing antiferrotoroidic stacking together with a field-orderable $_6$72-axis dipolar component (Jenkins et al., 8 Jun 2026).

Magnetometry is consistent with this toroidal ground state. For Dy$_6$73-2, Curie–Weiss temperatures are $_6$74 for $_6$75, $_6$76 for $_6$77, and $_6$78 for powder; $_6$79 at $_6$80 is $_6$81 and $_6$82 for $_6$83 and $_6$84, respectively; and at $_6$85 the isothermal magnetization is $_6$86 per trimer for $_6$87 and $_6$88 for $_6$89. Powder samples show hysteresis at $_6$90, whereas single crystals do not, consistent with cancellation between clockwise and counterclockwise toroidal domains within the unit cell (Jenkins et al., 8 Jun 2026).

The Fe$_6$91Dy$_6$92 study addresses a different experimental bottleneck: direct detection of the toroidal degree of freedom itself. Its calculated splittings under realistic magnetic-field curls and its sizable finite-temperature toroidal susceptibility are presented as signatures amenable to direct observation, in contrast with Dy$_6$93, where the smaller molecular radius strongly limits detectability (Soncini et al., 5 Sep 2025).

6. Preparation, manipulation, and design principles

A major obstacle in SMT research has been the selective preparation of one member of a degenerate toroidal doublet. The microwave-control study proposes a realistic pulsed-EPR protocol that avoids nanometer-scale static field gradients. In Dy$_6$94, an in-plane static field

$_6$95

moves the toroidal doublet $_6$96 above a singly degenerate magnetic ground state $_6$97. Because linearly polarized microwave radiation perpendicular to the triangle plane flips one Dy pseudospin through the small transverse components $_6$98, the transition $_6$99 is resonant while the transition to $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$00 requires two flips and is off-resonant. A single resonant $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$01 pulse at $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$02 therefore selectively prepares one toroidal state, with pulse duration

$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$03

This is a genuinely coherent preparation mechanism, not a thermally biased one (Hymas et al., 11 Apr 2025).

In AlDy$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$04, preparation of a ferrotoroidic state requires three Dy flips. With

$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$05

and a weak canting asymmetry $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$06, three sequential resonant pulses at

$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$07

selectively transfer population from $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$08 to $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$09. The toroidal polarization is monitored through

$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$10

A simultaneous-pulse variant, using three pulses over $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$11 with $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$12 and $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$13, yields $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$14 population in $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$15 in the simulations. Dissipative dynamics treated with a Redfield-type master equation show robustness over the literature-estimated relaxation range $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$16–$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$17, while very large relaxation rates drive population back toward ground or metastable states (Hymas et al., 11 Apr 2025).

A complementary control route is direct coupling to $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$18. For Fe$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$19Dy$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$20, Maxwell–Ampère,

$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$21

suggests two experimental realizations: local current injection perpendicular to the molecular ring, for example by STM tips, and displacement-current curls generated by spatially focused, shaped femtosecond laser pulses. For a Ti:sapphire pulse at $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$22, $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$23, $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$24 spot, and peak power $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$25, the estimated peak field is $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$26, giving

$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$27

The corresponding toroidal splittings are then

$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$28

with the larger wheel predicted to be readily observable at liquid-helium temperatures. Proposed signatures include pump–probe changes in magnetization, shifts or enhancements in low-temperature specific heat, and AC susceptibility under modulated $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$29 (Soncini et al., 5 Sep 2025).

Several design rules recur across these studies. Larger rings amplify $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$30 and hence $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$31; high-moment ions such as Dy$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$32 maximize $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$33; and tangential alignment of Dy easy axes around a ring or triangle stabilizes the vortex texture. Mixed $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$34–$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$35 architectures can use Fe$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$36 spins to shape the spectrum while keeping Fe–Fe antiferromagnetism weak enough not to suppress toroidicity. Large Kramers-doublet gaps, such as the $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$37–$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$38 range in Fe$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$39Dy$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$40, help justify Ising projections at liquid-helium temperatures. In Dy$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$41-2, the $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$42-OH proton is chemically decisive: removing it and replacing $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$43-OH by $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$44-O$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$45 reorients the anisotropy to $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$46, $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$47, destroying toroidal alignment (Jenkins et al., 8 Jun 2026).

These properties motivate qubit-oriented interpretations. Toroidal doublets have zero total magnetic dipole moment and are first-order insensitive to uniform magnetic fields, so they are less sensitive to stray-field noise than conventional dipolar states. The microwave study presents coherent $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$48-pulse-like transfers between toroidal configurations, while the Fe$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$49Dy$M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$50 analysis argues that protection against uniform $M=\mathrm{Al}^{3+},\mathrm{Cr}^{3+}$51 and direct addressability by field curls make larger SMTs attractive for robust spin states. A plausible implication is that future progress will depend on combining chemically engineered toroidicity, direct toroidal probes, and coherent control protocols in the same molecular platform (Hymas et al., 11 Apr 2025).