Source author record

Salvatore Torquato

Salvatore Torquato appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

73works
20topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

73 published item(s)

preprint2023arXiv

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

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.

preprint2023arXiv

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

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)

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

Extraordinary Disordered Hyperuniform Multifunctional Composites

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

preprint2022arXiv

Local Order Metrics for Two-Phase Media Across Length Scales

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.

preprint2022arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

Characterizing the Hyperuniformity of Ordered and Disordered Two-Phase Media

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.

preprint2020arXiv

Cloaking the Underlying Long-Range Order of Randomly Perturbed Lattices

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.

preprint2020arXiv

Multifunctional Composites for Elastic and Electromagnetic Wave Propagation

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.

preprint2020arXiv

Optimized Large Hyperuniform Binary Colloidal Suspensions in Two Dimensions

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.

preprint2020arXiv

Predicting permeability via statistical learning on higher-order microstructural information

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2019arXiv

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

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.

preprint2016arXiv

Characterization of Maximally Random Jammed Sphere Packings: II. Correlation Functions and Density Fluctuations

In the first paper of this series, we introduced Voronoi correlation functions to characterize the structure of maximally random jammed (MRJ) sphere packings across length scales. In the present paper, we determine a variety of correlation functions that can be rigorously related to effective physical properties of MRJ sphere packings and compare them to the corresponding statistical descriptors for overlapping spheres and equilibrium hard-sphere systems. Such structural descriptors arise in rigorous bounds and formulas for effective transport properties, diffusion and reactions constants, elastic moduli, and electromagnetic characteristics. First, we calculate the two-point, surface-void, and surface-surface correlation functions, for which we derive explicit analytical formulas for finite hard-sphere packings. We show analytically how the contacts between spheres in the MRJ packings translate into distinct functional behaviors of these two-point correlation functions that do not arise in the other two models examined here. Then, we show how the spectral density distinguishes the MRJ packings from the other disordered systems in that the spectral density vanishes in the limit of infinite wavelengths. These packings are hyperuniform, which means that density fluctuations on large length scales are anomalously suppressed. Moreover, we study and compute exclusion probabilities and pore size distributions as well as local density fluctuations. We conjecture that for general disordered hard-sphere packings, a central limit theorem holds for the number of points within an spherical observation window. Our analysis links problems of interest in material science, chemistry, physics, and mathematics. In the third paper, we will evaluate bounds and estimates of a host of different physical properties of the MRJ sphere packings based on the structural characteristics analyzed in this paper.

preprint2016arXiv

Critical slowing down and hyperuniformity on approach to jamming

Hyperuniformity characterizes a state of matter that is poised at a critical point at which density or volume-fraction fluctuations are anomalously suppressed at infinite wavelengths. Recently, much attention has been given to the link between strict jamming and hyperuniformity in frictionless hard-particle packings. Doing so requires one to study very large packings, which can be difficult to jam properly. We modify the rigorous linear programming method of Donev et al. [J. Comp. Phys. 197, 139 (2004)] in order to test for jamming in putatively jammed packings of hard-disks in two dimensions. We find that various standard packing protocols struggle to reliably create packings that are jammed for even modest system sizes; importantly, these packings appear to be jammed by conventional tests. We present evidence that suggests that deviations from hyperuniformity in putative maximally random jammed (MRJ) packings can in part be explained by a shortcoming in generating exactly-jammed configurations due to a type of "critical slowing down" as the necessary rearrangements become difficult to realize by numerical protocols. Additionally, various protocols are able to produce packings exhibiting hyperuniformity to different extents, but this is because certain protocols are better able to approach exactly-jammed configurations. Nonetheless, while one should not generally expect exact hyperuniformity for disordered packings with rattlers, we find that when jamming is ensured, our packings are very nearly hyperuniform, and deviations from hyperuniformity correlate with an inability to ensure jamming, suggesting that strict jamming and hyperuniformity are indeed linked. This raises the possibility that the ideal MRJ packings have no rattlers. Our work provides the impetus for the development of packing algorithms that produce large disordered strictly jammed packings that are rattler-free.

preprint2016arXiv

Disordered Hyperuniform Heterogeneous Materials

Disordered hyperuniform many-body systems are distinguishable states of matter that lie between a crystal and liquid: they are like perfect crystals in the way they suppress large-scale density fluctuations and yet are like liquids or glasses in that they are statistically isotropic with no Bragg peaks. These systems play a vital role in a number of fundamental and applied problems: glass formation, jamming, rigidity, photonic and electronic band structure, localization of waves and excitations, self-organization, fluid dynamics, quantum systems, and pure mathematics. systems. Here, we derive new rigorous criteria that disordered hyperuniform two-phase heterogeneous materials must obey and explore their consequences. Two-phase heterogeneous media are ubiquitous, examples include composites and porous media, biological media, foams, polymer blends, granular media, cellular solids, and colloids. We rigorously establish the requirements for sphere packings to be "multihyperuniform." We apply realizability conditions for an autocovariance function and its associated spectral density of a two-phase medium, and then incorporate hyperuniformity as a constraint in order to derive new conditions. We show that some functional forms can immediately be eliminated from consideration and identify other forms that are allowable. Specific examples and counterexamples are described. Contact is made with well-known microstructural models as well as irregular phase-separation and Turing-type patterns. We also ascertain a family of spectral densities that are realizable by disordered hyperuniform two-phase media in any space dimension, and present explicit constructions. These studies provide insight into the nature of disordered hyperuniformity in the context of heterogeneous materials and have implications for the design of such novel amorphous materials.

preprint2016arXiv

Hyperuniformity and its Generalizations

Disordered many-particle hyperuniform systems are exotic amorphous states characterized by anomalous suppression of large-scale density fluctuations. Here we substantially broaden the hyperuniformity concept along four different directions. This includes generalizations to treat fluctuations in the interfacial area in heterogeneous media and surface-area driven evolving microstructures, random scalar fields, divergence-free random vector fields, as well as statistically anisotropic many-particle systems and two-phase media. Interfacial-area fluctuations play a major role in characterizing the microstructure of two-phase systems , physical properties that intimately depend on the geometry of the interface, and evolving two-phase microstructures that depend on interfacial energies (e.g., spinodal decomposition). In the instances of divergence-free random vector fields and statistically anisotropic structures, we show that the standard definition of hyperuniformity must be generalized such that it accounts for the dependence of the relevant spectral functions on the direction in which the origin in Fourier space (nonanalyticities at the origin). Using this analysis, we place some well-known energy spectra from the theory of isotropic turbulence in the context of this generalization of hyperuniformity. We show that there exist many-particle ground-state configurations in which directional hyperuniformity imparts exotic anisotropic physical properties (e.g., elastic, optical and acoustic characteristics) to these states of matter. Such tunablity could have technological relevance for manipulating light and sound waves in ways heretofore not thought possible. We show that disordered many-particle systems that respond to external fields (e.g., magnetic and electric fields) are a natural class of materials to look for directional hyperuniformity.

preprint2016arXiv

Static Structural Signatures of Nearly Jammed Disordered and Ordered Hard-Sphere Packings: Direct Correlation Function

Dynamical signatures are known to precede jamming in hard-particle systems, but static structural signatures have proven more elusive. The observation that compressing hard-particle packings towards jamming causes growing hyperuniformity has paved the way for the analysis of jamming as an "inverted critical point" in which the direct correlation function $c(r)$ diverges. We establish quantitative relationships between various singularities in $c(r)$ and the total correlation function $h(r)$ that provide a concrete means of identifying features that must be expressed in $c(r)$ if one hopes to reproduce details in the pair correlation function accurately. We also analyze systems of three-dimensional monodisperse hard-spheres of diameter $D$ as they approach ordered and disordered jammed configurations. For the latter, we use the Lubachevsky-Stillinger (LS) and Torquato-Jiao (TJ) packing algorithms, which both generate disordered packings, but can show perceptible structural differences. We identify a short-ranged scaling $c(r) \propto -1/r$ as $r \rightarrow 0$ and show that this, along with the developing delta function at $c(D)$, is a consequence of the growing long-rangedness in $c(r)$. Near the freezing density, we identify qualitative differences in the structure factor $S(k)$ as well as $c(r)$ between TJ- and LS-generated configurations and link them to differences in the protocols' packing dynamics. Configurations from both algorithms have structure factors that approach zero in the low-wavenumber limit as jamming is approached and are shown to exhibit a corresponding power-law decay in $c(r)$ for large $r$ as a consequence. Our work advances the notion that static signatures are exhibited by hard-particle packings as they approach jamming and underscores the utility of the direct correlation function as a means of monitoring for an incipient rigid network.

preprint2015arXiv

Characterization of maximally random jammed sphere packings: Voronoi correlation functions

We characterize the structure of maximally random jammed (MRJ) sphere packings by computing the Minkowski functionals (volume, surface area, and integrated mean curvature) of their associated Voronoi cells. The probability distribution functions of these functionals of Voronoi cells in MRJ sphere packings are qualitatively similar to those of an equilibrium hard-sphere liquid and partly even to the uncorrelated Poisson point process, implying that such local statistics are relatively structurally insensitive. This is not surprising because the Minkowski functionals of a single Voronoi cell incorporate only local information and are insensitive to global structural information. To improve upon this, we introduce descriptors that incorporate nonlocal information via the correlation functions of the Minkowski functionals of two cells at a given distance as well as certain cell-cell probability density functions. We evaluate these higher-order functions for our MRJ packings as well as equilibrium hard spheres and the Poisson point process. We find strong anticorrelations in the Voronoi volumes for the hyperuniform MRJ packings, consistent with previous findings for other pair correlations [A. Donev et al., Phys. Rev. Lett. 95, 090604 (2005)], indicating that large-scale volume fluctuations are suppressed by accompanying large Voronoi cells with small cells, and vice versa. In contrast to the aforementioned local Voronoi statistics, the correlation functions of the Voronoi cells qualitatively distinguish the structure of MRJ sphere packings (prototypical glasses) from that of the correlated equilibrium hard-sphere liquids. Moreover, while we did not find any perfect icosahedra (the locally densest possible structure in which a central sphere contacts 12 neighbors) in the MRJ packings, a preliminary Voronoi topology analysis indicates the presence of strongly distorted icosahedra.

preprint2015arXiv

Confined disordered strictly jammed binary sphere packings

Disordered jammed packings under confinement have received considerably less attention than their \textit{bulk} counterparts and yet arise in a variety of practical situations. In this work, we study binary sphere packings that are confined between two parallel hard planes, and generalize the Torquato-Jiao (TJ) sequential linear programming algorithm [Phys. Rev. E {\bf 82}, 061302 (2010)] to obtain putative maximally random jammed (MRJ) packings that are exactly isostatic with high fidelity over a large range of plane separation distances $H$, small to large sphere radius ratio $α$ and small sphere relative concentration $x$. We find that packing characteristics can be substantially different from their bulk analogs, which is due to what we term "confinement frustration". Rattlers in confined packings are generally more prevalent than those in their bulk counterparts. We observe that packing fraction, rattler fraction and degree of disorder of MRJ packings generally increase with $H$, though exceptions exist. Discontinuities in the packing characteristics as $H$ varies in the vicinity of certain values of $H$ are due to associated discontinuous transitions between different jammed states. We also apply the local volume-fraction variance $σ_τ^2(R)$ to characterize confined packings and find that these packings possess essentially the same level of hyperuniformity as their bulk counterparts. Our findings are generally relevant to confined packings that arise in biology (e.g., structural color in birds and insects) and may have implications for the creation of high-density powders and improved battery designs.

preprint2015arXiv

Ground states of stealthy hyperuniform potentials: I. Entropically favored configurations

Systems of particles interacting with "stealthy" pair potentials have been shown to possess infinitely degenerate disordered hyperuniform classical ground states with novel physical properties. Previous attempts to sample the infinitely degenerate ground states used energy minimization techniques, introducing algorithmic dependence that is artificial in nature. Recently, an ensemble theory of stealthy hyperuniform ground states was formulated to predict the structure and thermodynamics that was shown to be in excellent agreement with corresponding computer simulation results in the canonical ensemble (in the zero-temperature limit). In this paper, we provide details and justifications of the simulation procedure, which involves performing molecular dynamics simulations at sufficiently low temperatures and minimizing the energy of the snapshots for both the high-density disordered regime, where the theory applies, as well as lower densities. We also use numerical simulations to extend our study to the lower-density regime. We report results for the pair correlation functions, structure factors, and Voronoi cell statistics. In the high-density regime, we verify the theoretical ansatz that stealthy disordered ground states behave like "pseudo" disordered equilibrium hard-sphere systems in Fourier space. These results show that as the density decreases from the high-density limit, the disordered ground states in the canonical ensemble are characterized by an increasing degree of short-range order and eventually the system undergoes a phase transition to crystalline ground states. We also provide numerical evidence suggesting that different forms of stealthy pair potentials produce the same ground-state ensemble in the zero-temperature limit. Our techniques may be applied to sample this limit of the canonical ensemble of other potentials with highly degenerate ground states.

preprint2015arXiv

Ground states of stealthy hyperuniform potentials. II. Stacked-slider phases

Stealthy potentials, a family of long-range isotropic pair potentials, produce infinitely degenerate disordered ground states at high densities and crystalline ground states at low densities in d-dimensional Euclidean space R^d. In the previous paper in this series, we numerically studied the entropically favored ground states in the canonical ensemble in the zero-temperature limit across the first three Euclidean space dimensions. In this paper, we investigate using both numerical and theoretical techniques metastable stacked-slider phases, which are part of the ground-state manifold of stealthy potentials at densities in which crystal ground states are favored entropically. Our numerical results enable us to devise analytical models of this phase in two, three, and higher dimensions. Utilizing this model, we estimated the size of the feasible region in configuration space of the stacked-slider phase, finding it to be smaller than that of crystal structures in the infinite-system-size limit, which is consistent with our recent previous work. In two dimensions, we also determine exact expressions for the pair correlation function and structure factor of the analytical model of stacked-slider phases and analyze the connectedness of the ground-state manifold of stealthy potentials in this density regime. We demonstrate that stacked-slider phases are distinguishable states of matter; they are nonperiodic, statistically anisotropic structures that possess long-range orientational order but have zero shear modulus. We outline some possible future avenues of research to elucidate our understanding of this unusual phase of matter.

preprint2015arXiv

Inverse design of disordered stealthy hyperuniform spin chains

Positioned between crystalline solids and liquids, disordered many-particle systems which are stealthy and hyperuniform represent new states of matter that are endowed with novel physical and thermodynamic properties. Such stealthy and hyperuniform states are unique in that they are transparent to radiation for a range of wavenumbers around the origin. In this work, we employ recently developed inverse statistical-mechanical methods, which seek to obtain the optimal set of interactions that will spontaneously produce a targeted structure or configuration as a unique ground state, to investigate the spin-spin interaction potentials required to stabilize disordered stealthy hyperuniform one-dimensional (1D) Ising-like spin chains. By performing an exhaustive search over the spin configurations that can be enumerated on periodic 1D integer lattices containing $N=2,3,\ldots,36$ sites, we were able to identify and structurally characterize \textit{all} stealthy hyperuniform spin chains in this range of system sizes. Within this pool of stealthy hyperuniform spin configurations, we then utilized such inverse optimization techniques to demonstrate that stealthy hyperuniform spin chains can be realized as either unique or degenerate disordered ground states of radial long-ranged (relative to the spin chain length) spin-spin interactions. Such exotic ground states are distinctly different from spin glasses in both their inherent structural properties and the nature of the spin-spin interactions required to stabilize them. As such, the implications and significance of the existence of such disordered stealthy hyperuniform ground state spin systems warrants further study, including whether their bulk physical properties and excited states, like their many-particle system counterparts, are singularly remarkable, and can be experimentally realized.

preprint2014arXiv

A Cellular Automaton Model for Tumor Dormancy: Emergence of a Proliferative Switch

Malignant cancers that lead to fatal outcomes for patients may remain dormant for very long periods of time. Although individual mechanisms such as cellular dormancy, angiogenic dormancy and immunosurveillance have been proposed, a comprehensive understanding of cancer dormancy and the "switch" from a dormant to a proliferative state still needs to be strengthened from both a basic and clinical point of view. Computational modeling enables one to explore a variety of scenarios for possible but realistic microscopic dormancy mechanisms and their predicted outcomes. The aim of this paper is to devise such a predictive computational model of dormancy with an emergent "switch" behavior. Specifically, we generalize a previous cellular automaton (CA) model for proliferative growth of solid tumor that now incorporates a variety of cell-level tumor-host interactions and different mechanisms for tumor dormancy, for example the effects of the immune system. Our new CA rules induce a natural "competition" between the tumor and tumor suppression factors in the microenvironment. This competition either results in a "stalemate" for a period of time in which the tumor either eventually wins (spontaneously emerges) or is eradicated; or it leads to a situation in which the tumor is eradicated before such a "stalemate" could ever develop. We also predict that if the number of actively dividing cells within the proliferative rim of the tumor reaches a critical, yet low level, the dormant tumor has a high probability to resume rapid growth. Our findings may shed light on the fundamental understanding of cancer dormancy.

preprint2014arXiv

Avian photoreceptor patterns represent a disordered hyperuniform solution to a multiscale packing problem

Optimal spatial sampling of light rigorously requires that identical photoreceptors be arranged in perfectly regular arrays in two dimensions. Examples of such perfect arrays in nature include the compound eyes of insects and the nearly crystalline photoreceptor patterns of some fish and reptiles. Birds are highly visual animals with five different cone photoreceptor subtypes, yet their photoreceptor patterns are not perfectly regular. By analyzing the chicken cone photoreceptor system consisting of five different cell types using a variety of sensitive microstructural descriptors, we find that the disordered photoreceptor patterns are ``hyperuniform'' (exhibiting vanishing infinite-wavelength density fluctuations), a property that had heretofore been identified in a unique subset of physical systems, but had never been observed in any living organism. Remarkably, the photoreceptor patterns of both the total population and the individual cell types are simultaneously hyperuniform. We term such patterns ``multi-hyperuniform'' because multiple distinct subsets of the overall point pattern are themselves hyperuniform. We have devised a unique multiscale cell packing model in two dimensions that suggests that photoreceptor types interact with both short- and long-ranged repulsive forces and that the resultant competition between the types gives rise to the aforementioned singular spatial features characterizing the system, including multi-hyperuniformity.

preprint2014arXiv

Dense periodic packings of tori

Dense packings of nonoverlapping bodies in three-dimensional Euclidean space are useful models of the structure of a variety of many-particle systems that arise in the physical and biological sciences. Here we investigate the packing behavior of congruent ring tori, which are multiply connected nonconvex bodies of genus 1, as well as horn and spindle tori. We analytically construct a family of dense periodic packings of unlinked tori guided by the organizing principles originally devised for simply connected solid bodies [Torquato and Jiao, PRE 86, 011102 (2012)]. We find that the horn tori as well as certain spindle and ring tori can achieve a packing density higher than the densest known packing of both sphere and ellipsoids. In addition, we study dense packings of cluster of pair-linked ring tori (i.e., Hopf links).

preprint2014arXiv

Diagnosing hyperuniformity in two-dimensional disordered jammed-packings of soft spheres

Hyperuniformity characterizes a state of matter for which density fluctuations diminish towards zero at the largest length scales. However, the task of determining whether or not an experimental system is hyperuniform is experimentally challenging due to finite-resolution, noise and sample-size effects that influence characterization measurements. Here we explore these issues, employing video optical microscopy to study hyperuniformity phenomena in disordered two-dimensional jammed packings of soft spheres. Using a combination of experiment and simulation we characterize the detrimental effects of particle polydispersity, image noise, and finite-size effects on the assignment of hyperuniformity, and we develop a methodology that permits improved diagnosis of hyperuniformity from real-space measurements. The key to this improvement is a simple packing reconstruction algorithm that incorporates particle polydispersity to minimize free volume. In addition, simulations show that hyperuniformity can be ascertained more accurately in direct space than in reciprocal space as a result of finite sample-size. Finally, experimental colloidal packings of soft polymeric spheres are shown to be hyperuniform.

preprint2014arXiv

Equilibrium Phase Behavior and Maximally Random Jammed State of Truncated Tetrahedra

Systems of hard nonspherical particles exhibit a variety of stable phases with different degrees of translational and orientational order, including isotropic liquid, solid crystal, rotator and a variety of liquid crystal phases. In this paper, we employ a Monte Carlo implementation of the adaptive-shrinking-cell (ASC) numerical scheme and free-energy calculations to ascertain with high precision the equilibrium phase behavior of systems of congruent Archimedean truncated tetrahedra over the entire range of possible densities up to the maximal nearly space-filling density. In particular, we find that the system undergoes two first-order phase transitions as the density increases: first a liquid-solid transition and then a solid-solid transition. The isotropic liquid phase coexists with the Conway-Torquato (CT) crystal phase at intermediate densities. At higher densities, we find that the CT phase undergoes another first-order phase transition to one associated with the densest-known crystal. We find no evidence for stable rotator (or plastic) or nematic phases. We also generate the maximally random jammed (MRJ) packings of truncated tetrahedra, which may be regarded to be the glassy end state of a rapid compression of the liquid. We find that such MRJ packings are hyperuniform with an average packing fraction of 0.770, which is considerably larger than the corresponding value for identical spheres (about 0.64). We conclude with some simple observations concerning what types of phase transitions might be expected in general hard-particle systems based on the particle shape and which would be good glass formers.

preprint2014arXiv

Jammed lattice sphere packings

We generate and study an ensemble of isostatic jammed hard-sphere lattices. These lattices are obtained by compression of a periodic system with an adaptive unit cell containing a single sphere until the point of mechanical stability. We present detailed numerical data about the densities, pair correlations, force distributions, and structure factors of such lattices. We show that this model retains many of the crucial structural features of the classical hard-sphere model and propose it as a model for the jamming and glass transitions that enables exploration of much higher dimensions than are usually accessible.

preprint2014arXiv

Marginal stability in jammed packings: quasicontacts and weak contacts

Maximally random jammed (MRJ) sphere packing is a prototypical example of a system naturally poised at the margin between underconstraint and overconstraint. This marginal stability has traditionally been understood in terms of isostaticity, the equality of the number of mechanical contacts and the number of degrees of freedom. Quasicontacts, pairs of spheres on the verge of coming in contact, are irrelevant for static stability, but they come into play when considering dynamic stability, as does the distribution of contact forces. We show that the effects of marginal dynamic stability, as manifested in the distributions of quasicontacts and weak contacts, are consequential and nontrivial. We study these ideas first in the context of MRJ packing of d-dimensional spheres, where we show that the abundance of quasicontacts grows at a faster rate than that of contacts. We reexamine a calculation of Jin et al. (Phys. Rev. E 82, 051126, 2010), where quasicontacts were originally neglected, and we explore the effect of their inclusion in the calculation. This analysis yields an estimate of the asymptotic behavior of the packing density in high dimensions. We argue that this estimate should be reinterpreted as a lower bound. The latter part of the paper is devoted to Bravais lattice packings that possess the minimum number of contacts to maintain mechanical stability. We show that quasicontacts play an even more important role in these packings. We also show that jammed lattices are a useful setting for studying the Edwards ensemble, which weights each mechanically stable configuration equally and does not account for dynamics. This ansatz fails to predict the power-law distribution of near-zero contact forces, $P(f)\sim f^θ$.

preprint2014arXiv

Maximally dense packings of two-dimensional convex and concave noncircular particles

Dense packings of hard particles have important applications in many fields, including condensed matter physics, discrete geometry and cell biology. In this paper, we employ a stochastic search implementation of the Torquato-Jiao Adaptive-Shrinking-Cell optimization scheme [Nature 460, 876 (2009)] to find maximally dense particle packings in d-dimensional Euclidean space $\mathbb{R}^d$. While the original implementation was designed to study spheres and convex polyhedra in $d\ge3$, our implementation focuses on $d=2$ and extends the algorithm to include both concave polygons and certain convex or concave non-polygonal particle shapes. We verify the robustness of this packing protocol by successfully reproducing the known putative optimal packings of congruent copies of regular pentagons and octagons, then employ it to suggest dense packing arrangements of congruent copies of certain families of concave crosses, convex and concave curved triangles (incorporating shapes resembling the Mercedes-Benz logo), and "moon-like" shapes. Analytical constructions are determined subsequently to obtain the densest known packings of these particle shapes. For the examples considered, we find that the densest packings of both convex and concave particles with central symmetry are achieved by their corresponding optimal Bravais lattice packings; for particles lacking central symmetry, the densest packings obtained are non-lattice periodic packings, which are consistent with recently-proposed general organizing principles for hard particles. Moreover, we find that the densest known packings of certain curved triangles are periodic with a four-particle basis, and we find that the densest known periodic packings of certain moon-like shapes possess no inherent symmetries. Our work adds to the growing evidence that particle shape can be used as a tuning parameter to achieve a diversity of packing structures.

preprint2013arXiv

An Efficient Linear Programming Algorithm to Generate the Densest Lattice Sphere Packings

Finding the densest sphere packing in $d$-dimensional Euclidean space $\mathbb{R}^d$ is an outstanding fundamental problem with relevance in many fields, including the ground states of molecular systems, colloidal crystal structures, coding theory, discrete geometry, number theory, and biological systems. Numerically generating the densest sphere packings becomes very challenging in high dimensions due to an exponentially increasing number of possible sphere contacts and sphere configurations, even for the restricted problem of finding the densest lattice sphere packings. In this paper, we apply the Torquato-Jiao packing algorithm, which is a method based on solving a sequence of linear programs, to robustly reproduce the densest known lattice sphere packings for dimensions 2 through 19. We show that the TJ algorithm is appreciably more efficient at solving these problems than previously published methods. Indeed, in some dimensions, the former procedure can be as much as three orders of magnitude faster at finding the optimal solutions than earlier ones. We also study the suboptimal local density-maxima solutions (inherent structures or "extreme" lattices) to gain insight about the nature of the topography of the "density" landscape.

preprint2013arXiv

Effect of Dimensionality on the Percolation Threshold of Overlapping Nonspherical Hyperparticles

A set of lower bounds on the continuum percolation threshold $η_c$ of overlapping convex hyperparticles of general nonspherical (anisotropic) shape with a specified orientational probability distribution in $d$-dimensional Euclidean space have been derived [S. Torquato, J. Chem. Phys. {\bf 136}, 054106 (2012)]. The simplest of these lower bounds is given by $η_c \ge v/v_{ex}$, where $v_{ex}$ is the $d$-dimensional exclusion volume of a hyperparticle and $v$ is its $d$-dimensional volume. In order to study the effect of dimensionality on the threshold $η_c$ of overlapping nonspherical convex hyperparticles with random orientations here, we obtain a scaling relation for $η_c$ that is based on this lower bound and a conjecture that hyperspheres provide the highest threshold among all convex hyperparticle shapes for any $d$. This scaling relation exploits the principle that low-dimensional continuum percolation behavior encodes high-dimensional information. We derive a formula for the exclusion volume $v_{ex}$ of a hyperparticle in terms of its $d$-dimensional volume $v$, surface area $s$ and {\it radius of mean curvature} ${\bar R}$ (or, equivalently, {\it mean width}). These basic geometrical properties are computed for a wide variety of nonspherical hyperparticle shapes with random orientations across all dimensions, including, among other shapes, various polygons for $d=2$, Platonic solids, spherocylinders, parallepipeds and zero-volume plates for $d=3$ and their appropriate generalizations for $d \ge 4$. We then compute the lower bound and scaling relation for $η_c$ for this comprehensive set of continuum percolation models across dimensions. We show that the scaling relation provides accurate {\it upper-bound} estimates of the threshold $η_c$ across dimensions and becomes increasingly accurate as the space $d$ increases.

preprint2013arXiv

Effect of Dimensionality on the Percolation Thresholds of Various $d$-Dimensional Lattices

We show analytically that the $[0,1]$, $[1,1]$ and $[2,1]$ Pad{é} approximants of the mean cluster number $S(p)$ for site and bond percolation on general $d$-dimensional lattices are upper bounds on this quantity in any Euclidean dimension $d$, where $p$ is the occupation probability. These results lead to certain lower bounds on the percolation threshold $p_c$ that become progressively tighter as $d$ increases and asymptotically exact as $d$ becomes large. These lower-bound estimates depend on the structure of the $d$-dimensional lattice and whether site or bond percolation is being considered. We obtain explicit bounds on $p_c$ for both site and bond percolation on five different lattices: $d$-dimensional generalizations of the simple-cubic, body-centered-cubic and face-centered-cubic Bravais lattices as well as the $d$-dimensional generalizations of the diamond and kagom{é} (or pyrochlore) non-Bravais lattices. These analytical estimates are used to assess available simulation results across dimensions (up through $d=13$ in some cases). It is noteworthy that the tightest lower bound provides reasonable estimates of $p_c$ in relatively low dimensions and becomes increasingly accurate as $d$ grows. We also derive high-dimensional asymptotic expansions for $p_c$ for the ten percolation problems and compare them to the Bethe-lattice approximation. Finally, we remark on the radius of convergence of the series expansion of $S$ in powers of $p$ as the dimension grows.

preprint2013arXiv

Extreme lattices: symmetries and decorrelations

We study statistical and structural properties of extreme lattices, which are the local minima in the density landscape of lattice sphere packings, in $d$-dimensional Euclidean space $\mathbb{R}^d$. Specifically, we ascertain the distributions of densities, kissing numbers and numbers of symmetries of the packings across a wide range of dimensions using the stochastic Voronoi algorithm. The degree to which the packings decorrelate as well as the correlations between the density maxima as the space dimension increases is also investigated. We find that the extreme lattices decorrelate with increasing dimension, the least symmetric lattices decorrelate faster. The extreme lattices in a fixed dimension of space $d$ ($d\geq 8$) are dominated by typical lattices that have similar packing properties, like packing densities and kissing numbers, while the best and the worst packers are in the long tails of the distribution of the extreme lattices.

preprint2013arXiv

Isotropic Band Gaps and Freeform Waveguides Observed in Hyperuniform Disordered Photonic Solids

Recently, disordered photonic media and random textured surfaces have attracted increasing attention as strong light diffusers with broadband and wide-angle properties. We report the first experimental realization of an isotropic complete photonic band gap (PBG) in a two-dimensional (2D) disordered dielectric structure. This structure is designed by a constrained-optimization method, which combines advantages of both isotropy due to disorder and controlled scattering properties due to low density fluctuations (hyperuniformity) and uniform local topology. Our experiments use a modular design composed of Al2O3 walls and cylinders arranged in a hyperuniform disordered network. We observe a complete PBG in the microwave region, in good agreement with theoretical simulations, and show that the intrinsic isotropy of this novel class of PBG materials enables remarkable design freedom, including the realization of waveguides with arbitrary bending angles impossible in photonic crystals. This first experimental verification of a complete PBG and realization of functional defects in this new class of materials demonstrates their potential as building blocks for precise manipulation of photons in planar optical micro-circuits and has implications for disordered acoustic and electronic bandgap materials.

preprint2013arXiv

Optical cavities and waveguides in hyperuniform disordered photonic solids

Using finite difference time domain and band structure computer simulations, we show that it is possible to construct optical cavities and waveguide architectures in hyperuniform disordered photonic solids that are unattainable in photonic crystals. The cavity modes can be classified according to the symmetry (monopole, dipole, quadrupole,etc.) of the confined electromagnetic wave pattern. Owing to the isotropy of the band gaps characteristic of hyperuniform disordered solids, high-quality waveguides with freeform geometries (e.g., arbitrary bending angles) can be constructed that have no analogue in periodic or quasiperiodic solids. These capabilities have implications for many photonic applications.

preprint2013arXiv

Photonic band gap in isotropic hyperuniform disordered solids with low dielectric contrast

We report the first experimental demonstration of a TE-polarization photonic band gap (PBG) in a 2D isotropic hyperuniform disordered solid (HUDS) made of dielectric media with a index contrast of 1.6:1, very low for PBG formation. The solid is composed of a connected network of dielectric walls enclosing air-filled cells. Direct comparison with photonic crystals and quasicrystals permitted us to investigate band-gap properties as a function of increasing rotational isotropy. We present results from numerical simulations proving that the PBG observed experimentally for HUDS at low index contrast has zero density of states. The PBG is associated with the energy difference between complementary resonant modes above and below the gap, with the field predominantly concentrated in the air or in the dielectric. The intrinsic isotropy of HUDS may offer unprecedented flexibilities and freedom in applications (i. e. defect architecture design) not limited by crystalline symmetries.

preprint2012arXiv

Densest binary sphere packings

The densest binary sphere packings in the alpha-x plane of small to large sphere radius ratio alpha and small sphere relative concentration x have historically been very difficult to determine. Previous research had led to the prediction that these packings were composed of a few known "alloy" phases including, for example, the AlB2 (hexagonal omega), HgBr2, and AuTe2 structures, and to XYn structures composed of close-packed large spheres with small spheres (in a number ratio of n to 1) in the interstices, e.g., the NaCl packing for n = 1. However, utilizing an implementation of the Torquato-Jiao sphere-packing algorithm [S. Torquato and Y. Jiao, Phys. Rev. E 82, 061302 (2010)], we have discovered that many more structures appear in the densest packings. For example, while all previously known densest structures were composed of spheres in small to large number ratios of one to one, two to one, and very recently three to one, we have identified densest structures with number ratios of seven to three and five to two. In a recent work [A. B. Hopkins, Y. Jiao, F. H. Stillinger, and S. Torquato, Phys. Rev. Lett. 107, 125501 (2011)], we summarized these findings. In this work, we present the structures of the densest-known packings and provide details about their characteristics. Our findings demonstrate that a broad array of different densest mechanically stable structures consisting of only two types of components can form without any consideration of attractive or anisotropic interactions. In addition, the novel structures that we have identified may correspond to currently unidentified stable phases of certain binary atomic and molecular systems, particularly at high temperatures and pressures.

preprint2012arXiv

Density of States for a Specified Correlation Function and the Energy Landscape

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

preprint2012arXiv

Designed Diamond Ground State via Optimized Isotropic Monotonic Pair Potentials

We apply inverse statistical-mechanical methods to find a simple family of optimized isotropic, monotonic pair potentials, under certain constraints, whose ground states for a wide range of pressures is the diamond crystal. These constraints include desirable phonon spectra and the widest possible pressure range for stability. We also ascertain the ground-state phase diagram for a specific optimized potential to show that other crystal structures arise for other pressures. Cooling disordered configurations interacting with our optimized potential to absolute zero frequently leads to the desired diamond crystal ground state, revealing that the capture basin for the global energy minimum is large and broad relative to the local energy minima basins.

preprint2012arXiv

Effect of Dimensionality on the Continuum Percolation of Overlapping Hyperspheres and Hypercubes: II. Simulation Results and Analyses

In the first paper of this series [S. Torquato, J. Chem. Phys. {\bf 136}, 054106 (2012)], analytical results concerning the continuum percolation of overlapping hyperparticles in $d$-dimensional Euclidean space $\mathbb{R}^d$ were obtained, including lower bounds on the percolation threshold. In the present investigation, we provide additional analytical results for certain cluster statistics, such as the concentration of $k$-mers and related quantities, and obtain an upper bound on the percolation threshold $η_c$. We utilize the tightest lower bound obtained in the first paper to formulate an efficient simulation method, called the {\it rescaled-particle} algorithm, to estimate continuum percolation properties across many space dimensions with heretofore unattained accuracy. This simulation procedure is applied to compute the threshold $η_c$ and associated mean number of overlaps per particle ${\cal N}_c$ for both overlapping hyperspheres and oriented hypercubes for $ 3 \le d \le 11$. These simulations results are compared to corresponding upper and lower bounds on these percolation properties. We find that the bounds converge to one another as the space dimension increases, but the lower bound provides an excellent estimate of $η_c$ and ${\cal N}_c$, even for relatively low dimensions. We confirm a prediction of the first paper in this series that low-dimensional percolation properties encode high-dimensional information. We also show that the concentration of monomers dominate over concentration values for higher-order clusters (dimers, trimers, etc.) as the space dimension becomes large. Finally, we provide accurate analytical estimates of the pair connectedness function and blocking function at their contact values for any $d$ as a function of density.

preprint2012arXiv

Exotic Ground States of Directional Pair Potentials via Collective-Density Variables

Collective-density variables have proved to be a useful tool in the prediction and manipulation of how spatial patterns form in the classical many-body problem. Previous work has employed properties of collective-density variables along with a robust numerical optimization technique to find the classical ground states of many-particle systems subject to radial pair potentials in one, two and three dimensions. That work led to the identification of ordered and disordered classical ground states. In this paper, we extend these collective-coordinate studies by investigating the ground states of directional pair potentials in two dimensions. Our study focuses on directional potentials whose Fourier representations are non-zero on compact sets that are symmetric with respect to the origin and zero everywhere else. We choose to focus on one representative set which has exotic ground-state properties: two circles whose centers are separated by some fixed distance. We obtain ground states for this "two-circle" potential that display large void regions in the disordered regime. As more degrees of freedom are constrained the ground states exhibit a collapse of dimensionality characterized by the emergence of filamentary structures and linear chains. This collapse of dimensionality has not been observed before in related studies.

preprint2012arXiv

Improved reconstructions of random media using dilation and erosion processes

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

preprint2012arXiv

Microstructural Degeneracy associated with a Two-Point Correlation Function and its Information Content

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

preprint2012arXiv

Nonequilibrium static diverging length scales on approaching a prototypical model glassy state

Maximally random jammed states of hard spheres are prototypical glasses. We study the small wavenumber $k$ behavior of the structure factor $S(k)$ of overcompressed million-sphere packings as a function of density up to the jammed state. We find both a precursor to the glassy jammed state evident long before the jamming density is reached and an associated growing length scale, extracted from the volume integral of the direct correlation function $c(r)$, which diverges at the "critical" jammed state. We also define a nonequilibrium index $X$ and use it to demonstrate that the packings studied are intrinsically nonequilibrium in nature well before the critical point is reached.

preprint2012arXiv

Nonequilibrium static growing length scales in supercooled liquids on approaching the glass transition

The small wavenumber $k$ behavior of the structure factor $S(k)$ of overcompressed amorphous hard-sphere configurations was previously studied for a wide range of densities up to the maximally random jammed state, which can be viewed as a prototypical glassy state [A. Hopkins, F. H. Stillinger and S. Torquato, Phys. Rev. E, 86, 021505 (2012)]. It was found that a precursor to the glassy jammed state was evident long before the jamming density was reached as measured by a growing nonequilibrium length scale extracted from the volume integral of the direct correlation function $c(r)$, which becomes long-ranged as the critical jammed state is reached. The present study extends that work by investigating via computer simulations two different atomic models: the single-component Z2 Dzugutov potential in three dimensions and the binary-mixture Kob-Andersen potential in two dimensions. Consistent with the aforementioned hard-sphere study, we demonstrate that for both models a signature of the glass transition is apparent well before the transition temperature is reached as measured by the length scale determined from from the volume integral of the direct correlation function in the single-component case and a generalized direct correlation function in the binary-mixture case. The latter quantity is obtained from a generalized Orstein-Zernike integral equation for a certain decoration of the atomic point configuration. We also show that these growing length scales, which are a consequence of the long-range nature of the direct correlation functions, are intrinsically nonequilibrium in nature as determined by an index $X$ that is a measure of deviation from thermal equilibrium. It is also demonstrated that this nonequilibrium index, which increases upon supercooling, is correlated with a characteristic relaxation time scale.

preprint2012arXiv

Organizing Principles for Dense Packings of Nonspherical Hard Particles: Not All Shapes Are Created Equal

We have recently devised organizing principles to obtain maximally dense packings of the Platonic and Archimedean solids, and certain smoothly-shaped convex nonspherical particles [Torquato and Jiao, Phys. Rev. E 81, 041310 (2010)]. Here we generalize them in order to guide one to ascertain the densest packings of other convex nonspherical particles as well as concave shapes. Our generalized organizing principles are explicitly stated as four distinct propositions. We apply and test all of these organizing principles to the most comprehensive set of both convex and concave particle shapes to date, including Catalan solids, prisms, antiprisms, cylinders, dimers of spheres and various concave polyhedra. We demonstrate that all of the densest known packings associated with this wide spectrum of nonspherical particles are consistent with our propositions. Among other applications, our general organizing principles enable us to construct analytically the densest known packings of certain convex nonspherical particles, including spherocylinders, "lens-shaped" particles, square pyramids and rhombic pyramids. Moreover, we show how to apply these principles to infer the high-density equilibrium crystalline phases of hard convex and concave particles. We also discuss the unique packing attributes of maximally random jammed packings of nonspherical particles.

preprint2012arXiv

Quantitative Characterization of the Microstructure and Transport Properties of Biopolymer Networks

Biopolymer networks are of fundamental importance to many biological processes in normal and tumorous tissues. In this paper, we employ the panoply of theoretical and simulation techniques developed for characterizing heterogeneous materials to quantify the microstructure and effective diffusive transport properties (diffusion coefficient $D_e$ and mean survival time $τ$) of collagen type I networks at various collagen concentrations. In particular, we compute the pore-size probability density function $P(δ)$ for the networks and present a variety of analytical estimates of the effective diffusion coefficient $D_e$ for finite-sized diffusing particles. The Hashin-Strikman upper bound on the effective diffusion coefficient $D_e$ and the pore-size lower bound on the mean survival time $τ$ are used as benchmarks to test our analytical approximations and numerical results. Moreover, we generalize the efficient first-passage-time techniques for Brownian-motion simulations in suspensions of spheres to the case of fiber networks and compute the associated effective diffusion coefficient $D_e$ as well as the mean survival time $τ$, which is related to nuclear magnetic resonance (NMR) relaxation times. Specifically, the Torquato approximation provides the most accurate estimates of $D_e$ for all collagen concentrations among all of the analytical approximations we consider. We formulate a universal curve for $τ$ for the networks at different collagen concentrations. We apply rigorous cross-property relations to estimate the effective bulk modulus of collagen networks from a knowledge of the effective diffusion coefficient computed here.

preprint2012arXiv

Rigidity of spherical codes

A packing of spherical caps on the surface of a sphere (that is, a spherical code) is called rigid or jammed if it is isolated within the space of packings. In other words, aside from applying a global isometry, the packing cannot be deformed. In this paper, we systematically study the rigidity of spherical codes, particularly kissing configurations. One surprise is that the kissing configuration of the Coxeter-Todd lattice is not jammed, despite being locally jammed (each individual cap is held in place if its neighbors are fixed); in this respect, the Coxeter-Todd lattice is analogous to the face-centered cubic lattice in three dimensions. By contrast, we find that many other packings have jammed kissing configurations, including the Barnes-Wall lattice and all of the best kissing configurations known in four through twelve dimensions. Jamming seems to become much less common for large kissing configurations in higher dimensions, and in particular it fails for the best kissing configurations known in 25 through 31 dimensions. Motivated by this phenomenon, we find new kissing configurations in these dimensions, which improve on the records set in 1982 by the laminated lattices.

preprint2011arXiv

High-dimensional generalizations of the kagome and diamond crystals and the decorrelation principle for periodic sphere packings

In this paper, we introduce constructions of the high-dimensional generalizations of the kagome and diamond crystals. The two-dimensional kagome crystal and its three-dimensional counterpart, the pyrochlore crystal, have been extensively studied in the context of geometric frustration in antiferromagnetic materials. Similarly, the polymorphs of elemental carbon include the diamond crystal and the corresponding two-dimensional honeycomb structure, adopted by graphene. The kagome crystal in d Euclidean dimensions consists of vertex-sharing d-dimensional simplices in which all of the points are topologically equivalent. The d-dimensional generalization of the diamond crystal can then be obtained from the centroids of each of the simplices, and we show that this natural construction of the diamond crystal is distinct from the D_d^+ family of crystals for all dimensions d\neq 3. We analyze the structural properties of these high-dimensional crystals, including the packing densities, coordination numbers, void exclusion probability functions, covering radii, and quantizer errors. Our results demonstrate that the so-called decorrelation principle, which formally states that unconstrained correlations vanish in asymptotically high dimensions, remarkably applies to the case of periodic point patterns with inherent long-range order. We show that the decorrelation principle is already exhibited in periodic crystals in low dimensions via a "smoothed" pair correlation function obtained by convolution with a Gaussian kernel. These observations support the universality of the decorrelation principle for any point pattern in high dimensions, whether disordered or not. This universal property in turn suggests that the best conjectural lower bound on the maximal sphere-packing density in high Euclidean dimensions derived by Torquato and Stillinger [Expt. Math. 15, 307 (2006)] is in fact optimal.

preprint2011arXiv

Hyperuniform long-range correlations are a signature of disordered jammed hard-particle packings

We show that quasi-long-range (QLR) pair correlations that decay asymptotically with scaling $r^{-(d+1)}$ in $d$-dimensional Euclidean space $\mathbb{R}^d$, trademarks of certain quantum systems and cosmological structures, are a universal signature of maximally random jammed (MRJ) hard-particle packings. We introduce a novel hyperuniformity descriptor in MRJ packings by studying local-volume-fraction fluctuations and show that infinite-wavelength fluctuations vanish even for packings with size- and shape-distributions. Special void statistics induce hyperuniformity and QLR pair correlations.

preprint2011arXiv

Hyperuniformity, quasi-long-range correlations, and void-space constraints in maximally random jammed particle packings. I. Polydisperse spheres

Hyperuniform many-particle distributions possess a local number variance that grows more slowly than the volume of an observation window, implying that the local density is effectively homogeneous beyond a few characteristic length scales. Previous work on maximally random strictly jammed sphere packings in three dimensions has shown that these systems are hyperuniform and possess unusual quasi-long-range pair correlations, resulting in anomalous logarithmic growth in the number variance. However, recent work on maximally random jammed sphere packings with a size distribution has suggested that such quasi-long-range correlations and hyperuniformity are not universal among jammed hard-particle systems. In this paper we show that such systems are indeed hyperuniform with signature quasi-long-range correlations by characterizing the more general local-volume-fraction fluctuations. We argue that the regularity of the void space induced by the constraints of saturation and strict jamming overcomes the local inhomogeneity of the disk centers to induce hyperuniformity in the medium with a linear small-wavenumber nonanalytic behavior in the spectral density, resulting in quasi-long-range spatial correlations. A numerical and analytical analysis of the pore-size distribution for a binary MRJ system in addition to a local characterization of the n-particle loops governing the void space surrounding the inclusions is presented in support of our argument. This paper is the first part of a series of two papers considering the relationships among hyperuniformity, jamming, and regularity of the void space in hard-particle packings.

preprint2011arXiv

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

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

preprint2011arXiv

Inherent Structures for Soft Long-Range Interactions in Two-Dimensional Many-Particle Systems

We generate inherent structures, local potential-energy minima, of the "$k$-space overlap potential" in two-dimensional many-particle systems using a cooling and quenching simulation technique. The ground states associated with the $k$-space overlap potential are stealthy ({\it i.e.,} completely suppress single scattering of radiation for a range of wavelengths) and hyperuniform ({\it i.e.,} infinite wavelength density fluctuations vanish). However, we show via quantitative metrics that the inherent structures exhibit a range of stealthiness and hyperuniformity depending on the fraction of degrees of freedom that are constrained. Inherent structures in two dimensions typically contain five-particle rings, wavy grain boundaries, and vacancy-interstitial defects. The structural and thermodynamic properties of inherent structures are relatively insensitive to the temperature from which they are sampled, signifying that the energy landscape is relatively flat and devoid of deep wells. Using the nudged-elastic-band algorithm, we construct paths from ground-state configurations to inherent structures and identify the transition points between them. In addition, we use point patterns generated from a random sequential addition (RSA) of hard disks, which are nearly stealthy, and examine the particle rearrangements necessary to make the configurations absolutely stealthy. We introduce a configurational proximity metric to show that only small local, but collective, particle rearrangements are needed to drive initial RSA configurations to stealthy disordered ground states. These results lead to a more complete understanding of the unusual behaviors exhibited by the family of "collective-coordinate" potentials to which the $k$-space overlap potential belongs.

preprint2011arXiv

Phase diagram and structural diversity of the densest binary sphere packings

The densest binary sphere packings have historically been very difficult to determine. The only rigorously known packings in the alpha-x plane of sphere radius ratio alpha and relative concentration x are at the Kepler limit alpha = 1, where packings are monodisperse. Utilizing an implementation of the Torquato-Jiao sphere-packing algorithm [S. Torquato and Y. Jiao, Phys. Rev. E 82, 061302 (2010)], we present the most comprehensive determination to date of the phase diagram in (alpha,x) for the densest binary sphere packings. Unexpectedly, we find many distinct new densest packings.

preprint2010arXiv

Anomalous local coordination, density fluctuations, and void statistics in disordered hyperuniform many-particle ground states

We provide numerical constructions of one-dimensional hyperuniform many-particle distributions that exhibit unusual clustering and asymptotic local number density fluctuations growing more slowly than the volume of an observation window but faster than the surface area. By targeting a specified form of the structure factor at small wavenumbers using collective density variables, we are able to tailor the form of asymptotic local density fluctuations while simultaneously measuring the effect of imposing weak and strong constraints on the available degrees of freedom within the system. Even in one dimension, the long-range effective interactions induce clustering and nontrivial phase transitions in the resulting ground-state configurations. We provide an analytical connection between the fraction of contrained degrees of freedom within the system and the disorder-order phase transition for a class of target structure factors by examining the realizability of the constrained contribution to the pair correlation function. Our results explicitly demonstrate that disordered hyperuniform many-particle ground states, and therefore also point distributions, with substantial clustering can be constructed. We directly relate the local coordination structure of our point patterns to the distribution of the void space external to the particles, and we provide a scaling argument for the configurational entropy of the systems when only a small fraction of the degrees of freedom are constrained. By emphasizing the intimate connection between geometrical constraints on the particle distribution and structural regularity, our work has direct implications for higher-dimensional systems, including an understanding of the appearance of hyperuniformity and quasi-long-range pair correlations in maximally random strictly jammed packings of hard spheres.

preprint2010arXiv

Complete Band Gaps in 2D Photonic Quasicrystals

We introduce a novel optimization method to design the first examples of photonic quasicrystals with substantial, complete photonic band gaps (PBGs): that is, a range of frequencies over which electromagnetic wave propagation is forbidden for all directions and polarizations. The method can be applied to photonic quasicrystals with arbitrary rotational symmetry; here, we illustrate the results for 5- and 8-fold symmetric quasicrystals. The optimized band gaps are highly isotropic, which may offer advantages over photonic crystals for certain applications.

preprint2010arXiv

Densest local packing diversity. II. Application to three dimensions

The densest local packings of N three-dimensional identical nonoverlapping spheres within a radius Rmin(N) of a fixed central sphere of the same size are obtained for selected values of N up to N = 1054. In the predecessor to this paper [A.B. Hopkins, F.H. Stillinger and S. Torquato, Phys. Rev. E 81 041305 (2010)], we described our method for finding the putative densest packings of N spheres in d-dimensional Euclidean space Rd and presented those packings in R2 for values of N up to N = 348. We analyze the properties and characteristics of the densest local packings in R3 and employ knowledge of the Rmin(N), using methods applicable in any d, to construct both a realizability condition for pair correlation functions of sphere packings and an upper bound on the maximal density of infinite sphere packings. In R3, we find wide variability in the densest local packings, including a multitude of packing symmetries such as perfect tetrahedral and imperfect icosahedral symmetry. We compare the densest local packings of N spheres near a central sphere to minimal-energy configurations of N+1 points interacting with short-range repulsive and long-range attractive pair potentials, e.g., 12-6 Lennard-Jones, and find that they are in general completely different, a result that has possible implications for nucleation theory. We also compare the densest local packings to finite subsets of stacking variants of the densest infinite packings in R3 (the Barlow packings) and find that the densest local packings are almost always most similar, as measured by a similarity metric, to the subsets of Barlow packings with the smallest number of coordination shells measured about a single central sphere, e.g., a subset of the FCC Barlow packing. We additionally observe that the densest local packings are dominated by the spheres arranged with centers at precisely distance Rmin(N) from the fixed sphere's center.

preprint2010arXiv

Designer disordered materials with large complete photonic band gaps

We present designs of 2D isotropic, disordered photonic materials of arbitrary size with complete band gaps blocking all directions and polarizations. The designs with the largest gaps are obtained by a constrained optimization method that starts from a hyperuniform disordered point pattern, an array of points whose number variance within a spherical sampling window grows more slowly than the volume. We argue that hyperuniformity, combined with uniform local topology and short-range geometric order, can explain how complete photonic band gaps are possible without long-range translational order. We note the ramifications for electronic and phononic band gaps in disordered materials.

preprint2010arXiv

Duality Relations for the Classical Ground States of Soft-Matter Systems

Bounded interactions are particularly important in soft-matter systems, such as colloids, microemulsions, and polymers. We derive new duality relations for a class of soft potentials, including three-body and higher-order functions, that can be applied to ordered and disordered classical ground states. These duality relations link the energy of configurations associated with a real-space potential to the corresponding energy of the dual (Fourier-transformed) potential. We apply the duality relations by demonstrating how information about the classical ground states of short-ranged potentials can be used to draw new conclusions about the ground states of long-ranged potentials and vice versa. The duality relations also lead to bounds on the T=0 system energies in density intervals of phase coexistence. Additionally, we identify classes of "self-similar" potentials, for which one can relate low- and high-density ground-state energies. We analyze the ground state configurations and thermodynamic properties of a one-dimensional system previously thought to exhibit an infinite number of structural phase transitions and comment on the known ground states of purely repulsive monotonic potentials in the context of our duality relations.

preprint2010arXiv

Effects of Random Link Removal on the Photonic Band Gaps of Honeycomb Networks

We explore the effects of random link removal on the photonic band gaps of honeycomb networks. Missing or incomplete links are expected to be common in practical realizations of this class of connected network structures due to unavoidable flaws in the fabrication process. We focus on the collapse of the photonic band gap due to the defects induced by the link removal. We show that the photonic band gap is quite robust against this type of random decimation and survives even when almost 58% of the network links are removed.

preprint2010arXiv

Jammed Hard-Particle Packings: From Kepler to Bernal and Beyond

This review describes the diversity of jammed configurations attainable by frictionless convex nonoverlapping (hard) particles in Euclidean spaces and for that purpose it stresses individual-packing geometric analysis. A fundamental feature of that diversity is the necessity to classify individual jammed configurations according to whether they are locally, collectively, or strictly jammed. Each of these categories contains a multitude of jammed configurations spanning a wide and (in the large system limit) continuous range of intensive properties, including packing fraction $ϕ$, mean contact number $Z$, and several scalar order metrics. Application of these analytical tools to spheres in three dimensions (an analog to the venerable Ising model) covers a myriad of jammed states, including maximally dense packings (as Kepler conjectured), low-density strictly-jammed tunneled crystals, and a substantial family of amorphous packings. With respect to the last of these, the current approach displaces the historically prominent but ambiguous idea of ``random close packing" (RCP) with the precise concept of ``maximally random jamming" (MRJ). This review also covers recent advances in understanding jammed packings of polydisperse sphere mixtures, as well as convex nonspherical particles, e.g., ellipsoids, ``superballs", and polyhedra. Because of their relevance to error-correcting codes and information theory, sphere packings in high-dimensional Euclidean spaces have been included as well. We also make some remarks about packings in (curved) non-Euclidean spaces. In closing this review, several basic open questions for future research to consider have been identified.

preprint2010arXiv

Novel Ground-State Crystals with Controlled Vacancy Concentrations: From Kagomé to Honeycomb to Stripes

We introduce a one-parameter family, $0 \leq H \leq 1$, of pair potential functions with a single relative energy minimum that stabilize a range of vacancy-riddled crystals as ground states. The "quintic potential" is a short-ranged, nonnegative pair potential with a single local minimum of height $H$ at unit distance and vanishes cubically at a distance of $\rt$. We have developed this potential to produce ground states with the symmetry of the triangular lattice while favoring the presence of vacancies. After an exhaustive search using various optimization and simulation methods, we believe that we have determined the ground states for all pressures, densities, and $0 \leq H \leq 1$. For specific areas below $3\rt/2$, the ground states of the "quintic potential" include high-density and low-density triangular lattices, kagomé and honeycomb crystals, and stripes. We find that these ground states are mechanically stable but are difficult to self-assemble in computer simulations without defects. For specific areas above $3\rt/2$, these systems have a ground-state phase diagram that corresponds to hard disks with radius $\rt$. For the special case of H=0, a broad range of ground states is available. Analysis of this case suggests that among many ground states, a high-density triangular lattice, low-density triangular lattice, and striped phases have the highest entropy for certain densities. The simplicity of this potential makes it an attractive candidate for experimental realization with application to the development of novel colloidal crystals or photonic materials.

preprint2010arXiv

Optimized Monotonic Convex Pair Potentials Stabilize Low-Coordinated Crystals

We have previously used inverse statistical-mechanical methods to optimize isotropic pair interactions with multiple extrema to yield low-coordinated crystal classical ground states (e.g., honeycomb and diamond structures) in d-dimensional Euclidean space R^d. Here we demonstrate the counterintuitive result that no extrema are required to produce such low-coordinated classical ground states. Specifically, we show that monotonic convex pair potentials can be optimized to yield classical ground states that are the square and honeycomb crystals in R^2 over a non-zero number density range. Such interactions may be feasible to achieve experimentally using colloids and polymers.

preprint2010arXiv

Phase Behavior of Colloidal Superballs: Shape Interpolation from Spheres to Cubes

The phase behavior of hard superballs is examined using molecular dynamics within a deformable periodic simulation box. A superball's interior is defined by the inequality $|x|^{2q} + |y|^{2q} + |z|^{2q} \leq 1$, which provides a versatile family of convex particles ($q \geq 0.5$) with cube-like and octahedron-like shapes as well as concave particles ($q < 0.5$) with octahedron-like shapes. Here, we consider the convex case with a deformation parameter q between the sphere point (q = 1) and the cube (q = 1). We find that the asphericity plays a significant role in the extent of cubatic ordering of both the liquid and crystal phases. Calculation of the first few virial coefficients shows that superballs that are visually similar to cubes can have low-density equations of state closer to spheres than to cubes. Dense liquids of superballs display cubatic orientational order that extends over several particle lengths only for large q. Along the ordered, high-density equation of state, superballs with 1 < q < 3 exhibit clear evidence of a phase transition from a crystal state to a state with reduced long-ranged orientational order upon the reduction of density. For $q \geq 3$, long-ranged orientational order persists until the melting transition. The width of coexistence region between the liquid and ordered, high-density phase decreases with q up to q = 4.0. The structures of the high-density phases are examined using certain order parameters, distribution functions, and orientational correlation functions. We also find that a fixed simulation cell induces artificial phase transitions that are out of equilibrium. Current fabrication techniques allow for the synthesis of colloidal superballs, and thus the phase behavior of such systems can be investigated experimentally.

preprint2010arXiv

Toward an Ising Model of Cancer and Beyond

Theoretical and computational tools that can be used in the clinic to predict neoplastic progression and propose individualized optimal treatment strategies to control cancer growth is desired. To develop such a predictive model, one must account for the complex mechanisms involved in tumor growth. Here we review resarch work that we have done toward the development of an "Ising model" of cancer. The review begins with a description of a minimalist four-dimensional (three in space and one in time) cellular automaton (CA) model of cancer in which healthy cells transition between states (proliferative, hypoxic, and necrotic) according to simple local rules and their present states, which can viewed as a stripped-down Ising model of cancer. This model is applied to model the growth of glioblastoma multiforme, the most malignant of brain cancers. This is followed by a discussion of the extension of the model to study the effect on the tumor dynamics and geometry of a mutated subpopulation. A discussion of how tumor growth is affected by chemotherapeutic treatment is then described. How angiogenesis as well as the heterogeneous and confined environment in which a tumor grows is incorporated in the CA model is discussed. The characterization of the level of organization of the invasive network around a solid tumor using spanning trees is subsequently described. Then, we describe open problems and future promising avenues for future research, including the need to develop better molecular-based models that incorporate the true heterogeneous environment over wide range of length and time scales (via imaging data), cell motility, oncogenes, tumor suppressor genes and cell-cell communication. The need to bring to bear the powerful machinery of the theory of heterogeneous media to better understand the behavior of cancer in its microenvironment is presented.

preprint2009arXiv

Hyperuniformity in point patterns and two-phase random heterogeneous media

Hyperuniform point patterns are characterized by vanishing infinite wavelength density fluctuations and encompass all crystal structures, certain quasi-periodic systems, and special disordered point patterns. This article generalizes the notion of hyperuniformity to include also two-phase random heterogeneous media. Hyperuniform random media do not possess infinite-wavelength volume fraction fluctuations, implying that the variance in the local volume fraction in an observation window decays faster than the reciprocal window volume as the window size increases. For microstructures of impenetrable and penetrable spheres, we derive an upper bound on the asymptotic coefficient governing local volume fraction fluctuations in terms of the corresponding quantity describing the variance in the local number density (i.e., number variance). Extensive calculations of the asymptotic number variance coefficients are performed for a number of disordered (e.g., quasiperiodic tilings, classical stealth disordered ground states, and certain determinantal point processes), quasicrystal, and ordered (e.g., Bravais and non-Bravais periodic systems) hyperuniform structures in various Euclidean space dimensions, and our results are consistent with a quantitative order metric characterizing the degree of hyperuniformity. We also present corresponding estimates for the asymptotic local volume fraction coefficients for several lattice families. Our results have interesting implications for a certain problem in number theory.

preprint2004arXiv

Superdense Crystal Packings of Ellipsoids

Particle packing problems have fascinated people since the dawn of civilization, and continue to intrigue mathematicians and scientists. Resurgent interest has been spurred by the recent proof of Kepler's conjecture: the face-centered cubic lattice provides the densest packing of equal spheres with a packing fraction $ϕ\approx0.7405$ \cite{Kepler_Hales}. Here we report on the densest known packings of congruent ellipsoids. The family of new packings are crystal (periodic) arrangements of nearly spherically-shaped ellipsoids, and always surpass the densest lattice packing. A remarkable maximum density of $ϕ\approx0.7707$ is achieved for both prolate and oblate ellipsoids with aspect ratios of $\sqrt{3}$ and $1/\sqrt{3}$, respectively, and each ellipsoid has 14 touching neighbors. Present results do not exclude the possibility that even denser crystal packings of ellipsoids could be found, and that a corresponding Kepler-like conjecture could be formulated for ellipsoids.

preprint2004arXiv

Voronoi and Voids Statistics for Super-homogeneous Point Processes

We study the Voronoi and void statistics of super-homogeneous (or hyperuniform) point patterns in which the infinite-wavelength density fluctuations vanish. Super-homogeneous or hyperuniform point patterns arise in one-component plasmas, primordial density fluctuations in the Universe, and in jammed hard-particle packings. We specifically analyze a certain one-dimensional model by studying size fluctuations and correlations of the associated Voronoi cells. We derive exact results for the complete joint statistics of the size of two Voronoi cells. We also provide a sum rule that the correlation matrix for the Voronoi cells must obey in any space dimension. In contrast to the conventional picture of super-homogeneous systems, we show that infinitely large Voronoi cells or voids can exist in super-homogeneous point processes in any dimension. We also present two heuristic conditions to identify and classify any super-homogeneous point process in terms of the asymptotic behavior of the void size distribution.