Ferronics: Dynamic Ferroic Order and Applications

• Ferronics is the study of dynamic ferroic order, focusing on collective excitations called ferrons that are electric analogues of magnons.
• Experimental investigations in materials like NbOX2 demonstrate room-temperature THz emission with high quality factors and controllable polarization transport.
• The field supports diverse applications including RF tunable devices, reconfigurable electronics, and spintronic memories through engineered ferroic functionalities.

Ferronics denotes a developing body of research in which ferroic order is treated as an active dynamical resource rather than a static material property. In its narrowest and now most technically developed sense, ferronics is the study and technology of collective excitations of electric polarization in ferroelectrics—quasiparticles called ferrons—which are the electric analogues of magnons; in broader usage, the term also covers RF electronics built on ferroic-material functionalities, ferroelectric domain-programmed electronic structures, and strain- or carrier-controlled magnetic semiconductor devices (Tang et al., 2022, Bauer et al., 2023, Monroy-Villa et al., 25 Jun 2025). Across these usages, the common theme is the direct generation, transport, modulation, and readout of ferroic order parameters by light, electric fields, magnetic fields, strain, or thermal gradients.

1. Definitions and semantic scope

In the ferron-centered literature, ferrons are the elementary excitations of ferroelectric order that carry a finite electric dipole. In a ferroelectric, the order parameter is the spontaneous polarization $P$, and a ferron is a coherent, propagating or localized oscillation of that order. In displacive ferroelectrics, ferrons emerge from the concerted action of anharmonicity and broken inversion symmetry; the relevant fluctuations are longitudinal with respect to the equilibrium polarization $P_0$, and their dipole can be written as $$\delta p_{\mathbf q} = -\partial(\hbar \omega_{\mathbf q})/\partial E$$ (Tang et al., 2022, Bauer et al., 2023).

A second, explicitly transport-oriented definition describes ferronics as the transport and manipulation of electric polarization by lattice excitations in ferroelectrics, in direct analogy to magnonics and phononics. Within this formulation, ferrons transport polarization and heat, support polarization and heat currents, and admit linear-response coefficients such as ferron conductivity, ferroelectric Seebeck coefficients, and polarization Peltier coefficients (Bauer et al., 2021, Lopez et al., 29 Jun 2026).

The term is also used more broadly. One line of work defines ferronics as RF electronics built on ferroic-material functionalities, exemplified by voltage-tunable ferroelectric reactive networks for superconducting RF cavities. Another uses “Ferrotronic” or “Graphene Ferronics” for electronic effects arising from engineered ferroelectric domain polarization patterns in a substrate. A further extension treats ferronics as a device-level paradigm exploiting coupling among ferroelectricity, ferromagnetism, antiferromagnetism, and ferroelasticity for sensing, memory, and logic (Monroy-Villa et al., 25 Jun 2025, Wan et al., 2021, Liu et al., 2019). This semantic plurality is itself a feature of the field: ferronics names both a specific quasiparticle program and a broader ferroic-device agenda.

2. Theoretical foundations

The long-wavelength theory of ferrons is usually formulated with Landau–Ginzburg–Devonshire free energies. A representative ferroelectric free energy is

$$F(P)=\frac{\alpha}{2}P^2+\frac{\beta}{4}P^4+\frac{\gamma_6}{6}P^6+\dots-E(t)P,$$

or, in spatially resolved form,

$$F[\mathbf P]=\int d^3r\left(\tfrac{\alpha}{2}|\mathbf P|^2+\tfrac{\beta}{4}|\mathbf P|^4+\tfrac{\kappa}{2}|\nabla \mathbf P|^2-\mathbf E\cdot \mathbf P\right).$$

Near the ordered state, Landau–Khalatnikov or Landau–Khalatnikov–Tani dynamics gives

$$\Gamma \frac{dP}{dt}=-\frac{\partial F}{\partial P}+E(t), \qquad m_p\partial_t^2 P+\gamma \partial_t P=-\frac{\delta F}{\delta P},$$

while coupled phonon–polarization descriptions introduce a soft optical phonon coordinate $u$ with $P \approx P_0+\lambda u$ (Zhang et al., 7 Sep 2025, Tang et al., 2022).

Linearization around $P_0$ yields a gapped ferron dispersion. In one common notation,

$$\omega^2(\mathbf q)=\omega_0^2+c^2q^2,\qquad \omega_0^2=\frac{A}{m_p},\qquad c^2=\frac{\kappa}{m_p},$$

with $P_0$0. In surface geometries, long-range dipolar interactions generate distinct surface ferron branches; the high-frequency branch is pushed toward the surface ionic plasma frequency, while the low-frequency branch is highly anisotropic and remains in the THz band (Zhou et al., 2022).

The transport formulation parallels magnonics. Ferron and heat currents can be written as

$P_0$1

and the polarization accumulation obeys a diffusion equation with diffusion length $P_0$2 (Bauer et al., 2021, Bauer et al., 2023). In ultrafast THz-emission problems, the same order parameter enters radiative electrodynamics through

$P_0$3

which is central to ferron-based THz generation (Zhang et al., 7 Sep 2025).

3. Experimental establishment of ferrons

Room-temperature experimental ferronics was established most clearly in layered van der Waals niobium oxide dihalides $P_0$4 with $P_0$5. In these materials, each 2D layer consists of $P_0$6 octahedra; the polar axis is along $P_0$7, the $P_0$8-axis is nonpolar, and spontaneous polarization arises from Nb displacement toward a bridging O, producing two distinct Nb–O bond lengths $P_0$9. Raman-active and infrared-active optical phonons are stronger along the $$\delta p_{\mathbf q} = -\partial(\hbar \omega_{\mathbf q})/\partial E$$0 axis than along the $$\delta p_{\mathbf q} = -\partial(\hbar \omega_{\mathbf q})/\partial E$$1 axis, consistent with ferroelectric coupling (Zhang et al., 7 Sep 2025).

Ultrafast experiments on exfoliated NbOX$$\delta p_{\mathbf q} = -\partial(\hbar \omega_{\mathbf q})/\partial E$$2 nanofilms revealed intense, narrowband THz emission with multiple ferron modes at room temperature. For NbOI$$\delta p_{\mathbf q} = -\partial(\hbar \omega_{\mathbf q})/\partial E$$3, distinct peaks were observed at $$\delta p_{\mathbf q} = -\partial(\hbar \omega_{\mathbf q})/\partial E$$4 THz and $$\delta p_{\mathbf q} = -\partial(\hbar \omega_{\mathbf q})/\partial E$$5 THz, with additional modes near $$\delta p_{\mathbf q} = -\partial(\hbar \omega_{\mathbf q})/\partial E$$6 THz; NbOBr$$\delta p_{\mathbf q} = -\partial(\hbar \omega_{\mathbf q})/\partial E$$7 showed a $$\delta p_{\mathbf q} = -\partial(\hbar \omega_{\mathbf q})/\partial E$$8 THz peak and NbOCl$$\delta p_{\mathbf q} = -\partial(\hbar \omega_{\mathbf q})/\partial E$$9 a $$F(P)=\frac{\alpha}{2}P^2+\frac{\beta}{4}P^4+\frac{\gamma_6}{6}P^6+\dots-E(t)P,$$0 THz peak. Reported linewidths and quality factors include $$F(P)=\frac{\alpha}{2}P^2+\frac{\beta}{4}P^4+\frac{\gamma_6}{6}P^6+\dots-E(t)P,$$1 at $$F(P)=\frac{\alpha}{2}P^2+\frac{\beta}{4}P^4+\frac{\gamma_6}{6}P^6+\dots-E(t)P,$$2 THz in NbOI$$F(P)=\frac{\alpha}{2}P^2+\frac{\beta}{4}P^4+\frac{\gamma_6}{6}P^6+\dots-E(t)P,$$3, $$F(P)=\frac{\alpha}{2}P^2+\frac{\beta}{4}P^4+\frac{\gamma_6}{6}P^6+\dots-E(t)P,$$4 at $$F(P)=\frac{\alpha}{2}P^2+\frac{\beta}{4}P^4+\frac{\gamma_6}{6}P^6+\dots-E(t)P,$$5 THz in NbOBr$$F(P)=\frac{\alpha}{2}P^2+\frac{\beta}{4}P^4+\frac{\gamma_6}{6}P^6+\dots-E(t)P,$$6, $$F(P)=\frac{\alpha}{2}P^2+\frac{\beta}{4}P^4+\frac{\gamma_6}{6}P^6+\dots-E(t)P,$$7 at $$F(P)=\frac{\alpha}{2}P^2+\frac{\beta}{4}P^4+\frac{\gamma_6}{6}P^6+\dots-E(t)P,$$8 THz in NbOCl$$F(P)=\frac{\alpha}{2}P^2+\frac{\beta}{4}P^4+\frac{\gamma_6}{6}P^6+\dots-E(t)P,$$9, and a low-frequency NbOI$$F[\mathbf P]=\int d^3r\left(\tfrac{\alpha}{2}|\mathbf P|^2+\tfrac{\beta}{4}|\mathbf P|^4+\tfrac{\kappa}{2}|\nabla \mathbf P|^2-\mathbf E\cdot \mathbf P\right).$$0 ferron at $$F[\mathbf P]=\int d^3r\left(\tfrac{\alpha}{2}|\mathbf P|^2+\tfrac{\beta}{4}|\mathbf P|^4+\tfrac{\kappa}{2}|\nabla \mathbf P|^2-\mathbf E\cdot \mathbf P\right).$$1 THz with FWHM $$F[\mathbf P]=\int d^3r\left(\tfrac{\alpha}{2}|\mathbf P|^2+\tfrac{\beta}{4}|\mathbf P|^4+\tfrac{\kappa}{2}|\nabla \mathbf P|^2-\mathbf E\cdot \mathbf P\right).$$2 GHz and $$F[\mathbf P]=\int d^3r\left(\tfrac{\alpha}{2}|\mathbf P|^2+\tfrac{\beta}{4}|\mathbf P|^4+\tfrac{\kappa}{2}|\nabla \mathbf P|^2-\mathbf E\cdot \mathbf P\right).$$3. The time-domain signal comprises an initial single-cycle THz pulse followed by long-lived periodic oscillations, with damped oscillations persisting for up to $$F[\mathbf P]=\int d^3r\left(\tfrac{\alpha}{2}|\mathbf P|^2+\tfrac{\beta}{4}|\mathbf P|^4+\tfrac{\kappa}{2}|\nabla \mathbf P|^2-\mathbf E\cdot \mathbf P\right).$$4 ps and a fitted decay time $$F[\mathbf P]=\int d^3r\left(\tfrac{\alpha}{2}|\mathbf P|^2+\tfrac{\beta}{4}|\mathbf P|^4+\tfrac{\kappa}{2}|\nabla \mathbf P|^2-\mathbf E\cdot \mathbf P\right).$$5 ps in GaP detection. The maximum detected THz field from NbOI$$F[\mathbf P]=\int d^3r\left(\tfrac{\alpha}{2}|\mathbf P|^2+\tfrac{\beta}{4}|\mathbf P|^4+\tfrac{\kappa}{2}|\nabla \mathbf P|^2-\mathbf E\cdot \mathbf P\right).$$6 was $$F[\mathbf P]=\int d^3r\left(\tfrac{\alpha}{2}|\mathbf P|^2+\tfrac{\beta}{4}|\mathbf P|^4+\tfrac{\kappa}{2}|\nabla \mathbf P|^2-\mathbf E\cdot \mathbf P\right).$$7 kV/cm at room temperature, and the per-unit-thickness radiation efficiency was reported as up to five orders of magnitude greater than state-of-the-art semiconductor THz emitters (Zhang et al., 7 Sep 2025).

Independent evidence for coherent ferrons came from transient reflectance microscopy and stroboSCAT in NbOI$$F[\mathbf P]=\int d^3r\left(\tfrac{\alpha}{2}|\mathbf P|^2+\tfrac{\beta}{4}|\mathbf P|^4+\tfrac{\kappa}{2}|\nabla \mathbf P|^2-\mathbf E\cdot \mathbf P\right).$$8. A coherent mode at $$F[\mathbf P]=\int d^3r\left(\tfrac{\alpha}{2}|\mathbf P|^2+\tfrac{\beta}{4}|\mathbf P|^4+\tfrac{\kappa}{2}|\nabla \mathbf P|^2-\mathbf E\cdot \mathbf P\right).$$9 THz, matching the dominant infrared-active transverse-optical phonon along the polar axis, was launched by femtosecond optical pulses and propagated along the $$\Gamma \frac{dP}{dt}=-\frac{\partial F}{\partial P}+E(t), \qquad m_p\partial_t^2 P+\gamma \partial_t P=-\frac{\delta F}{\delta P},$$0-axis at extremely hypersonic velocities exceeding $$\Gamma \frac{dP}{dt}=-\frac{\partial F}{\partial P}+E(t), \qquad m_p\partial_t^2 P+\gamma \partial_t P=-\frac{\delta F}{\delta P},$$1 m/s. Reported group velocities ranged from $$\Gamma \frac{dP}{dt}=-\frac{\partial F}{\partial P}+E(t), \qquad m_p\partial_t^2 P+\gamma \partial_t P=-\frac{\delta F}{\delta P},$$2 km/s in $$\Gamma \frac{dP}{dt}=-\frac{\partial F}{\partial P}+E(t), \qquad m_p\partial_t^2 P+\gamma \partial_t P=-\frac{\delta F}{\delta P},$$3 nm flakes to $$\Gamma \frac{dP}{dt}=-\frac{\partial F}{\partial P}+E(t), \qquad m_p\partial_t^2 P+\gamma \partial_t P=-\frac{\delta F}{\delta P},$$4 km/s in $$\Gamma \frac{dP}{dt}=-\frac{\partial F}{\partial P}+E(t), \qquad m_p\partial_t^2 P+\gamma \partial_t P=-\frac{\delta F}{\delta P},$$5 nm flakes, with $$\Gamma \frac{dP}{dt}=-\frac{\partial F}{\partial P}+E(t), \qquad m_p\partial_t^2 P+\gamma \partial_t P=-\frac{\delta F}{\delta P},$$6 km/s measured in a $$\Gamma \frac{dP}{dt}=-\frac{\partial F}{\partial P}+E(t), \qquad m_p\partial_t^2 P+\gamma \partial_t P=-\frac{\delta F}{\delta P},$$7 nm flake. No measurable propagation was found along the nonpolar $$\Gamma \frac{dP}{dt}=-\frac{\partial F}{\partial P}+E(t), \qquad m_p\partial_t^2 P+\gamma \partial_t P=-\frac{\delta F}{\delta P},$$8-axis. The same work emphasized that these velocities are incompatible with bare optical phonons near $$\Gamma \frac{dP}{dt}=-\frac{\partial F}{\partial P}+E(t), \qquad m_p\partial_t^2 P+\gamma \partial_t P=-\frac{\delta F}{\delta P},$$9 and that the narrow-band THz emission at $u$0 is incompatible with a conventional phonon-polariton interpretation (Choe et al., 28 May 2025).

Electrical control completed the experimental case. In NbOI$u$1 films on interdigital electrodes with $u$2m gaps, fields of $u$3 kV/cm applied along the $u$4-axis reversed the polarity of ferron oscillations without changing their frequencies. The emitted THz amplitude versus applied field showed a pronounced hysteresis loop, and the simultaneously measured current showed similar hysteresis, establishing non-volatile ferron control through ferroelectric domain switching (Zhang et al., 7 Sep 2025).

4. Nonlinear, surface, and transverse transport regimes

Ferronics rapidly expanded beyond linear THz emission. In NbOI$u$5, resonant THz excitation of a $u$6 THz ferron coherently upconverted energy to a $u$7 THz optical phonon. Two-dimensional THz spectroscopy showed an off-diagonal cross peak linking $u$8 THz to $u$9 THz, establishing coherent intermode coupling. First-principles calculations and analytical modeling identified a dominant cubic coupling

$P \approx P_0+\lambda u$0

with fitted $P \approx P_0+\lambda u$1 meV $P \approx P_0+\lambda u$2 amu$P \approx P_0+\lambda u$3. Electric-field switching between $P \approx P_0+\lambda u$4 kV/cm reversed the phase of both the ferron and the upconverted phonon and produced hysteresis in the complex $P \approx P_0+\lambda u$5 THz amplitude, demonstrating electrically programmable nonlinear phononics (Subedi et al., 19 Mar 2026).

Thermal transport generated a distinct transverse effect. The ferron Hall effect was predicted for ferroelectrics in which a longitudinal thermal gradient drives polarization-carrying lattice excitations and a magnetic field produces Hall-type deflection. In a local-mode model for BaTiO$P \approx P_0+\lambda u$6,

$P \approx P_0+\lambda u$7

and the stochastic dynamics includes an antisymmetric gyroscopic term $P \approx P_0+\lambda u$8. The field-odd transverse redistribution of longitudinal kinetic energy generates a field-odd edge polarization accumulation,

$P \approx P_0+\lambda u$9

with a computed edge accumulation $P_0$0 in BaTiO$P_0$1 (Lopez et al., 29 Jun 2026).

Surface ferronics introduces a further mode class. For a uniaxial ferroelectric slab with polarization parallel to the surface, long-range dipolar interactions produce surface ferron branches with dispersions

$P_0$2

The lower branch is highly anisotropic and linearly polarized transverse to its in-plane wavevector, while its evanescent stray field in vacuum is circularly polarized with momentum-locked chirality. Focused THz illumination is predicted to excite directional ferron beams and optical routing in ferroelectric devices (Zhou et al., 2022).

A more speculative extension is the antiferron. Using a generalized LGD framework, dynamically stabilized excitations around an inverted polarization configuration were proposed under a high-frequency drive $P_0$3. Time averaging renormalizes the curvature to

$P_0$4

and metastability requires $P_0$5. These negative-energy low-$P_0$6 modes remain theoretical, but they define a Floquet-like branch of ferronics in which polarization channels are created and erased by the drive itself (Galvez-Poblete et al., 12 Aug 2025).

5. Broader ferronics device landscape

Outside ferron-specific THz physics, ferronics has also been used for devices in which ferroic materials provide high-speed tuning, non-volatile patterning, or cross-coupled switching. A prominent RF example is the Ferroelectric Fast Reactive Tuner for 1.3 GHz TESLA-type superconducting cavities at MESA. This FE-FRT uses a BaTiO$P_0$7/SrTiO$P_0$8–Mg ferroelectric element whose permittivity is voltage-tunable with measured bulk response $P_0$9 ns and cavity-integrated response $$\omega^2(\mathbf q)=\omega_0^2+c^2q^2,\qquad \omega_0^2=\frac{A}{m_p},\qquad c^2=\frac{\kappa}{m_p},$$0 ns at 400 MHz. For microphonics compensation, the design provides a tuning range of $$\omega^2(\mathbf q)=\omega_0^2+c^2q^2,\qquad \omega_0^2=\frac{A}{m_p},\qquad c^2=\frac{\kappa}{m_p},$$1 Hz and reduces the maximal forward RF power from $$\omega^2(\mathbf q)=\omega_0^2+c^2q^2,\qquad \omega_0^2=\frac{A}{m_p},\qquad c^2=\frac{\kappa}{m_p},$$2 W without an FE-FRT to $$\omega^2(\mathbf q)=\omega_0^2+c^2q^2,\qquad \omega_0^2=\frac{A}{m_p},\qquad c^2=\frac{\kappa}{m_p},$$3 W with the FE-FRT, a $$\omega^2(\mathbf q)=\omega_0^2+c^2q^2,\qquad \omega_0^2=\frac{A}{m_p},\qquad c^2=\frac{\kappa}{m_p},$$4 reduction (Monroy-Villa et al., 25 Jun 2025).

A distinct electronic use of the term appears in graphene/ferroelectric superlattices. A periodically poled PZT substrate writes a one-dimensional square-wave electrostatic potential into monolayer graphene, with domain periods $$\omega^2(\mathbf q)=\omega_0^2+c^2q^2,\qquad \omega_0^2=\frac{A}{m_p},\qquad c^2=\frac{\kappa}{m_p},$$5–$$\omega^2(\mathbf q)=\omega_0^2+c^2q^2,\qquad \omega_0^2=\frac{A}{m_p},\qquad c^2=\frac{\kappa}{m_p},$$6 nm and KPFM-measured barrier heights $$\omega^2(\mathbf q)=\omega_0^2+c^2q^2,\qquad \omega_0^2=\frac{A}{m_p},\qquad c^2=\frac{\kappa}{m_p},$$7–$$\omega^2(\mathbf q)=\omega_0^2+c^2q^2,\qquad \omega_0^2=\frac{A}{m_p},\qquad c^2=\frac{\kappa}{m_p},$$8 meV. In transport, the device shows a flattened conductance region of width $$\omega^2(\mathbf q)=\omega_0^2+c^2q^2,\qquad \omega_0^2=\frac{A}{m_p},\qquad c^2=\frac{\kappa}{m_p},$$9 meV around the Dirac point, with low-temperature miniband openings and closures near the superlattice Brillouin-zone boundaries. This work presents ferronics as nonvolatile substrate-programmed electronic band engineering (Wan et al., 2021).

Multiferroic and spintronic usages are equally explicit. In a ferromagnet/insulator/ferroelectric nanometer multilayer, field-induced strain couples a Stoner–Wohlfarth ferromagnet to a Landau ferroelectric, enabling multiple nonvolatile magnetization–polarization states. In PZTFWx thin films, modest magnetic fields $P_0$00 T were reported to drive the polarization from $P_0$01 C/m$P_0$02 to zero at room temperature, supporting a three-state logic picture $P_0$03. Antiferromagnetic piezospintronics further places ferronics within strain-controlled antiferromagnetic memory, with room-temperature MnPt/PMN–PT devices switched by $P_0$04 kV/cm and $P_0$05 kV/cm gate fields and a piezo-strain-controlled antiferromagnetic tunnel junction exhibiting $P_0$06 TAMR (Zhang, 2012, Kumar et al., 2010, Liu et al., 2019).

Ferromagnetic semiconductor work extends the same logic to carrier-controlled magnetic order. The first intrinsic n-type electron-induced ferromagnetic semiconductor, (In,Fe)As, established independent tuning of Fe content and electron density in a III–V host, while the later narrow-gap system (In,Fe)Sb reached $P_0$07 K at $P_0$08 and produced a room-temperature anomalous-Hall sensor with $P_0$09 mV·(mT·V)$P_0$10, surpassing the best commercial InSb benchmark cited in that study. In this broader vocabulary, ferronics denotes electronics whose operational state is ferromagnetic order integrated into semiconductor platforms (Hai et al., 2011, Tu et al., 2017, Hasan, 2024).

6. Applications, open problems, and conceptual tensions

Ferronics now supports a wide application map. In the ferron-centered branch, narrowband room-temperature THz sources, THz high-$P_0$11 spectroscopy, ferron-based modulators and detectors, resonant control of lattice and molecular rotations or vibrations, ferronic polariton THz lasers, ultrafast electronics, photonics, quantum interconnects, and next-generation wireless communication are all stated targets. In the broader ferroic-device branch, the applications include microphonics suppression in SRF cavities, reconfigurable graphene miniband electronics, multiferroic RAM and three-state logic, magnetic-field-insensitive antiferromagnetic memory, and Hall sensing based on room-temperature ferromagnetic semiconductors (Zhang et al., 7 Sep 2025, Monroy-Villa et al., 25 Jun 2025, Wan et al., 2021).

Several misconceptions are already identifiable. Ferrons are not merely ordinary optical phonons: the NbOI$P_0$12 experiments reported hypersonic propagation at velocities exceeding $P_0$13 m/s and a ferroelectric-order-dependent second-order THz-emission mechanism, explicitly distinguishing the observed modes from bare TO phonons and from conventional phonon polaritons (Choe et al., 28 May 2025). Nor is ferronics terminologically uniform: some papers restrict it to polarization-carrying collective modes in ferroelectrics, while others use it for ferroic-enabled RF, electronic, or spintronic architectures (Bauer et al., 2023, Monroy-Villa et al., 25 Jun 2025).

The open problems are correspondingly diverse. In ferron spectroscopy and transport, repeatedly stated questions include direct measurement of $P_0$14 beyond the near-$P_0$15 regime, coherent transport and interference conditions, disorder and domain-wall effects on $P_0$16, coupling to excitons, magnons, polaritons, strain, and defects, and the limits of radiation efficiency and absolute power (Zhang et al., 7 Sep 2025, Lopez et al., 29 Jun 2026). Surface ferrons and antiferrons remain primarily theoretical; the former await systematic near-field observation and routing experiments, while the latter require experimentally realistic high-frequency drive protocols and damping analyses (Zhou et al., 2022, Galvez-Poblete et al., 12 Aug 2025). In device-oriented ferronics, fatigue, leakage, field generation, thermal management, domain stability, and large-area reproducibility remain the dominant engineering constraints (Kumar et al., 2010, Liu et al., 2019).

A plausible implication is that ferronics is currently best understood not as a single stabilized subdiscipline, but as a converging set of research programs organized around ferroic order as a controllable dynamical variable. The ferron program has supplied the field with its most specific quasiparticle content and some of its strongest recent experimental results, especially in layered ferroelectrics. The broader device program shows that the same organizing principle also scales naturally into RF tuning, multiferroic memory, programmable electronic structure, and semiconductor spintronics.