Reuse & Permissions

It is not necessary to obtain permission to reuse this article or its components as it is available under the terms of the Creative Commons Attribution 4.0 International license. This license permits unrestricted use, distribution, and reproduction in any medium, provided attribution to the author(s) and the published article's title, journal citation, and DOI are maintained. Please note that some figures may have been included with permission from other third parties. It is your responsibility to obtain the proper permission from the rights holder directly for these figures.

Export citation

Export citation

Choose format for download:

Download Citation
  • Open Access

Global thermodynamic manifold for conservative control of stochastic systems

Jordan R. Sawchuk* and David A. Sivak†

  • *Contact author: jordan_sawchuk@sfu.ca
  • †Contact author: dsivak@sfu.ca

Phys. Rev. Research 8, 013004 – Published 5 January, 2026

DOI: https://doi.org/10.1103/j59j-q88v

Abstract

Optimal control of stochastic systems plays a central role in nonequilibrium physics, with applications in the study of biological molecular motors and the design of single-molecule experiments. While exact analytic solutions to optimization problems are rare, under slow driving conditions, the problem can be reformulated geometrically solely in terms of equilibrium properties. In this framework, minimum-work protocols are geodesics on a thermodynamic manifold whose metric is a generalized friction tensor. Here, we introduce a foundation for this friction-tensor formalism for conservatively driven systems on either discrete- or continuous-state spaces. Under complete control of the potential energy, a global thermodynamic manifold (on which points are identified with instantaneous energy landscapes) has as its metric a full-control friction tensor. Arbitrary partial-control friction tensors arise naturally as inherited metrics on submanifolds of this global manifold. We leverage a simple mathematical relationship between system dynamics and the geometry of the global manifold to derive expressions for the generalized friction and linear-response excess work. We show that the friction tensor, usually defined as an infinite-time integral, may be expressed as a product of straightforwardly computed matrices. We furthermore show that the linear-response excess power is decomposable into relaxation-time-scaled projections of the control velocity onto eigenmodes of the rate matrix or Fokker-Planck operator. This intuitive result provides a systematic means of interpreting otherwise-mysterious efficient control strategies, unveiling how multiparameter optimal control takes advantage of multiexponential relaxation dynamics to reduce dissipation relative to single-parameter control. We demonstrate the computational and conceptual advances in three illustrative examples.

View figure in article

Physics Subject Headings (PhySH)

Article Text

References (73)

  1. B. Andresen, Current trends in finite-time thermodynamics, Angew. Chem. Int. Ed. 50, 2690 (2011).
  2. U. Seifert, Stochastic thermodynamics, fluctuation theorems, and molecular machines, Rep. Prog. Phys. 75, 126001 (2012).
  3. A. I. Brown and D. A. Sivak, Theory of nonequilibrium free energy transduction by molecular machines, Chem. Rev. 120, 434 (2020).
  4. S. Blaber and D. A. Sivak, Optimal control in stochastic thermodynamics, J. Phys. Commun. 7, 033001 (2023).
  5. D. Guéry-Odelin, C. Jarzynski, C. A. Plata, A. Prados, and E. Trizac, Driving rapidly while remaining in control: Classical shortcuts from Hamiltonian to stochastic dynamics, Rep. Prog. Phys. 86, 035902 (2023).
  6. T. Schmiedl and U. Seifert, Optimal finite-time processes in stochastic thermodynamics, Phys. Rev. Lett. 98, 108301 (2007).
  7. A. Gomez-Marin, T. Schmiedl, and U. Seifert, Optimal protocols for minimal work processes in underdamped stochastic thermodynamics, J. Chem. Phys. 129, 024114 (2008).
  8. M. Esposito, R. Kawai, K. Lindenberg, and C. Broeck, Finite time thermodynamics for a single level quantum dot, Europhys. Lett. 89, 20003 (2010).
  9. E. Aurell, K. Gawędzki, C. Mejía-Monasterio, R. Mohayaee, and P. Muratore-Ginanneschi, Refined second law of thermodynamics for fast random processes, J. Stat. Phys. 147, 487 (2012).
  10. C. A. Plata, D. Guéry-Odelin, E. Trizac, and A. Prados, Optimal work in a harmonic trap with bounded stiffness, Phys. Rev. E 99, 012140 (2019).
  11. A. Zhong and M. R. DeWeese, Limited-control optimal protocols arbitrarily far from equilibrium, Phys. Rev. E 106, 044135 (2022).
  12. J. Yan, H. Touchette, and G. M. Rotskoff, Learning nonequilibrium control forces to characterize dynamical phase transitions, Phys. Rev. E 105, 024115 (2022).
  13. M. C. Engel, J. A. Smith, and M. P. Brenner, Optimal control of nonequilibrium systems through automatic differentiation, Phys. Rev. X 13, 041032 (2023).
  14. S. Blaber, M. D. Louwerse, and D. A. Sivak, Steps minimize dissipation in rapidly driven stochastic systems, Phys. Rev. E 104, L022101 (2021).
  15. D. A. Sivak and G. E. Crooks, Thermodynamic metrics and optimal paths, Phys. Rev. Lett. 108, 190602 (2012).
  16. M. V. S. Bonanča and S. Deffner, Optimal driving of isothermal processes close to equilibrium, J. Chem. Phys. 140, 244119 (2014).
  17. D. A. Sivak and G. E. Crooks, Thermodynamic geometry of minimum-dissipation driven barrier crossing, Phys. Rev. E 94, 052106 (2016).
  18. J. N. E. Lucero, A. Mehdizadeh, and D. A. Sivak, Optimal control of rotary motors, Phys. Rev. E 99, 012119 (2019).
  19. 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).
  20. S. Tafoya, S. J. Large, S. Liu, C. Bustamante, and D. A. Sivak, Using a system’s equilibrium behavior to reduce its energy dissipation in nonequilibrium processes, Proc. Natl. Acad. Sci. USA 116, 5920 (2019).
  21. S. Blaber and D. A. Sivak, Optimal control of protein copy number, Phys. Rev. E 101, 022118 (2020).
  22. J.-F. Chen, R.-X. Zhai, C. P. Sun, and H. Dong, Geodesic lower bound of the energy consumption to achieve membrane separation within finite time, PRX Energy 2, 033003 (2023).
  23. S. Deffner, Kibble-Zurek scaling of the irreversible entropy production, Phys. Rev. E 96, 052125 (2017).
  24. P. R. Zulkowski and M. R. DeWeese, Optimal control of overdamped systems, Phys. Rev. E 92, 032117 (2015).
  25. V. Lindahl, J. Lidmar, and B. Hess, Riemann metric approach to optimal sampling of multidimensional free-energy landscapes, Phys. Rev. E 98, 023312 (2018).
  26. M. Lundborg, J. Lidmar, and B. Hess, On the path to optimal alchemistry, Protein J. 42, 477 (2023).
  27. A. Zhong and M. DeWeese, Beyond linear response: Equivalence between thermodynamic geometry and optimal transport, Phys. Rev. Lett. 133, 057102 (2024).
  28. W. Cheney, Analysis for Applied Mathematicians (Springer-Verlag New York, New York, 2001).
  29. B. Remlein and U. Seifert, Optimality of nonconservative driving for finite-time processes with discrete states, Phys. Rev. E 103, L050105 (2021).
  30. The technical condition (the Bakry-Emery criterion) is that V be λ-convex [31]. A quadratic potential in Rd is an example of a λ-convex potential.
  31. G. A. Pavliotis, Stochastic Processes and Applications: Diffusion Processes, the Fokker-Planck and Langevin Equations (Springer Nature, New York, 2014), Vol. 60.
  32. P. Salamon and A. Nitzan, Finite time optimizations of a Newton’s law Carnot cycle, J. Chem. Phys. 74, 3546 (1981).
  33. K. Sekimoto, Langevin equation and thermodynamics, Prog. Theor. Phys. Suppl. 130, 17 (1998).
  34. D. Chandler, Introduction to Modern Statistical Mechanics (Oxford University Press, Oxford, 1987).
  35. V. V. Prasolov, Differential Geometry (Springer, Cham, Switzerland, 2022).
  36. P. R. Zulkowski and M. R. DeWeese, Optimal finite-time erasure of a classical bit, Phys. Rev. E 89, 052140 (2014).
  37. A. G. Frim and M. R. DeWeese, Optimal finite-time Brownian Carnot engine, Phys. Rev. E 105, L052103 (2022).
  38. G. M. Rotskoff and G. E. Crooks, Optimal control in nonequilibrium systems: Dynamic Riemannian geometry of the Ising model, Phys. Rev. E 92, 060102(R) (2015).
  39. M. D. Louwerse and D. A. Sivak, Multidimensional minimum-work control of a 2D Ising model, J. Chem. Phys. 156, 194108 (2022).
  40. 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).
  41. While this section is intended to be self-contained, the essential geometric ideas will be unfamiliar to many readers. For a thorough overview of differential geometry including submanifolds and degenerate metrics with notation that will be familiar to readers versed in general relativity, see Ref. [72]. For a more modern treatment of the subject, see Ref. [73].
  42. Taking the inner product of both sides of Eq. (27) with the constant function 1 and applying (1,δp)=0 leads to the same conclusion. Since ζ is symmetric, in this context the arbitrariness of the potential-energy zero point and conservation of probability are equivalent. This highlights a duality between energy and probability, manifest in the geometry: The linear-response relation [Eq. (27)] states that δp is the vector dual to dV/dt. This is also reflected in the conjugate forces to the state energies [Eq. (22)].
  43. J. Liphardt, S. Dumont, S. B. Smith, I. Tinoco, and C. Bustamante, Equilibrium information from nonequilibrium measurements in an experimental test of Jarzynski’s equality, Science 296, 1832 (2002).
  44. D. Collin, F. Ritort, C. Jarzynski, S. B. Smith, I. Tinoco, and C. Bustamante, Verification of the Crooks fluctuation theorem and recovery of RNA folding free energies, Nature (London) 437, 231 (2005).
  45. G. E. Crooks, On the Drazin inverse of the rate matrix, 2018, https://threeplusone.com/pubs/on_drazin.pdf.
  46. D. Mandal and C. Jarzynski, Analysis of slow transitions between nonequilibrium steady states, J. Stat. Mech. (2016) 063204.
  47. M. Scandi and M. Perarnau-Llobet, Thermodynamic length in open quantum systems, Quantum 3, 197 (2019).
  48. H. J. D. Miller, M. Scandi, J. Anders, and M. Perarnau-Llobet, Work fluctuations in slow processes: Quantum signatures and optimal control, Phys. Rev. Lett. 123, 230603 (2019).
  49. D. D. Yao, First-passage-time moments of Markov processes, J. Appl. Probab. 22, 939 (1985).
  50. P. Coolen-Schrijner and E. A. van Doorn, The deviation matrix of a continuous-time Markov chain, Prob. Eng. Inf. Sci. 16, 351 (2002).
  51. G. Wang, Y. Wei, and S. Qiao, Generalized Inverses: Theory and Computations, Developments in Mathematics (Springer, Singapore, 2018), Vol. 53.
  52. N. G. Van Kampen, The master equation, in Stochastic Processes in Physics and Chemistry, 3rd ed. (Elsevier, Amsterdam, 2007).
  53. Property (36b) can be re-written as (AD)2A=AD by the commutative property (36c). Thus, Ax=0⇒ADx=0. An identical argument holds for right multiplication.
  54. H. Risken, in The Fokker-Planck Equation: Methods of Solution and Applications, 2nd ed., edited by T. Frank (Springer, Berlin, 1996).
  55. F. Covito, F. G. Eich, R. Tuovinen, M. A. Sentef, and A. Rubio, Transient charge and energy flow in the wide-band limit, J. Chem. Theory Comput. 14, 2495 (2018).
  56. G. M. Rotskoff, G. E. Crooks, and E. Vanden-Eijnden, Geometric approach to optimal nonequilibrium control: Minimizing dissipation in nanomagnetic spin systems, Phys. Rev. E 95, 012148 (2017).
  57. M. D. Louwerse and D. A. Sivak, Connections between efficient control and spontaneous transitions in an Ising model, Phys. Rev. E 106, 064124 (2022).
  58. W. Huang, D. A. Kitchaev, S. T. Dacek, Z. Rong, A. Urban, S. Cao, C. Luo, and G. Ceder, Finding and proving the exact ground state of a generalized Ising model by convex optimization and MAX-SAT, Phys. Rev. B 94, 134424 (2016).
  59. M. Mézard and A. Montanari, Spin glasses, in Information, Physics, and Computation, edited by M. Mézard and A. Montanari (Oxford University Press, Oxford, 2009).
  60. A. Decelle, C. Furtlehner, A. D. J. N. Gómez, and B. Seoane, Inferring effective couplings with restricted Boltzmann machines, SciPost Phys. 16, 095 (2024).
  61. N. Bulso and Y. Roudi, Restricted Boltzmann machines as models of interacting variables, Neural Comput. 33, 2646 (2021).
  62. R. J. Glauber, Time-dependent statistics of the Ising model, J. Math. Phys 4, 294 (1963).
  63. A. Mariz and H. Herrmann, Mixing heat-bath and Glauber dynamics: Damage spreading in the Ising model, J. Phys. A: Math. Gen. 22, L1081 (1989).
  64. E. Aurell, G. Del Ferraro, E. Domínguez, and R. Mulet, Cavity master equation for the continuous time dynamics of discrete-spin models, Phys. Rev. E 95, 052119 (2017).
  65. A. M. Mariz, F. D. Nobre, and C. Tsallis, Generalized single-spin-flip dynamics for the Ising model and thermodynamic properties, Phys. Rev. B 49, 3576 (1994).
  66. S. Blaber and D. A. Sivak, Efficient two-dimensional control of barrier crossing, Europhys. Lett. 139, 17001 (2022).
  67. B. Roux, Transition rate theory, spectral analysis, and reactive paths, J. Chem. Phys. 156, 134111 (2022).
  68. alliancecan.ca.
  69. J. Sawchuk, 4 spin system protocols and spectra [Data set] Zenodo (2025), https://doi.org/10.5281/zenodo.15678284.
  70. H. Yanai, K. Takeuchi, and Y. Takane, Projection Matrices, Generalized Inverse Matrices, and Singular Value Decomposition (Springer, New York, NY, 2011).
  71. L. N. Trefethen and I. D. Bau, Numerical Linear Algebra (Society for Industrial and Applied Mathematics, Philadelphia, PA, 1997).
  72. L. P. Eisenhart, Riemannian Geometry (Princeton University Press, Princeton, NJ, 1926).
  73. B.-Y. Chen, Riemannian submanifolds, in Handbook of Differential Geometry, edited by F. J. Dillen and L. C. Verstraelen, Vol. 1 (North-Holland, Amsterdam, 2000), Chap. 3, pp. 187–418.

Outline

Information

Sign In to Your Journals Account

Filter

Filter

Article Lookup

Enter a citation