Calculation of the three-particle contribution to the configurational entropy for two different models of amorphous Si

Equilibrium states corresponding to targeted hyperuniform nonequilibrium pair statistics

The Zhang-Torquato conjecture [G. Zhang and S. Torquato, Phys. Rev. E, 2020, 101, 032124.] states that any realizable pair correlation function g₂(r) or structure factor S(k) of a translationally invariant nonequilibrium system can be attained by an equilibrium ensemble involving only (up to) effective two-body interactions. To further test and study this conjecture, we consider two singular nonequilibrium models of recent interest that also have the exotic hyperuniformity property: a 2D "perfect glass" and a 3D critical absorbing-state model. We find that each nonequilibrium target can be achieved accurately by equilibrium states with effective one- and two-body potentials, lending further support to the conjecture. To characterize the structural degeneracy of such a nonequilibrium-equilibrium correspondence, we compute higher-order statistics for both models, as well as those for a hyperuniform 3D uniformly randomized lattice (URL), whose higher-order statistics can be very precisely ascertained. Interestingly, we find that the differences in the higher-order statistics between nonequilibrium and equilibrium systems with matching pair statistics, as measured by the "hole" probability distribution, provide measures of the degree to which a system is out of equilibrium. We show that all three systems studied possess the bounded-hole property and that holes near the maximum hole size in the URL are much rarer than those in the underlying simple cubic lattice. Remarkably, upon quenching, the effective potentials for all three systems possess local energy minima (i.e., inherent structures) with stronger forms of hyperuniformity compared to their target counterparts. Our methods are expected to facilitate the self-assembly of tunable hyperuniform soft-matter systems.

Precise determination of pair interactions from pair statistics of many-body systems in and out of equilibrium

The determination of the pair potential v(r) that accurately yields an equilibrium state at positive temperature T with a prescribed pair correlation function g_{2}(r) or corresponding structure factor S(k) in d-dimensional Euclidean space R^{d} is an outstanding inverse statistical mechanics problem with far-reaching implications. Recently, Zhang and Torquato [Phys. Rev. E 101, 032124 (2020)2470-004510.1103/PhysRevE.101.032124] conjectured that any realizable g_{2}(r) or S(k) corresponding to a translationally invariant nonequilibrium system can be attained by a classical equilibrium ensemble involving only (up to) effective pair interactions. Testing this conjecture for nonequilibrium systems as well as for nontrivial equilibrium states requires improved inverse methodologies. We have devised an optimization algorithm to precisely determine effective pair potentials that correspond to pair statistics of general translationally invariant disordered many-body equilibrium or nonequilibrium systems at positive temperatures. This methodology utilizes a parameterized family of pointwise basis functions for the potential function whose initial form is informed by small-, intermediate- and large-distance behaviors dictated by statistical-mechanical theory. Subsequently, a nonlinear optimization technique is utilized to minimize an objective function that incorporates both the target pair correlation function g_{2}(r) and structure factor S(k) so that the small intermediate- and large-distance correlations are very accurately captured. To illustrate the versatility and power of our methodology, we accurately determine the effective pair interactions of the following four diverse target systems: (1) Lennard-Jones system in the vicinity of its critical point, (2) liquid under the Dzugutov potential, (3) nonequilibrium random sequential addition packing, and (4) a nonequilibrium hyperuniform "cloaked" uniformly randomized lattice. We found that the optimized pair potentials generate corresponding pair statistics that accurately match their corresponding targets with total L_{2}-norm errors that are an order of magnitude smaller than that of previous methods. The results of our investigation lend further support to the Zhang-Torquato conjecture. Furthermore, our algorithm will enable one to probe systems with identical pair statistics but different higher-body statistics, which will shed light on the well-known degeneracy problem of statistical mechanics.

Probing the triplet correlation function in liquid water by experiments and molecular simulations

Three-body information of liquid water is extracted using X-ray diffraction experiment as well as in molecular simulations via isothermal pressure derivative of structure factor term.

Triplet correlation functions in liquid water

Triplet correlations have been shown to play a crucial role in the transformation of simple liquids to anomalous tetrahedral fluids [M. Singh, D. Dhabal, A. H. Nguyen, V. Molinero, and C. Chakravarty, Phys. Rev. Lett. 112, 147801 (2014)]. Here we examine triplet correlation functions for water, arguably the most important tetrahedral liquid, under ambient conditions, using configurational ensembles derived from molecular dynamics (MD) simulations and reverse Monte Carlo (RMC) datasets fitted to experimental scattering data. Four different RMC data sets with widely varying hydrogen-bond topologies fitted to neutron and x-ray scattering data are considered [K. T. Wikfeldt, M. Leetmaa, M. P. Ljungberg, A. Nilsson, and L. G. M. Pettersson, J. Phys. Chem. B 113, 6246 (2009)]. Molecular dynamics simulations are performed for two rigid-body effective pair potentials (SPC/E and TIP4P/2005) and the monatomic water (mW) model. Triplet correlation functions are compared with other structural measures for tetrahedrality, such as the O-O-O angular distribution function and the local tetrahedral order distributions. In contrast to the pair correlation functions, which are identical for all the RMC ensembles, the O-O-O triplet correlation function can discriminate between ensembles with different degrees of tetrahedral network formation with the maximally symmetric, tetrahedral SYM dataset displaying distinct signatures of tetrahedrality similar to those obtained from atomistic simulations of the SPC/E model. Triplet correlations from the RMC datasets conform closely to the Kirkwood superposition approximation, while those from MD simulations show deviations within the first two neighbour shells. The possibilities for experimental estimation of triplet correlations of water and other tetrahedral liquids are discussed.