---
title: 'Molar Toroidal Susceptibility: Theory & Applications'
url: https://www.emergentmind.com/topics/molar-toroidal-susceptibility
type: topic
---

# Molar Toroidal Susceptibility: Theory & Applications

Searching arXiv for the cited literature and topic scope.
arXiv search: molar toroidal susceptibility
Molar toroidal susceptibility is the thermodynamic linear-response coefficient that quantifies how a bulk ensemble develops a toroidal polarization per mole under a field conjugate to a toroidal moment. In the strictest formulation presently available in the cited literature, it is introduced for the Fe\(_{10}\)Dy\(_{10}\) molecular ring as a response to a magnetic field with a small non-vanishing curl, in direct analogy with the Van Vleck magnetic susceptibility [2509.05424]. Across the broader toroidal-multipole literature, the same quantity is often implicit rather than explicitly named: it is reconstructed from free-energy couplings, symmetry-allowed magnetoelectric tensors, nonlinear transport, or microscopic toroidal operators in localized-spin, plaquette, cluster, and atomic-scale models [1006.5199], [2208.01754], [1903.04708], [2403.09492].

## 1. Definition, conjugate fields, and thermodynamic status

In the Fe\(_{10}\)Dy\(_{10}\) work, the toroidal susceptibility tensor is defined through the Helmholtz free energy \(F\) as
\[
\xi_{\alpha \beta} = - \left. \frac{\partial^2 F}{\partial \left( \nabla \times \mathbf{B} \right)_{\alpha} \partial \left( \nabla \times \mathbf{B} \right)_{\beta} } \right|_{\nabla \times \mathbf{B} = 0},
\]
with \(\nabla\times\mathbf B\) treated as the external variable conjugate to the toroidal moment operator \(\boldsymbol{\tau}\) [2509.05424]. The corresponding toroidal polarization satisfies
\[
P_{\tau,\alpha} = -\frac{\partial F}{\partial (\nabla\times\mathbf B)_\alpha},
\]
so that, in linear response, \(P_{\tau,\alpha} = \sum_\beta \xi_{\alpha\beta}(\nabla\times\mathbf B)_\beta\) [2509.05424].

A more general symmetry framework distinguishes electric, magnetic, electric-toroidal, and magnetic-toroidal multipoles. Within that classification, the magnetic toroidal dipole \(\mathbf T\) is time-odd and parity-odd, while the electric toroidal dipole \(\mathbf G\) is time-even and parity-even; their conjugated fields differ accordingly. For rank-1 multipoles, \(\mathbf T\) is conjugated to \(\nabla\times\mathbf B\) or \(\partial_t \mathbf E\), whereas \(\mathbf G\) is conjugated to \(\nabla\times\mathbf E\) or \(\partial_t \mathbf B\) [2403.09492]. This point is fundamental: a toroidal susceptibility is not, in general, the response to a uniform \(\mathbf B\) or \(\mathbf E\), even though effective low-field proxies can emerge through symmetry-allowed couplings.

A closely related construction appears in the two-dimensional plaquette model of toroidal order, where the magnetoelectric free energy is written as
\[
{\cal F}_{\rm ME}
= T_z \left\{
A(H_x E_y - H_y E_x) + B H_z E_z + C(H_x E_x + H_y E_y)
\right\},
\]
and the effective conjugate field is therefore
\[
f_z(\mathbf E,\mathbf H)
=
A(H_x E_y - H_y E_x) + B H_z E_z + C(H_x E_x + H_y E_y).
\]
In that formulation, the toroidal susceptibility is
\[
\chi^{(T)}_{zz}
=
\left.\frac{\partial T_z}{\partial f_z}\right|_{f_z=0}.
\]
This establishes a second, equally legitimate thermodynamic definition: the conjugate field may be a symmetry-allowed bilinear combination of \(\mathbf E\) and \(\mathbf H\), rather than a bare curl field [1006.5199].

## 2. Microscopic realizations of toroidal moments

The microscopic object whose molar response is being measured is model dependent. In cluster and spin-texture language, toroidal dipoles are vortex-like moments of dipole distributions, schematically \(\mathbf T \propto \sum_j \mathbf R_j \times \mathbf M_j\) for magnetic dipoles and \(\mathbf G \propto \sum_j \mathbf R_j \times \mathbf Q_j\) for electric dipoles [2403.09492]. In a localized-spin cluster, this reduces to the familiar discrete spin form used for trimers and molecular rings.

In the 3-spin toroidal plaquette model, the plaquette toroidicity is defined so that the leading magnetostatic interaction between two plaquettes is proportional to \(T_A T_B/R^7\). The explicit scalar toroidal strength is
\[
T_i = \frac{3}{2} m_i r^2 = 3 Q_i \chi_i r^3,
\]
where the sign is an Ising-like variable set by chirality [1006.5199]. That definition is deliberately adapted to the plaquette interaction and is not universal, but it provides a concrete microscopic toroidal variable whose susceptibility can be defined thermodynamically.

In Ce\(_3\)TiBi\(_5\), the relevant object is a magnetic toroidal dipole generated by staggered antiferromagnetic order on locally noncentrosymmetric zigzag chains. The local toroidal component is written as
\[
T_i^y \propto (\mathbf V_i \times \mathbf m_i)^y,
\]
with staggered local crystal fields \(\mathbf V_{\rm A}=-\mathbf V_{\rm B}\parallel \hat{\mathbf x}\) and staggered moments \(\mathbf m_{\rm A}=-\mathbf m_{\rm B}\parallel \hat{\mathbf z}\), producing a ferroic in-plane toroidal dipole. At the unit-cell level, these chain toroidal moments form cluster multipoles \(T_x^{(\mathrm c)}\), \(T_y^{(\mathrm c)}\), and \(T_{3b}^{(\mathrm c)}\), and the partial-disordered toroidal phase is identified with an in-plane cluster toroidal dipole [2208.01754].

At the atomic scale, the two-orbital \(d\)-\(f\) model with odd-parity hybridization yields toroidal operators directly in orbital-spin space. The in-plane and out-of-plane magnetic toroidal dipoles are
\[
(\hat T_x,\hat T_y)=(-\sigma_y\tau_x,\sigma_x\tau_x),\qquad
\hat T_z=-\sigma_0\tau_y,
\]
with \(T_x\) and \(T_y\) spin dependent and \(T_z\) purely orbital [1903.04708]. In that setting, toroidal susceptibility refers to the field response of expectation values \(\langle \hat T_\mu\rangle\), not to a geometric spin texture alone.

## 3. From microscopic response to a molar quantity

The transition from a microscopic toroidal response to a molar susceptibility is conceptually straightforward but convention dependent. The core step is to identify the toroidal moment per microscopic unit—plaquette, site, unit cell, or molecule—and then multiply the corresponding single-unit response by the number of such units per mole.

For the plaquette model, the molar toroidal susceptibility is written as
\[
\chi^{(T)}_{\rm mol}
=
n_p N_A
\left.\frac{\partial \langle T_i\rangle}{\partial f_z}\right|_{f_z=0},
\]
where \(n_p\) is the number of toroidal plaquettes per formula unit and \(N_A\) is Avogadro’s number [1006.5199]. The same paper notes that the precise units depend on how \(T\) and the conjugate field \(f_z\) are defined.

For Ce\(_3\)TiBi\(_5\), where the order parameter is a cluster toroidal dipole \(T_\gamma^{\rm cell}\), the corresponding molar quantity is
\[
\chi^{T,\mathrm{mol}}_{\gamma\gamma}
=
\frac{N_A}{n_{\mathrm{Ce}}}
\left.
\frac{\partial T_\gamma^{\mathrm{cell}}}{\partial \mathcal B_\gamma}
\right|_{\mathcal B_\gamma=0},
\]
with \(n_{\mathrm{Ce}}=3\) Ce ions per formula unit if the susceptibility is normalized per mole of Ce, or \(n_{\mathrm{f.u.}}=1\) if it is normalized per mole of formula units [2208.01754]. This makes explicit that “molar toroidal susceptibility” is not a unique number until the normalization convention is stated.

A complementary route starts from ordinary magnetic susceptibility tensors of localized centers. For a cluster of spins at positions \(\mathbf r_i\) with local linear responses \(\mathbf m_i=\boldsymbol{\chi}^{(i)}\mathbf H\), the toroidal moment is
\[
\mathbf T = \frac{1}{2}\sum_i \mathbf r_i \times \boldsymbol{\chi}^{(i)}\mathbf H,
\]
and the toroidal susceptibility tensor becomes
\[
\chi^T_{\alpha\beta}
=
\frac{1}{2}
\sum_i
\epsilon_{\alpha\gamma\delta}\,
r_{i,\gamma}\,
\chi^{(i)}_{\delta\beta}.
\]
Once the unit-cell geometry is known, this tensor can be converted to a molar toroidal susceptibility by the usual scaling from unit-cell to per-mole quantities [1312.2401].

The Fe\(_{10}\)Dy\(_{10}\) paper uses the same logic at the molecular level. It computes \(\xi_{\alpha\beta}\) for a single ring and then interprets it as the molar toroidal susceptibility by the usual scaling with Avogadro’s number when estimating experimental signals [2509.05424].

## 4. Magnetoelectric, transport, and optical proxies

In much of the literature, molar toroidal susceptibility is not measured directly; it is inferred from symmetry-related response functions. The most common proxy is the linear magnetoelectric tensor. In the ferrotoroidal phase of the plaquette model, the allowed linear magnetoelectric components follow directly from the free energy \( {\cal F}_{\rm ME}\), and they vanish when the uniform toroidal moment \(T_z\) vanishes [1006.5199]. This makes the slope of toroidal-order-dependent magnetoelectric coefficients a practical surrogate for toroidal susceptibility.

Ce\(_3\)TiBi\(_5\) provides a metallic realization of this logic. The paper does not introduce an explicit tensor named toroidal susceptibility, but it identifies toroidal-origin response functions in the partial-disordered toroidal phase. The linear magnetoelectric effect is written as \(M_\nu=\alpha_{\mu;\nu}E_\mu\), with symmetry-allowed components in the PD-MT state given by \(\alpha_{y;z}\propto \alpha_z\) and \(\alpha_{z;y}\propto \alpha_z^2\). Nonlinear transport coefficients,
\[
J_\mu=\sigma_{\mu;\nu\eta}E_\nu E_\eta,\qquad
J_\mu^{(\mathrm s)}=\sigma_{\mu;\nu\eta}^{(\mathrm s)}E_\nu E_\eta,
\]
also become finite only in phases with a net toroidal dipole. The paper therefore treats toroidal susceptibility as something that must be inferred from magnetoelectric and nonlinear transport responses, not read literally from the Hamiltonian [2208.01754].

The atomic-scale \(d\)-\(f\) theory makes this relation even more explicit. The symmetry form of the field-driven magnetoelectric tensor is
\[
\hat{\alpha}^{(\mathrm E)}
\sim
\begin{pmatrix}
M_u & T_z & 0\\
-T_z & M_u & 0\\
T_y & -T_x & 0
\end{pmatrix},
\]
so the antisymmetric in-plane part tracks \(T_z\), while the \(zx\) and \(zy\) components track in-plane toroidal dipoles \(T_x\) and \(T_y\) [1903.04708]. Within that model, odd-parity hybridization enhances the magnetoelectric effect for in-plane magnetic toroidal dipoles and suppresses it for the out-of-plane one. A plausible implication is that the associated molar toroidal susceptibilities inherit the same enhancement and suppression trends.

Optical activity provides a further, nonlocal proxy. In the metamaterial study of toroidal optical activity, the toroidal dipole moment
\[
\mathbf T=\frac{1}{10c}\int \left[(\mathbf r\cdot\mathbf j)\mathbf r-2r^2\mathbf j\right]\,d^3r
\]
must be included to reproduce the observed circular dichroism at one resonance, and the analysis attributes the effect to toroidal dipole plus electric quadrupole contributions rather than to the usual electric and magnetic dipoles [1508.06192]. This suggests a frequency-dependent toroidal susceptibility in the optical channel, although the paper does not normalize it in molar form.

## 5. Material realizations and characteristic temperature or field dependence

The first explicit molar toroidal susceptibility in the cited corpus is the Fe\(_{10}\)Dy\(_{10}\) ring. There the molecular toroidal moment operator is
\[
\boldsymbol{\tau}=\sum_i \mathbf r_i \times \mathbf M_i,
\]
and the Hamiltonian contains an explicit toroidal coupling
\[
H_{\rm Tor}=\boldsymbol{\tau}\cdot(\nabla\times\mathbf B).
\]
The calculated \(\xi_{\alpha\alpha}(T)T\) tends to a finite constant as \(T\to 0\), which the authors use to identify a non-zero toroidal moment in the ground state. The maximal toroidal moment is estimated as \(|\boldsymbol{\tau}_{\mathrm{Fe}_{10}\mathrm{Dy}_{10}}|\approx 1500\,\mu_B\cdot\text{\AA}\), compared with \(|\boldsymbol{\tau}_{\mathrm{Dy}_3}|\approx 60\,\mu_B\cdot\text{\AA}\) for Dy\(_3\), and for an estimated curl \(|\nabla\times\mathbf B|\sim 5\times 10^{-4}\,\mathrm{T}\,\text{\AA}^{-1}\) the corresponding toroidal splitting is \(\approx 1.5\,\mathrm{cm}^{-1}\) rather than \(\approx 0.06\,\mathrm{cm}^{-1}\) [2509.05424]. In this literature, that scale difference is the most direct sense in which a “giant” molar toroidal susceptibility is claimed.

In the two-dimensional plaquette model, the relevant critical behavior is Ising-like. The gauge-toroid phase is characterized by ordering of the composite variable \(Y_i=T_i\sin(3\phi_i)\), whereas true ferrotoroidal order appears only when the gauge symmetry is broken at a lower temperature. The paper states that the toroidal susceptibility should diverge or show a strong peak at the ferrotoroidal transition; within Landau mean-field theory \(\chi^{(T)}\propto 1/(T-T_c)\), while in the 2D Ising universality class one expects \(\chi^{(T)}\sim |T-T_c|^{-\gamma}\) with \(\gamma=7/4\) [1006.5199].

In BaCoSiO\(_4\), the paper does not define toroidal susceptibility explicitly, but it presents a field-tunable toroidal order parameter. The ferritoroidal state in zero field carries a net \(+1\mathbf t\) or \(-1\mathbf t\) per macroscopic domain, while a field \(H\parallel c\) drives a transition at \(\mu_0H_c\sim 1.2\,\mathrm T\) to a ferrotoroidal \(+3\mathbf t\) or \(-3\mathbf t\) state [2103.01360]. In the paper’s own reconstruction, one can naturally define
\[
\chi_T = \frac{\partial T_z}{\partial H_z},
\]
or, in molar form, \( \chi_T^{(\mathrm{mol})}=N_A\,\partial T_z^{(\mathrm{f.u.})}/\partial H_z\). Under that definition, the first-order toroidal transition yields a large anomaly in toroidal susceptibility [2103.01360].

Ce\(_3\)TiBi\(_5\) contributes a different hallmark: an unusual anisotropic magnetic susceptibility in the partial-disordered toroidal regime. At \(J_2/J_1=-0.8\) and \(J_3/J_1=1.2\), \(\chi^{zz}\) shows a cusp at \(T_N/J_1=3\), whereas \(\chi^{xx}\) does not show a cusp and continues to increase below \(T_N\), with \(\chi^{xx}>\chi^{zz}\) across the full temperature range [2208.01754]. The paper explicitly interprets this as an indirect fingerprint of a soft toroidal order parameter, although it does not compute \(\chi^T\) itself.

## 6. Conceptual limits, common confusions, and current directions

A first recurrent confusion is to identify toroidal susceptibility with ordinary magnetic susceptibility. The literature does not support that identification. In the strict thermodynamic construction, the magnetic toroidal dipole is conjugated to \(\nabla\times\mathbf B\) or \(\partial_t\mathbf E\), or to symmetry-allowed field bilinears such as \(f_z(\mathbf E,\mathbf H)\); a uniform magnetic field becomes relevant only indirectly, through composite couplings or through specific lattice symmetries [2403.09492], [1006.5199], [2509.05424].

A second confusion is to assume that toroidal order is equivalent to a net magnetic dipole or to a spontaneous polarization. The plaquette model explicitly states that no phase permits a spontaneous polarization under its assumed mirror symmetry, even in the ferrotoroidal phase [1006.5199]. Conversely, Ce\(_3\)TiBi\(_5\) realizes an in-plane toroidal dipole generated by staggered antiferromagnetic order in a globally centrosymmetric but locally noncentrosymmetric lattice, and its partial-disordered toroidal phase contains one nonmagnetic chain rather than a conventional uniform magnetization [2208.01754].

A third issue is methodological. Many papers provide all ingredients needed to define molar toroidal susceptibility but do not actually tabulate it. The atomic-scale visualization of toroidal order in Dy\(_3\) shows how polarized neutron diffraction yields local Dy\(^{3+}\) susceptibility tensors, variable-field neutron diffraction resolves field-induced toroidal ordering, and ab initio calculations recover the toroidal ground state, yet the paper states that toroidal susceptibility is not explicitly defined or numerically tabulated there [2606.10077]. Similarly, the coherent-control work on Dy-based molecular toroics computes time-dependent \(\langle \tau_z(t)\rangle\) under pulsed microwave fields and argues that the ensuing magneto-electric properties can be used as a read-out mechanism, but it does not introduce a closed-form molar susceptibility [2504.08701].

The present state of the subject is therefore bifurcated. On one side, a fully thermodynamic molar toroidal susceptibility has been written down and evaluated for a specific molecular ring [2509.05424]. On the other, a larger body of work treats toroidal susceptibility as a derived or inferred quantity, reconstructed from toroidal order parameters, symmetry-allowed free-energy terms, local susceptibility tensors, or cross-correlated responses [2208.01754], [1903.04708], [1006.5199], [2403.09492]. This suggests that, in current usage, “molar toroidal susceptibility” is best understood not as a single universally standardized observable, but as a family of response coefficients whose common structure is the derivative of a molar toroidal polarization with respect to the correctly specified conjugate field.

Source: https://www.emergentmind.com/topics/molar-toroidal-susceptibility