Precise algorithms to compute surface correlation functions of two-phase heterogeneous media and their applications

The quantitative characterization of the microstructure of random heterogeneous media in $d$-dimensional Euclidean space $\mathbb{R}^d$ via a variety of $n$-point correlation functions is of great importance, since the respective infinite set determines the effective physical properties of the media. In particular, surface-surface $F_{ss}$ and surface-void $F_{sv}$ correlation functions (obtainable from radiation scattering experiments) contain crucial interfacial information that enables one to estimate transport properties of the media (e.g., the mean survival time and fluid permeability) and complements the information content of the conventional two-point correlation function. However, the current technical difficulty involved in sampling surface correlation functions has been a stumbling block in their widespread use. We first present a concise derivation of the small-$r$ behaviors of these functions, which are linked to the \textit{mean curvature} of the system. Then we demonstrate that one can reduce the computational complexity of the problem by extracting the necessary interfacial information from a cut of the system with an infinitely long line. Accordingly, we devise algorithms based on this idea and test them for two-phase media in continuous and discrete spaces. Specifically for the exact benchmark model of overlapping spheres, we find excellent agreement between numerical and exact results. We compute surface correlation functions and corresponding local surface-area variances for a variety of other model microstructures, including hard spheres in equilibrium, decorated "stealthy" patterns, as well as snapshots of evolving pattern formation processes (e.g., spinodal decomposition). It is demonstrated that the precise determination of surface correlation functions provides a powerful means to characterize a wide class of complex multiphase microstructures

Random Scalar Fields and Hyperuniformity

Disordered many-particle hyperuniform systems are exotic amorphous states of matter that lie between crystals and liquids. Hyperuniform systems have attracted recent attention because they are endowed with novel transport and optical properties. Recently, the hyperuniformity concept has been generalized to characterize scalar fields, two-phase media and random vector fields. In this paper, we devise methods to explicitly construct hyperuniform scalar fields. We investigate explicitly spatial patterns generated from Gaussian random fields, which have been used to model the microwave background radiation and heterogeneous materials, the Cahn-Hilliard equation for spinodal decomposition, and Swift-Hohenberg equations that have been used to model emergent pattern formation, including Rayleigh-B{\' e}nard convection. We show that the Gaussian random scalar fields can be constructed to be hyperuniform. We also numerically study the time evolution of spinodal decomposition patterns and demonstrate that these patterns are hyperuniform in the scaling regime. Moreover, we find that labyrinth-like patterns generated by the Swift-Hohenberg equation are effectively hyperuniform. We show that thresholding a hyperuniform Gaussian random field to produce a two-phase random medium tends to destroy the hyperuniformity of the progenitor scalar field. We then propose guidelines to achieve effectively hyperuniform two-phase media derived from thresholded non-Gaussian fields. Our investigation paves the way for new research directions to characterize the large-structure spatial patterns that arise in physics, chemistry, biology and ecology. Moreover, our theoretical results are expected to guide experimentalists to synthesize new classes of hyperuniform materials with novel physical properties via coarsening processes and using state-of-the-art techniques, such as stereolithography and 3D printing.Comment: 16 pages, 18 figure