Geometrical ambiguity of pair statistics. II. Heterogeneous media

Hierarchical reconstruction of three-dimensional porous media from a single two-dimensional image with multiscale entropy statistics

Despite the development of 3D imaging technology, the reconstruction of three-dimensional (3D) microstructure from a single two-dimensional (2D) image is still a prominent problem. In this paper, we propose a hierarchical reconstruction method based on simulated annealing, which is named hierarchical simulated annealing method (HSA), with the multiscale entropy statistics as the morphological information descriptor to reconstruct its corresponding three-dimensional (3D) microstructure from a single two-dimensional (2D) image. Both hierarchical simulated annealing (HSA) method and simulated annealing (SA) method are used to perform on the 2D and 3D microstructure reconstruction from a single 2D image, where the two-point cluster function and the standard two-point correlation function are used as the measurement metrics for the reconstructed 2D and 3D structures. From the 2D reconstructions, it can be seen that all the reconstructions of HSA method and SA method not only captures the similar morphological information with the original images, but also have a good agreement with the target microstructures in two-point cluster function. For the reconstructed 3D microstructures, the comparison of two-point correlation function shows that both HSA method and SA method can effectively reconstruct its 3D microstructure and the comparison of the reconstruction time between HSA method and SA method shows that the reconstruction speed of HSA method is an order of magnitude faster than that of SA method.

Fast reconstruction of multiphase microstructures based on statistical descriptors

In this paper, we propose a hierarchical simulated annealing of erosion method (HSAE) to improve the computational efficiency of multiphase microstructure reconstruction, whose computational efficiency can be improved by an order of magnitude. Reconstruction of the two-dimensional (2D) and three-dimensional (3D) multiphase microstructures (pore, grain, and clay) based on simulated annealing (SA) and HSAE are performed. In the reconstruction of multiphase microstructure with HSAE and SA, three independent two-point correlation functions are chosen as the morphological information descriptors. The two-point cluster function which contains significant high-order statistical information is used to verify the reconstruction results. From the analysis of 2D reconstruction, it can find that the proposed HSAE technique not only improves the quality of reconstruction, but also improves the computational efficiency. The reconstructions of our proposed method are still imperfect. This is because the used two-point correlation functions contain insufficient information. For the 3D reconstruction, the two-point correlation functions of the 3D generation are in excellent agreement with those of the original 2D image, which illustrates that our proposed method is effective for the reconstruction of 3D microstructure. The comparison of the energy vs computational time between the SA and HSAE methods shows that our presented method is an order of magnitude faster than the SA method. That is because only some of the pixels in the overall hierarchy need to be considered for sampling.

Quantifying microstructural evolution via time-dependent reduced-dimension metrics based on hierarchical n-point polytope functions

We devise reduced-dimension metrics for effectively measuring the distance between two points (i.e., microstructures) in the microstructure space and quantifying the pathway associated with microstructural evolution, based on a recently introduced set of hierarchical n-point polytope functions P_{n}. The P_{n} functions provide the probability of finding particular n-point configurations associated with regular n polytopes in the material system, and are a special subset of the standard n-point correlation functions S_{n} that effectively decompose the structural features in the system into regular polyhedral basis with different symmetries. The nth order metric Ω_{n} is defined as the L_{1} norm associated with the P_{n} functions of two distinct microstructures. By choosing a reference initial state (i.e., a microstructure associated with t_{0}=0), the Ω_{n}(t) metrics quantify the evolution of distinct polyhedral symmetries and can in principle capture emerging polyhedral symmetries that are not apparent in the initial state. To demonstrate their utility, we apply the Ω_{n} metrics to a two-dimensional binary system undergoing spinodal decomposition to extract the phase separation dynamics via the temporal scaling behavior of the corresponding Ω_{n}(t), which reveals mechanisms governing the evolution. Moreover, we employ Ω_{n}(t) to analyze pattern evolution during vapor deposition of phase-separating alloy films with different surface contact angles, which exhibit rich evolution dynamics including both unstable and oscillating patterns. The Ω_{n} metrics have potential applications in establishing quantitative processing-structure-property relationships, as well as real-time processing control and optimization of complex heterogeneous material systems.

Understanding degeneracy of two-point correlation functions via Debye random media

It is well known that the degeneracy of two-phase microstructures with the same volume fraction and two-point correlation function S_{2}(r) is generally infinite. To elucidate the degeneracy problem explicitly, we examine Debye random media, which are entirely defined by a purely exponentially decaying two-point correlation function S_{2}(r). In this work, we consider three different classes of Debye random media. First, we generate the "most probable" class using the Yeong-Torquato construction algorithm [Yeong and Torquato, Phys. Rev. E 57, 495 (1998)1063-651X10.1103/PhysRevE.57.495]. A second class of Debye random media is obtained by demonstrating that the corresponding two-point correlation functions are effectively realized in the first three space dimensions by certain models of overlapping, polydisperse spheres. A third class is obtained by using the Yeong-Torquato algorithm to construct Debye random media that are constrained to have an unusual prescribed pore-size probability density function. We structurally discriminate these three classes of Debye random media from one another by ascertaining their other statistical descriptors, including the pore-size, surface correlation, chord-length probability density, and lineal-path functions. We also compare and contrast the percolation thresholds as well as the diffusion and fluid transport properties of these degenerate Debye random media. We find that these three classes of Debye random media are generally distinguished by the aforementioned descriptors, and their microstructures are also visually distinct from one another. Our work further confirms the well-known fact that scattering information is insufficient to determine the effective physical properties of two-phase media. Additionally, our findings demonstrate the importance of the other two-point descriptors considered here in the design of materials with a spectrum of physical properties.

On the generation of periodic discrete structures with identical two-point correlation

Strategies for the generation of periodic discrete structures with identical two-point correlation-called 2PC-equivalent-are developed. It is shown that starting from a set of 2PC-equivalent root structures, 2PC-equivalent child structures of arbitrary resolution and number of phases (e.g. material phases) can be generated based on phase extension through trivial embeddings, kernel-based extension and phase coalescence. Proofs are provided by means of discrete Fourier transform theory. A Python 3 implementation is offered for reproduction of examples and future applications.

Probing information content of hierarchical n-point polytope functions for quantifying and reconstructing disordered systems

Disordered systems are ubiquitous in physical, biological, and material sciences. Examples include liquid and glassy states of condensed matter, colloids, granular materials, porous media, composites, alloys, packings of cells in avian retina, and tumor spheroids, to name but a few. A comprehensive understanding of such disordered systems requires, as the first step, systematic quantification, modeling, and representation of the underlying complex configurations and microstructure, which is generally very challenging to achieve. Recently, we introduced a set of hierarchical statistical microstructural descriptors, i.e., the "n-point polytope functions" P_{n}, which are derived from the standard n-point correlation functions S_{n}, and successively included higher-order n-point statistics of the morphological features of interest in a concise, explainable, and expressive manner. Here we investigate the information content of the P_{n} functions via optimization-based realization rendering. This is achieved by successively incorporating higher-order P_{n} functions up to n=8 and quantitatively assessing the accuracy of the reconstructed systems via unconstrained statistical morphological descriptors (e.g., the lineal-path function). We examine a wide spectrum of representative random systems with distinct geometrical and topological features. We find that, generally, successively incorporating higher-order P_{n} functions and, thus, the higher-order morphological information encoded in these descriptors leads to superior accuracy of the reconstructions. However, incorporating more P_{n} functions into the reconstruction also significantly increases the complexity and roughness of the associated energy landscape for the underlying stochastic optimization, making it difficult to convergence numerically.

Three-dimensional pore network characterization of reconstructed extracellular matrix

The extracellular matrix (ECM) has a fiber network that provides physical scaffolds to cells and plays important roles by regulating cellular functions. Some previous works characterized the mechanical and geometrical properties of the ECM fiber network using reconstituted collagen-I. However, the characterization of the porous structure of reconstituted collagen-I has been limited to the pore diameter measurement, and pore network extraction has not been applied to reconstituted collagen-I despite the importance of pore interconnectivity. Here, we aim to show the importance of characterizing the pore network of reconstituted collagen-I by comparing the pore networks of structures that have different fiber alignments. We show that the fiber alignment significantly changes the pore throat area but not the pore diameter. Also, we demonstrate that larger pore throats are directed in the direction of the fiber alignment, which may help in understanding the enhanced cell migration when fibers are aligned.

Inferring low-dimensional microstructure representations using convolutional neural networks

We apply recent advances in machine learning and computer vision to a central problem in materials informatics: the statistical representation of microstructural images. We use activations in a pretrained convolutional neural network to provide a high-dimensional characterization of a set of synthetic microstructural images. Next, we use manifold learning to obtain a low-dimensional embedding of this statistical characterization. We show that the low-dimensional embedding extracts the parameters used to generate the images. According to a variety of metrics, the convolutional neural network method yields dramatically better embeddings than the analogous method derived from two-point correlations alone.

Dynamic reconstruction of heterogeneous materials and microstructure evolution

Reconstructing heterogeneous materials from limited structural information has been a topic that attracts extensive research efforts and still poses many challenges. The Yeong-Torquato procedure is one of the most popular reconstruction techniques, in which the material reconstruction problem based on a set of spatial correlation functions is formulated as a constrained energy minimization (optimization) problem and solved using simulated annealing [Yeong and Torquato, Phys. Rev. E 57, 495 (1998)]. The standard two-point correlation function S2 has been widely used in reconstructions, but can also lead to large structural degeneracy for certain nearly percolating systems. To improve reconstruction accuracy and reduce structural degeneracy, one can successively incorporate additional morphological information (e.g., nonconventional or higher-order correlation functions), which amounts to reshaping the energy landscape to create a deep (local) energy minimum. In this paper, we present a dynamic reconstruction procedure that allows one to use a series of auxiliary S2 to achieve the same level of accuracy as those incorporating additional nonconventional correlation functions. In particular, instead of randomly sampling the microstructure space as in the simulated annealing scheme, our procedure utilizes a series of auxiliary microstructures that mimic a physical structural evolution process (e.g., grain growth). This amounts to constructing a series auxiliary energy landscapes that bias the convergence of the reconstruction to a favored (local) energy minimum. Moreover, our dynamic procedure can be naturally applied to reconstruct an actual microstructure evolution process. In contrast to commonly used evolution reconstruction approaches that separately generate individual static configurations, our procedure continuously evolves a single microstructure according to a time-dependent correlation function. The utility of our procedure is illustrated by successfully reconstructing nearly percolating hard-sphere packings and particle-reinforced composites as well as the coarsening process in a binary metallic alloy.

Universal spatial correlation functions for describing and reconstructing soil microstructure

Structural features of porous materials such as soil define the majority of its physical properties, including water infiltration and redistribution, multi-phase flow (e.g. simultaneous water/air flow, or gas exchange between biologically active soil root zone and atmosphere) and solute transport. To characterize soil microstructure, conventional soil science uses such metrics as pore size and pore-size distributions and thin section-derived morphological indicators. However, these descriptors provide only limited amount of information about the complex arrangement of soil structure and have limited capability to reconstruct structural features or predict physical properties. We introduce three different spatial correlation functions as a comprehensive tool to characterize soil microstructure: 1) two-point probability functions, 2) linear functions, and 3) two-point cluster functions. This novel approach was tested on thin-sections (2.21×2.21 cm2) representing eight soils with different pore space configurations. The two-point probability and linear correlation functions were subsequently used as a part of simulated annealing optimization procedures to reconstruct soil structure. Comparison of original and reconstructed images was based on morphological characteristics, cluster correlation functions, total number of pores and pore-size distribution. Results showed excellent agreement for soils with isolated pores, but relatively poor correspondence for soils exhibiting dual-porosity features (i.e. superposition of pores and micro-cracks). Insufficient information content in the correlation function sets used for reconstruction may have contributed to the observed discrepancies. Improved reconstructions may be obtained by adding cluster and other correlation functions into reconstruction sets. Correlation functions and the associated stochastic reconstruction algorithms introduced here are universally applicable in soil science, such as for soil classification, pore-scale modelling of soil properties, soil degradation monitoring, and description of spatial dynamics of soil microbial activity.

Ab-initio reconstruction of complex Euclidean networks in two dimensions

Reconstruction of complex structures is an inverse problem arising in virtually all areas of science and technology, from protein structure determination to bulk heterostructure solar cells and the structure of nanoparticles. We cast this problem as a complex network problem where the edges in a network have weights equal to the Euclidean distance between their endpoints. We present a method for reconstruction of the locations of the nodes of the network given only the edge weights of the Euclidean network. The theoretical foundations of the method are based on rigidity theory, which enables derivation of a polynomial bound on its efficiency. An efficient implementation of the method is discussed and timing results indicate that the run time of the algorithm is polynomial in the number of nodes in the network. We have reconstructed Euclidean networks of about 1000 nodes in approximately 24 h on a desktop computer using this implementation. We also reconstruct Euclidean networks corresponding to polymer chains in two dimensions and planar graphene nanoparticles. We have also modified our base algorithm so that it can successfully solve random point sets when the input data are less precise.

Stable-phase method for hierarchical annealing in the reconstruction of porous media images

In this paper, we introduce a stable-phase approach for hierarchical annealing which addresses the very large computational costs associated with simulated annealing for the reconstruction of large-scale binary porous media images. Our presented method, which uses the two-point correlation function as the morphological descriptor, involves the reconstruction of three-phase and two-phase structures. We consider reconstructing the three-phase structures based on standard annealing and the two-phase structures based on standard and hierarchical annealings. From the result of the two-dimensional (2D) reconstruction, we find that the 2D generation does not fully capture the morphological information of the original image, even though the two-point correlation function of the reconstruction is in excellent agreement with that of the reference image. For the reconstructed three-dimensional (3D) microstructure, we calculate its permeability and compare it to that of the reference 3D microstructure. The result indicates that the reconstructed structure has a lower degree of connectedness than that of the actual sandstone. We also compare the computation time of our presented method to that of the standard annealing, which shows that our presented method of orders of magnitude improves the convergence rate. That is because only a small part of the pixels in the overall hierarchy need to be considered for sampling by the annealer.

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.

Improved reconstructions of random media using dilation and erosion processes

By using the most sensitive two-point correlation functions introduced to date, we reconstruct the microstructures of two-phase random media with heretofore unattained accuracy. Such media arise in a host of contexts, including porous and composite media, ecological structures, biological media, and astrophysical structures. The aforementioned correlation functions are special cases of the so-called canonical n-point correlation function H(n) and generalize the ones that have been recently employed to advance our ability to reconstruct complex microstructures [Y. Jiao, F. H. Stillinger, and S. Torquato, Proc. Natl. Acad. Sci. 106, 17634 (2009)]. The use of these generalized correlation functions is tantamount to dilating or eroding a reference phase of the target medium and incorporating the additional topological information of the modified media, thereby providing more accurate reconstructions of percolating, filamentary, and other topologically complex microstructures. We apply our methods to a multiply connected "donut" medium and a dilute distribution of "cracks" (a set of essentially zero measure), demonstrating improved accuracy in both cases with implications for higher-dimensional and biconnected two-phase systems. The high information content of the generalized two-point correlation functions suggests that it would be profitable to explore their use to characterize the structural and physical properties of not only random media, but also molecular systems, including structural glasses.

Hyperuniformity, quasi-long-range correlations, and void-space constraints in maximally random jammed particle packings. II. Anisotropy in particle shape

We extend the results from the first part of this series of two papers by examining hyperuniformity in heterogeneous media composed of impenetrable anisotropic inclusions. Specifically, we consider maximally random jammed (MRJ) packings of hard ellipses and superdisks and show that these systems both possess vanishing infinite-wavelength local-volume-fraction fluctuations and quasi-long-range pair correlations scaling as r(-(d+1)) in d Euclidean dimensions. Our results suggest a strong generalization of a conjecture by Torquato and Stillinger [Phys. Rev. E 68, 041113 (2003)], namely, that all strictly jammed saturated packings of hard particles, including those with size and shape distributions, are hyperuniform with signature quasi-long-range correlations. We show that our arguments concerning the constrained distribution of the void space in MRJ packings directly extend to hard-ellipse and superdisk packings, thereby providing a direct structural explanation for the appearance of hyperuniformity and quasi-long-range correlations in these systems. Additionally, we examine general heterogeneous media with anisotropic inclusions and show unexpectedly that one can decorate a periodic point pattern to obtain a hard-particle system that is not hyperuniform with respect to local-volume-fraction fluctuations. This apparent discrepancy can also be rationalized by appealing to the irregular distribution of the void space arising from the anisotropic shapes of the particles. Our work suggests the intriguing possibility that the MRJ states of hard particles share certain universal features independent of the local properties of the packings, including the packing fraction and average contact number per particle.