Ferro-Rotational Systems Overview

• Ferro-rotational systems are ordered states where materials display a uniform sense of rotation defined by axial order parameters such as n, T, or G.
• They are characterized through advanced techniques like EQ-SHG and RA-SHG that reveal broken mirror symmetry, domain formation, and weak first-order phase transitions.
• This concept extends to nuclear structure as ferro-deformation, where reinforced axial quadrupole deformations quantified by R values highlight cross-disciplinary applications.

Ferro-rotational (FR) systems are ordered states in which a material, lattice, or finite many-body system acquires a uniform sense of rotation rather than a net electric polarization or magnetization. In condensed-matter usage, FR order is also described as “ferro-axial” or “2D chirality,” and is encoded by axial quantities such as $\mathbf n$, $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$, $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$, or $G=\int \mathbf r\times \mathbf p(\mathbf r)\,d^3r$. In nuclear-structure usage, a related “ferro-deformation” denotes a reinforced axial quadrupole deformation diagnosed by $R\equiv E(4_1^+)/E(2_1^+)\simeq 3.3$. Across these literatures, FR behavior is associated with broken mirror symmetry, domain formation, sharp phase transitions or crossovers, and response functions that are not captured by the standard ferroelectric or ferromagnetic paradigms (Jin et al., 2019, Guo et al., 2022, Moon, 2016, Moon, 2016).

1. Definitions, order parameters, and symmetry class

In the condensed-matter ferroic taxonomy, ferro-rotational order is the vector-type order whose order parameter is an axial vector that describes a uniform sense of rotation of structural units such as octahedra. One explicit classification writes ferroelectric $P$ as TR$(+)$, SI$(-)$; ferromagnetic $M$ as TR$(-)$, SI$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$0; ferro-toroidal $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$1 as TR$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$2, SI$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$3; and ferro-rotational $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$4 as TR$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$5, SI$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$6. In this formulation, $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$7 transforms as the $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$8 irreducible representation of the parent $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$9 point group, and FR order completed the roster of vector ferroics when it was directly observed in RbFe(MoO$\Phi=\sum_i \mathbf r_i\times \mathbf P_i$0)$\Phi=\sum_i \mathbf r_i\times \mathbf P_i$1 (Jin et al., 2019).

Several later condensed-matter formulations use closely related axial variables. In NiTiO$\Phi=\sum_i \mathbf r_i\times \mathbf P_i$2, the FR order parameter is written as $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$3, with inversion-even and time-even transformation properties. In Ni$\Phi=\sum_i \mathbf r_i\times \mathbf P_i$4TeO$\Phi=\sum_i \mathbf r_i\times \mathbf P_i$5, the axial vector $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$6 is likewise inversion-even and mirror-odd, and its existence is presented as the prerequisite background that binds polar and chiral orders. In MnTiO$\Phi=\sum_i \mathbf r_i\times \mathbf P_i$7, the corresponding pseudovector $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$8 is also inversion-even and time-even, while the ordered state preserves $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$9 and $G=\int \mathbf r\times \mathbf p(\mathbf r)\,d^3r$0 but breaks the vertical mirrors containing the FR axis (Guo et al., 2022, Zhang et al., 10 Sep 2025, Zhang et al., 2024).

A distinct but related formulation appears in the electrotoroidicity literature, where FR order is described as a uniform sense of rotation of local electric dipoles within each unit cell, and the natural order parameter is the electrotoroidal moment

$G=\int \mathbf r\times \mathbf p(\mathbf r)\,d^3r$1

That work also states $G=\int \mathbf r\times \mathbf p(\mathbf r)\,d^3r$2 and $G=\int \mathbf r\times \mathbf p(\mathbf r)\,d^3r$3, while describing FR materials as preserving $G=\int \mathbf r\times \mathbf p(\mathbf r)\,d^3r$4 and $G=\int \mathbf r\times \mathbf p(\mathbf r)\,d^3r$5 and breaking all mirror planes and two-fold axes perpendicular to the rotation axis. This indicates that recent FR literature uses distinct but overlapping conventions for the axial variable that encodes the ordered state (Du et al., 1 Oct 2025).

2. Landau and Ginzburg–Landau descriptions

A minimal Landau–Ginzburg description near the FR transition in RbFe(MoO$G=\int \mathbf r\times \mathbf p(\mathbf r)\,d^3r$6)$G=\int \mathbf r\times \mathbf p(\mathbf r)\,d^3r$7 is

$G=\int \mathbf r\times \mathbf p(\mathbf r)\,d^3r$8

Here $G=\int \mathbf r\times \mathbf p(\mathbf r)\,d^3r$9, $R\equiv E(4_1^+)/E(2_1^+)\simeq 3.3$0, and $R\equiv E(4_1^+)/E(2_1^+)\simeq 3.3$1 yield a weak first-order transition, while $R\equiv E(4_1^+)/E(2_1^+)\simeq 3.3$2 is the conjugate field. Minimization gives a discontinuous onset of the FR amplitude below $R\equiv E(4_1^+)/E(2_1^+)\simeq 3.3$3, consistent with the observed jump of the EQ-SHG component $R\equiv E(4_1^+)/E(2_1^+)\simeq 3.3$4 at $R\equiv E(4_1^+)/E(2_1^+)\simeq 3.3$5 K (Jin et al., 2019).

When FR order is coupled to magnetization, the electrotoroidicity framework uses a Landau expansion

$R\equiv E(4_1^+)/E(2_1^+)\simeq 3.3$6

The trilinear term $R\equiv E(4_1^+)/E(2_1^+)\simeq 3.3$7 is the lowest-order coupling that is linear in $R\equiv E(4_1^+)/E(2_1^+)\simeq 3.3$8, even under $R\equiv E(4_1^+)/E(2_1^+)\simeq 3.3$9, and selects the Hall-like geometry $P$0. For a spontaneous $P$1, minimization yields a susceptibility tensor of the form

$P$2

so the off-diagonal response is odd in the FR chirality. The same analysis notes that perfect $P$3 symmetry can forbid the linear term, whereas strain or lower symmetry can restore a nonzero $P$4 (Du et al., 1 Oct 2025).

In materials with intertwined orders, the FR background is embedded in a larger Ginzburg–Landau functional. For Ni$P$5TeO$P$6,

$P$7

with $P$8, $P$9, $(+)$0, and $(+)$1. The trilinear term $(+)$2 enforces the “closed set” of intertwined orders: whenever any two are nonzero, the third is induced. In the rigid-background limit $(+)$3, the trilinear term becomes an effective bilinear $(+)$4, and the sign of $(+)$5 locks the signs of $(+)$6 and $(+)$7 within a given FR background (Zhang et al., 10 Sep 2025).

3. Optical detection, symmetry resolution, and domain imaging

Because FR order is invariant under both spatial inversion and time reversal in much of the condensed-matter literature, conventional linear probes are often ineffective. Electric-quadrupole second-harmonic generation (EQ-SHG) has therefore become a principal symmetry-selective probe. In centrosymmetric media, the leading bulk SHG contribution is

$(+)$8

and rotational-anisotropy SHG (RA-SHG) uses the angular dependence of the resulting intensity to discriminate FR symmetry sectors and domain states (Jin et al., 2019).

In RbFe(MoO$(+)$9)$(-)$0, RA-SHG directly coupled to the centrosymmetric FR order through the EQ channel. Above $(-)$1 K, the pattern showed pure threefold symmetry locked to the mirror planes of $(-)$2; below $(-)$3, two domain types $(-)$4 and $(-)$5 with opposite ferro-rotational vectors appeared, and the observed $(-)$6 became a weighted sum of the two opposite-domain patterns. Immediately below $(-)$7, the domain weight $(-)$8; upon cooling, $(-)$9, consistent with domain nucleation, growth, and coarsening. The jump of $M$0 at $M$1, the spike of $M$2, and the abrupt onset of the rotation angle $M$3 established the transition as weak first order (Jin et al., 2019).

In 1T-TaS$M$4, EQ RA-SHG established the ferro-rotational nature of the commensurate charge-density-wave phase. Below $M$5 K, the crossed-channel intensity takes the symmetry-allowed form

$M$6

where $M$7 is the quadrupole amplitude and $M$8 is the FR domain angle. The onset of long-range FR order is tracked by a kink in $M$9 and a jump of $(-)$0 by a few degrees at the transition (Luo et al., 2021).

NiTiO$(-)$1 provided the first direct real-space visualization of FR domains and walls by scanning EQ-SHG microscopy. Above $(-)$2 K, the crystal is $(-)$3 with point group $(-)$4; below $(-)$5, ordering of Ni and Ti layers lowers the symmetry to $(-)$6 with point group $(-)$7, producing two FR domains related by the broken vertical mirrors. In a representative map, the domain population was $(-)$8, $(-)$9, the lateral domain-size distribution peaked in the $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$00–$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$01 range, and domain walls appeared as dark lines with a $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$02 SHG suppression and width $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$03. Local RA-SHG measured on the walls yielded a bow-tie pattern consistent with restoration of the mirror symmetry and a wall point group $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$04, proving the walls to be nonpolar (Guo et al., 2022).

Ni$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$05TeO$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$06 extended this methodology to a multimodal setting. RA-SHG resolved the point-group content of polarity and the FR background, polarization-resolved transmission circular birefringence (tCB) mapped chirality without admixture of $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$07, and scanning SHG microscopy revealed bright domain walls associated with emergent in-plane polarization. Correlated SHG and tCB linecuts showed anti-correlated enhancement of in-plane $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$08 and suppression of $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$09 at the walls (Zhang et al., 10 Sep 2025).

System	FR signature	Primary probe
RbFe(MoO$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$10)$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$11	weak first-order FR transition at $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$12 K	EQ RA-SHG
1T-TaS$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$13	six-lobe $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$14 below $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$15 K	EQ RA-SHG
NiTiO$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$16	dark, nonpolar walls with restored $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$17 symmetry	scanning EQ-SHG
Ni$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$18TeO$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$19	interlocked polar, chiral, and FR domains	RA-SHG, tCB, scanning SHG

4. Domain walls, switching, and dynamical control

Electrical switching of FR domains was demonstrated in nano-thick 1T-TaS$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$20, where the relevant axial order parameter $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$21 is even under both $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$22 and $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$23. Symmetry forbids a bilinear $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$24 coupling, so a uniform in-plane electric field cannot directly select one FR state over the other. The switching mechanism instead acts through domain walls: local inversion breaking in TaS$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$25 octahedra generates local dipoles $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$26Å, and a coarse-grained coupling of the form

$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$27

drives wall motion without coupling the field directly to the bulk $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$28. In the NCCDW and CCDW phases, the two FR states are the $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$29 and $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$30 star-of-David tilings rotated by $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$31 and $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$32, respectively. Bias-cooling from $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$33 K to $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$34 K under $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$35 selected a monodomain state according to the sign of $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$36, while isothermal switching at $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$37–$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$38 K occurred once $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$39 exceeded a critical $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$40. At $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$41 K, $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$42 V for $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$43, corresponding to $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$44 V/cm. The switched state was nonvolatile, rectangular hysteresis was observed, $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$45 cycles at $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$46 K showed no fatigue, and pulsed switching in a $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$47 nm device gave an average wall velocity $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$48 m/s (Liu et al., 2022).

The Ni$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$49TeO$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$50 domain-wall theory made the FR background explicit. For a rigid $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$51, the out-of-plane polarization follows a kink profile,

$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$52

the induced chirality follows

$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$53

and second-derivative couplings generate localized in-plane polarization components

$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$54

Because $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$55 is perpendicular to the wall and $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$56 is parallel to the wall, the result is a mixed Néel–Bloch wall rather than a purely Néel or purely Bloch object (Zhang et al., 10 Sep 2025).

NiTiO$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$57 shows a different wall phenotype. There the local restoration of $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$58 symmetry, together with the suppression rather than enhancement of SHG intensity at the wall, establishes that FR domain walls can be nonpolar despite the broken mirror symmetry of the adjacent domains. This sharply contrasts with ferroelectric walls, where a strong electric-dipole SHG contribution would normally brighten the wall (Guo et al., 2022).

A mechanically mediated route to symmetry circumvention was analyzed for a composite multiferroic torsional oscillator. A single-domain ferromagnetic particle with built-in electric polarization, attached to a torsional cantilever, experiences an effective Larmor field in the rotating frame,

$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$59

so an applied electric field can switch the magnetic moment through the chain $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$60 mechanical torque $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$61 rotation $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$62. In this model, $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$63 magnetization switching follows a soliton-like trajectory, parameter windows for switching were obtained in $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$64 and $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$65 phase diagrams, and the device-scale estimates were GHz torsional modes, $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$66–$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$67 V drive voltages, switching times of order $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$68 ns or less, and dissipation on the order of $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$69–$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$70 J (Chudnovsky et al., 2015).

5. Electrotoroidicity, nonlinear optics, orbital currents, and ferro-rotational phonons

The most extensive recent generalization of FR phenomenology is electrotoroidicity: a family of transverse electromagnetic responses that arise from spontaneous electrotoroidal moments in FR materials even when $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$71 and $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$72 remain unbroken. In doped ilmenite Fe$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$73Ti$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$74O$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$75, magnetic-force microscopy revealed a reduced diagonal susceptibility at FR domain walls, parameterized as

$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$76

with $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$77–$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$78 and $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$79 in the reported crystals. The parent $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$80 structure exhibits FR order below $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$81 K, partial Fe doping raises the ferrimagnetic Curie temperature to $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$82 K while preserving the FR space group $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$83, the spontaneous $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$84 from first-principles modeling is on the order of $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$85Å per cell, and MFM under an in-plane $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$86 kOe field yielded an effective $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$87 contrast corresponding to $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$88 emu/mol. Density-functional calculations under $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$89 in-plane strain predicted a linear $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$90 of order $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$91–$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$92 emu/mol for hole doping near $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$93 eV, with sign reversal between CW and CCW FR domains (Du et al., 1 Oct 2025).

MnTiO$\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$94 exhibited directionally asymmetric nonlinear optics without breaking inversion or time reversal. In point group $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$95, the FR order allows two independent magnetic-dipole SHG tensor elements, $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$96 and $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$97, with $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$98 changing sign under enantiomorph reversal. In circular-polarization SHG, the conversion asymmetry depends on interference between these tensor elements. At $\mathbf T=\sum_i \mathbf r_i\times \mathbf p_i$99 nm and $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$00 K, $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$01 and $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$02, giving maximal contrast: $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$03 in one orientation, and the contrast reversed upon flipping the crystal by $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$04 about an in-plane axis. The normalized contrast

$\Phi=\sum_i \mathbf r_i\times \mathbf P_i$05

reached $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$06, whereas above $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$07 K the asymmetry vanished as $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$08 (Zhang et al., 2024).

FR order can also generate purely orbital transport channels. In TiAu$\Phi=\sum_i \mathbf r_i\times \mathbf P_i$09, the static structural rotation is associated with an electric hexadecapole moment

$\Phi=\sum_i \mathbf r_i\times \mathbf P_i$10

which hybridizes $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$11 and $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$12 orbitals in a minimal tight-binding model. Under an applied electric field, this FR-induced hybridization generates longitudinal and unconventional Hall orbital currents. First-principles calculations for space group $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$13 TiAu$\Phi=\sum_i \mathbf r_i\times \mathbf P_i$14, where Au squares are rotated by $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$15, gave at the Fermi level a conventional Hall conductivity $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$16, together with rotation-induced $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$17 and $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$18 for $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$19, and $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$20 and $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$21 for $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$22 (Jo et al., 7 May 2025).

An additional dynamical layer appears in MnTiO$\Phi=\sum_i \mathbf r_i\times \mathbf P_i$23 through ferro-rotational phonons. Resonant inelastic X-ray scattering identified circularly polarized $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$24 phonons, with the lower-energy $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$25 meV mode corresponding to bond-bending octahedral rotation motion in TiO$\Phi=\sum_i \mathbf r_i\times \mathbf P_i$26 cages. The coupling

$\Phi=\sum_i \mathbf r_i\times \mathbf P_i$27

links the FR order parameter to phonon angular momentum. Experimentally, the $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$28 meV mode displayed a large $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$29 circular dichroism for $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$30, the sign of $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$31 flipped upon reversing $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$32, no dichroism appeared for large in-plane momentum, and the effect persisted above and below the Mn Néel temperature $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$33 K, establishing that the response is tied to FR order rather than magnetism. The standing-wave condensate picture proposed in that work treats the static FR order as emerging from degenerate circular phonons at $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$34 with zero net momentum but finite $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$35 (Huang et al., 26 Dec 2025).

6. Nuclear ferro-deformation and ferro-rotational analogies

In nuclear structure, the FR label is used for “ferro-deformation”: a reinforced axial quadrupole deformation in even–even nuclei, explicitly compared to spontaneous magnetization in a ferromagnet. Its principal diagnostic is

$\Phi=\sum_i \mathbf r_i\times \mathbf P_i$36

with the empirical correspondence $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$37 for spherical single-particle structure, $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$38 for vibrational nuclei, and $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$39 for an axial rotor. The microscopic mechanism is formulated in terms of proton and neutron pseudo-shells constructed by coherent mixing of spherical subshells. In the $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$40, $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$41 region, representative pseudo-shells include $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$42, $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$43, $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$44, $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$45, and $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$46; half-filling of the proton $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$47 at $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$48 and neutron $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$49 at $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$50 coincides with maximal $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$51. Systematics of $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$52 and $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$53 reveal an abrupt shape phase transition between $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$54 and $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$55 for $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$56, a peak $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$57 when $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$58 or $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$59 correlates with $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$60, and a deformation plateau over $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$61 and $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$62, centered at $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$63 or $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$64 and $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$65 or $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$66. The proposed microscopic driver is a strong isospin-dependent spin–orbit interaction between neutrons in $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$67 and protons in $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$68, encoded schematically in

$\Phi=\sum_i \mathbf r_i\times \mathbf P_i$69

At the borders $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$70 or $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$71 with $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$72 and $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$73, the literature expects shape coexistence between a strong axial-rotor branch and a near-spherical vibrational branch, specifically in $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$74, $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$75, $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$76, $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$77, $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$78, and $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$79 (Moon, 2016).

A second critical region appears for $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$80 and $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$81. There, ferro-deformation arises suddenly at $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$82 or $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$83, $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$84, and approaches $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$85 centered at $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$86, $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$87. The relevant pseudo-shells are the proton $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$88, formed from $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$89, and the neutron $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$90, formed from $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$91. When both are half-filled, the collective quadrupole moment is maximized. The associated shape coexistence is predicted around $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$92 in $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$93, $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$94, $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$95, $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$96, $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$97, and $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$98, and the driving interaction is identified as the strong neutron–proton spin–orbit coupling between neutron $\Phi=\sum_i \mathbf r_i\times \mathbf P_i$99 and proton $G=\int \mathbf r\times \mathbf p(\mathbf r)\,d^3r$00 orbitals with the same orbital angular momentum $G=\int \mathbf r\times \mathbf p(\mathbf r)\,d^3r$01 (Moon, 2016).

Taken together, the nuclear papers preserve the core “ferro” analogy while shifting the object of ordering from a crystallographic rotation pattern to a saturated collective deformation. This suggests that FR nomenclature now spans two technically distinct domains: hidden axial order in solids, and half-filled pseudo-shell reinforcement of rotational collectivity in finite nuclei.