Non-equilibrium Trajectory Sampling (NETS) method for generating free-energy landscapes and steady-state distributions

Jonathan K. Whitmer,Akwasi Nana Prempeh Ansah-Antwi,Pedro H. Amorim Valenca,Mohsen Farshad
DOI: https://doi.org/10.26434/chemrxiv-2023-n1r3c
2023-12-06
Abstract:This study presents a novel method for constructing free-energy profiles and steady-state distributions from either equilibrium or non-equilibrium trajectories along a defined reaction coordinate. The method works by tracking the final states of a swarm of short simulations launched from different initial conditions with no prior knowledge of the free energy landscape. Subsequently, this trajectory information is used to build a transition matrix whose primary eigenvector captures the steady-state occupation probability for each value of the reaction coordinate, yielding the free energy profile in equilibrium. This innovative method holds potential for many new materials and engineering applications where it is desired to know the free energy of rate-limiting configurations as may be relevant for transport processes (in, e.g. battery electrolytes and nano-filtration membranes), complexation (in, e.g., self-assembly and ligand-binding interactions) or in tuning properties such as adsorption. We illustrate the effectiveness of the method by capturing the free energy associated with a one-dimensional barrier potential modeling a separation membrane, and the particle distribution associated with thermophoresis under a temperature gradient. Further extensions and applications of the method are also discussed.
Chemistry
What problem does this paper attempt to address?
This paper presents a new method called Non-Equilibrium Trajectory Sampling (NETS) for constructing free energy landscapes and steady-state distributions from isothermal or non-isothermal trajectories. This method tracks the final states of a group of short simulations starting from different initial conditions, without the need for prior knowledge of the free energy landscape. In this way, a transition matrix can be constructed, whose dominant eigenvectors capture the stationary occupancy probabilities along the reaction coordinate, thereby yielding the free energy profile at constant temperature. NETS has the potential to understand the free energy in materials and engineering applications, such as the rate-limiting configuration of free energy in battery electrolytes and nanofiltration membranes for processes such as transport, self-assembly, and ligand binding. The paper demonstrates the effectiveness of this method through a one-dimensional potential barrier model and the distribution of thermally migrating particles under a temperature gradient, and discusses the extension and application of the method. This research aims to address how to effectively compute and understand free energy in non-equilibrium systems, in order to promote the development of new materials and technologies.