Counting metastable states within the adsorption/desorption hysteresis loop: A molecular simulation study of confinement in heterogeneous pores

Microstructural degeneracy associated with a two-point correlation function and its information content

A two-point correlation function provides a crucial yet an incomplete characterization of a microstructure because distinctly different microstructures may have the same correlation function. In an earlier Letter [Gommes, Jiao, and Torquato, Phys. Rev. Lett. 108, 080601 (2012)], we addressed the microstructural degeneracy question: What is the number of microstructures compatible with a specified correlation function? We computed this degeneracy, i.e., configurational entropy, in the framework of reconstruction methods, which enabled us to map the problem to the determination of ground-state degeneracies. Here, we provide a more comprehensive presentation of the methodology and analyses, as well as additional results. Since the configuration space of a reconstruction problem is a hypercube on which a Hamming distance is defined, we can calculate analytically the energy profile of any reconstruction problem, corresponding to the average energy of all microstructures at a given Hamming distance from a ground state. The steepness of the energy profile is a measure of the roughness of the energy landscape associated with the reconstruction problem, which can be used as a proxy for the ground-state degeneracy. The relationship between this roughness metric and the ground-state degeneracy is calibrated using a Monte Carlo algorithm for determining the ground-state degeneracy of a variety of microstructures, including realizations of hard disks and Poisson point processes at various densities as well as those with known degeneracies (e.g., single disks of various sizes and a particular crystalline microstructure). We show that our results can be expressed in terms of the information content of the two-point correlation functions. From this perspective, the a priori condition for a reconstruction to be accurate is that the information content, expressed in bits, should be comparable to the number of pixels in the unknown microstructure. We provide a formula to calculate the information content of any two-point correlation function, which makes our results broadly applicable to any field in which correlation functions are employed.

Density of States for a specified correlation function and the energy landscape

The degeneracy of two-phase disordered microstructures consistent with a specified correlation function is analyzed by mapping it to a ground-state degeneracy. We determine for the first time the associated density of states via a Monte Carlo algorithm. Our results are explained in terms of the roughness of an energy landscape, defined on a hypercubic configuration space. The use of a Hamming distance in this space enables us to define a roughness metric, which is calculated from the correlation function alone and related quantitatively to the structural degeneracy. This relation is validated for a wide variety of disordered structures.

Stability intervals of metastable states in hysteretic systems

Hysteresis in disordered systems originates in a plethora of metastable states. Previous works focused on their distribution inside the hysteresis. In contrast, an analysis of their range of metastability is proposed. This model, designed to catch the main features of fluid adsorption in porous materials, shows strong evidence that, in the thermodynamic limit, despite that metastable states of finite range can be found, they are exponentially dominated by those infinitely localized states.

Numerical characterization of the density of metastable states within the hysteresis loop in disordered systems

An improved approach is proposed to analyze the density of metastable states within any hysteresis loop, such as those observed in magnetic materials or for adsorption in porous materials. Except for a few analytically tractable models, most calculations have to be performed numerically on finite systems. The main points to be addressed thus concern the average over various material samples (the so-called realizations of the disorder), and the finite size analysis to estimate the thermodynamic limit. As an improvement of previously existing methods, it is proposed to introduce the Fourier transform of the density of metastable states (characteristic function). Its logarithm is shown to be additive and can straightforwardly be averaged over disorder. This procedure leads to a new definition of the complexity in finite size, giving the usual quenched complexity in the thermodynamic limit, while being better suited to performing finite size analysis. The calculations are illustrated on a molecular simulation based model for a simple fluid adsorbed in heterogeneous siliceous tubular pores mimicking mesoporous materials like MCM-41 or porous silicon. This approach is expected to be of general interest for hysteresis phenomena, including magnetic materials.