1
|
Sammüller F, Hermann S, Schmidt M. Why neural functionals suit statistical mechanics. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2024; 36:243002. [PMID: 38467072 DOI: 10.1088/1361-648x/ad326f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/29/2023] [Accepted: 03/11/2024] [Indexed: 03/13/2024]
Abstract
We describe recent progress in the statistical mechanical description of many-body systems via machine learning combined with concepts from density functional theory and many-body simulations. We argue that the neural functional theory by Sammülleret al(2023Proc. Natl Acad. Sci.120e2312484120) gives a functional representation of direct correlations and of thermodynamics that allows for thorough quality control and consistency checking of the involved methods of artificial intelligence. Addressing a prototypical system we here present a pedagogical application to hard core particle in one spatial dimension, where Percus' exact solution for the free energy functional provides an unambiguous reference. A corresponding standalone numerical tutorial that demonstrates the neural functional concepts together with the underlying fundamentals of Monte Carlo simulations, classical density functional theory, machine learning, and differential programming is available online athttps://github.com/sfalmo/NeuralDFT-Tutorial.
Collapse
Affiliation(s)
- Florian Sammüller
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| | - Sophie Hermann
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| | - Matthias Schmidt
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| |
Collapse
|
2
|
Eckert T, Stuhlmüller NCX, Sammüller F, Schmidt M. Local measures of fluctuations in inhomogeneous liquids: statistical mechanics and illustrative applications. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2023; 35:425102. [PMID: 37414000 DOI: 10.1088/1361-648x/ace50c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/12/2023] [Accepted: 07/06/2023] [Indexed: 07/08/2023]
Abstract
We show in detail how three one-body fluctuation profiles, namely the local compressibility, the local thermal susceptibility, and the reduced density, can be obtained from a statistical mechanical many-body description of classical particle-based systems. We present several different and equivalent routes to the definition of each fluctuation profile, facilitating their explicit numerical calculation in inhomogeneous equilibrium systems. This underlying framework is used for the derivation of further properties such as hard wall contact theorems and novel types of inhomogeneous one-body Ornstein-Zernike equations. The practical accessibility of all three fluctuation profiles is exemplified by grand canonical Monte Carlo simulations that we present for hard sphere, Gaussian core and Lennard-Jones fluids in confinement.
Collapse
Affiliation(s)
- Tobias Eckert
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95440 Bayreuth, Germany
| | - Nico C X Stuhlmüller
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95440 Bayreuth, Germany
| | - Florian Sammüller
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95440 Bayreuth, Germany
| | - Matthias Schmidt
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95440 Bayreuth, Germany
| |
Collapse
|
3
|
Sammüller F, Hermann S, de Las Heras D, Schmidt M. Noether-Constrained Correlations in Equilibrium Liquids. PHYSICAL REVIEW LETTERS 2023; 130:268203. [PMID: 37450808 DOI: 10.1103/physrevlett.130.268203] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/26/2023] [Accepted: 06/12/2023] [Indexed: 07/18/2023]
Abstract
Liquid structure carries deep imprints of an inherent thermal invariance against a spatial transformation of the underlying classical many-body Hamiltonian. At first order in the transformation field Noether's theorem yields the local force balance. Three distinct two-body correlation functions emerge at second order, namely the standard two-body density, the localized force-force correlation function, and the localized force gradient. An exact Noether sum rule interrelates these correlators. Simulations of Lennard-Jones, Yukawa, soft-sphere dipolar, Stockmayer, Gay-Berne and Weeks-Chandler-Andersen liquids, of monatomic water and of a colloidal gel former demonstrate the fundamental role in the characterization of spatial structure.
Collapse
Affiliation(s)
- Florian Sammüller
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| | - Sophie Hermann
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| | - Daniel de Las Heras
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| | - Matthias Schmidt
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| |
Collapse
|
4
|
de Las Heras D, Zimmermann T, Sammüller F, Hermann S, Schmidt M. Perspective: How to overcome dynamical density functional theory. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2023; 35:271501. [PMID: 37023762 DOI: 10.1088/1361-648x/accb33] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/28/2023] [Accepted: 04/06/2023] [Indexed: 06/19/2023]
Abstract
We argue in favour of developing a comprehensive dynamical theory for rationalizing, predicting, designing, and machine learning nonequilibrium phenomena that occur in soft matter. To give guidance for navigating the theoretical and practical challenges that lie ahead, we discuss and exemplify the limitations of dynamical density functional theory (DDFT). Instead of the implied adiabatic sequence of equilibrium states that this approach provides as a makeshift for the true time evolution, we posit that the pending theoretical tasks lie in developing a systematic understanding of the dynamical functional relationships that govern the genuine nonequilibrium physics. While static density functional theory gives a comprehensive account of the equilibrium properties of many-body systems, we argue that power functional theory is the only present contender to shed similar insights into nonequilibrium dynamics, including the recognition and implementation of exact sum rules that result from the Noether theorem. As a demonstration of the power functional point of view, we consider an idealized steady sedimentation flow of the three-dimensional Lennard-Jones fluid and machine-learn the kinematic map from the mean motion to the internal force field. The trained model is capable of both predicting and designing the steady state dynamics universally for various target density modulations. This demonstrates the significant potential of using such techniques in nonequilibrium many-body physics and overcomes both the conceptual constraints of DDFT as well as the limited availability of its analytical functional approximations.
Collapse
Affiliation(s)
- Daniel de Las Heras
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| | - Toni Zimmermann
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| | - Florian Sammüller
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| | - Sophie Hermann
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| | - Matthias Schmidt
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| |
Collapse
|
5
|
Renner J, Schmidt M, de Las Heras D. Reduced-variance orientational distribution functions from torque sampling. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2023; 35:235901. [PMID: 36974000 DOI: 10.1088/1361-648x/acc522] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/15/2022] [Accepted: 03/08/2023] [Indexed: 06/18/2023]
Abstract
We introduce a method to sample the orientational distribution function in computer simulations. The method is based on the exact torque balance equation for classical many-body systems of interacting anisotropic particles in equilibrium. Instead of the traditional counting of events, we reconstruct the orientational distribution function via an orientational integral of the torque acting on the particles. We test the torque sampling method in two- and three-dimensions, using both Langevin dynamics and overdamped Brownian dynamics, and with two interparticle interaction potentials. In all cases the torque sampling method produces profiles of the orientational distribution function with better accuracy than those obtained with the traditional counting method. The accuracy of the torque sampling method is independent of the bin size, and hence it is possible to resolve the orientational distribution function with arbitrarily small angular resolutions.
Collapse
Affiliation(s)
- Johannes Renner
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95440 Bayreuth, Germany
| | - Matthias Schmidt
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95440 Bayreuth, Germany
| | - Daniel de Las Heras
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95440 Bayreuth, Germany
| |
Collapse
|
6
|
Sammüller F, Hermann S, Schmidt M. Comparative study of force-based classical density functional theory. Phys Rev E 2023; 107:034109. [PMID: 37072997 DOI: 10.1103/physreve.107.034109] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2022] [Accepted: 02/17/2023] [Indexed: 04/20/2023]
Abstract
We reexamine results obtained with the recently proposed density functional theory framework based on forces (force-DFT) [S. M. Tschopp et al., Phys. Rev. E 106, 014115 (2022)2470-004510.1103/PhysRevE.106.014115]. We compare inhomogeneous density profiles for hard sphere fluids to results from both standard density functional theory and from computer simulations. Test situations include the equilibrium hard sphere fluid adsorbed against a planar hard wall and the dynamical relaxation of hard spheres in a switched harmonic potential. The comparison to grand canonical Monte Carlo simulation profiles shows that equilibrium force-DFT alone does not improve upon results obtained with the standard Rosenfeld functional. Similar behavior holds for the relaxation dynamics, where we use our event-driven Brownian dynamics data as benchmark. Based on an appropriate linear combination of standard and force-DFT results, we investigate a simple hybrid scheme which rectifies these deficiencies in both the equilibrium and the dynamical case. We explicitly demonstrate that although the hybrid method is based on the original Rosenfeld fundamental measure functional, its performance is comparable to that of the more advanced White Bear theory.
Collapse
Affiliation(s)
- Florian Sammüller
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| | - Sophie Hermann
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| | - Matthias Schmidt
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| |
Collapse
|
7
|
Tschopp SM, Brader JM. First-principles superadiabatic theory for the dynamics of inhomogeneous fluids. J Chem Phys 2022; 157:234108. [PMID: 36550050 DOI: 10.1063/5.0131441] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022] Open
Abstract
For classical many-body systems subject to Brownian dynamics, we develop a superadiabatic dynamical density functional theory (DDFT) for the description of inhomogeneous fluids out-of-equilibrium. By explicitly incorporating the dynamics of the inhomogeneous two-body correlation functions, we obtain superadiabatic forces directly from the microscopic interparticle interactions. We demonstrate the importance of these nonequilibrium forces for an accurate description of the one-body density by numerical implementation of our theory for three-dimensional hard-spheres in a time-dependent planar potential. The relaxation of the one-body density in superadiabatic-DDFT is found to be slower than that predicted by standard adiabatic DDFT and significantly improves the agreement with Brownian dynamics simulation data. We attribute this improved performance to the correct treatment of structural relaxation within the superadiabatic-DDFT. Our approach provides fundamental insight into the underlying structure of dynamical density functional theories and makes possible the study of situations for which standard approaches fail.
Collapse
Affiliation(s)
- S M Tschopp
- Department of Physics, University of Fribourg, CH-1700 Fribourg, Switzerland
| | - J M Brader
- Department of Physics, University of Fribourg, CH-1700 Fribourg, Switzerland
| |
Collapse
|
8
|
Te Vrugt M, Wittkowski R. Perspective: New directions in dynamical density functional theory. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2022; 35:041501. [PMID: 35917827 DOI: 10.1088/1361-648x/ac8633] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/23/2022] [Accepted: 08/01/2022] [Indexed: 06/15/2023]
Abstract
Classical dynamical density functional theory (DDFT) has become one of the central modeling approaches in nonequilibrium soft matter physics. Recent years have seen the emergence of novel and interesting fields of application for DDFT. In particular, there has been a remarkable growth in the amount of work related to chemistry. Moreover, DDFT has stimulated research on other theories such as phase field crystal models and power functional theory. In this perspective, we summarize the latest developments in the field of DDFT and discuss a variety of possible directions for future research.
Collapse
Affiliation(s)
- Michael Te Vrugt
- Institut für Theoretische Physik, Center for Soft Nanoscience, Westfälische Wilhelms-Universität Münster, 48149 Münster, Germany
| | - Raphael Wittkowski
- Institut für Theoretische Physik, Center for Soft Nanoscience, Westfälische Wilhelms-Universität Münster, 48149 Münster, Germany
| |
Collapse
|
9
|
Hermann S, Schmidt M. Why Noether's theorem applies to statistical mechanics. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2022; 34:213001. [PMID: 35255482 DOI: 10.1088/1361-648x/ac5b47] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/08/2021] [Accepted: 03/07/2022] [Indexed: 06/14/2023]
Abstract
Noether's theorem is familiar to most physicists due its fundamental role in linking the existence of conservation laws to the underlying symmetries of a physical system. Typically the systems are described in the particle-based context of classical mechanics or on the basis of field theory. We have recently shown (2021Commun. Phys.4176) that Noether's reasoning also applies to thermal systems, where fluctuations are paramount and one aims for a statistical mechanical description. Here we give a pedagogical introduction based on the canonical ensemble and apply it explicitly to ideal sedimentation. The relevant mathematical objects, such as the free energy, are viewed as functionals. This vantage point allows for systematic functional differentiation and the resulting identities express properties of both macroscopic average forces and molecularly resolved correlations in many-body systems, both in and out-of-equilibrium, and for active Brownian particles. To provide further background, we briefly describe the variational principles of classical density functional theory, of power functional theory, and of classical mechanics.
Collapse
Affiliation(s)
- Sophie Hermann
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| | - Matthias Schmidt
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| |
Collapse
|