REVIEW 3 major objections 5 minor 33 references
Optimizing Metachronal Paddling with Reinforcement Learning at Low Reynolds Number
T0 review · 3 major / 5 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read The paper claims that for a rigid-paddle swimmer at zero Reynolds number, the back-to-front (antiplectic) wave is the most efficient learned stroke for two, three, and four paddle pairs, even though the fastest stroke changes with paddle…
desk verdict A clean RL study that recovers antiplectic metachrony as efficient, but the efficiency claim rests on a drag coefficient measured for one paddle posture and used for all strokes. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is a two-dimensional 'paddler': a long body with pairs of straight, inflexible paddles, each with one angular degree of freedom discretized into 11 states and a stepwise action moving it left, right, or not at all. The learning loop is tabular Q-learning with the update rule $Q(s,a)\leftarrow (1-\alpha)Q(s,a)+\alpha[r+\gamma\max_{a'}Q(s',a')]$, and the reward for each step is the swimmer's net forward displacement over that time interval. Hydrodynamic forces and swimming speed are computed by solving the Stokes equations with the method of regularized Stokeslets, subject to zero net force on the swimmer. Stroke performance is then compared with the efficiency ratio $\eta = \zeta U^2/P$, using a drag coefficient measured by towing the swimmer with paddles held perpendicular to the body. Since the limbs cannot bend, any swimming motion must come from rhythm alone, which is what lets the paper attribute the learned gaits purely to coordination timing.
What would settle it
Take two learned strokes at the same spacing and paddle count, measure the time-averaged drag force during each stroke by integrating the computed forces over one cycle, and recompute each efficiency with the measured drag value. If the front-to-back stroke's efficiency then reaches or exceeds the back-to-front stroke's, the paper's central conclusion fails for that configuration.
Extended reading notes
Core claim
Across the strokes the agent learns, the back-to-front wave is the most efficient coordination pattern for two, three, and four paddle pairs, while the fastest stroke is not universal: with two and four paddles a front-to-back, inward-tilted stroke is faster at wide spacings, and with three paddles the back-to-front stroke is fastest. The gaits are learned by Q-learning maximizing each move's net forward displacement, so efficiency is not part of the reward; it is measured afterward using the standard ratio $\eta = \zeta U^2/P$, where $U$ is the average swimming speed, $P$ the average mechanical power, and $\zeta$ the drag coefficient. Because the paddles are rigid and each has a single angular degree of freedom, the only time asymmetry available to the agent is the ordering and timing of paddle motions. The paper takes this result as evidence that the antiplectic rhythm seen in nature can be an efficiency-optimal solution even when the swimmer is free to choose any coordination.
Load-bearing premise
The paper assumes that the water resistance measured when the paddles are held straight out is close to the resistance the swimmer actually feels during each learned stroke; if that is wrong, the efficiency ranking could change.
Editorial extensions
If this is right
- In this model, an antiplectic back-to-front wave is the efficiency-optimal coordination for two, three, and four paddle pairs, so the biologically common rhythm can arise from a pure efficiency objective without any prescribed wave rule.
- Because the speed-optimal and efficiency-optimal strokes differ in most configurations, a microrobot designer must choose explicitly between maximizing distance per time and maximizing distance per energy; the same gait will not do both.
- The learned gait switches from back-to-front to front-to-back near a paddle spacing of about 2 body units for all paddle counts, identifying appendage spacing as a control parameter for gait selection.
- The back-to-front stroke is both fastest and most efficient at tight spacings, which is consistent with closely packed cilia arrays in biology and suggests a simple design rule: pack appendages tightly for efficient antiplectic swimming.
Reading between the lines
- The efficiency comparison could be sharpened by measuring each stroke's own time-averaged drag instead of using the single towed value; this is a direct, testable check of whether the efficiency ranking survives the paddles' actual tilted orientations.
- An agent trained to maximize efficiency directly, rather than speed, might settle on a different gait; the paper notes this requires changing the action space or reward, so an actor-critic or deep-Q variant is the natural next experiment.
- Because Stokes flow has no intrinsic length scale, the spacing threshold near 2 units should scale with paddle length and body geometry; varying those dimensions in the same model would test whether the threshold is a ratio rather than an absolute distance.
- The fact that speed and efficiency diverge most at wide spacings suggests that organisms with widely spaced limbs may face a sharper speed-efficiency trade-off than those with tight arrays; comparing learned gaits across more paddle counts would map this trade-off.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper uses tabular Q-learning to discover optimal paddling gaits for a two-dimensional swimmer at zero Reynolds number. The swimmer consists of an elongated body with two, three, or four pairs of rigid paddles, and the paddle angles are discretized into eleven states. The reward is the net forward displacement per action, so the learning agent optimizes swimming speed. For tight paddle spacings, the agent learns a back-to-front (antiplectic) metachronal stroke, while for wide spacings it learns front-to-back strokes in which paddles tilt inward. The speed-optimal stroke depends on paddle number, but the authors report that the back-to-front stroke is the most efficient across two, three, and four paddle pairs. The hydrodynamic problem is solved with regularized Stokeslets, and efficiency is computed with a Lighthill-type metric using a drag coefficient obtained by towing the swimmer with perpendicular paddles. A parameter study in Section 3.5 examines sensitivity to training length, discount factor, exploration rate, and learning rate.
Significance. If the results are robust, the paper provides a useful demonstration that model-free reinforcement learning can recover biologically observed antiplectic metachrony in a minimal rigid-paddle model, and it identifies a clear separation between speed-optimal and efficiency-optimal coordination. The qualitative comparison with the biological literature on ciliates and other metachronal swimmers is a strength, and the parameter study in Section 3.5 is a genuine robustness check. The hydrodynamic model and Q-learning setup are standard, and the finding that efficiency, rather than speed, selects the antiplectic wave across different paddle numbers is potentially valuable for designing low-Reynolds-number robotic swimmers.
major comments (3)
- [Section 3.4, Appendix B (Table B1)] The efficiency metric η = ζU²/P uses a single drag coefficient ζ per paddle number and spacing, obtained by towing the swimmer with all paddles fixed perpendicular to the body (θ=0). The learned strokes are not small perturbations of this reference: Table C2 reports mean paddle states of about ±2 (≈18°) for the back-to-front strokes and values up to ±3.72 (≈33°) for the front-to-back strokes. Because the drag of an inclined flat appendage in Stokes flow depends on orientation, the perpendicular-paddle ζ may not be representative of either stroke. The bias likely runs in the direction of overestimating the efficiency of the more tilted front-to-back strokes, which would make the main conclusion conservative, but the paper does not state this or test it. I recommend adding a sensitivity analysis that recomputes η using ζ estimated from a stroke-specific reference configuration (for example, the time-averaged or mean-tilt configuration) and reporting whether the efficiency ordering changes.
- [Section 3.5.1, Figure 9] The parameter study reveals a strong sensitivity to the discount factor. At paddle spacing 2 with the longest training loop, γ=0.99 converges to the front-to-back stroke in every trial, but γ=0.98 converges to a suboptimal back-to-front stroke in every trial, and γ=0.97 produces ineffective swimming in all trials. This indicates that the learned strokes are not robustly optimal with respect to small hyperparameter changes. Since the abstract and discussion repeatedly refer to 'optimal' strokes, the paper should either provide a formal optimality check (for example, policy iteration on the converged Q-table, multiple random seeds with the final chosen parameters, or a systematic grid over γ) or consistently soften the terminology to 'learned' or 'RL-selected' strokes. As written, the claim of optimality is stronger than the evidence.
- [Section 3.4, Figure 8] The headline efficiency claim rests on comparisons of point estimates from only five runs per spacing. The text states that for the two- and four-paddle swimmers the maximum front-to-back efficiency is 'comparable' to the minimum back-to-front efficiency, which suggests a small margin that could be comparable to run-to-run variation. No confidence intervals, standard deviations, or effect sizes are reported for η, and Figure 8 does not include error bars. Given that the central conclusion is an efficiency ranking, the authors should report the quantitative margins (for example, the maximum and minimum η values for each stroke type) and assess whether the ranking is robust to the observed variation across the five runs.
minor comments (5)
- [Section 2.2] The phrase '968 states' should be '968 state-action pairs': the calculation is 121 paddle configurations times 8 actions, which is the number of state-action pairs, not the number of states.
- [Appendix C, Table C2 caption] The word 'mantitude' should be 'magnitude' in the caption of Table C2.
- [Figure 9] The caption does not clearly explain the row and column layout of the histograms or what the numbers above the bars represent beyond 'mean swimming speed'; adding a sentence that defines the training-loop combinations for each panel would improve readability.
- [Section 3.5.1] The term 'suboptimal' is used without defining the comparison baseline; since the strokes are outputs of the same learning algorithm, the authors should specify that 'suboptimal' means slower than the best stroke found with other hyperparameter choices.
- [Data Availability] The code for the Q-learning algorithm and the regularized Stokeslet solver is not provided, which makes it difficult to reproduce the stochastic RL results; making the code and random seeds available would strengthen reproducibility.
Circularity Check
No circularity: the RL objective is displacement, efficiency is evaluated post hoc, and no load-bearing step reduces to a fitted input or self-citation.
full rationale
The paper's central claim—that a back-to-front metachronal stroke is the most efficient while the fastest stroke depends on paddle number—is not circular. The reinforcement learning reward is the net displacement of the swimmer (Eq. 3), so the learning objective is swimming speed; efficiency is computed afterward using the Lighthill metric η = ζU²/P. The drag coefficient ζ in Appendix B is a model parameter obtained by towing the same discretized swimmer with paddles fixed perpendicular to the body; it is not fitted to the efficiency results, and it is not defined in terms of the stroke being evaluated. The efficiency ranking therefore has independent content: it compares strokes that were selected for speed, not for efficiency. There are no self-citations by Bailey and Guy in the reference list, and the cited external results about antiplectic metachrony are used as context, not as premises that force the numerical outcome. The reader's concern that the fixed perpendicular-paddle drag coefficient may not represent the time-averaged drag of tilted learned strokes is a legitimate modeling-robustness question, but it is a correctness risk rather than circularity: no quantity in the derivation is equivalent by construction to the quantity it is supposed to predict. Accordingly, no circular step is identified.
Assumptions & free parameters
free parameters (7)
- regularization parameter epsilon =
0.05
- Stokeslet discretization spacing =
0.1
- paddle angle discretization =
11 states, pi/20 increments over [-pi/4, pi/4]
- time per action =
1 unit
- discount factor gamma =
0.99 (2 paddle), 0.999 (3,4 paddle)
- number of episodes and learning steps =
50 episodes x 50,000 steps (2); 500 x 500,000 (3,4)
- drag coefficient for efficiency =
listed per spacing in Table B1
assumptions (6)
- domain assumption The fluid is governed by the incompressible Stokes equations at zero Reynolds number, solved with regularized Stokeslets in free space.
- standard math Reciprocal motions produce zero net displacement (scallop theorem), so propulsion requires time-asymmetric coordination.
- domain assumption Top and bottom paddle pairs beat symmetrically, so the body does not rotate and motion is purely horizontal.
- domain assumption Tabular Q-learning with the chosen discretization, discount factor, and training length converges to the optimal policy for the discounted reward objective.
- domain assumption The drag coefficient measured with paddles perpendicular to the body is representative of the drag during the learned strokes.
- domain assumption The 2D model with rigid paddles captures the essential physics of metachronal paddling relevant to the conclusions.
Cite this review
Pith. "Pith review of Optimizing Metachronal Paddling with Reinforcement Learning at Low Reynolds Number." pith.science (2026). https://pith.science/paper/STB2WWS2
@misc{pith2026250718849,
author = {Pith},
title = {Pith review of: Optimizing Metachronal Paddling with Reinforcement Learning at Low Reynolds Number},
year = {2026},
howpublished = {\url{https://pith.science/paper/STB2WWS2}},
note = {Machine review of arXiv:2507.18849}
}
read the original abstract
Metachronal paddling is a swimming strategy in which an organism oscillates sets of adjacent limbs with a constant phase lag, propagating a metachronal wave through its limbs and propelling it forward. This limb coordination strategy is utilized by swimmers across a wide range of Reynolds numbers, which suggests that this metachronal rhythm was selected for its optimality of swimming performance. In this study, we apply reinforcement learning to a swimmer at zero Reynolds number and investigate whether the learning algorithm selects this metachronal rhythm, or if other coordination patterns emerge. We design the swimmer agent with an elongated body and pairs of straight, inflexible paddles placed along the body for various fixed paddle spacings. Based on paddle spacing, the swimmer agent learns qualitatively different coordination patterns. At tight spacings, a back-to-front metachronal wave-like stroke emerges which resembles the commonly observed biological rhythm, but at wide spacings, different limb coordinations are selected. Across all resulting strokes, the fastest stroke is dependent on the number of paddles, however, the most efficient stroke is a back-to-front wave-like stroke regardless of the number of paddles.
Reference graph
Works this paper leans on
-
[1]
Purcell, E.M.: Life at low reynolds number. In: Physics and Our World: Reissue of the Proceedings of a Symposium in Honor of Vic- tor F Weisskopf, pp. 47–67 (2014). World Scientific
work page 2014
-
[2]
Reports on progress in physics 78(5), 056601 (2015)
Elgeti, J., Winkler, R.G., Gompper, G.: Physics of microswimmers—single particle motion and collective behavior: a review. Reports on progress in physics 78(5), 056601 (2015)
work page 2015
-
[3]
: Metachronal motion across scales: current challenges and future directions
Byron, M.L., Murphy, D.W., Katija, K., Hoover, A.P., Daniels, J., Garayev, K., Tak- agi, D., Kanso, E., Gemmell, B.J., Ruszczyk, M., et al. : Metachronal motion across scales: current challenges and future directions. Inte- grative and comparative biology 61(5), 1674– 1688 (2021)
work page 2021
-
[4]
Physical Review E 95(5), 052412 (2017)
Ghorbani, A., Najafi, A.: Symplectic and antiplectic waves in an array of beating cilia attached to a closed body. Physical Review E 95(5), 052412 (2017)
work page 2017
-
[5]
Journal of Cell Science 3(32), 503–521 (1954)
Knight-Jones, E.: Relations between metachronism and the direction of ciliary beat in metazoa. Journal of Cell Science 3(32), 503–521 (1954)
work page 1954
-
[6]
Physics of Fluids 26(9) (2014)
Guo, H., Nawroth, J., Ding, Y., Kanso, E.: Cilia beating patterns are not hydrodynami- cally optimal. Physics of Fluids 26(9) (2014)
work page 2014
-
[7]
Marine biology 158, 2541–2554 (2011)
Murphy, D., Webster, D., Kawaguchi, S., King, R., Yen, J.: Metachronal swimming in antarctic krill: gait kinematics and sys- tem design. Marine biology 158, 2541–2554 (2011)
work page 2011
-
[8]
Integrative and Comparative Biology 61(5), 1631–1643 (2021)
Garayev, K., Murphy, D.W.: Metachronal swimming of mantis shrimp: kinematics and interpleopod vortex interactions. Integrative and Comparative Biology 61(5), 1631–1643 (2021)
work page 2021
Show all 33 references
-
[9]
Scientific reports 6(1), 30713 (2016)
Kim, S., Lee, S., Lee, J., Nelson, B.J., Zhang, L., Choi, H.: Fabrication and manipulation of ciliary microrobots with non-reciprocal mag- netic actuation. Scientific reports 6(1), 30713 (2016)
2016
-
[10]
Journal of Intelli- gent & Robotic Systems 63, 399–416 (2011)
Ghanbari, A., Bahrami, M.: A novel swim- ming microrobot based on artificial cilia for biomedical applications. Journal of Intelli- gent & Robotic Systems 63, 399–416 (2011)
2011
-
[11]
Ocean Engineering 297, 116963 (2024)
Chen, G., Xu, Y., Yang, X., Hu, H., Cheng, H., Zhu, L., Zhang, J., Shi, J., Chai, X.: Target tracking control of a bionic mantis shrimp robot with closed-loop central pattern generators. Ocean Engineering 297, 116963 (2024)
2024
-
[12]
Bioinspiration & Biomimetics 18(1), 16 011003 (2022)
Lim, S., Du, Y., Lee, Y., Panda, S.K., Tong, D., Jawed, M.K.: Fabrication, control, and modeling of robots inspired by flagella and cilia. Bioinspiration & Biomimetics 18(1), 16 011003 (2022)
2022
-
[13]
Pro- ceedings of the National Academy of Sciences 108(38), 15727–15732 (2011)
Osterman, N., Vilfan, A.: Finding the ciliary beating pattern with optimal efficiency. Pro- ceedings of the National Academy of Sciences 108(38), 15727–15732 (2011)
2011
-
[14]
Physical Review E 92(2), 023020 (2015)
Takagi, D.: Swimming with stiff legs at low reynolds number. Physical Review E 92(2), 023020 (2015)
2015
-
[15]
Computers & Fluids 225, 104973 (2021)
Garnier, P., Viquerat, J., Rabault, J., Larcher, A., Kuhnle, A., Hachem, E.: A review on deep reinforcement learning for fluid mechanics. Computers & Fluids 225, 104973 (2021)
2021
-
[16]
Physical Review Fluids 5(7), 074101 (2020)
Tsang, A.C.H., Tong, P.W., Nallan, S., Pak, O.S.: Self-learning how to swim at low reynolds number. Physical Review Fluids 5(7), 074101 (2020)
2020
-
[17]
Physics of Fluids 35(3) (2023)
Qin, K., Zou, Z., Zhu, L., Pak, O.S.: Rein- forcement learning of a multi-link swimmer at low reynolds numbers. Physics of Fluids 35(3) (2023)
2023
-
[18]
Physical Review Fluids 6(5), 050505 (2021)
Jiao, Y., Ling, F., Heydari, S., Heess, N., Merel, J., Kanso, E.: Learning to swim in potential flow. Physical Review Fluids 6(5), 050505 (2021)
2021
-
[19]
Sutton, R.S., Barto, A.G.: Reinforcement learning: An introduction (2018)
2018
-
[20]
Robotics and Autonomous Systems 175, 104659 (2024)
Jebellat, I., Jebellat, E., Amiri-Margavi, A., Vahidi-Moghaddam, A., Pishkenari, H.N.: A reinforcement learning approach to find optimal propulsion strategy for microrobots swimming at low reynolds number. Robotics and Autonomous Systems 175, 104659 (2024)
2024
-
[21]
Advances in Neural Information Processing Systems 36 (2024)
Zhang, S., Li, H., Wang, M., Liu, M., Chen, P.-Y., Lu, S., Liu, S., Murugesan, K., Chaudhury, S.: On the convergence and sam- ple complexity analysis of deep q-networks with epsilon-greedy exploration. Advances in Neural Information Processing Systems 36 (2024)
2024
-
[22]
SIAM Journal on Scientific Com- puting 23(4), 1204–1225 (2001)
Cortez, R.: The method of regularized stokeslets. SIAM Journal on Scientific Com- puting 23(4), 1204–1225 (2001)
2001
-
[23]
SIAM, Philadelphia (1975)
Lighthill, J.: Mathematical Biofluiddynamics. SIAM, Philadelphia (1975)
1975
-
[24]
Communications Physics 5(1), 158 (2022)
Zou, Z., Liu, Y., Young, Y.-N., Pak, O.S., Tsang, A.C.: Gait switching and targeted navigation of microswimmers via deep rein- forcement learning. Communications Physics 5(1), 158 (2022)
2022
-
[25]
arXiv preprint arXiv:2209.10935 (2022)
Zhang, J., Zhou, L., Cao, B.: Learning swim- ming via deep reinforcement learning. arXiv preprint arXiv:2209.10935 (2022)
2022 arXiv
-
[26]
arXiv preprint arXiv:2506.00084 (2025)
Lai, Y., Heydari, S., Pak, O.S., Man, Y.: Navigation of a three-link microswimmer via deep reinforcement learning. arXiv preprint arXiv:2506.00084 (2025)
2025 arXiv
-
[27]
Proceedings of the National Academy of Sciences 117(48), 30201–30207 (2020)
Omori, T., Ito, H., Ishikawa, T.: Swim- ming microorganisms acquire optimal effi- ciency with multiple cilia. Proceedings of the National Academy of Sciences 117(48), 30201–30207 (2020)
2020
-
[28]
Journal of Fluid Mechan- ics 55(1), 1–23 (1972)
Blake, J.: A model for the micro-structure in ciliated organisms. Journal of Fluid Mechan- ics 55(1), 1–23 (1972)
1972
-
[29]
arXiv preprint arXiv:1701.07274 (2017)
Li, Y.: Deep reinforcement learning: An overview. arXiv preprint arXiv:1701.07274 (2017)
2017 arXiv
-
[30]
arXiv preprint arXiv:2409.08815 (2024)
Chen, Y., Yang, Y.: Deep reinforcement learning for tracking a moving target in jellyfish-like swimming. arXiv preprint arXiv:2409.08815 (2024)
2024 arXiv
-
[31]
Proceedings of the National Academy of Sciences 115(23), 5849–5854 (2018)
Verma, S., Novati, G., Koumoutsakos, P.: Efficient collective swimming by harnessing vortices through deep reinforcement learn- ing. Proceedings of the National Academy of Sciences 115(23), 5849–5854 (2018)
2018
-
[32]
arXiv preprint arXiv:2405.11457 (2024) 17
Jiao, Y., Ling, F., Heydari, S., Heess, N., Merel, J., Kanso, E.: Deep dive into model- free reinforcement learning for biological and robotic systems: theory and practice. arXiv preprint arXiv:2405.11457 (2024) 17
2024 arXiv
-
[33]
Scientific Reports 13(1), 9397 (2023) 18
Liu, Y., Zou, Z., Pak, O.S., Tsang, A.C.: Learning to cooperate for low-reynolds- number swimming: a model problem for gait coordination. Scientific Reports 13(1), 9397 (2023) 18
2023
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.