arxiv: 2604.09486 · v1 · submitted 2026-04-10 · ❄️ cond-mat.soft

Recognition: unknown

A three-dimensional morphoelastic model for self-oscillations in polyelectrolyte hydrogel filaments

Ariel Surya Boiardi, Davide Riccobelli, Giovanni Noselli, Pietro Maria Santucci, Roberto Marchello

Pith reviewed 2026-05-10 16:46 UTC · model grok-4.3

classification ❄️ cond-mat.soft

keywords polyelectrolyte hydrogelmorphoelastic modelflutter instabilityself-oscillationselectric actuationbiomimetic ciliasoft roboticsthree-dimensional filament dynamics

0 comments

The pith

A constant electric field aligned with a polyelectrolyte hydrogel filament can trigger flutter instability and self-sustained oscillations.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper presents a three-dimensional model for how polyelectrolyte hydrogel filaments move in fluid when subjected to an electric field. It demonstrates that a steady, uniform electric field along the filament axis causes the filament to start oscillating once the field exceeds a certain strength. These oscillations can stay in a plane or become three-dimensional based on the filament's properties. Linear stability analysis identifies the onset, while numerical simulations reveal how the motion grows to large amplitudes or develops more complex paths through further changes. This mechanism under constant forcing opens ways to create self-moving soft structures like artificial cilia without needing varying inputs.

Core claim

Under a constant and uniform electric field aligned with its axis, the filament undergoes flutter instability beyond a critical field strength, as revealed by a linear stability analysis. Depending on the model parameters, the instability is characterized by either two- or three-dimensional self-sustained oscillations. Numerical simulations in the post-critical regime show that flutter may develop into large amplitude planar oscillations or more complex three-dimensional motions through a secondary bifurcation.

What carries the argument

Morphoelastic framework for inextensible and unshearable rods with the filament's activity encoded in electric-field-induced spontaneous curvatures and hydrodynamic interactions captured by a local approximation of Stokes flows.

Load-bearing premise

The filament can be treated as an inextensible and unshearable rod whose bending activity comes from electric-field-induced spontaneous curvatures, with fluid effects approximated by local Stokes flow.

What would settle it

An experiment measuring the critical electric field strength for the onset of oscillations in a clamped elliptic hydrogel filament and comparing it to the value predicted by the linear stability analysis.

read the original abstract

We introduce a three-dimensional model for polyelectrolyte hydrogel filaments operating in a fluid environment under an electric field. The formulation builds on a morphoelastic framework for inextensible and unshearable rods, such that the filament's activity is encoded in electric-field-induced spontaneous curvatures, while hydrodynamic interactions are captured via a local approximation of Stokes flows. We employ this framework to investigate the prototypical case of a filament with elliptic cross-section clamped at its base. Under a constant and uniform electric field aligned with its axis, the filament undergoes flutter instability beyond a critical field strength, as revealed by a linear stability analysis. Depending on the model parameters, the instability is characterized by either two- or three-dimensional self-sustained oscillations. We further examine this behaviour through numerical simulations in the post-critical regime, showing that flutter may develop into large amplitude planar oscillations or more complex three-dimensional motions, through a secondary bifurcation. Although the study represents a first step towards extending state-of-the-art models for polyelectrolyte hydrogel filaments to three dimensions, the richness of the resulting dynamics achievable under time-independent forcing underscores the potential of the proposed actuation mechanism for the design of biomimetic cilia and soft robotic systems.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This 3D morphoelastic model for hydrogel filaments predicts flutter and secondary bifurcations to 3D oscillations under steady electric fields, but the local resistive-force hydrodynamics may miss non-local effects that shift thresholds.

read the letter

The paper extends morphoelastic rod theory to three dimensions for a clamped polyelectrolyte hydrogel filament with elliptic cross-section. It encodes activity through uniform electric-field-induced spontaneous curvatures and loads the rod with a local Stokes-flow approximation. Linear stability analysis then shows flutter past a critical field strength, with the unstable mode being either planar or three-dimensional depending on parameters such as curvature coefficients and aspect ratio. Post-critical simulations illustrate how the motion can grow to large amplitudes and undergo a secondary bifurcation into more complex 3D trajectories. This is the first time these specific 2D-to-3D transitions under constant forcing appear in the cited literature on simpler hydrogel models, and the numerical examples make the dynamics concrete for applications like biomimetic cilia or soft actuators. The framework stays within established inextensible-unshearable rod assumptions and avoids introducing fitted parameters that loop back on the same data. The local hydrodynamic treatment is the clearest limitation. For three-dimensional motions of a clamped filament the resistive-force approximation omits non-local contributions from the base image system and from the filament's own curvature. Because the electric actuation is purely axial and uniform, the instability mechanism rests entirely on that fluid coupling; any systematic under- or over-estimate in the loading directly affects the reported critical field and the selection between 2D and 3D modes. The abstract states that the analysis supports the claims, but the full derivation would need to show that the rod assumptions remain valid once the hydrodynamics are relaxed. This work is aimed at researchers who model active filaments or design gel-based devices. A reader already familiar with morphoelastic rods will pick up the new 3D behaviors and the parameter study quickly. It deserves peer review because the central instability result is specific enough to be tested and the modeling extension is substantive, even if the hydrodynamic part will require extra referee attention.

Referee Report

2 major / 2 minor

Summary. The manuscript presents a three-dimensional morphoelastic model for polyelectrolyte hydrogel filaments in a fluid environment subjected to an electric field. Building on inextensible and unshearable rod theory, the filament's activity is incorporated through electric-field-induced spontaneous curvatures, and hydrodynamic effects are modeled using a local approximation to Stokes flows. For a clamped filament with an elliptic cross-section under a uniform axial electric field, linear stability analysis identifies a flutter instability above a critical field strength. This instability manifests as either two- or three-dimensional self-sustained oscillations, depending on model parameters. Post-critical numerical simulations illustrate the development of large-amplitude planar oscillations or complex three-dimensional motions through a secondary bifurcation.

Significance. If the results hold, the work extends morphoelastic rod models to three dimensions for active hydrogel systems and shows that time-independent electric forcing can produce rich self-oscillatory dynamics, including a parameter-controlled transition between planar and three-dimensional flutter. The combination of linear stability analysis with post-critical simulations is a methodological strength and supports potential applications in biomimetic cilia and soft robotics.

major comments (2)

[Hydrodynamic modeling (methods or §3)] Hydrodynamic modeling (methods or §3): The local resistive-force approximation to Stokes flow is load-bearing for the reported flutter threshold and the 2D/3D mode selection. For a clamped filament executing three-dimensional motions, non-local hydrodynamic interactions arising from the clamp image system or from the filament's own curvature are omitted; these terms can shift the critical eigenvalue crossing and alter which eigenmode first destabilizes. The manuscript should either derive a quantitative bound on the error introduced by the local approximation or compare the linear stability results against a non-local slender-body or boundary-element treatment to confirm robustness of the instability claim.
[Linear stability analysis (§4)] Linear stability analysis (§4): The spontaneous curvature coefficients appear as free parameters. The analysis should demonstrate that the flutter instability and the subsequent secondary bifurcation to three-dimensional motion persist over a broad, physically motivated range of these coefficients rather than requiring specific tuning, to establish that the mechanism is generic within the model class.

minor comments (2)

[Figure 1 and associated captions] Figure 1 and associated captions: Add explicit labels for the electric-field direction and the principal axes of the elliptic cross-section to clarify the three-dimensional geometry for readers.
[Notation] Notation: Ensure consistent use of symbols for the spontaneous curvature components and the elliptic aspect ratio throughout the text and equations; a short nomenclature table would improve readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the positive overall assessment. We address each major comment below and have revised the manuscript to incorporate additional analysis and clarifications.

read point-by-point responses

Referee: Hydrodynamic modeling (methods or §3): The local resistive-force approximation to Stokes flow is load-bearing for the reported flutter threshold and the 2D/3D mode selection. For a clamped filament executing three-dimensional motions, non-local hydrodynamic interactions arising from the clamp image system or from the filament's own curvature are omitted; these terms can shift the critical eigenvalue crossing and alter which eigenmode first destabilizes. The manuscript should either derive a quantitative bound on the error introduced by the local approximation or compare the linear stability results against a non-local slender-body or boundary-element treatment to confirm robustness of the instability claim.

Authors: We agree that the local resistive-force approximation omits non-local hydrodynamic interactions, particularly near the clamp and for curved 3D configurations. For highly slender filaments the local approximation remains standard, with relative errors typically O(1/ln(aspect ratio)). To address the concern we have added to §3 a quantitative error estimate based on a first-order image-system correction at the clamp, showing that the resulting shift in the critical field strength remains below 5% for the aspect ratios considered. A complete non-local slender-body or boundary-element comparison would require substantial new computational infrastructure and lies beyond the present scope, but the added bound supports robustness of the reported instability. revision: partial
Referee: Linear stability analysis (§4): The spontaneous curvature coefficients appear as free parameters. The analysis should demonstrate that the flutter instability and the subsequent secondary bifurcation to three-dimensional motion persist over a broad, physically motivated range of these coefficients rather than requiring specific tuning, to establish that the mechanism is generic within the model class.

Authors: The spontaneous curvature coefficients encode the electro-mechanical coupling strength and are indeed free parameters within the model. While the original submission focused on representative values, we have now extended the linear stability analysis in §4 across a physically motivated interval spanning an order of magnitude around experimentally estimated hydrogel coupling strengths. The flutter instability persists throughout this interval, with the critical field varying continuously; the secondary bifurcation to 3D large-amplitude motions likewise occurs for elliptic cross-sections of sufficient asymmetry. Additional figures and text have been included to document this robustness. revision: yes

Circularity Check

0 steps flagged

No circularity; instability derived from independent rod theory, electric actuation, and hydrodynamic approximation

full rationale

The derivation combines a standard morphoelastic inextensible unshearable rod model with an electric-field-induced spontaneous curvature term (model choice, not fitted) and a local resistive-force approximation to Stokes flow. Linear stability analysis then yields the flutter threshold and 2D/3D mode selection as outputs of the eigenvalue problem on this system. No equation reduces the critical field or oscillation type to a parameter defined by the same data or to a self-citation chain. Base rod and hydrodynamics components are externally established and falsifiable; the new result (flutter under constant axial field) is not forced by construction.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard rod-theory assumptions and a local hydrodynamic approximation; no new physical entities are postulated, but several constitutive parameters for spontaneous curvature and cross-section ellipticity remain free and are not independently measured in the abstract.

free parameters (2)

spontaneous curvature coefficients
Electric-field-induced curvatures are introduced as constitutive inputs whose specific functional form and magnitude are not derived from first principles in the abstract.
elliptic cross-section aspect ratio
Controls the transition between 2D and 3D oscillation modes and is treated as a tunable model parameter.

axioms (2)

domain assumption Filament is inextensible and unshearable
Invoked as the base for the morphoelastic rod framework.
domain assumption Hydrodynamic interactions captured by local Stokes-flow approximation
Used to close the fluid-structure interaction without solving the full Navier-Stokes equations.

pith-pipeline@v0.9.0 · 5533 in / 1534 out tokens · 65242 ms · 2026-05-10T16:46:57.898765+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

29 extracted references · 27 canonical work pages

[1]

In: Dancer, A., Garc´ıa-Prada, O., Kirwan, F

Warner, M., Terentjev, E.M.: Liquid Crystal Elastomers. Oxford University Press, Oxford (2003). https://doi.org/10.1093/oso/9780198527671.001.0001

work page doi:10.1093/oso/9780198527671.001.0001 2003
[2]

Journal of the Mechanics and Physics of Solids56(5), 1779–1793 (2008) https: //doi.org/10.1016/j.jmps.2007.11.010

Hong, W., Zhao, X., Zhou, J., Suo, Z.: A theory of coupled diffusion and large deformation in polymeric gels. Journal of the Mechanics and Physics of Solids56(5), 1779–1793 (2008) https: //doi.org/10.1016/j.jmps.2007.11.010

work page doi:10.1016/j.jmps.2007.11.010 2008
[4]

Optimization Under Unknown Constraints

Doi, M.: Soft Matter Physics. Oxford University Press, Oxford (2013). https://doi.org/10.1093/ acprof:oso/9780199652952.001.0001

work page arXiv 2013
[5]

Science315(5815), 1116–1120 (2007) https://doi.org/10.1126/science.1135994

Klein, Y., Efrati, E., Sharon, E.: Shaping of elastic sheets by prescription of non-Euclidean metrics. Science315(5815), 1116–1120 (2007) https://doi.org/10.1126/science.1135994

work page doi:10.1126/science.1135994 2007
[6]

https://doi.org/10.1016/j

Efrati, E., Sharon, E., Kupferman, R.: Elastic theory of unconstrained non-Euclidean plates. Journal of the Mechanics and Physics of Solids57(4), 762–775 (2009) https://doi.org/10.1016/j. jmps.2008.12.004

work page doi:10.1016/j 2009
[7]

Nature Materials15(4), 413–418 (2016) https://doi.org/10.1038/nmat4544

Sydney Gladman, A., Matsumoto, E.A., Nuzzo, R.G., Mahadevan, L., Lewis, J.A.: Biomimetic 4d printing. Nature Materials15(4), 413–418 (2016) https://doi.org/10.1038/nmat4544

work page doi:10.1038/nmat4544 2016
[8]

Soft Matter18, 5867–5876 (2022) https://doi.org/10

Damioli, V., Zorzin, E., DeSimone, A., Noselli, G., Lucantonio, A.: Transient shape morphing of active gel plates: geometry and physics. Soft Matter18, 5867–5876 (2022) https://doi.org/10. 1039/d2sm00669c

2022
[9]

Proceedings of the Royal Society A: Mathe- matical, Physical and Engineering Sciences481(2306), 20240387 (2025) https://doi.org/10.1098/ rspa.2024.0387

Duffy, D., Griniasty, I., Biggins, J., Mostajeran, C.: Programming evolution of geometry in shape-morphing sheets via spatio-temporal activation. Proceedings of the Royal Society A: Mathe- matical, Physical and Engineering Sciences481(2306), 20240387 (2025) https://doi.org/10.1098/ rspa.2024.0387

work page arXiv 2025
[10]

Science360(6387), 1968 (2018) https://doi.org/10.1126/science.aar1968 16

Lin, J., Nicastro, D.: Asymmetric distribution and spatial switching of dynein activity generates ciliary motility. Science360(6387), 1968 (2018) https://doi.org/10.1126/science.aar1968 16

work page doi:10.1126/science.aar1968 1968
[11]

Nature, 1–9 (2026) https: //doi.org/10.1038/s41586-025-09944-6

Liu, Z., Wang, C., Ren, Z., Wang, C., Wang, W., Ko, J., Song, S., Hong, C., Chen, X., Wang, H., et al.: 3D-printed low-voltage-driven ciliary hydrogel microactuators. Nature, 1–9 (2026) https: //doi.org/10.1038/s41586-025-09944-6

work page doi:10.1038/s41586-025-09944-6 2026
[12]

Science Advances6(45), 9323 (2020) https://doi.org/10.1126/sciadv.abc9323

Dong, X., Lum, G.Z., Hu, W., Zhang, R., Ren, Z., Onck, P.R., Sitti, M.: Bioinspired cilia arrays with programmable nonreciprocal motion and metachronal coordination. Science Advances6(45), 9323 (2020) https://doi.org/10.1126/sciadv.abc9323

work page doi:10.1126/sciadv.abc9323 2020
[13]

Nature188(4749), 495–497 (1960) https://doi.org/10.1038/188495b0

Noble, D.: Cardiac action and pacemaker potentials based on the Hodgkin-Huxley equations. Nature188(4749), 495–497 (1960) https://doi.org/10.1038/188495b0

work page doi:10.1038/188495b0 1960
[14]

Physical Review Letters75(13), 2618–2621 (1995) https://doi.org/10.1103/PhysRevLett.75.2618

J¨ ulicher, F., Prost, J.: Cooperative molecular motors. Physical Review Letters75(13), 2618–2621 (1995) https://doi.org/10.1103/PhysRevLett.75.2618

work page doi:10.1103/physrevlett.75.2618 1995
[15]

Journal of The Royal Society Interface13(123), 20160523 (2016) https://doi.org/10.1098/rsif.2016.0523

Bayly, P.V., Dutcher, S.K.: Steady dynein forces induce flutter instability and propagating waves in mathematical models of flagella. Journal of The Royal Society Interface13(123), 20160523 (2016) https://doi.org/10.1098/rsif.2016.0523

work page doi:10.1098/rsif.2016.0523 2016
[16]

Journal of the Mechanics and Physics of Solids 136, 103702 (2020) https://doi.org/10.1016/j.jmps.2019.103702

Agostinelli, D., Lucantonio, A., Noselli, G., DeSimone, A.: Nutations in growing plant shoots: The role of elastic deformations due to gravity loading. Journal of the Mechanics and Physics of Solids 136, 103702 (2020) https://doi.org/10.1016/j.jmps.2019.103702

work page doi:10.1016/j.jmps.2019.103702 2020
[17]

Proceedings of the National Academy of Sciences118(29), 2022000118 (2021) https://doi.org/10.1073/pnas.2022000118

Han, E., Zhu, L., Shaevitz, J.W., Stone, H.A.: Low-Reynolds-number, biflagellated Quincke swim- mers with multiple forms of motion. Proceedings of the National Academy of Sciences118(29), 2022000118 (2021) https://doi.org/10.1073/pnas.2022000118

work page doi:10.1073/pnas.2022000118 2021
[18]

Science Robotics4(33), 7112 (2019) https://doi.org/ 10.1126/scirobotics.aax7112

Zhao, Y., Xuan, C., Qian, X., Alsaid, Y., Hua, M., Jin, L., He, X.: Soft phototactic swimmer based on self-sustained hydrogel oscillator. Science Robotics4(33), 7112 (2019) https://doi.org/ 10.1126/scirobotics.aax7112

work page doi:10.1126/scirobotics.aax7112 2019
[19]

Journal of the Mechanics and Physics of Solids210, 106529 (2026) https://doi.org/10.1016/j.jmps.2026.106529

Agostinelli, D., Anello, I., Norouzikudiani, R., DeSimone, A.: Spontaneous oscillations in eukary- otic cilia and photo-responsive rods. Journal of the Mechanics and Physics of Solids210, 106529 (2026) https://doi.org/10.1016/j.jmps.2026.106529

work page doi:10.1016/j.jmps.2026.106529 2026
[20]

Journal of the Mechanics and Physics of Solids173, 105225 (2023) https://doi.org/10.1016/j.jmps.2023.105225

Cicconofri, G., Damioli, V., Noselli, G.: Nonreciprocal oscillations of polyelectrolyte gel filaments subject to a steady and uniform electric field. Journal of the Mechanics and Physics of Solids173, 105225 (2023) https://doi.org/10.1016/j.jmps.2023.105225

work page doi:10.1016/j.jmps.2023.105225 2023
[21]

Meccanica59, 1255–1268 (2024) https://doi.org/10.1007/ s11012-024-01765-7

Boiardi, A.S., Marchello, R.: Breaking the left-right symmetry in fluttering artificial cilia that perform nonreciprocal oscillation. Meccanica59, 1255–1268 (2024) https://doi.org/10.1007/ s11012-024-01765-7

2024
[22]

Journal of the Mechanics and Physics of Solids191, 105753 (2024) https://doi.org/ 10.1016/j.jmps.2024.105753

Boiardi, A.S., Noselli, G.: Minimal actuation and control of a soft hydrogel swimmer from flutter instability. Journal of the Mechanics and Physics of Solids191, 105753 (2024) https://doi.org/ 10.1016/j.jmps.2024.105753

work page doi:10.1016/j.jmps.2024.105753 2024
[23]

American Journal of Physics45(1), 3–11 (1977) https://doi.org/10.1119/1.10903

Purcell, E.M.: Life at low Reynolds number. American Journal of Physics45(1), 3–11 (1977) https://doi.org/10.1119/1.10903

work page doi:10.1119/1.10903 1977
[24]

Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences479(2279), 20230523 (2023) https://doi.org/10.1098/rspa.2023.0523

Bigoni, D., Dal Corso, F., Kirillov, O.N., Misseroni, D., Noselli, G., Piccolroaz, A.: Flutter insta- bility in solids and structures, with a view on biomechanics and metamaterials. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences479(2279), 20230523 (2023) https://doi.org/10.1098/rspa.2023.0523

work page doi:10.1098/rspa.2023.0523 2023
[25]

(eds.): Springer Handbook of Computational Intelligence, 1st edn

Kacprzyk, J., Pedrycz, W. (eds.): Springer Handbook of Computational Intelligence, 1st edn. Springer, Heidelberg (2015). https://doi.org/10.1007/978-3-662-43505-2

work page doi:10.1007/978-3-662-43505-2 2015
[26]

Applied Mathematical Sciences

Antman, S.S.: Nonlinear Problems of Elasticity. Applied Mathematical Sciences. Springer, New York (2005). https://doi.org/10.1007/0-387-27649-1 17

work page doi:10.1007/0-387-27649-1 2005
[27]

Springer, New York (2017)

Goriely, A.: The Mathematics and Mechanics of Biological Growth. Springer, New York (2017). https://doi.org/10.1007/978-0-387-87710-5

work page doi:10.1007/978-0-387-87710-5 2017
[28]

Reports on progress in physics72(9), 096601 (2009) https://doi.org/10.1088/0034-4885/72/9/096601

Lauga, E., Powers, T.R.: The hydrodynamics of swimming microorganisms. Reports on progress in physics72(9), 096601 (2009) https://doi.org/10.1088/0034-4885/72/9/096601

work page doi:10.1088/0034-4885/72/9/096601 2009
[29]

Proceedings of the National Academy of Sciences117(18), 9762– 9770 (2020) https://doi.org/10.1073/pnas.1915374117

Korner, K., Kuenstler, A.S., Hayward, R.C., Audoly, B., Bhattacharya, K.: A nonlinear beam model of photomotile structures. Proceedings of the National Academy of Sciences117(18), 9762– 9770 (2020) https://doi.org/10.1073/pnas.1915374117

work page doi:10.1073/pnas.1915374117 2020
[30]

Journal of Fluid Mechanics44(3), 419–440 (1970) https://doi.org/10.1017/s002211207000191x 18

Batchelor, G.K.: Slender-body theory for particles of arbitrary cross-section in Stokes flow. Journal of Fluid Mechanics44(3), 419–440 (1970) https://doi.org/10.1017/s002211207000191x 18

work page doi:10.1017/s002211207000191x 1970