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

Machine learning approach to fast thermal equilibration

Diego Rengifo*

Gabriel Téllez†

  • *Contact author: d.rengifo@uniandes.edu.co
  • †Contact author: gtellez@uniandes.edu.co

Phys. Rev. E 111, 065311 – Published 23 June, 2025

DOI: https://doi.org/10.1103/q18f-yrm2

Abstract

We present a method to design driving protocols that achieve fast thermal equilibration of a system of interest using techniques inspired by machine learning training algorithms. For example, consider a Brownian particle manipulated by optical tweezers. The force on the particle can be controlled and adjusted over time, resulting in a driving protocol that transitions the particle from an initial state to a final state. Once the driving protocol has been completed, the system requires additional time to relax to thermal equilibrium. Designing driving protocols that bypass the relaxation period is of interest so that, at the end of the protocol, the system is either in thermal equilibrium or very close to it. Several studies have addressed this problem through reverse engineering methods, which involve prescribing a specific evolution for the probability density function of the system and then deducing the corresponding form of the driving protocol potential. Here we propose a new method that can be applied to more complex systems where reverse engineering is not feasible. We simulate the evolution of a large ensemble of trajectories while tracking the gradients with respect to a parametrization of the driving protocol. The final probability density function is compared to the target equilibrium one. Using machine learning libraries, the gradients are computed via backpropagation and the protocol is iteratively adjusted until the optimal protocol is achieved. We demonstrate the effectiveness of our approach with several examples.

View figure in article

Physics Subject Headings (PhySH)

Article Text

References (45)

  1. J. Bechhoefer, Control Theory for Physicists (Cambridge University Press, Cambridge, UK, 2021).
  2. G. M. Wang, E. M. Sevick, E. Mittag, D. J. Searles, and D. J. Evans, Phys. Rev. Lett. 89, 050601 (2002).
  3. C. Gonzalez-Ballestero, M. Aspelmeyer, L. Novotny, R. Quidant, and O. Romero-Isart, Science 374, eabg3027 (2021) .
  4. D. Raynal, T. de Guillebon, D. Guéry-Odelin, E. Trizac, J.-S. Lauret, and L. Rondin, Phys. Rev. Lett. 131, 087101 (2023).
  5. I. Martínez, A. Petrosyan, D. Guéry-Odelin, E. Trizac, and S. Ciliberto, Nat. Phys. 12, 843 (2016).
  6. G. Li, H. T. Quan, and Z. C. Tu, Phys. Rev. E 96, 012144 (2017).
  7. D. Guéry-Odelin, C. Jarzynski, C. A. Plata, A. Prados, and E. Trizac, Rep. Prog. Phys. 86, 035902 (2023).
  8. S. Dago, B. Besga, R. Mothe, D. Guéry-Odelin, E. Trizac, A. Petrosyan, L. Bellon, and S. Ciliberto, SciPost Phys. 9, 064 (2020).
  9. A. G. Baydin, B. A. Pearlmutter, A. A. Radul, and J. M. Siskind, J. Mach. Learn. Res. 18, 1 (2018).
  10. J. Sohl-Dickstein, E. Weiss, N. Maheswaranathan, and S. Ganguli, in Proceedings of the 32nd International Conference on Machine Learning, Proceedings of Machine Learning Research Vol. 37, edited by F. Bach and D. Blei (PMLR, Lille, France, 2015), pp. 2256–2265.
  11. J. Ho, A. Jain, and P. Abbeel, in 34th Conference on Neural Information Processing Systems (NeurIPS 2020) (2020) arXiv:2006.11239.
  12. M. C. Engel, J. A. Smith, and M. P. Brenner, Phys. Rev. X 13, 041032 (2023).
  13. C. Casert and S. Whitelam, Nat. Commun. 15, 9128 (2024).
  14. A. Karpathy, The spelled-out intro to neural networks and backpropagation: Building micrograd, https://www.youtube.com/watch?v=VMj-3S1tku0.
  15. Pytorch, https://pytorch.org/.
  16. Tensorflow, https://www.tensorflow.org/.
  17. Jax: High-performance array computing, https://jax.readthedocs.io/.
  18. K. Hornik, M. Stinchcombe, and H. White, Neural Networks 2, 359 (1989).
  19. S. Whitelam and I. Tamblyn, Phys. Rev. E 101, 052604 (2020).
  20. C. P. Goodrich, E. M. King, S. S. Schoenholz, E. D. Cubuk, and M. P. Brenner, Proc. Natl. Acad. Sci. USA 118, e2024083118 (2021).
  21. S. Whitelam, Phys. Rev. X 13, 021005 (2023).
  22. S. Whitelam, Phys. Rev. E 108, 044138 (2023).
  23. N. Barros, S. Whitelam, S. Ciliberto, and L. Bellon, Phys. Rev. E 111, 044114 (2025).
  24. S. A. M. Loos, S. Monter, F. Ginot, and C. Bechinger, Phys. Rev. X 14, 021032 (2024).
  25. G. Biroli and M. Mézard, J. Stat. Mech. (2023) 093402.
  26. G. Biroli, T. Bonnaire, V. de Bortoli, and M. Mézard, Nat. Commun. 15, 9957 (2024).
  27. S. Whitelam, V. Selin, I. Benlolo, C. Casert, and I. Tamblyn, Mach. Learn.: Sci. Technol. 3, 045026 (2022).
  28. E. Aurell, C. Mejía-Monasterio, and P. Muratore-Ginanneschi, Phys. Rev. Lett. 106, 250601 (2011).
  29. C. Villani, Topics in Optimal Transportation, Graduate Studies in Mathematics (American Mathematical Society, Providence, RI, 2003).
  30. J. Sanders, M. Baldovin, and P. Muratore-Ginanneschi, Phys. Rev. E 111, 034127 (2025).
  31. J. Klinger and G. M. Rotskoff, Phys. Rev. E 111, 014114 (2025).
  32. P. Muratore-Ginanneschi and K. Schwieger, Phys. Rev. E 90, 060102(R) (2014).
  33. R. T. Q. Chen, Y. Rubanova, J. Bettencourt, and D. Duvenaud, Neural ordinary differential equations, arXiv:1806.07366.
  34. L. S. Pontryagin, V. G. Boltyanskii, R. V. Gamkrelidze, and E. F. Mishchenko, The Mathematical Theory of Optimal Processes (Interscience Publishers, Geneva, 1962).
  35. X. Li, T.-K. L. Wong, R. T. Q. Chen, and D. Duvenaud, Scalable gradients for stochastic differential equations, arXiv:2001.01328.
  36. D. Rengifo and G. Téllez, Phys. Rev. E 109, 014145 (2024).
  37. L. Peliti and S. Pigolotti, Stochastic Thermodynamics: An Introduction (Princeton University Press, Princeton, NJ, 2021).
  38. T. Schmiedl and U. Seifert, Phys. Rev. Lett. 98, 108301 (2007).
  39. T. Schmiedl and U. Seifert, Europhys. Lett. 81, 20003 (2008).
  40. K. Proesmans, J. Ehrich, and J. Bechhoefer, Phys. Rev. Lett. 125, 100602 (2020).
  41. S. Dago, J. Pereda, N. Barros, S. Ciliberto, and L. Bellon, Phys. Rev. Lett. 126, 170601 (2021).
  42. D. P. Kingma and J. Ba, Adam: A method for stochastic optimization, arXiv:1412.6980.
  43. A. Gomez-Marin, T. Schmiedl, and U. Seifert, J. Chem. Phys. 129, 024114 (2008).
  44. Langesim optim library, https://github.com/GabrielTellez/langesim_optim.
  45. Langesim optim documentation, https://gabrieltellez.github.io/langesim_optim.

Outline

Information

Sign In to Your Journals Account

Filter

Filter

Article Lookup

Enter a citation