Showing 1–50 of 189 results for author: Torquato, S

Hole Statistics of Equilibrium 2D and 3D Hard-Sphere Crystals

Authors: Haina Wang, David A. Huse, Salvatore Torquato

Abstract: The probability of finding a spherical "hole" of a given radius $r$ contains crucial structural information about many-body systems. Such hole statistics, including the void conditional nearest-neighbor probability functions $G_V(r)$, have been well studied for hard-sphere fluids in $d$-dimensional Euclidean space $\mathbb{R}^d$. However, little is known about these functions for hard-sphere cryst… ▽ More The probability of finding a spherical "hole" of a given radius $r$ contains crucial structural information about many-body systems. Such hole statistics, including the void conditional nearest-neighbor probability functions $G_V(r)$, have been well studied for hard-sphere fluids in $d$-dimensional Euclidean space $\mathbb{R}^d$. However, little is known about these functions for hard-sphere crystals for values of $r$ beyond the hard-sphere diameter, as large holes are extremely rare in crystal phases. To overcome these computational challenges, we introduce a biased-sampling scheme that accurately determines hole statistics for equilibrium hard spheres on ranges of $r$ that far extend those that could be previously explored. We discover that $G_V(r)$ in crystal and hexatic states exhibits oscillations whose amplitudes increase rapidly with the packing fraction, which stands in contrast to $G_V(r)$ for fluid states that is monotonic in $r$. The oscillations in $G_V(r)$ for 2D crystals are strongly correlated with the local orientational order metric in the vicinity of the holes, and variations in $G_V(r)$ for 3D states indicate a transition between tetrahedral and octahedral holes, demonstrating the power of $G_V(r)$ as a probe of local coordination geometry. To further study the statistics of interparticle spacing in hard-sphere systems, we compute the local packing fraction distribution $f(φ_l)$ of Delaunay cells, and find that for $d\leq 3$, the excess kurtosis of $f(φ_l)$ switches sign at a certain transitional global packing fraction. Our investigation facilitates the study of structural and bulk properties of materials that involve the creation of rare large holes, such as the solubility of alloys. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Effective medium properties of stealthy hyperuniform photonic structures using multiscale physics-informed neural networks

Authors: Roberto Riganti, Yilin Zhu, Wei Cai, Salvatore Torquato, Luca Dal Negro

Abstract: In this article, we employ multiscale physics-informed neural networks (MscalePINNs) for the inverse retrieval of the effective permittivity and homogenization of finite-size photonic media with stealthy hyperuniform (SHU) disordered geometries. Specifically, we show that MscalePINNs are capable of capturing the fast spatial variations of complex fields scattered by arrays of dielectric nanocylind… ▽ More In this article, we employ multiscale physics-informed neural networks (MscalePINNs) for the inverse retrieval of the effective permittivity and homogenization of finite-size photonic media with stealthy hyperuniform (SHU) disordered geometries. Specifically, we show that MscalePINNs are capable of capturing the fast spatial variations of complex fields scattered by arrays of dielectric nanocylinders arranged according to isotropic SHU point patterns, thus enabling a systematic methodology to inverse retrieve their effective dielectric profiles. Our approach extends the recently developed high-frequency homogenization theory of hyperuniform media and retrieves more general permittivity profiles for applications-relevant finite-size SHU systems, unveiling unique features related to their isotropic nature. In particular, we demonstrate the existence of a transparency region beyond the long-wavelength approximation, enabling effective and isotropic homogenization even without disorder-averaging, in contrast to the case of uncorrelated Poisson random patterns. We believe that the multiscale network approach introduced here enables the efficient inverse design of general effective media and finite-size metamaterials with isotropic electromagnetic responses beyond the limitations of traditional homogenization theories. △ Less

Submitted 14 May, 2024; v1 submitted 13 May, 2024; originally announced May 2024.

Hyperuniformity Classes of Quasiperiodic Tilings via Diffusion Spreadability

Authors: Adam Hitin-Bialus, Charles Emmett Maher, Paul J. Steinhardt, Salvatore Torquato

Abstract: Hyperuniform point patterns can be classified by the hyperuniformity scaling exponent $α> 0$, that characterizes the power-law scaling behavior of the structure factor $S(\mathbf{k})$ as a function of wavenumber $k\equiv|\mathbf{k}|$ in the vicinity of the origin, e.g., $S(\mathbf{k})\sim|\mathbf{k}|^α$ in cases where $S(\mathbf{k})$ varies continuously with $k$ as $k\rightarrow0$. In this paper,… ▽ More Hyperuniform point patterns can be classified by the hyperuniformity scaling exponent $α> 0$, that characterizes the power-law scaling behavior of the structure factor $S(\mathbf{k})$ as a function of wavenumber $k\equiv|\mathbf{k}|$ in the vicinity of the origin, e.g., $S(\mathbf{k})\sim|\mathbf{k}|^α$ in cases where $S(\mathbf{k})$ varies continuously with $k$ as $k\rightarrow0$. In this paper, we show that the spreadability is an effective method for determining $α$ for quasiperiodic systems where $S(\mathbf{k})$ is discontinuous and consists of a dense set of Bragg peaks. We first transform quasiperiodic and limit-periodic point patterns into two-phase media by mapping them onto packings of identical nonoverlapping disks, where space interior to the disks represents one phase and the space in exterior to them represents the second phase. We then compute the spectral density of the packings, and finally compute and fit the long-time behavior of their excess spreadabilities. Specifically, we show that the excess spreadability can be used to accurately extract $α$ for the 1D limit-periodic period doubling chain and the 1D quasicrystalline Fibonacci chain to within $0.02\%$ of the analytically known exact results. Moreover, we obtain a value of $α= 5.97\pm0.06$ for the 2D Penrose tiling, which had not been computed previously. We also show that one can truncate the small-$k$ region of the scattering information used to compute the spreadability and still obtain an accurate value of $α$. The methods described here offer a simple way to characterize the large-scale translational order present in quasicrystalline and limit-periodic media in any space dimension that are self-similar. Moreover, the scattering information extracted from these two-phase media encoded in the spectral density can be used to estimate their physical properties. (abridged) △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: 15 pages, 11 figures

Local number fluctuations in ordered and disordered phases of water across temperatures: Higher-order moments and degrees of tetrahedrality

Authors: Michael A. Klatt, Jaeuk Kim, Thomas E. Gartner III, Salvatore Torquato

Abstract: The isothermal compressibility (i.e., the asymptotic number variance) of equilibrium liquid water as a function of temperature is minimal near ambient conditions. This anomalous non-monotonic temperature dependence is due to a balance between thermal fluctuations and the formation of tetrahedral hydrogen-bond networks. Since tetrahedrality is a many-body property, it will also influence the higher… ▽ More The isothermal compressibility (i.e., the asymptotic number variance) of equilibrium liquid water as a function of temperature is minimal near ambient conditions. This anomalous non-monotonic temperature dependence is due to a balance between thermal fluctuations and the formation of tetrahedral hydrogen-bond networks. Since tetrahedrality is a many-body property, it will also influence the higher-order moments of density fluctuations, including the skewness and kurtosis. To gain a more complete picture, we examine these higher-order moments that encapsulate many-body correlations using a recently developed, advanced platform for local density fluctuations. We study an extensive set of simulated phases of water across a range of temperatures (80 K to 1600 K) with various degrees of tetrahedrality, including ice phases, equilibrium liquid water, supercritical water, and disordered nonequilibrium quenches. We find clear signatures of tetrahedrality in the higher-order moments, including the skewness and excess kurtosis, that scale for all cases with the degree of tetrahedrality. More importantly, this scaling behavior leads to non-monotonic temperature dependencies in the higher-order moments for both equilibrium and non-equilibrium phases. Specifically, at near-ambient conditions, the higher-order moments vanish most rapidly for large length scales, and the distribution quickly converges to a Gaussian in our metric. However, at non-ambient conditions, higher-order moments vanish more slowly and hence become more relevant especially for improving information-theoretic approximations of hydrophobic solubility. The temperature non-monotonicity that we observe in the full distribution across length-scales could shed light on water's nested anomalies, i.e., reveal new links between structural, dynamic, and thermodynamic anomalies. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: 14 pages, 8 figures

Microstructural and Transport Characteristics of Triply Periodic Bicontinuous Materials

Authors: Salvatore Torquato, Jaeuk Kim

Abstract: 3D bicontinuous two-phase materials are increasingly gaining interest because of their unique multifunctional characteristics and advancements in techniques to fabricate them. Due to their complex topological and structural properties, it still has been nontrivial to develop explicit microstructure-dependent formulas to predict accurately their physical properties. A primary goal of the present pa… ▽ More 3D bicontinuous two-phase materials are increasingly gaining interest because of their unique multifunctional characteristics and advancements in techniques to fabricate them. Due to their complex topological and structural properties, it still has been nontrivial to develop explicit microstructure-dependent formulas to predict accurately their physical properties. A primary goal of the present paper is to ascertain various microstructural and transport characteristics of five different models of triply periodic bicontinuous porous materials at a porosity $φ_1=1/2$: those in which the two-phase interfaces are the Schwarz P, Schwarz D and Schoen G minimal surfaces as well as two different pore-channel structures. We ascertain their spectral densities, pore-size distribution functions, local volume-fraction variances, and hyperuniformity order metrics and then use this information to estimate certain effective transport properties via closed-form microstructure-property formulas. Specifically, we estimate the recently introduced time-dependent diffusion spreadability exactly from the spectral density. Moreover, we accurately estimate the fluid permeability of such porous materials from the second moment of the pore-size function and the formation factor, a measure of the tortuosity of the pore space. We also rigorously bound the permeability from above using the spectral density. For the five models with identical cubic unit cells, we find that the permeability, inverse of the specific surface, hyperuniformity order metric, pore-size second moment and long-time spreadability behavior are all positively correlated and rank order the structures in exactly the same way. We also conjecture what structures maximize the fluid permeability for arbitrary porosities and show that this conjecture must be true in the extreme porosity limits by identifying the corresponding optimal structures. △ Less

Submitted 8 July, 2024; v1 submitted 25 April, 2024; originally announced May 2024.

Comments: 16 pages, 6 figures, and 5 tables

Journal ref: Acta Materialia 276 (2024) 120142

Ordered and disordered stealthy hyperuniform point patterns across spatial dimensions

Authors: Peter K. Morse, Paul J. Steinhardt, Salvatore Torquato

Abstract: In previous work [Phys. Rev. X 5, 021020 (2015)], it was shown that stealthy hyperuniform systems can be regarded as hard spheres in Fourier-space in the sense that the the structure factor is exactly zero in a spherical region around the origin in analogy with the pair-correlation function of real-space hard spheres. In this work, we exploit this correspondence to confirm that the densest Fourier… ▽ More In previous work [Phys. Rev. X 5, 021020 (2015)], it was shown that stealthy hyperuniform systems can be regarded as hard spheres in Fourier-space in the sense that the the structure factor is exactly zero in a spherical region around the origin in analogy with the pair-correlation function of real-space hard spheres. In this work, we exploit this correspondence to confirm that the densest Fourier-space hard-sphere system is that of a Bravais lattice. This is in contrast to real-space hard-spheres, whose densest configuration is conjectured to be disordered. We also extend the virial series previously suggested for disordered stealthy hyperuniform systems to higher dimensions in order to predict spatial decorrelation as function of dimension. This prediction is then borne out by numerical simulations of disordered stealthy hyperuniform ground states in dimensions $d=2$-$8$. △ Less

Submitted 28 April, 2024; v1 submitted 25 April, 2024; originally announced April 2024.

Comments: 9 pages; 6 figures

Designer Pair Statistics of Disordered Many-Particle Systems with Novel Properties

Authors: Haina Wang, Salvatore Torquato

Abstract: Knowledge of exact analytical functional forms for the pair correlation function $g_2(r)$ and its corresponding structure factor $S(k)$ of disordered many-particle systems is limited. For fundamental and practical reasons, it is highly desirable to add to the existing data base of analytical functional forms for such pair statistics. Here, we design a plethora of such pair functions in direct and… ▽ More Knowledge of exact analytical functional forms for the pair correlation function $g_2(r)$ and its corresponding structure factor $S(k)$ of disordered many-particle systems is limited. For fundamental and practical reasons, it is highly desirable to add to the existing data base of analytical functional forms for such pair statistics. Here, we design a plethora of such pair functions in direct and Fourier spaces across the first three Euclidean space dimensions that are realizable by diverse many-particle systems with varying degrees of correlated disorder across length scales, spanning a wide spectrum of hyperuniform, typical nonhyperuniform and antihyperuniform ones. This is accomplished by utilizing an efficient inverse algorithm that determines equilibrium states with up to pair interactions at positive temperature that precisely match targeted forms for both $g_2(r)$ and $S(k)$. Among other results, we realize an example with the strongest hyperuniform property among known positive-temperature equilibrium states, critical-point systems (implying unusual 1D systems with phase transitions) that are not in the Ising universality class, systems that attain self-similar pair statistics under Fourier transformation, and an experimentally feasible polymer model. We show that our pair functions enable one to achieve systems with a wide range of translational order and self-diffusion coefficients $\cal D$, which are inversely related to one another. One can design other realizable pair statistics via linear combinations of our functions or by applying our inverse procedure to other desirable functional forms. Our approach facilitates the inverse design of materials with desirable physical and chemical properties by tuning their pair statistics. △ Less

Submitted 29 December, 2023; originally announced January 2024.

Theoretical Prediction of the Effective Dynamic Dielectric Constant of Disordered Hyperuniform Anisotropic Composites Beyond the Long-Wavelength Regime

Authors: Jaeuk Kim, Salvatore Torquato

Abstract: Torquato and Kim [Phys. Rev. X 11, 296 021002 (2021)] derived exact nonlocal strong-contrast expansions of the effective dynamic dielectric constant tensor that treat general three-dimensional (3D) two-phase composites, which are valid well beyond the long-wavelength regime. Here, we demonstrate that truncating this general rapidly converging series at the two- and three-point levels is a powerful… ▽ More Torquato and Kim [Phys. Rev. X 11, 296 021002 (2021)] derived exact nonlocal strong-contrast expansions of the effective dynamic dielectric constant tensor that treat general three-dimensional (3D) two-phase composites, which are valid well beyond the long-wavelength regime. Here, we demonstrate that truncating this general rapidly converging series at the two- and three-point levels is a powerful theoretical tool for extracting accurate approximations suited for various microstructural symmetries. We derive such closed-form formulas applicable to transverse polarization in layered media and transverse magnetic polarization in transversely isotropic media, respectively. We use these formulas to estimate effective dielectric constant for models of 3D disordered hyperuniform layered and transversely isotropic media: nonstealthy hyperuniform and stealthy hyperuniform (SHU) media. In particular, we show that SHU media are perfectly transparent (trivially implying no Anderson localization, in principle) within finite wave number intervals through the third-order terms. For these two models, we validate that the second-order formulas, which depend on the spectral density, are already very accurate well beyond the long-wavelength regime by showing very good agreement with the finite-difference time-domain simulations. The high predictive power of the second-order formulas implies that higher-order contributions are negligibly small, and thus, it very accurately approximates multiple scattering effects. Therefore, there can be no Anderson localization in practice within the predicted perfect transparency interval in SHU media because the localization length should be very large compared to any practically large sample size. Our predictive theory provides a foundation for the inverse design of novel effective wave characteristics of disordered and statistically anisotropic structures. △ Less

Submitted 22 December, 2023; v1 submitted 8 December, 2023; originally announced December 2023.

Comments: 16 pages, 7 figures

Journal ref: Opt. Mater. Express 14, 194-215 (2024)

Hyperuniformity of Maximally Random Jammed Packings of Hyperspheres Across Spatial Dimensions

Authors: Charles Emmett Maher, Yang Jiao, Salvatore Torquato

Abstract: The maximally random jammed (MRJ) state is the most random configuration of strictly jammed (mechanically rigid) nonoverlapping objects. MRJ packings are hyperuniform, meaning their long-wavelength density fluctuations are anomalously suppressed compared to typical disordered systems, i.e., their structure factors $S(\mathbf{k})$ tend to zero as the wavenumber $|\mathbf{k}|$ tends to zero. Here, w… ▽ More The maximally random jammed (MRJ) state is the most random configuration of strictly jammed (mechanically rigid) nonoverlapping objects. MRJ packings are hyperuniform, meaning their long-wavelength density fluctuations are anomalously suppressed compared to typical disordered systems, i.e., their structure factors $S(\mathbf{k})$ tend to zero as the wavenumber $|\mathbf{k}|$ tends to zero. Here, we show that generating high-quality strictly jammed states for space dimensions $d = 3,4,$ and $5$ is of paramount importance in ensuring hyperuniformity and extracting precise values of the hyperuniformity exponent $α> 0$ for MRJ states, defined by the power-law behavior of $S(\mathbf{k})\sim|\mathbf{k}|^α$ in the limit $|\mathbf{k}|\rightarrow0$. Moreover, we show that for fixed $d$ it is more difficult to ensure jamming as the particle number $N$ increases, which results in packings that are nonhyperuniform. Free-volume theory arguments suggest that the ideal MRJ state does not contain rattlers, which act as defects in numerically generated packings. As $d$ increases, we find that the fraction of rattlers decreases substantially. Our analysis of the largest truly jammed packings suggests that the ideal MRJ packings for all dimensions $d\geq3$ are hyperuniform with $α= d - 2$, implying the packings become more hyperuniform as $d$ increases. The differences in $α$ between MRJ packings and recently proposed Manna-class random close packed (RCP) states, which were reported to have $α= 0.25$ in $d=3$ and be nonhyperuniform ($α= 0$) for $d = 4$ and $d = 5$, demonstrate the vivid distinctions between the large-scale structure of RCP and MRJ states in these dimensions. Our work clarifies the importance of the link between true jamming and hyperuniformity and motivates the development of an algorithm to produce rattler-free three-dimensional MRJ packings. △ Less

Submitted 8 December, 2023; v1 submitted 12 November, 2023; originally announced November 2023.

Comments: 15 pages, 8 figures

Journal ref: Physical Review E, 108, 064602 (2023)

Computational Design of Anisotropic Stealthy Hyperuniform Composites with Engineered Directional Scattering Properties

Authors: Wenlong Shi, David Keeney, Duyu Chen, Yang Jiao, Salvatore Torquato

Abstract: Disordered hyperuniform materials are an emerging class of exotic amorphous states of matter that endow them with singular physical properties. Here, we generalize the Fourier-space based numerical construction procedure for designing {\it isotropic} disordered hyperuniform two-phase heterogeneous materials (i.e., composites) developed by Chen and Torquato [Acta Mater. {\bf 142}, 152 (2018)] to {\… ▽ More Disordered hyperuniform materials are an emerging class of exotic amorphous states of matter that endow them with singular physical properties. Here, we generalize the Fourier-space based numerical construction procedure for designing {\it isotropic} disordered hyperuniform two-phase heterogeneous materials (i.e., composites) developed by Chen and Torquato [Acta Mater. {\bf 142}, 152 (2018)] to {\it anisotropic} microstructures by explicitly incorporating the {\it vector-dependent} spectral density function ${\tilde χ}_{_V}({\bf k})$ of {\it arbitrary form} that is realizable. We demonstrate the utility of the procedure by generating a wide spectrum of {\it anisotropic} stealthy hyperuniform (SHU) microstructures with ${\tilde χ}_{_V}({\bf k}) = 0$ for ${\bf k} \in Ω$. We show how different exclusion-region shapes with various discrete symmetries and varying size affect the resulting statistically anisotropic microstructures as a function of the and phase volume fraction. We find that, among other properties, the directional hyperuniform behaviors imposed by the shape asymmetry (or anisotropy) of certain exclusion regions give rise to distinct anisotropic structures and degree of uniformity in the distribution of the phases on intermediate and large length scales along different directions. Moreover, while the anisotropic exclusion regions impose strong constraints on the {\it global} symmetry of the resulting media, they can still possess almost isotropic {\it local} structures. Our construction algorithm enables one to control the statistical anisotropy of composite microstructures which is crucial to engineering directional optical, transport and mechanical properties of two-phase composite media. △ Less

Submitted 25 July, 2023; originally announced July 2023.

Comments: 22 figures, 24 pages

Effective Electromagnetic Wave Properties of Disordered Stealthy Hyperuniform Layered Media Beyond the Quasistatic Regime

Authors: Jaeuk Kim, Salvatore Torquato

Abstract: Disordered stealthy hyperuniform dielectric composites exhibit novel electromagnetic wave transport properties in two and three dimensions. Here, we carry out the first study of the electromagnetic properties of one-dimensional (1D) disordered stealthy hyperuniform layered media. From an exact nonlocal theory, we derive an approximation formula for the effective dynamic dielectric constant tensor… ▽ More Disordered stealthy hyperuniform dielectric composites exhibit novel electromagnetic wave transport properties in two and three dimensions. Here, we carry out the first study of the electromagnetic properties of one-dimensional (1D) disordered stealthy hyperuniform layered media. From an exact nonlocal theory, we derive an approximation formula for the effective dynamic dielectric constant tensor ${\boldsymbol \varepsilon}_e({\bf k}_q,ω)$ of general 1D media that is valid well beyond the quasistatic regime and apply it to 1D stealthy hyperuniform systems. We consider incident waves of transverse polarization, frequency $ω$, and wavenumber $k_q$. Our formula for ${\boldsymbol \varepsilon}_e({k}_q,ω)$, which is given in terms of the spectral density, leads to a closed-form relation for the transmittance $T$. Our theoretical predictions are in excellent agreement with finite-difference time-domain (FDTD) simulations. Stealthy hyperuniform layered media have perfect transparency intervals up to a finite wavenumber, implying no Anderson localization, but non-stealthy hyperuniform media are not perfectly transparent. Our predictive theory provides a new path for the inverse design of the wave characteristics of disordered layered media, which are readily fabricated, by engineering their spectral densities. △ Less

Submitted 31 July, 2023; v1 submitted 22 May, 2023; originally announced May 2023.

Comments: 8 pages, 6 figures

Journal ref: Optica Vol. 10, Issue 8, pp. 965-972 (2023)

Generating large disordered stealthy hyperuniform systems with ultra-high accuracy to determine their physical properties

Authors: Peter K. Morse, Jaeuk Kim, Paul J. Steinhardt, Salvatore Torquato

Abstract: Hyperuniform many-particle systems are characterized by a structure factor $S({\mathbf{k}})$ that is precisely zero as $|\mathbf{k}|\rightarrow0$; and stealthy hyperuniform systems have $S({\mathbf{k}})=0$ for the finite range $0 < |{\mathbf{k}}| \le K$, called the "exclusion region." Through a process of collective-coordinate optimization, energy-minimizing disordered stealthy hyperuniform system… ▽ More Hyperuniform many-particle systems are characterized by a structure factor $S({\mathbf{k}})$ that is precisely zero as $|\mathbf{k}|\rightarrow0$; and stealthy hyperuniform systems have $S({\mathbf{k}})=0$ for the finite range $0 < |{\mathbf{k}}| \le K$, called the "exclusion region." Through a process of collective-coordinate optimization, energy-minimizing disordered stealthy hyperuniform systems of moderate size have been made to high accuracy, and their novel physical properties have shown great promise. However, minimizing $S(\mathbf{k})$ in the exclusion region is computationally intensive as the system size becomes large. In this Letter, we present an improved methodology to generate such states using double-double precision calculations on GPUs that reduces the deviations from zero within the exclusion region by a factor of approximately $10^{30}$ for systems sizes more than an order of magnitude larger. We further show that this ultra-high accuracy is required to draw conclusions about their corresponding characteristics, such as the nature of the associated energy landscape and the presence or absence of Anderson localization, which might be masked, even when deviations are relatively small. △ Less

Submitted 18 April, 2023; originally announced April 2023.

Comments: 7 pages, 3 figures

Realizability of Iso-$g_2$ Processes via Effective Pair Interactions

Authors: Haina Wang, Frank H. Stillinger, Salvatore Torquato

Abstract: An outstanding problem in statistical mechanics is the determination of whether prescribed functional forms of the pair correlation function $g_2(r)$ [or equivalently, structure factor $S(k)$] at some number density $ρ$ can be achieved by $d$-dimensional many-body systems. The Zhang-Torquato conjecture states that any realizable set of pair statistics, whether from a nonequilibrium or equilibrium… ▽ More An outstanding problem in statistical mechanics is the determination of whether prescribed functional forms of the pair correlation function $g_2(r)$ [or equivalently, structure factor $S(k)$] at some number density $ρ$ can be achieved by $d$-dimensional many-body systems. The Zhang-Torquato conjecture states that any realizable set of pair statistics, whether from a nonequilibrium or equilibrium system, can be achieved by equilibrium systems involving up to two-body interactions. In light of this conjecture, we study the realizability problem of the nonequilibrium iso-$g_2$ process, i.e., the determination of density-dependent effective potentials that yield equilibrium states in which $g_2$ remains invariant for a positive range of densities. Using a precise inverse methodology that determines effective potentials that match hypothesized functional forms of $g_2(r)$ for all $r$ and $S(k)$ for all $k$, we show that the unit-step function $g_2$, which is the zero-density limit of the hard-sphere potential, is remarkably numerically realizable up to the packing fraction $φ=0.49$ for $d=1$. For $d=2$ and 3, it is realizable up to the maximum ``terminal'' packing fraction $φ_c=1/2^d$, at which the systems are hyperuniform, implying that the explicitly known necessary conditions for realizability are sufficient up through $φ_c$. For $φ$ near but below $φ_c$, the large-$r$ behaviors of the effective potentials are given exactly by the functional forms $\exp[-κ(φ) r]$ for $d=1$, $r^{-1/2}\exp[-κ(φ) r]$ for $d=2$, and $r^{-1}\exp[-κ(φ) r]$ (Yukawa form) for $d=3$, where $κ^{-1}(φ)$ is a screening length, and for $φ=φ_c$, the potentials at large $r$ are given by the pure Coulomb forms in the respective dimensions, as predicted by Torquato and Stillinger [\textit{Phys. Rev. E}, 68, 041113 1-25 (2003)]. △ Less

Submitted 24 November, 2022; originally announced November 2022.

Equilibrium States Corresponding to Targeted Hyperuniform Nonequilibrium Pair Statistics

Authors: Haina Wang, Salvatore Torquato

Abstract: The Zhang-Torquato conjecture [Phys. Rev. E 101, 032124 (2020)] states that any realizable pair correlation function $g_2({\bf r})$ or structure factor $S({\bf k})$ of a translationally invariant nonequilibrium system can be attained by an equilibrium ensemble involving only (up to) effective two-body interactions. To test this conjecture, we consider two singular nonequilibrium models of recent i… ▽ More The Zhang-Torquato conjecture [Phys. Rev. E 101, 032124 (2020)] states that any realizable pair correlation function $g_2({\bf r})$ or structure factor $S({\bf k})$ of a translationally invariant nonequilibrium system can be attained by an equilibrium ensemble involving only (up to) effective two-body interactions. To test 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 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, provides measures of the degree to which a system is out of equilibrium. We show that all three systems studied possess the \textit{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 with stronger forms of hyperuniformity compared to their target counterparts. Our work is expected to facilitate the self-assembly of tunable hyperuniform soft-matter systems. △ Less

Submitted 19 December, 2022; v1 submitted 12 October, 2022; originally announced October 2022.

Precise Determination of Pair Interactions from Pair Statistics of Many-Body Systems In and Out of Equilibrium

Authors: Salvatore Torquato, Haina Wang

Abstract: The determination of the pair potential $v({\bf r})$ that accurately yields an equilibrium state at positive temperature $T$ with a prescribed pair correlation function $g_2({\bf r})$ or corresponding structure factor $S({\bf k})$ in $d$-dimensional Euclidean space $\mathbb{R}^d$ is an outstanding inverse statistical mechanics problem with far-reaching implications. Recently, Zhang and Torquato co… ▽ More The determination of the pair potential $v({\bf r})$ that accurately yields an equilibrium state at positive temperature $T$ with a prescribed pair correlation function $g_2({\bf r})$ or corresponding structure factor $S({\bf k})$ in $d$-dimensional Euclidean space $\mathbb{R}^d$ is an outstanding inverse statistical mechanics problem with far-reaching implications. Recently, Zhang and Torquato conjectured that any realizable $g_2({\bf r})$ or $S({\bf 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 a novel optimization algorithm to find 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- 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({\bf r})$ and structure factor $S({\bf k})$ so that the small- 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. 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. △ Less

Submitted 22 September, 2022; v1 submitted 11 July, 2022; originally announced July 2022.

The Disordered Heterogeneous Universe: Galaxy Distribution and Clustering Across Length Scales

Authors: Oliver H. E. Philcox, Salvatore Torquato

Abstract: Studies of disordered heterogeneous media and galaxy cosmology share a common goal: analyzing the distribution of particles at `microscales' to predict physical properties at `macroscales', whether for a liquid, composite material, or entire Universe. The former theory provides an array of techniques to characterize a wide class of microstructures; in this work, we apply them to the distributions… ▽ More Studies of disordered heterogeneous media and galaxy cosmology share a common goal: analyzing the distribution of particles at `microscales' to predict physical properties at `macroscales', whether for a liquid, composite material, or entire Universe. The former theory provides an array of techniques to characterize a wide class of microstructures; in this work, we apply them to the distributions of galaxies. We focus on the lower-order correlation functions, `void' and `particle' nearest-neighbor functions, pair-connectedness functions, percolation properties, and a scalar order metric. Compared to homogeneous Poisson and typical disordered systems, the cosmological simulations exhibit enhanced large-scale clustering and longer tails in the nearest-neighbor functions, due to the presence of quasi-long-range correlations. On large scales, the system appears `hyperuniform', due to primordial density fluctuations, whilst on the smallest scales, the system becomes almost `antihyperuniform', and, via the order metric, is shown to be a highly correlated disordered system. Via a finite scaling analysis, we show that the percolation threshold of the galaxy catalogs is significantly lower than for Poisson realizations; this is consistent with the observation that the galaxy distribution contains larger voids. However, the two sets of simulations share a fractal dimension, implying that they lie in the same universality class. Finally, we consider the ability of large-scale clustering statistics to constrain cosmological parameters using simulation-based inference. Both the nearest-neighbor distribution and pair-connectedness function considerably tighten bounds on the amplitude of cosmological fluctuations at a level equivalent to observing twenty-five times more galaxies. These are a useful alternative to the three-particle correlation, and are computable in much reduced time. (Abridged) △ Less

Submitted 10 January, 2023; v1 submitted 1 July, 2022; originally announced July 2022.

Comments: 29 pages, accepted by Phys. Rev. X

Local Order Metrics for Two-Phase Media Across Length Scales

Authors: Salvatore Torquato, Murray Skolnick, Jaeuk Kim

Abstract: The capacity to devise order metrics for microstructures of multiphase heterogeneous media is a highly challenging task, given the richness of the possible geometries and topologies of the phases that can arise. This investigation initiates a program to formulate order metrics to characterize the degree of order/disorder of the microstructures of two-phase media in $d$-dimensional Euclidean space… ▽ More The capacity to devise order metrics for microstructures of multiphase heterogeneous media is a highly challenging task, given the richness of the possible geometries and topologies of the phases that can arise. This investigation initiates a program to formulate order metrics to characterize the degree of order/disorder of the microstructures of two-phase media in $d$-dimensional Euclidean space $\mathbb{R}^d$ across length scales. In particular, we propose the use of the local volume-fraction variance $σ^2_{_V}(R)$ associated with a spherical window of radius $R$ as an order metric. We determine $σ^2_{_V}(R)$ as a function of $R$ for 22 different models across the first three space dimensions, including both hyperuniform and nonhyperuniform systems with varying degrees of short- and long-range order. We find that the local volume-fraction variance as well as asymptotic coefficients and integral measures derived from it provide reasonably robust and sensitive order metrics to categorize disordered and ordered two-phase media across all length scales. △ Less

Submitted 27 April, 2022; originally announced April 2022.

Extraordinary Disordered Hyperuniform Multifunctional Composites

Authors: Salvatore Torquato

Abstract: A variety of performance demands are being placed on material systems, including desirable mechanical, thermal, electrical, optical, acoustic and flow properties. The purpose of the present article is to review the emerging field of disordered hyperuniform composites and their novel multifunctional characteristics. Disordered hyperuniform media are exotic amorphous states of matter that are charac… ▽ More A variety of performance demands are being placed on material systems, including desirable mechanical, thermal, electrical, optical, acoustic and flow properties. The purpose of the present article is to review the emerging field of disordered hyperuniform composites and their novel multifunctional characteristics. Disordered hyperuniform media are exotic amorphous states of matter that are characterized by an anomalous suppression of large-scale volume-fraction fluctuations compared to those in "garden-variety" disordered materials. Such unusual composites can have advantages over their periodic counterparts, such as unique or nearly optimal, direction-independent physical properties and robustness against defects. It will be shown that disordered hyperuniform composites and porous media can be endowed with a broad spectrum of extraordinary physical properties, including photonic, phononic, transport, chemical and mechanical characteristics that are only beginning to be discov △ Less

Submitted 24 April, 2022; originally announced April 2022.

Comments: 12 figures

Exclusion Volumes of Convex Bodies in High Space Dimensions: Applications to Virial Coefficients and Continuum Percolation

Authors: Salvatore Torquato, Yang Jiao

Abstract: Using the concepts of mixed volumes and quermassintegrals of convex geometry, we derive an exact formula for the exclusion volume for a general convex body that applies in any space dimension, including both the rotationally-averaged exclusion volume and with the same orientation. We show that the sphere minimizes the dimensionless exclusion volume $v_{ex}(K)/v(K)$ among all convex bodies, whether… ▽ More Using the concepts of mixed volumes and quermassintegrals of convex geometry, we derive an exact formula for the exclusion volume for a general convex body that applies in any space dimension, including both the rotationally-averaged exclusion volume and with the same orientation. We show that the sphere minimizes the dimensionless exclusion volume $v_{ex}(K)/v(K)$ among all convex bodies, whether randomly oriented or uniformly oriented, for any $d$, where $v(K)$ is the volume of $K$. When the bodies have the same orientation, the simplex maximizes the dimensionless exclusion volume for any $d$ with a large-$d$ asymptotic scaling behavior of $2^{2d}/d^{3/2}$, which is to be contrasted with the scaling of $2^d$ for the sphere. We present explicit formulas for quermassintegrals for many nonspherical convex bodies as well as as well as lower-dimensional bodies. These results are utilized to determine the rotationally-averaged exclusion volume for these shapes for dimensions 2 through 12. While the sphere is the shape possessing the minimal dimensionless exclusion volume, among the convex bodies considered that are sufficiently compact, the simplex possesses the maximal dimensionless exclusion volume with a scaling behavior of $2^{1.6618\ldots d}$. We also determine the corresponding second virial coefficient $B_2(K)$ of the aforementioned hard hyperparticles and compute estimates of the continuum percolation threshold $η_c$ derived previously by the authors. We conjecture that overlapping spheres possess the maximal value of $η_c$ among all identical nonzero-volume convex overlapping bodies for $d \ge 2$, randomly or uniformly oriented, and that, among all identical, oriented nonzero-volume convex bodies, overlapping simplices have the minimal value of $η_c$ for $d\ge 2$. △ Less

Submitted 21 September, 2022; v1 submitted 30 March, 2022; originally announced March 2022.

Journal ref: Journal of Statistical Mechanics: Theory and Experiment, (2022) 093404

Dynamic Measure of Hyperuniformity and Nonhyperuniformity in Heterogeneous Media via the Diffusion Spreadability

Authors: Haina Wang, Salvatore Torquato

Abstract: The recently developed concept of spreadability, $\mathcal{S}(t)$, provides a direct link between time-dependent diffusive transport and the microstructure of two-phase media across length scales. We explicitly compute $\mathcal{S}(t)$ for well-known two-dimensional and three-dimensional model structures, including nonhyperuniform fully penetrable spheres and equilibrium hard spheres, as well as h… ▽ More The recently developed concept of spreadability, $\mathcal{S}(t)$, provides a direct link between time-dependent diffusive transport and the microstructure of two-phase media across length scales. We explicitly compute $\mathcal{S}(t)$ for well-known two-dimensional and three-dimensional model structures, including nonhyperuniform fully penetrable spheres and equilibrium hard spheres, as well as hyperuniform sphere packings derived from perfect glasses, uniformly randomized lattices (URL), disordered stealthy point processes and Bravais lattices. We further confirm that the small-, intermediate- and long-time behaviors of $\mathcal{S}(t)$ sensitively capture the small-, intermediate- and large-scale characteristics of the models. In instances in which the spectral density $\tildeχ_{_V}(\mathbf{k})$ has a power-law form $B|\mathbf{k}|^α$ in the limit $|\mathbf{k}|\rightarrow 0$, the long-time spreadability provides a simple means to extract the value of the coefficients $α$ and $B$ that is robust against noise in $\tildeχ_{_V}(\mathbf{k})$ at small wavenumbers. Interestingly, the excess spreadability $\mathcal{S}(\infty)-\mathcal{S}(t)$ for URL packings has nearly exponential decay at small to intermediate $t$, but transforms to a power-law decay at large $t$. Our study of the aforementioned models enables us to devise an algorithm that accurately extracts large-scale behaviors from diffusion data alone. Lessons learned from such analyses of our models are used to determine accurately the large-scale structural characteristics of a sample Fontainebleau sandstone, which we show is nonhyperuniform. Our study demonstrates the practical applicability of the diffusion spreadability to extract crucial microstructural information from real data across length scales and provides a basis for the inverse design of materials with desirable time-dependent diffusion properties. △ Less

Submitted 14 February, 2022; v1 submitted 16 December, 2021; originally announced December 2021.

Comments: 21 pages, 11 figures

Characterization of Void Space, Large-Scale Structure, and Transport Properties of Maximally Random Jammed Packings of Superballs

Authors: Charles Emmett Maher, Frank H. Stillinger, Salvatore Torquato

Abstract: Dense, disordered packings of particles are useful models of low-temperature amorphous phases of matter, biological systems, granular media, and colloidal systems. The study of dense packings of nonspherical particles enables one to ascertain how rotational degrees of freedom affect packing behavior. Here, we study superballs, a large family of deformations of the sphere, defined in three dimensio… ▽ More Dense, disordered packings of particles are useful models of low-temperature amorphous phases of matter, biological systems, granular media, and colloidal systems. The study of dense packings of nonspherical particles enables one to ascertain how rotational degrees of freedom affect packing behavior. Here, we study superballs, a large family of deformations of the sphere, defined in three dimensions by $|x_1|^{2p}+|x_2|^{2p}+|x_3|^{2p}\leq1$, where $p>0$ is a deformation parameter indicating to what extent the shape deviates from a sphere. As $p$ increases from the sphere point ($p=1$), the superball attains cubic symmetry, and attains octahedral symmetry when $p<1$. Previous characterization of superball packings has shown that they have a maximally random jammed (MRJ) state, whose properties (e.g., packing fraction $φ$) vary nonanalytically as $p$ diverges from unity. Here, we use an event-driven molecular dynamics algorithm to produce MRJ superball packings. We characterize their large-scale structure by examining the small-$Q$ behavior of their structure factors and spectral densities, which indicate these packings are effectively hyperuniform. We show that the mean width $\bar{w}$ is a useful length scale to make distances dimensionless in order to compare systematically superballs of different shape. We also compute the complementary cumulative pore-size distribution $F(δ)$ and find the pore sizes tend to decrease as $|1-p|$ increases. From $F(δ)$, we estimate the fluid permeability, mean survival time, and principal diffusion relaxation time of the packings. Additionally, we compute the diffusion spreadability and find the long-time power-law scaling indicates these packings are hyperuniform. Our results can be used to help inform the design of granular materials with targeted densities and transport properties. △ Less

Submitted 6 May, 2022; v1 submitted 29 November, 2021; originally announced November 2021.

Comments: 16 Pages, 17 Figures. SI: 4 Pages, 6 Figures

Journal ref: Physical Review Materials, 6 025603 (2022)

Diffusion Spreadability as a Probe of the Microstructure of Complex Media Across Length Scales

Authors: Salvatore Torquato

Abstract: Consider the time-dependent problem of mass transfer of a solute between two phases and assume that the solute is initially distributed in one phase (phase 2) and absent from the other (phase 1). We desire the fraction of total solute present in phase 1 as a function of time, ${\cal S}(t)$, which we call the {\it spreadability}, since it is a measure of the spreadability of diffusion information a… ▽ More Consider the time-dependent problem of mass transfer of a solute between two phases and assume that the solute is initially distributed in one phase (phase 2) and absent from the other (phase 1). We desire the fraction of total solute present in phase 1 as a function of time, ${\cal S}(t)$, which we call the {\it spreadability}, since it is a measure of the spreadability of diffusion information as a function of time. We derive exact direct-space formulas for ${\cal S}(t)$ in any Euclidean space dimension $d$ in terms of the autocovariance function as well as corresponding Fourier representations of ${\cal S}(t)$ in terms of the spectral density. We derive closed-form general formulas for the short- and long-time behaviors of the spreadability in terms of crucial small- and large-scale microstructural information, respectively. The long-time behavior of ${\cal S}(t)$ enables one to distinguish the entire spectrum of microstructures that span from hyperuniform to nonhyperuniform media. For hyperuniform media, disordered or not, we show that the "excess" spreadability, ${\cal S}(\infty)-{\cal S}(t)$, decays to its long-time behavior exponentially faster than that of any nonhyperuniform two-phase medium, the "slowest" being antihyperuniform media. The stealthy hyperuniform class is characterized by an excess spreadability with the fastest decay rate among all translationally invariant microstructures. Moreover, we establish a remarkable connection between the spreadability and an outstanding problem in discrete geometry, namely, microstructures with "fast" spreadabilities are also those that can be derived from efficient "coverings" of space. We also identify heretofore unnoticed remarkable links between the spreadability ${\cal S}(t)$ and NMR pulsed field gradient spin-echo amplitude as well as diffusion MRI measurements. △ Less

Submitted 19 September, 2021; originally announced September 2021.

Comments: 18 pages, 10 figures

Sub-diffusive wave transport and weak localization transition in three-dimensional stealthy hyperuniform disordered systems

Authors: Fabrizio Sgrignuoli, Salvatore Torquato, Luca Dal Negro

Abstract: The purpose of this work is to understand the fundamental connection between structural correlations and light localization in three-dimensional (3D) open scattering systems of finite size. We numerically investigate the transport of vector electromagnetic waves scattered by resonant electric dipoles spatially arranged in 3D space by stealthy hyperuniform disordered point patterns. Three-dimension… ▽ More The purpose of this work is to understand the fundamental connection between structural correlations and light localization in three-dimensional (3D) open scattering systems of finite size. We numerically investigate the transport of vector electromagnetic waves scattered by resonant electric dipoles spatially arranged in 3D space by stealthy hyperuniform disordered point patterns. Three-dimensional stealthy hyperuniform disordered systems (3D-SHDS) are engineered with different structural correlation properties determined by their degree of stealthiness $χ$. Such fine control of exotic states of amorphous matter enables the systematic design of optical media that interpolate in a tunable fashion between uncorrelated random structures and crystalline materials. By solving the electromagnetic multiple scattering problem using Green's matrix spectral method, we establish a transport phase diagram that demonstrates a distinctive transition from a diffusive to a weak localization regime beyond a critical scattering density that depends on $χ$. The transition is characterized by studying the Thouless number and the spectral statistics of the scattering resonances. In particular, by tuning the $χ$ parameter, we demonstrate large spectral gaps and suppressed sub-radiant proximity resonances, facilitating light localization. Moreover, consistently with previous studies, our results show a region of the transport phase diagram where the investigated scattering systems become transparent. Our work provides a systematic description of the transport and weak localization properties of light in stealthy hyperuniform structures and motivates the engineering of novel photonic systems with enhanced light-matter interactions for applications to both classical and quantum devices. △ Less

Submitted 21 February, 2022; v1 submitted 8 September, 2021; originally announced September 2021.

Comments: 12 pages, 7 figures

Understanding Degeneracy of Two-Point Correlation Functions via Debye Random Media

Authors: Murray Skolnick, Salvatore Torquato

Abstract: It is well-known that the degeneracy of two-phase microstructures with the same volume fraction and two-point correlation function $S_2(\mathbf{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 clas… ▽ More It is well-known that the degeneracy of two-phase microstructures with the same volume fraction and two-point correlation function $S_2(\mathbf{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)]. 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. △ Less

Submitted 14 October, 2021; v1 submitted 27 July, 2021; originally announced July 2021.

Journal ref: Phys. Rev. E 104, 045306 (2021)

Manifestations of metastable criticality in the long-range structure of model water glasses

Authors: Thomas E. Gartner III, Salvatore Torquato, Roberto Car, Pablo G. Debenedetti

Abstract: Much attention has been devoted to water's metastable phase behavior, including polyamorphism (multiple amorphous solid phases), and the hypothesized liquid-liquid transition and associated critical point. However, the possible relationship between these phenomena remains incompletely understood. Using molecular dynamics simulations of the realistic TIP4P/2005 model, we found a striking signature… ▽ More Much attention has been devoted to water's metastable phase behavior, including polyamorphism (multiple amorphous solid phases), and the hypothesized liquid-liquid transition and associated critical point. However, the possible relationship between these phenomena remains incompletely understood. Using molecular dynamics simulations of the realistic TIP4P/2005 model, we found a striking signature of the liquid-liquid critical point in the structure of water glasses, manifested as a pronounced increase in long-range density fluctuations at pressures proximate to the critical pressure. By contrast, these signatures were absent in glasses of two model systems that lack a critical point. We also characterized the departure from equilibrium upon vitrification via the non-equilibrium index; water-like systems exhibited a strong pressure dependence in this metric, whereas simple liquids did not. These results reflect a surprising relationship between the metastable equilibrium phenomenon of liquid-liquid criticality and the non-equilibrium structure of glassy water, with implications for our understanding of water phase behavior and glass physics. Our calculations suggest a possible experimental route to probing the existence of the liquid-liquid transition in water and other fluids. △ Less

Submitted 7 June, 2021; originally announced June 2021.

Comments: 9 pages, 5 figures. Published in Nature Communications (2021)

Journal ref: Nat Commun 12, 3398 (2021)

Critical pore radius and transport properties of disordered hard- and overlapping-sphere models

Authors: Michael A. Klatt, Robert M. Ziff, Salvatore Torquato

Abstract: Descriptors that characterize the geometry and topology of the pore space of porous media are intimately linked to their transport properties. We quantify such descriptors, including pore-size functions and the critical pore radius $δ_c$, for four different models: maximally random jammed sphere packings, overlapping spheres, equilibrium hard spheres, and inherent structures of the quantizer energ… ▽ More Descriptors that characterize the geometry and topology of the pore space of porous media are intimately linked to their transport properties. We quantify such descriptors, including pore-size functions and the critical pore radius $δ_c$, for four different models: maximally random jammed sphere packings, overlapping spheres, equilibrium hard spheres, and inherent structures of the quantizer energy. For precise estimates of the percolation thresholds, we use a strict relation of the void percolation around sphere configurations to weighted bond percolation on the corresponding Voronoi networks. We use the Newman-Ziff algorithm to determine the percolation threshold using universal properties of the cluster size distribution. Often, $δ_c$ is used as the key characteristic length scale that determines the fluid permeability $k$. A recent study [Torquato. Adv. Wat. Resour. 140, 103565 (2020)] suggested for porous media with a well-connected pore space an alternative estimate of $k$ based on the second moment of the pore size $\langleδ^2\rangle$. Here, we confirm that, for all porosities and all models considered, $δ_c^2$ is to a good approximation proportional to $\langleδ^2\rangle$. However, unlike $\langleδ^2\rangle$, the permeability estimate based on $δ_c^2$ does not predict the correct ranking of $k$ for our models. Thus, we confirm $\langleδ^2\rangle$ to be a promising candidate for convenient and reliable estimates of $k$ for porous media with a well-connected pore space. Moreover, we compare the fluid permeability of our models with varying degrees of order, as measured by the $τ$ order metric. We find that (effectively) hyperuniform models tend to have lower values of $k$ than their nonhyperuniform counterparts. Our findings could facilitate the design of porous media with desirable transport properties via targeted pore statistics. △ Less

Submitted 26 July, 2021; v1 submitted 19 April, 2021; originally announced April 2021.

MSC Class: 82B43; 82D30; 74A40;

Journal ref: Phys. Rev. E 104, 014127 (2021)

Structural Characterization of Many-Particle Systems on Approach to Hyperuniform States

Authors: Salvatore Torquato

Abstract: We explore quantitative descriptors that herald when a many-particle system in $d$-dimensional Euclidean space $\mathbb{R}^d$ approaches a hyperuniform state as a function of the relevant control parameter. We establish quantitative criteria to ascertain the extent of hyperuniform and nonhyperuniform distance-scaling regimes n terms of the ratio $B/A$, where $A$ is "volume" coefficient and $B$ is… ▽ More We explore quantitative descriptors that herald when a many-particle system in $d$-dimensional Euclidean space $\mathbb{R}^d$ approaches a hyperuniform state as a function of the relevant control parameter. We establish quantitative criteria to ascertain the extent of hyperuniform and nonhyperuniform distance-scaling regimes n terms of the ratio $B/A$, where $A$ is "volume" coefficient and $B$ is "surface-area" coefficient associated with the local number variance $σ^2(R)$ for a spherical window of radius $R$. To complement the known direct-space representation of the coefficient $B$ in terms of the total correlation function $h({\bf r})$, we derive its corresponding Fourier representation in terms of the structure factor $S({\bf k})$, which is especially useful when scattering information is available experimentally or theoretically. We show that the free-volume theory of the pressure of equilibrium packings of identical hard spheres that approach a strictly jammed state either along the stable crystal or metastable disordered branch dictates that such end states be exactly hyperuniform. Using the ratio $B/A$, the hyperuniformity index $H$ and the direct-correlation function length scale $ξ_c$, we study three different exactly solvable models as a function of the relevant control parameter, either density or temperature, with end states that are perfectly hyperuniform. We analyze equilibrium hard rods and "sticky" hard-sphere systems in arbitrary space dimension $d$ as a function of density. We also examine low-temperature excited states of many-particle systems interacting with "stealthy" long-ranged pair interactions as the temperature tends to zero. The capacity to identify hyperuniform scaling regimes should be particularly useful in analyzing experimentally- or computationally-generated samples that are necessarily of finite size. △ Less

Submitted 30 March, 2021; v1 submitted 27 March, 2021; originally announced March 2021.

Comments: 19 pages, 14 figures

Journal ref: Phys. Rev. E 103, 052126 (2021)

Kinetic Frustration Effects on Dense Two-Dimensional Packings of Convex Particles and Their Structural Characteristics

Authors: Charles Emmett Maher, Frank H. Stillinger, Salvatore Torquato

Abstract: The study of hard-particle packings is of fundamental importance in physics, chemistry, cell biology, and discrete geometry. Much of the previous work on hard-particle packings concerns their densest possible arrangements. By contrast, we examine kinetic effects inevitably present in both numerical and experimental packing protocols. Specifically, we determine how changing the compression/shear ra… ▽ More The study of hard-particle packings is of fundamental importance in physics, chemistry, cell biology, and discrete geometry. Much of the previous work on hard-particle packings concerns their densest possible arrangements. By contrast, we examine kinetic effects inevitably present in both numerical and experimental packing protocols. Specifically, we determine how changing the compression/shear rate of a two-dimensional packing of noncircular particles causes it to deviate from its densest possible configuration, which is always periodic. The adaptive shrinking cell (ASC) optimization scheme maximizes the packing fraction of a hard-particle packing by first applying random translations and rotations to the particles and then isotropically compressing and shearing the simulation box repeatedly until a possibly jammed state is reached. We use a stochastic implementation of the ASC optimization scheme to mimic different effective time scales by varying the number of particle moves between compressions/shears. We generate dense, effectively jammed, monodisperse, two-dimensional packings of obtuse scalene triangle, rhombus, curved triangle, lens, and "ice cream cone" (a semicircle grafted onto an isosceles triangle) shaped particles, with a wide range of packing fractions and degrees of order. To quantify these kinetic effects, we introduce the kinetic frustration index $K$, which measures the deviation of a packing from its maximum possible packing fraction. To investigate how kinetics affect short- and long-range ordering in these packings, we compute their spectral densities and characterize their contact networks. We find that kinetic effects are most significant when the particles have greater asphericity, less curvature, and less rotational symmetry. This work may be relevant to the design of laboratory packing protocols. △ Less

Submitted 10 March, 2021; originally announced March 2021.

Comments: 25 pages, 18 figures

Effective Elastic Wave Characteristics of Composite Media

Authors: J. Kim, S. Torquato

Abstract: We derive exact expressions for effective elastodynamic properties of two-phase composites in the long-wavelength (quasistatic) regime via homogenized constitutive relations that are local in space. This is accomplished by extending the "strong-contrast" expansion formalism that was previously applied to the static problem. These strong-contrast expansions explicitly incorporate complete microstru… ▽ More We derive exact expressions for effective elastodynamic properties of two-phase composites in the long-wavelength (quasistatic) regime via homogenized constitutive relations that are local in space. This is accomplished by extending the "strong-contrast" expansion formalism that was previously applied to the static problem. These strong-contrast expansions explicitly incorporate complete microstructural information of the composite via an infinite set of $n$-point correlation functions. Utilizing the rapid-convergence properties of these series expansions (even for extreme contrast ratios), we extract accurate approximations that depend on the microstructure via the spectral density, which is easy to compute or measure for any composite. We also investigate the predictive power of modifications of such approximation formulas postulated elsewhere [J. Kim and S. Torquato, Proc. Nat. Acad. Sci. ${\bf 117}$, 8764 (2020)] to extend their applicability ${\it beyond ~the ~quasistatic ~regime.}$ The accuracy of these nonlocal microstructure-dependent approximations is validated by comparison to full-waveform simulation results for certain models of dispersions. We apply our formulas to a variety of models of nonhyperuniform and hyperuniform disordered composites. We demonstrate that hyperuniform systems are less lossy than their nonhyperuniform counterparts in the quasistatic regime, and stealthy hyperuniform media can be perfectly transparent for a wide range of wavenumber. Finally, we discuss how to utilize our approximations for engineering composites with prescribed elastic wave characteristics. △ Less

Submitted 30 December, 2020; originally announced December 2020.

Comments: 30 pages, 7 figures, 1 Supplementary Material. arXiv admin note: text overlap with arXiv:2007.00701

Quantum Phase Transitions in Long-Range Interacting Hyperuniform Spin Chains in a Transverse Field

Authors: Amartya Bose, Salvatore Torquato

Abstract: Hyperuniform states of matter are characterized by anomalous suppression of long-wavelength density fluctuations. While most of interesting cases of disordered hyperuniformity are provided by complex many-body systems like liquids or amorphous solids, classical spin chains with certain long-range interactions have been shown to demonstrate the same phenomenon. It is well-known that the transverse… ▽ More Hyperuniform states of matter are characterized by anomalous suppression of long-wavelength density fluctuations. While most of interesting cases of disordered hyperuniformity are provided by complex many-body systems like liquids or amorphous solids, classical spin chains with certain long-range interactions have been shown to demonstrate the same phenomenon. It is well-known that the transverse field Ising model shows a quantum phase transition (QPT) at zero temperature. Under the quantum effects of a transverse magnetic field, classical hyperuniform spin chains are expected to lose their hyperuniformity. High-precision simulations of these cases are complicated because of the presence of highly nontrivial long-range interactions. We perform extensive analysis of these systems using density matrix renormalization group to study the possibilities of phase transitions and the mechanism by which they lose hyperuniformity. We discover first-order QPTs in the hyperuniform spin chains. An interesting feature of the phase transitions in these disordered hyperuniform spin chains is that, depending on the parameter values, the presence of transverse magnetic field may remarkably lead to increase in the order of the ground state as measured by the "$τ$ order metric," even if hyperuniformity is lost. Therefore, it would be possible to design materials to target specific novel quantum behaviors in the presence of a transverse magnetic field. Our numerical investigations suggest that these spin chains can show no more than two QPTs. We further analyze the long-range interacting spin chains via the Jordan-Wigner mapping, showing that under the pairwise interacting approximation and a mean-field treatment, there can be at most two QPTs. Based on these numerical and theoretical explorations, we conjecture that these spin chains can show a maximum of two QPTs at zero temperature. △ Less

Submitted 14 December, 2020; v1 submitted 11 December, 2020; originally announced December 2020.

Comments: SM was missing in the previous version. Added here

Journal ref: Phys. Rev. B 103, 014118 (2021)

Local Number Fluctuations in Hyperuniform and Nonhyperuniform Systems: Higher-Order Moments and Distribution Functions

Authors: Salvatore Torquato, Jaeuk Kim, Michael A. Klatt

Abstract: The local number variance associated with a spherical sampling window of radius $R$ enables a classification of many-particle systems in $d$-dimensional Euclidean space according to the degree to which large-scale density fluctuations are suppressed, resulting in a demarcation between hyperuniform and nonhyperuniform phyla. To better characterize density fluctuations, we carry out an extensive stu… ▽ More The local number variance associated with a spherical sampling window of radius $R$ enables a classification of many-particle systems in $d$-dimensional Euclidean space according to the degree to which large-scale density fluctuations are suppressed, resulting in a demarcation between hyperuniform and nonhyperuniform phyla. To better characterize density fluctuations, we carry out an extensive study of higher-order moments, including the skewness $γ_1(R)$, excess kurtosis $γ_2(R)$ and the corresponding probability distribution function $P[N(R)]$ of a large family of models across the first three space dimensions, including both hyperuniform and nonhyperuniform models. Specifically, we derive explicit integral expressions for $γ_1(R)$ and $γ_2(R)$ involving up to three- and four-body correlation functions, respectively. We also derive rigorous bounds on $γ_1(R)$, $γ_2(R)$ and $P[N(R)]$. High-quality simulation data for these quantities are generated for each model. We also ascertain the proximity of $P[N(R)]$ to the normal distribution via a novel Gaussian distance metric $l_2(R)$. Among all models, the convergence to a central limit theorem (CLT) is generally fastest for the disordered hyperuniform processes. The convergence to a CLT is slower for standard nonhyperuniform models, and slowest for the antihyperuniform model studied here. We prove that one-dimensional hyperuniform systems of class I or any $d$-dimensional lattice cannot obey a CLT. Remarkably, we discovered that the gamma distribution provides a good approximation to $P[N(R)]$ for all models that obey a CLT, enabling us to estimate the large-$R$ scalings of $γ_1(R)$, $γ_2(R)$ and $l_2(R)$. For any $d$-dimensional model that "decorrelates" or "correlates" with $d$, we elucidate why $P[N(R)]$ increasingly moves toward or away from Gaussian-like behavior, respectively. △ Less

Submitted 7 December, 2020; v1 submitted 3 December, 2020; originally announced December 2020.

Comments: 23 pages, 8 figures

Journal ref: Phys. Rev. X 11, 021028 (2021)

Generation and Structural Characterization of Debye Random Media

Authors: Zheng Ma, Salvatore Torquato

Abstract: In their seminal paper on scattering by an inhomogeneous solid, Debye and coworkers proposed a simple exponentially decaying function for the two-point correlation function of an idealized class of two-phase random media. Such {\it Debye random media}, which have been shown to be realizable, are singularly distinct from all other models of two-phase media in that they are entirely defined by their… ▽ More In their seminal paper on scattering by an inhomogeneous solid, Debye and coworkers proposed a simple exponentially decaying function for the two-point correlation function of an idealized class of two-phase random media. Such {\it Debye random media}, which have been shown to be realizable, are singularly distinct from all other models of two-phase media in that they are entirely defined by their one- and two-point correlation functions. To our knowledge, there has been no determination of other microstructural descriptors of Debye random media. In this paper, we generate Debye random media in two dimensions using an accelerated Yeong-Torquato construction algorithm. We then ascertain microstructural descriptors of the constructed media, including their surface correlation functions, pore-size distributions, lineal-path function, and chord-length probability density function. Accurate semi-analytic and empirical formulas for these descriptors are devised. We compare our results for Debye random media to those of other popular models (overlapping disks and equilibrium hard disks), and find that the former model possesses a wider spectrum of hole sizes, including a substantial fraction of large holes. Our algorithm can be applied to generate other models defined by their two-point correlation functions, and their other microstructural descriptors can be determined and analyzed by the procedures laid out here. △ Less

Submitted 20 October, 2020; originally announced October 2020.

Journal ref: Physical Review E 102 (4), 043310 (2020)

Predicting permeability via statistical learning on higher-order microstructural information

Authors: Magnus Röding, Zheng Ma, Salvatore Torquato

Abstract: Quantitative structure-property relationships are crucial for the understanding and prediction of the physical properties of complex materials. For fluid flow in porous materials, characterizing the geometry of the pore microstructure facilitates prediction of permeability, a key property that has been extensively studied in material science, geophysics and chemical engineering. In this work, we s… ▽ More Quantitative structure-property relationships are crucial for the understanding and prediction of the physical properties of complex materials. For fluid flow in porous materials, characterizing the geometry of the pore microstructure facilitates prediction of permeability, a key property that has been extensively studied in material science, geophysics and chemical engineering. In this work, we study the predictability of different structural descriptors via both linear regressions and neural networks. A large data set of 30,000 virtual, porous microstructures of different types is created for this end. We compute permeabilities of these structures using the lattice Boltzmann method, and characterize the pore space geometry using one-point correlation functions (porosity, specific surface), two-point surface-surface, surface-void, and void-void correlation functions, as well as the geodesic tortuosity as an implicit descriptor. Then, we study the prediction of the permeability using different combinations of these descriptors. We obtain significant improvements of performance when compared to a Kozeny-Carman regression with only lowest-order descriptors (porosity and specific surface). We find that combining all three two-point correlation functions and tortuosity provides the best prediction of permeability, with the void-void correlation function being the most informative individual descriptor. Moreover, the combination of porosity, specific surface, and geodesic tortuosity provides very good predictive performance. This shows that higher-order correlation functions are extremely useful for forming a general model for predicting physical properties of complex materials. Additionally, our results suggest that neural networks are superior to the more conventional regression methods for establishing quantitative structure-property relationships. △ Less

Submitted 7 September, 2020; originally announced September 2020.

Journal ref: Sci Rep 10, 15239 (2020)

Nearest-Neighbor Functions for Disordered Stealthy Hyperuniform Many-Particle Systems

Authors: Timothy M. Middlemas, Salvatore Torquato

Abstract: Disordered stealthy many-particle systems in $\mathbb{R}^d$ are exotic states of matter that suppress single scattering events for a finite range of wavenumbers around the origin in reciprocal space. We derive analytical formulas for the nearest-neighbor functions of disordered stealthy systems. First, we analyze asymptotic small-$r$ approximations and bounding expressions of the nearest-neighbor… ▽ More Disordered stealthy many-particle systems in $\mathbb{R}^d$ are exotic states of matter that suppress single scattering events for a finite range of wavenumbers around the origin in reciprocal space. We derive analytical formulas for the nearest-neighbor functions of disordered stealthy systems. First, we analyze asymptotic small-$r$ approximations and bounding expressions of the nearest-neighbor functions based on the pseudo-hard-sphere ansatz. We then determine how many of the standard $n$-point correlation functions are needed to determine the nearest neighbor functions, and find that a finite number suffice. Via theoretical and computational methods, we compare the large-$r$ behavior of these functions for disordered stealthy systems to those belonging to crystalline lattices. Such ordered and disordered stealthy systems have bounded hole sizes. However, we find that the approach to the critical-hole size can be quantitatively different. We argue that the probability of finding a hole close to the critical-hole size should decrease as a power law with an exponent only dependent on the space dimension $d$ for ordered systems, but that this probability decays asymptotically faster for disordered systems. This implies that holes close to the critical-hole size are rarer in disordered systems. The rarity of observing large holes in disordered systems creates substantial numerical difficulties in sampling the nearest neighbor distributions near the critical-hole size. This motivates both the need for new computational methods for efficient sampling and the development of novel theoretical methods. We also devise a simple analytical formula that accurately describes these systems in the underconstrained regime for all $r$. These results provide a foundation for the analytical description of the nearest-neighbor functions of stealthy systems in the disordered, underconstrained regime. △ Less

Submitted 7 September, 2020; v1 submitted 31 August, 2020; originally announced September 2020.

Comments: 41 pages, 17 figures

Characterizing the Hyperuniformity of Ordered and Disordered Two-Phase Media

Authors: Jaeuk Kim, Salvatore Torquato

Abstract: The hyperuniformity concept provides a unified means to classify all perfect crystals, perfect quasicrystals, and exotic amorphous states of matter according to their capacity to suppress large-scale density fluctuations. While the classification of hyperuniform point configurations has received considerable attention, much less is known about the classification of hyperuniform heterogeneous two-p… ▽ More The hyperuniformity concept provides a unified means to classify all perfect crystals, perfect quasicrystals, and exotic amorphous states of matter according to their capacity to suppress large-scale density fluctuations. While the classification of hyperuniform point configurations has received considerable attention, much less is known about the classification of hyperuniform heterogeneous two-phase media, which include composites, porous media, foams, cellular solids, colloidal suspensions and polymer blends. The purpose of this article is to begin such a program for certain two-dimensional models of hyperuniform two-phase media by ascertaining their local volume-fraction variances $σ^2_{_V}(R)$ and the associated hyperuniformity order metrics $\overline{B}_V$. This is a highly challenging task because the geometries and topologies of the phases are generally much richer and more complex than point-configuration arrangements and one must ascertain a broadly applicable length scale to make key quantities dimensionless. Therefore, we focus on a certain class of two-dimensional periodic cellular networks, periodic and disordered/irregular packings, some of which maximize their effective transport and elastic properties. Among the cellular networks considered, the honeycomb networks have a minimal value of $\overline{B}_V$ across all volume fractions. On the other hand, among all packings considered, the triangular-lattice packings have the smallest values of $\overline{B}_V$ for the possible range of volume fractions. Among all structures studied here, the triangular-lattice packing has the minimal order metric for almost all volume fractions. Our study provides a theoretical foundation for the establishment of hyperuniformity order metrics for general two-phase media and a basis to discover new hyperuniform two-phase systems by inverse design procedures. △ Less

Submitted 16 August, 2020; originally announced August 2020.

Comments: 13 pages, 8 Figures, and 4 tables

Journal ref: Phys. Rev. E 103, 012123 (2021)

Nonlocal Effective Electromagnetic Wave Characteristics of Composite Media: Beyond the Quasistatic Regime

Authors: Salvatore Torquato, Jaeuk Kim

Abstract: We derive exact nonlocal expressions for the effective dielectric constant tensor ${\boldsymbol \varepsilon}_e({\bf k}_I, ω)$ of disordered two-phase composites and metamaterials from first principles. This formalism extends the long-wavelength limitations of conventional homogenization estimates of ${\boldsymbol \varepsilon}_e({\bf k}_I, ω)$ for arbitrary microstructures so that it can capture sp… ▽ More We derive exact nonlocal expressions for the effective dielectric constant tensor ${\boldsymbol \varepsilon}_e({\bf k}_I, ω)$ of disordered two-phase composites and metamaterials from first principles. This formalism extends the long-wavelength limitations of conventional homogenization estimates of ${\boldsymbol \varepsilon}_e({\bf k}_I, ω)$ for arbitrary microstructures so that it can capture spatial dispersion well beyond the quasistatic regime (where $ω$ and ${\bf k}_I$ are frequency and wavevector of the incident radiation). This is done by deriving nonlocal strong-contrast expansions that exactly account for multiple scattering for the range of wavenumbers for which our extended homogenization theory applies, i.e., $0 \le |{\bf k}_I| \ell \lesssim 1$ (where $\ell$ is a characteristic heterogeneity length scale). Due to the fast-convergence properties of such expansions, their lower-order truncations yield accurate closed-form approximate formulas for ${\varepsilon}_e({\bf k}_I,ω)$ that incorporate microstructural information via the spectral density, which is easy to compute for any composite. The accuracy of these microstructure-dependent approximations is validated by comparison to full-waveform simulation methods for both 2D and 3D ordered and disordered models of composite media. Thus, our closed-form formulas enable one to predict accurately and efficiently the effective wave characteristics well beyond the quasistatic regime without having to perform full-blown simulations. Among other results, we show that certain disordered hyperuniform particulate composites exhibit novel wave characteristics. Our results demonstrate that one can design the effective wave characteristics of a disordered composite by engineering the microstructure to possess tailored spatial correlations at prescribed length scales. △ Less

Submitted 16 March, 2021; v1 submitted 1 July, 2020; originally announced July 2020.

Comments: 26 pages, 13 figures (revised after the review process)

Journal ref: Phys. Rev. X 11, 021002 (2021)

Optimized Large Hyperuniform Binary Colloidal Suspensions in Two Dimensions

Authors: Zheng Ma, Enrique Lomba, Salvatore Torquato

Abstract: The creation of disordered hyperuniform materials with potentially extraordinary optical properties requires a capacity to synthesize large samples that are effectively hyperuniform down to the nanoscale. Motivated by this challenge, we propose a fabrication protocol using binary superparamagnetic colloidal particles confined in a 2D plane. The strong and long-ranged dipolar interaction induced by… ▽ More The creation of disordered hyperuniform materials with potentially extraordinary optical properties requires a capacity to synthesize large samples that are effectively hyperuniform down to the nanoscale. Motivated by this challenge, we propose a fabrication protocol using binary superparamagnetic colloidal particles confined in a 2D plane. The strong and long-ranged dipolar interaction induced by a tunable magnetic field is free from screening effects that attenuates long-ranged electrostatic interactions in charged colloidal systems. Specifically, we find a family of optimal size ratios that makes the two-phase system effectively hyperuniform. We show that hyperuniformity is a general consequence of low isothermal compressibilities, which makes our protocol suitable to systems with other long-ranged soft interactions, dimensionalities and/or polydispersity. Our methodology paves the way to synthesize large photonic hyperuniform materials that function in the visible to infrared range and hence may accelerate the discovery of novel photonic materials. △ Less

Submitted 19 May, 2020; originally announced May 2020.

Journal ref: Phys. Rev. Lett. 125, 068002 (2020)

A minimal statistical-mechanical model for multihyperuniform patterns in avian retina

Authors: Enrique Lomba, Jean-Jacques Weis, Leandro Guisández, Salvatore Torquato

Abstract: Birds are known for their extremely acute sense of vision. The very peculiar structural distribution of five different types of cones in the retina underlies this exquisite ability to sample light. It was recently found that each cone population as well as their total population display a disordered pattern in which long wave-length density fluctuations vanish. This property, known as hyperuniform… ▽ More Birds are known for their extremely acute sense of vision. The very peculiar structural distribution of five different types of cones in the retina underlies this exquisite ability to sample light. It was recently found that each cone population as well as their total population display a disordered pattern in which long wave-length density fluctuations vanish. This property, known as hyperuniformity is also present in perfect crystals. In situations like the avian retina in which both the global structure and that of each component display hyperuniformity, the system is said to be multi-hyperuniform. In this work, we aim at devising a minimal statistical-mechanical model that can reproduce the main features of the spatial distribution of photoreceptors in avian retina, namely the presence of disorder, multi-hyperuniformity and local hetero-coordination. This last feature is key to avoid local clustering of the same type of photoreceptors, an undesirable feature for the efficient sampling of light. For this purpose we formulate a simple model that definitively exhibits the required structural properties, namely an equimolar three-component mixture (one component to sample each primary color, red, green, and blue) of non-additive hard disks to which a long-range logarithmic repulsion is added between like particles. A Voronoi analysis of our idealized system of photoreceptors shows that the space-filling Voronoi polygons interestingly display a rather uniform area distribution, symmetrically centered around that of a regular lattice, a structural property also found in human retina. Disordered multi-hyperuniformity offers an alternative to generate photoreceptor patterns with minimal long-range concentration and density fluctuations. This is the key to overcome the difficulties in devising an efficient visual system in which crystal-like order is absent. △ Less

Submitted 4 May, 2020; originally announced May 2020.

Journal ref: Phys. Rev. E 102, 012134 (2020)

Realizable hyperuniform and nonhyperuniform particle configurations with targeted spectral functions via effective pair interactions

Authors: Ge Zhang, Salvatore Torquato

Abstract: The capacity to identify realizable many-body configurations associated with targeted functional forms for the pair correlation function $g_2(r)$ or its corresponding structure factor $S(k)$ is of great fundamental and practical importance. A major aim of this paper is to expand our theoretical knowledge of the class of pair correlation functions or structure factors that are realizable by classic… ▽ More The capacity to identify realizable many-body configurations associated with targeted functional forms for the pair correlation function $g_2(r)$ or its corresponding structure factor $S(k)$ is of great fundamental and practical importance. A major aim of this paper is to expand our theoretical knowledge of the class of pair correlation functions or structure factors that are realizable by classical disordered ensembles of particle configurations. We first introduce a theoretical formalism that provides a means to draw classical particle configurations from canonical ensembles with certain pairwise-additive potentials that could correspond to targeted analytical functional forms for the structure factor. This formulation enables us to devise an improved algorithm to construct systematically canonical-ensemble particle configurations with such targeted pair statistics. We test the algorithm by targeting several different structure factors across dimensions that are known to be realizable and one hyperuniform target that is known to be nontrivially unrealizable. Subsequently, we also target several families of structure-factor functions that meet the known necessary realizability conditions but were heretofore not known to be realizable. Moreover, we also explore unusual nonhyperuniform targets, including "hyposurficial" and "anti-hyperuniform" examples. In all of these instances, the targeted structure factors were achieved with high accuracy, suggesting that they are indeed realizable by equilibrium configurations with pairwise interactions at positive temperatures. Our findings lead us to the conjecture that any realizable structure factor corresponding to either a translationally invariant equilibrium or nonequilibrium system can be attained by an equilibrium ensemble involving only effective pair interactions. △ Less

Submitted 5 April, 2020; originally announced April 2020.

Journal ref: Phys. Rev. E 101, 032124 (2020)

Cloaking the Underlying Long-Range Order of Randomly Perturbed Lattices

Authors: Michael A. Klatt, Jaeuk Kim, Salvatore Torquato

Abstract: Random, uncorrelated displacements of particles on a lattice preserve the hyperuniformity of the original lattice, that is, normalized density fluctuations vanish in the limit of infinite wavelengths. In addition to a diffuse contribution, the scattering intensity from the the resulting point pattern typically inherits the Bragg peaks (long-range order) of the original lattice. Here we demonstrate… ▽ More Random, uncorrelated displacements of particles on a lattice preserve the hyperuniformity of the original lattice, that is, normalized density fluctuations vanish in the limit of infinite wavelengths. In addition to a diffuse contribution, the scattering intensity from the the resulting point pattern typically inherits the Bragg peaks (long-range order) of the original lattice. Here we demonstrate how these Bragg peaks can be hidden in the effective diffraction pattern of independent and identically distributed perturbations. All Bragg peaks vanish if and only if the sum of all probability densities of the positions of the shifted lattice points is a constant at all positions. The underlying long-range order is then 'cloaked' in the sense that it cannot be reconstructed from the pair correlation function alone. On the one hand, density fluctuations increase monotonically with the strength of perturbations $a$, as measured by the hyperuniformity order metric $\overlineΛ$. On the other hand, the disappearance and reemergence of long-range order, depending on whether the system is cloaked or not as the perturbation strength increases, is manifestly captured by the $τ$ order metric. Therefore, while the perturbation strength $a$ may seem to be a natural choice for an order metric of perturbed lattices, the $τ$ order metric is a superior choice. It is noteworthy that cloaked perturbed lattices allow one to easily simulate very large samples (with at least $10^6$ particles) of disordered hyperuniform point patterns without Bragg peaks. △ Less

Submitted 5 April, 2020; v1 submitted 22 January, 2020; originally announced January 2020.

Comments: 10 pages, 9 figures

Journal ref: Phys. Rev. E 101, 032118 (2020)

Hard convex lens-shaped particles: Characterization of dense disordered packings

Authors: Giorgio Cinacchi, Salvatore Torquato

Abstract: Among the family of hard convex lens-shaped particles (lenses), the one with aspect ratio equal to 2/3 is `optimal' in the sense that the maximally random jammed (MRJ) packings of such lenses achieve the highest packing fraction $φ_{\rm MRJ} \simeq 0.73$. This value is only a few percent lower than $φ_{\rm DKP} = 0.76210\dots$, the packing fraction of the corresponding densest-known crystalline (d… ▽ More Among the family of hard convex lens-shaped particles (lenses), the one with aspect ratio equal to 2/3 is `optimal' in the sense that the maximally random jammed (MRJ) packings of such lenses achieve the highest packing fraction $φ_{\rm MRJ} \simeq 0.73$. This value is only a few percent lower than $φ_{\rm DKP} = 0.76210\dots$, the packing fraction of the corresponding densest-known crystalline (degenerate) packings. By exploiting the appreciably reduced propensity that a system of such optimal lenses has to positionally and orientationally order, disordered packings of them are progressively generated by a Monte Carlo method-based procedure from the dilute equilibrium isotropic fluid phase to the dense nonequilibrium MRJ state. This allows one to closely monitor how the (micro)structure of these packings changes in the process of formation of the MRJ state. The gradual changes undergone by the many structural descriptors calculated can coherently and consistently be traced back to the gradual increase in contacts between the hard particles until the isostatic mean value of 10 contact neighbors per lens is reached at the effectively hyperuniform MRJ state. Compared to the MRJ state of hard spheres, the MRJ state of such optimal lenses is denser (less porous), more disordered, and rattler-free. This set of characteristics makes them good glass formers. It is possible that this conclusion may also hold for other hard convex uniaxial particles with a correspondingly similar aspect ratio, be they oblate or prolate, and that, by using suitable biaxial variants of them, that set of characteristics might further improve. △ Less

Submitted 3 September, 2019; originally announced September 2019.

Journal ref: Phys. Rev. E 100, 062902 (2019)

Multifunctional Composites for Elastic and Electromagnetic Wave Propagation

Authors: Jaeuk Kim, Salvatore Torquato

Abstract: Composites are ideally suited to achieve desirable multifunctional effective properties since the best properties of different materials can be judiciously combined with designed microstructures. Here we establish cross-property relations for two-phase composite media that link effective elastic and electromagnetic wave characteristics to one another, including the respective effective wave speeds… ▽ More Composites are ideally suited to achieve desirable multifunctional effective properties since the best properties of different materials can be judiciously combined with designed microstructures. Here we establish cross-property relations for two-phase composite media that link effective elastic and electromagnetic wave characteristics to one another, including the respective effective wave speeds and attenuation coefficients, which facilitate multifunctional material design. This is achieved by deriving accurate formulas for the effective electromagnetic and elastodynamic properties that depend on the wavelengths of the incident waves and the microstructure via the spectral density. Our formulas enable us to explore the wave characteristics of a broad class of disordered microstructures because they apply, unlike conventional formulas, for a wide range of incident wavelengths, i.e., well beyond the long-wavelength regime. This capability enables us to study the dynamic properties of exotic disordered ``hyperuniform'' composites that can have advantages over crystalline ones, such as nearly optimal, direction-independent properties and robustness against defects. We specifically show that disordered ``stealthy'' hyperuniform microstructures exhibit novel wave characteristics, e.g., low-pass filters that transmit waves ``isotropically'' up to a finite wavenumber. Our cross-property relations for the effective wave characteristics can be applied to design multifunctional composites via inverse techniques. Design examples include structural components that require high stiffness and electromagnetic absorption, heat-sinks for CPUs, and sound-absorbing housings for motors that have to efficiently emit thermal radiation and suppress mechanical vibrations, and nondestructive evaluation of the elastic moduli of materials from the effective dielectric response. △ Less

Submitted 11 May, 2020; v1 submitted 19 August, 2019; originally announced August 2019.

Comments: 11 pages and 5 figures; one supplementary information; considerable revisions are made

Journal ref: PNAS 117 (16), 8764-8774 (2020)

Methodology to Construct Large Realizations of Perfectly Hyperuniform Disordered Packings

Authors: Jaeuk Kim, Salvatore Torquato

Abstract: Disordered hyperuniform packings are unusual amorphous states of two-phase materials that are endowed with exotic physical properties. Such hyperuniform systems are characterized by an anomalous suppression of volume-fraction fluctuations at infinitely long-wavelengths, compared to ordinary disordered materials. While there has been growing interest in such singular states of amorphous matter, a m… ▽ More Disordered hyperuniform packings are unusual amorphous states of two-phase materials that are endowed with exotic physical properties. Such hyperuniform systems are characterized by an anomalous suppression of volume-fraction fluctuations at infinitely long-wavelengths, compared to ordinary disordered materials. While there has been growing interest in such singular states of amorphous matter, a major obstacle has been an inability to produce large hyperuniform samples due to practical limitations of conventional methods. To overcome these limitations, we introduce a general theoretical methodology to construct perfectly hyperuniform packings in $d$-dimensional Euclidean space. Specifically, beginning with an initial general tessellation of space that meets a "bounded-cell" condition, hard particles are placed inside each cell such that the packing fraction within the cell becomes identical to the global packing fraction. We prove that the constructed packings are perfectly hyperuniform in the infinite-sample-size limit and their hyperuniformity is independent of particle shapes, positions, and numbers per cell. We employ two distinct types of initial tessellations: Voronoi as well as sphere tessellations. Beginning with Voronoi tessellations, we show that our algorithm can remarkably convert extremely large nonhyperuniform packings into hyperuniform ones in $\mathbb{R}^2$ and $\mathbb{R}^3$. Implementing our theoretical methodology on sphere tessellations, we establish the hyperuniformity of the classical Hashin-Shtrikman multiscale coated-spheres structures that possess optimal effective transport and elastic properties. In addition, the tunability of our methodology offers promise for discovery of novel disordered hyperuniform two-phase materials. △ Less

Submitted 1 June, 2019; v1 submitted 28 January, 2019; originally announced January 2019.

Comments: 22 pages, 13 figures, and supplementary material (6 pages + 2 figures + 9 tables)

Journal ref: Phys. Rev. E 99, 052141 (2019)

New Tessellation-Based Procedure to Design Perfectly Hyperuniform Disordered Dispersions for Materials Discovery

Authors: Jaeuk Kim, Salvatore Torquato

Abstract: Disordered hyperuniform dispersions are exotic amorphous two-phase materials characterized by an anomalous suppression of long-wavelength volume-fraction fluctuations, endowing them with novel physical properties. While such unusual materials have received considerable attention, a stumbling block has been an inability to create large samples that are truly hyperuniform due to current computationa… ▽ More Disordered hyperuniform dispersions are exotic amorphous two-phase materials characterized by an anomalous suppression of long-wavelength volume-fraction fluctuations, endowing them with novel physical properties. While such unusual materials have received considerable attention, a stumbling block has been an inability to create large samples that are truly hyperuniform due to current computational and experimental limitations. To overcome such barriers, we introduce a new, simple construction procedure that guarantees perfect hyperuniformity for very large sample sizes. It involves tessellating space into cells and then inserting a particle into each cell such that the local-cell particle packing fractions are identical to the global packing fraction. We analytically prove that such dispersions are perfectly hyperuniform in the infinite-sample-size limit. Our methodology enables a remarkable mapping that converts a very large nonhyperuniform disordered dispersion into a perfectly hyperuniform one, which we numerically demonstrate in two and three dimensions. A similar analysis also establishes the hyperuniformity of the famous Hashin-Shtrikman multiscale dispersions, which possess optimal transport and elastic properties. Our hyperuniform designs can be readily fabricated using modern photolithographic and 3D printing technologies. The exploration of the enormous class of hyperuniform dispersions that can be designed and tuned by our tessellation-based methodology paves the way for accelerating the discovery of novel hyperuniform materials. △ Less

Submitted 9 March, 2019; v1 submitted 17 January, 2019; originally announced January 2019.

Comments: 11 pages, 6 figures, 1 supplementary material, and published in Acta Materialia (with several corrections.)

Journal ref: Acta Materialia 168 (2019) 143-151

Hyperuniformity on spherical surfaces

Authors: Ariel G. Meyra, Guillermo J. Zarragoicoechea, Alberto. L. Maltz, Enrique Lomba, Salvatore Torquato

Abstract: In this work we present a study on the characterization of ordered and disordered hyperuniform point distributions on spherical surfaces. In spite of the extensive literature on disordered hyperuniform systems in Euclidean geometries, to date few works have dealt with the problem of hyperuniformity in curved spaces. As a matter of fact, some systems that display disordered hyperuniformity, like th… ▽ More In this work we present a study on the characterization of ordered and disordered hyperuniform point distributions on spherical surfaces. In spite of the extensive literature on disordered hyperuniform systems in Euclidean geometries, to date few works have dealt with the problem of hyperuniformity in curved spaces. As a matter of fact, some systems that display disordered hyperuniformity, like the space distribution of photoreceptors in avian retina, actually occur on curved surfaces. Here we will focus on the local particle number variance and its dependence on the size of the sampling window (which we take to be a spherical cap) for regular and uniform point distributions, as well as for equilibrium configurations of fluid particles interacting through Lennard-Jones, dipole-dipole and charge-charge potentials. We will show how the scaling of the local number variance enables the characterization of hyperuniform point patterns also on spherical surfaces. △ Less

Submitted 11 December, 2018; originally announced December 2018.

Journal ref: Phys. Rev. E 100, 022107 (2019)

Hyperuniformity Order Metric of Barlow Packings

Authors: Timothy M. Middlemas, Frank H. Stillinger, Salvatore Torquato

Abstract: The concept of hyperuniformity has been a useful tool in the study of large-scale density fluctuations in systems ranging across the natural and mathematical sciences. One can rank a large class of hyperuniform systems by their ability to suppress long-range density fluctuations through the use of a hyperuniformity order metric $\barΛ$. We apply this order metric to the Barlow packings, which are… ▽ More The concept of hyperuniformity has been a useful tool in the study of large-scale density fluctuations in systems ranging across the natural and mathematical sciences. One can rank a large class of hyperuniform systems by their ability to suppress long-range density fluctuations through the use of a hyperuniformity order metric $\barΛ$. We apply this order metric to the Barlow packings, which are the infinitely degenerate densest packings of identical rigid spheres that are distinguished by their stacking geometries and include the commonly known fcc lattice and hcp crystal. The "stealthy stacking" theorem implies that these packings are all stealthy hyperuniform, a strong type of hyperuniformity which involves the suppression of scattering up to a wavevector $K$. We describe the geometry of three classes of Barlow packings, two disordered classes and small-period packings. In addition, we compute a lower bound on $K$ for all Barlow packings. We compute $\barΛ$ for the aforementioned three classes of Barlow packings and find that to a very good approximation, it is linear in the fraction of fcc-like clusters, taking values between those of least-ordered hcp and most-ordered fcc. This implies that the $\barΛ$ of all Barlow packings is primarily controlled by the local cluster geometry. These results indicate the special nature of anisotropic stacking disorder, which provides impetus for future research on the development of anisotropic order metrics and hyperuniformity properties. △ Less

Submitted 5 December, 2018; v1 submitted 23 October, 2018; originally announced October 2018.

Comments: 13 pages, 7 figures

Journal ref: Phys. Rev. E 99, 022111 (2019)

Hyperuniformity of generalized random organization models

Authors: Zheng Ma, Salvatore Torquato

Abstract: Studies of random organization models of monodisperse spherical particles have shown that a hyperuniform state is achievable when the system goes through an absorbing phase transition to a critical state. Here we investigate to what extent hyperuniformity is preserved when the model is generalized to particles with a size distribution and/or nonspherical shapes. We begin by examining binary disks… ▽ More Studies of random organization models of monodisperse spherical particles have shown that a hyperuniform state is achievable when the system goes through an absorbing phase transition to a critical state. Here we investigate to what extent hyperuniformity is preserved when the model is generalized to particles with a size distribution and/or nonspherical shapes. We begin by examining binary disks in two dimensions and demonstrate that their critical states are hyperuniform as two-phase media, but not hyperuniform nor multihyperuniform as point patterns formed by the particle centroids. We further confirm the generality of our findings by studying particles with a continuous size distribution. Finally, to study the effect of rotational degrees of freedom, we extend our model to noncircular particles, namely, hard rectangles with various aspect ratios, including the hard-needle limit. Although these systems exhibit only short-range orientational order, hyperuniformity is still preserved. Our analysis reveals that the redistribution of the "mass" of the particles rather than the particle centroids is central to this dynamical process. The consideration of the "active volume fraction" of generalized random organization models may help to resolve which universality class they belong to and hence may lead to a deeper theoretical understanding of absorbing-state models. Our results suggest that general particle systems subject to random organization can be a robust way to fabricate a wide class of hyperuniform states of matter by tuning the structures via different particle-size and -shape distributions. This in turn potentially enables the creation of multifunctional hyperuniform materials with desirable optical, transport, and mechanical properties. △ Less

Submitted 13 February, 2019; v1 submitted 14 September, 2018; originally announced September 2018.

Journal ref: Phys. Rev. E 99, 022115 (2019)

Multifunctional Hyperuniform Cellular Networks: Optimality, Anisotropy and Disorder

Authors: Salvatore Torquato, Duyu Chen

Abstract: Disordered hyperuniform heterogeneous materials are new, exotic amorphous states of matter that behave like crystals in the manner in which they suppress volume-fraction fluctuations at large length scales, and yet are statistically isotropic with no Bragg peaks. It has recently been shown that disordered hyperuniform dielectric two-dimensional cellular network solids possess complete photonic ban… ▽ More Disordered hyperuniform heterogeneous materials are new, exotic amorphous states of matter that behave like crystals in the manner in which they suppress volume-fraction fluctuations at large length scales, and yet are statistically isotropic with no Bragg peaks. It has recently been shown that disordered hyperuniform dielectric two-dimensional cellular network solids possess complete photonic band gaps comparable in size to photonic crystals, while at the same time maintaining statistical isotropy, enabling waveguide geometries not possible with photonic crystals. Motivated by these developments, we explore other functionalities of various two-dimensional ordered and disordered hyperuniform cellular networks, including their effective thermal or electrical conductivities and elastic moduli. We establish the multifunctionality of a class of such low-density networks by demonstrating that they maximize or virtually maximize the effective conductivities and elastic moduli. This is accomplished using the machinery of homogenization theory, including optimal bounds and cross-property bounds, and statistical mechanics. We have identified ordered and disordered hyperuniform low-weight cellular networks that are multifunctional with respect to transport (e.g., heat dissipation and fluid transport), mechanical and electromagnetic properties, which can be readily fabricated using 3D printing and lithographic technologies. △ Less

Submitted 4 August, 2018; originally announced August 2018.

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

Authors: Zheng Ma, Salvatore Torquato

Abstract: 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 (obtainabl… ▽ More 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. △ Less

Submitted 31 July, 2018; originally announced August 2018.

Journal ref: Phys. Rev. E 98, 013307, 2018

Inversion Problems for Fourier Transforms of Particle Distributions

Authors: Jaeuk Kim, Ge Zhang, Frank H. Stillinger, Salvatore Torquato

Abstract: Collective coordinates in a many-particle system are complex Fourier components of the particle density, and often provide useful physical insights. However, given collective coordinates, it is desirable to infer the particle coordinates via inverse transformations. In principle, a sufficiently large set of collective coordinates are equivalent to particle coordinates, but the nonlinear relation b… ▽ More Collective coordinates in a many-particle system are complex Fourier components of the particle density, and often provide useful physical insights. However, given collective coordinates, it is desirable to infer the particle coordinates via inverse transformations. In principle, a sufficiently large set of collective coordinates are equivalent to particle coordinates, but the nonlinear relation between collective and particle coordinates makes the inversion procedure highly nontrivial. Given a "target" configuration in one-dimensional Euclidean space, we investigate the minimal set of its collective coordinates that can be uniquely inverted into particle coordinates. For this purpose, we treat a finite number $M$ of the real and/or the imaginary parts of collective coordinates of the target configuration as constraints, and then reconstruct "solution" configurations whose collective coordinates satisfy these constraints. Both theoretical and numerical investigations reveal that the number of numerically distinct solutions depends sensitively on the chosen collective-coordinate constraints and target configurations. From detailed analysis, we conclude that collective coordinates at the $\lceil\frac{N}{2}\rceil$ smallest wavevectors is the minimal set of constraints for unique inversion, where $\lceil{\cdot}\rceil$ represents the ceiling function. This result provides useful groundwork to the inverse transform of collective coordinates in higher-dimensional systems. △ Less

Submitted 5 July, 2018; originally announced July 2018.

Comments: 21 pages (including references), 10 figures