🔬 condensed matter

How Does Hard Magnetic Soft Matter Deform?

This paper resolves conflicting models regarding the deformation of hard magnetic soft materials by constructing counterexamples and establishing a first-principles framework that unifies these approaches into a single optimal closed-form law for rational device design.

Original authors: Daniel Katusele, Carmel Majidi, Liping Liu, Pradeep Sharma, Kaushik Dayal

Published 2026-10-05

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Daniel Katusele, Carmel Majidi, Liping Liu, Pradeep Sharma, Kaushik Dayal

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a material that can bend, twist, and fold on command, not because it is powered by wires or batteries, but because it reacts to a magnetic field. This is the promise of hard-magnetic soft materials, a class of substances that combine the flexibility of rubber with the persistent magnetic strength of rare-earth magnets. Inside these soft blocks are tiny, rigid magnetic particles that hold a fixed magnetic direction, much like a compass needle that never forgets which way is north. When an external magnetic field is applied, these embedded particles try to align with it, exerting a twisting force on the surrounding soft rubber. This internal tug-of-war causes the entire material to contort into complex shapes, enabling the creation of untethered robots, flexible medical tools, and adaptive lenses that can change shape without any physical connection to a power source.

For engineers to design these devices, they must be able to predict exactly how the material will move. This requires understanding a specific mechanical puzzle: when the soft rubber stretches or shears, how do the rigid magnetic particles inside rotate? The answer is not obvious. If the particles rotate in one way, the device might bend as intended; if they rotate differently, the device might do nothing at all. For years, scientists have relied on different mathematical rules to guess this rotation, but these rules often disagreed with one another, leading to conflicting predictions about how the material would behave under stress. Some rules suggested the particles would simply follow the overall stretch of the material, while others argued they would only rotate based on the direction of the stretch, ignoring the shape of the particles themselves.

In a new study published in PNAS Nexus, researchers have resolved this long-standing debate by building a new framework from the ground up. Instead of guessing how the particles move, they calculated the most energetically efficient way for a rigid particle to rotate inside a stretching soft matrix. They treated the problem as a search for the path of least resistance, balancing the energy required to stretch the rubber against the energy needed to rotate the particle. By using a method that interpolates between the rigid motion of the particle and the smooth stretching of the surrounding material, they derived a single, clear rule that describes exactly how the particle turns. This new rule accounts for the shape of the particle and its initial angle, factors that previous models had either ignored or handled incorrectly.

The researchers tested their new rule against computer simulations that modeled the material in extreme detail, breaking it down into thousands of tiny points to see exactly how it deformed. They found that their new rule matched these detailed simulations with high accuracy, whereas the older, competing rules often failed. In particular, the study showed that the popular idea that particles only rotate based on the overall direction of the stretch is fundamentally flawed. The researchers demonstrated that if a particle is shaped like an ellipse rather than a circle, and it is tilted at an angle, it will rotate significantly even when the material is simply being pulled in a straight line. The older models, which assumed the particle would not rotate in such a scenario, predicted that the material would remain static, missing the entire point of how these magnetic devices work.

One of the most significant findings is that the conflicting models used by scientists in the past are not actually wrong in every situation; they are just special, limited versions of this new, more general rule. When the magnetic particles are perfectly round, or when the material is stretched in a very specific way, the new rule simplifies down to match the older predictions. However, for the vast majority of real-world shapes and angles, the older rules break down. The new approach reveals that the shape of the particle is just as important as the force applied to it. A long, thin particle will behave very differently from a round one, and the new formula captures this difference precisely. This means that engineers can now design soft robots with much greater confidence, knowing exactly how the internal magnetic particles will orient themselves as the device moves.

The study also highlighted a critical failure in the older models when applied to real-world scenarios. In a test where a magnetic field was applied to a material with tilted, elliptical particles, the older models predicted that the material would not move at all, effectively suggesting that the magnetic force had no effect on the shape. The new model, however, correctly predicted that the material would deform and rotate, aligning with the physical reality that magnetic forces should cause motion. This distinction is vital for the future of soft robotics. If a designer uses a model that predicts no movement, they might build a device that never activates. By using the new, more accurate rule, designers can create devices that reliably convert magnetic energy into complex, controlled motion, opening the door to more sophisticated and capable soft machines for medicine and exploration.

Ultimately, this work provides a unified way to understand how hard-magnetic soft matter deforms. It replaces a collection of conflicting guesses with a single, physics-based law that works across a wide range of shapes and conditions. The researchers showed that by focusing on the energy balance between the stretching rubber and the rotating particle, they could derive a closed-form expression that predicts the motion of the magnetic inclusions. This expression is simple enough to be used in design software yet accurate enough to match detailed computer simulations. The result is a clearer path forward for creating the next generation of wireless, programmable soft devices, ensuring that the materials we build will behave exactly as we intend when the magnetic field is turned on.

Technical Summary: How Does Hard Magnetic Soft Matter Deform?

Problem Statement
Hard magnetic soft materials (HMSMs), consisting of elastomers seeded with high-coercivity magnetic particles, are pivotal for untethered soft robotics and biomedical devices. Their functionality relies on the conversion of magnetic energy into complex mechanical motion. However, a fundamental controversy exists regarding the micromechanism governing how the remanent magnetization of embedded rigid particles transforms under finite deformation.

Current modeling approaches suffer from conflicting predictions:

Material Mapping (e.g., Zhao et al.): Assumes magnetization is convected by the full deformation gradient ( $F$ ). This predicts that magnetization magnitude and direction change with stretching, potentially leading to non-physical predictions where remanent magnetization depends on stretch in the absence of an external field.
Polar-Decomposition Mapping (e.g., Danas & Reis, Yan et al.): Assumes magnetization rotates only according to the rotational component ( $R_{PD}$ ) of the polar decomposition of $F$ , preserving magnitude. While this addresses the stretching issue, the authors demonstrate it fails to capture geometry-induced reorientation. Specifically, for anisotropic (e.g., elliptical) inclusions under pure stretch, this model predicts zero rotation because the polar rotation of pure stretch is the identity, contradicting physical expectations that elliptical particles should rotate to align with the extensional axis.

The central unresolved question is: What is the correct constitutive mapping that relates remanent magnetization to the deformed configuration, preserving magnitude while accurately capturing rotation across the full spectrum of deformations and inclusion geometries?

Methodology
The authors propose a first-principles, energetic micromechanical framework to derive an optimal deformation mapping.

Variational Formulation: The macroscopic energy is posed as a variational problem over admissible fields to construct a computable upper bound. The matrix is modeled as a near-incompressible Neo-Hookean solid, and the inclusion is treated as rigid.
Interpolation Ansatz: A trial deformation field is constructed that interpolates between the imposed homogeneous boundary deformation ( $y = Fx$ ) and the rigid-body motion of the inclusion ( $y = R_0x + c_0$ ). This interpolation occurs in a transition zone defined by a level-set function.
Energy Minimization: The rotation of the inclusion ( $R_0$ ) is determined by minimizing the upper-bound energy with respect to the free parameters (rotation and translation).
Derivation: The framework is specialized to a 2D plane-strain setting with a single rigid elliptical inclusion. The authors derive a closed-form expression for the inclusion rotation ( $\theta_{MIA}$ ) as a function of the macroscopic deformation gradient, the inclusion aspect ratio ( $\eta$ ), and the initial orientation ( $\alpha$ ).
Validation: The derived "Micromechanical Interpolation Approach" (MIA) mapping is compared against:
- Existing models (Material Mapping and Polar-Decomposition Mapping).
- Full-field numerical simulations using finite element analysis (FEniCSx) for single inclusions under simple shear and uniaxial tension.

Key Contributions and Results

Closed-Form Optimal Mapping: The paper derives an explicit law (Eq. 16) for the rotation of a rigid inclusion. Unlike previous models, this mapping depends on the inclusion's aspect ratio and initial orientation, allowing it to capture geometry-induced reorientation even under pure stretch.
Resolution of Model Conflicts: The study demonstrates that the seemingly conflicting prior models are actually special limit cases of the general MIA framework:
- The Polar-Decomposition model is recovered as the limit for circular inclusions ( $\eta \to 1$ ).
- The Material Mapping model is recovered as the limit for long, thin inclusions ( $\eta \to \infty$ ).
Counterexamples to Existing Models: The authors provide counterexamples showing that polar-decomposition-based models fail to predict rotation for inclined anisotropic inclusions under pure stretch, whereas the MIA model correctly predicts stretch- and geometry-dependent rotation.
Numerical Agreement: Under simple shear and tension, the MIA predictions show superior agreement with finite-element simulations compared to the Material Mapping and Polar-Decomposition models, particularly for elliptical inclusions.
Magneto-Mechanical Coupling: The authors illustrate that the failure of the Polar-Decomposition model to predict rotation in anisotropic cases leads to a "vacuous" prediction of zero magneto-mechanical coupling in simple device-scale examples. In contrast, the MIA framework correctly predicts induced deformation driven by magnetic torque, resolving the physical inconsistency of prior models.

Significance
The paper claims to resolve the debate on the fundamental micromechanisms of HMSMs by deriving an analytic, energy-based law that reconciles competing hypotheses. The resulting closed-form expressions provide a unified view of deformation in magnetic soft matter, showing that previous models are simply different limits of a general approach.

The significance of this work lies in its ability to:

Provide a systematic and predictive foundation for the rational design of magnetically actuated soft devices.
Offer a unified framework that accommodates inclusion shape and initial orientation, essential for designing programmable actuation in soft robotics and multifunctional materials.
Correctly predict geometry-induced rotations that are critical for accurate device-scale modeling, avoiding the "phantom device" problem where models predict no motion despite magnetic stimulation.

The authors note that the current derivation is restricted to dilute, single-inclusion systems in 2D plane strain. Future work is suggested to extend this to 3D ellipsoidal inclusions and to incorporate explicit magnetostatic interactions and field-gradient forces for concentrated, device-scale settings.

More like this