Point processes and stochastic displacement fields

Hyperuniformity classes of quasiperiodic tilings via diffusion spreadability

Hyperuniform point patterns can be classified by the hyperuniformity scaling exponent α>0, that characterizes the power-law scaling behavior of the structure factor S(k) as a function of wave number k≡|k| in the vicinity of the origin, e.g., S(k)∼|k|^{α} in cases where S(k) varies continuously with k as k→0. In this paper, we show that the spreadability is an effective method for determining α for quasiperiodic systems where S(k) is discontinuous and consists of a dense set of Bragg peaks. It has been shown in [Phys. Rev. E 104, 054102 (2021)10.1103/PhysRevE.104.054102] that, for media with finite α, the long-time behavior of the excess spreadability S(∞)-S(t) can be fit to a power law of the form t^{-(d-α)/2}, where d is the space dimension, to accurately extract α for the continuous case. We first transform quasiperiodic and limit-periodic point patterns into two-phase media by mapping them onto packings of identical nonoverlapping disks, where space interior to the disks represents one phase and the space in exterior to them represents the second phase. We then compute the spectral density χ[over ̃]_{_{V}}(k) of the packings, and finally compute and fit the long-time behavior of their excess spreadabilities. Specifically, we show that the excess spreadability can be used to accurately extract α for the one-dimensional (1D) limit-periodic period-doubling chain (α=1) and the 1D quasicrystalline Fibonacci chain (α=3) to within 0.02% of the analytically known exact results. Moreover, we obtain a value of α=5.97±0.06 for the two-dimensional Penrose tiling and present plausible theoretical arguments strongly suggesting that α is exactly equal to six. We also show that, due to the self-similarity of the structures examined here, one can truncate the small-k region of the scattering information used to compute the spreadability and obtain an accurate value of α, with a small deviation from the untruncated case that decreases as the system size increases. This strongly suggests that one can obtain a good estimate of α for an infinite self-similar quasicrystal from a modestly sized finite sample. The methods described here offer a simple and general procedure to characterize accurately the large-scale translational order present in quasicrystalline and limit-periodic media in any space dimension that are self-similar. Moreover, the scattering information extracted from these two-phase media encoded in χ[over ̃]_{_{V}}(k), can be used to estimate their physical properties, such as their effective dynamic dielectric constants, effective dynamic elastic constants, and fluid permeabilities.

Evidence of scale-free clusters of vegetation in tropical rainforests

Tropical rainforests exhibit a rich repertoire of spatial patterns emerging from the intricate relationship between the microscopic interaction between species. In particular, the distribution of vegetation clusters can shed much light on the underlying process that regulates the ecosystem. Analyzing the distribution of vegetation clusters at different resolution scales, we show the first robust evidence of scale-invariant clusters of vegetation, suggesting the coexistence of multiple intertwined scales in the collective dynamics of tropical rainforests. We use field data and computational simulations to confirm our hypothesis, proposing a predictor that could be particularly interesting to monitor the ecological resilience of the world's "green lungs."

Observation of magnetic structural universality and jamming transition with NMR

Nuclear magnetic resonance (NMR) has been instrumental in deciphering the structure of proteins. Here we show that transverse NMR relaxation, through its time-dependent relaxation rate, is distinctly sensitive to the structure of complex materials or biological tissues at the mesoscopic scale, from micrometers to tens of micrometers. Based on the ideas of universality, we show analytically and numerically that the time-dependent transverse relaxation rate approaches its long-time limit in a power-law fashion, with the dynamical exponent reflecting the universality class of mesoscopic magnetic structure. The spectral line shape acquires the corresponding non-analytic power law singularity at zero frequency. We experimentally detect the change in the dynamical exponent as a result of the transition into maximally random jammed state characterized by hyperuniform correlations. The relation between relaxational dynamics and magnetic structure opens the way for noninvasive characterization of porous media, complex materials and biological tissues.

Non-hyperuniform metastable states around a disordered hyperuniform state of densely packed spheres: stochastic density functional theory at strong coupling

The disordered and hyperuniform structures of densely packed spheres near and at jamming are characterized by vanishing of long-wavelength density fluctuations, or equivalently by long-range power-law decay of the direct correlation function (DCF). We focus on previous simulation results that exhibit the degradation of hyperuniformity in jammed structures while maintaining the long-range nature of the DCF to a certain length scale. Here we demonstrate that the field-theoretic formulation of stochastic density functional theory is relevant to explore the degradation mechanism. The strong-coupling expansion method of stochastic density functional theory is developed to obtain the metastable chemical potential considering the intermittent fluctuations in dense packings. The metastable chemical potential yields the analytical form of the metastable DCF that has a short-range cutoff inside the sphere while retaining the long-range power-law behavior. It is confirmed that the metastable DCF provides the zero-wavevector limit of the structure factor in quantitative agreement with the previous simulation results of degraded hyperuniformity. We can also predict the emergence of soft modes localized at the particle scale by plugging this metastable DCF into the linearized Dean-Kawasaki equation, a stochastic density functional equation.

Structural characterization of many-particle systems on approach to hyperuniform states

The study of hyperuniform states of matter is an emerging multidisciplinary field, impinging on topics in the physical sciences, mathematics, and biology. The focus of this work is the exploration of quantitative descriptors that herald when a many-particle system in d-dimensional Euclidean space R^{d} approaches a hyperuniform state as a function of the relevant control parameter. We establish quantitative criteria to ascertain the extent of hyperuniform and nonhyperuniform distance-scaling regimes as well as the crossover point between them in terms of the "volume" coefficient A and "surface-area" coefficient B associated with the local number variance σ^{2}(R) for a spherical window of radius R. The larger the ratio B/A, the larger the hyperuniform scaling regime, which becomes of infinite extent in the limit B/A→∞. To complement the known direct-space representation of the coefficient B in terms of the total correlation function h(r), we derive its corresponding Fourier representation in terms of the structure factor S(k), which is especially useful when scattering information is available experimentally or theoretically. We also demonstrate that the free-volume theory of the pressure of equilibrium packings of identical hard spheres that approach a strictly jammed state either along the stable crystal or metastable disordered branch dictates that such end states be exactly hyperuniform. Using the ratio B/A, as well as other diagnostic measures of hyperuniformity, including the hyperuniformity index H and the direct-correlation function length scale ξ_{c}, we study three different exactly solvable models as a function of the relevant control parameter, either density or temperature, with end states that are perfectly hyperuniform. Specifically, we analyze equilibrium systems of hard rods and "sticky" hard-sphere systems in arbitrary space dimension d as a function of density. We also examine low-temperature excited states of many-particle systems interacting with "stealthy" long-ranged pair interactions as the temperature tends to zero, where the ground states are disordered, hyperuniform, and infinitely degenerate. We demonstrate that our various diagnostic hyperuniformity measures are positively correlated with one another. The same diagnostic measures can be used to detect the degree to which imperfections in nearly hyperuniform systems cause deviations from perfect hyperuniformity. Moreover, the capacity to identify hyperuniform scaling regimes should be particularly useful in analyzing experimentally or computationally generated samples that are necessarily of finite size.

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 two-phase heterogeneous 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 σ_{_{V}}^{2}(R) and the associated hyperuniformity order metrics B[over ¯]_{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 purposely restrict ourselves to a certain class of two-dimensional periodic cellular networks as well as periodic and disordered or irregular packings of circular disks, some of which maximize their effective transport and elastic properties. Among the cellular networks considered, the honeycomb networks have minimal values of the hyperuniformity order metrics B[over ¯]_{V} across all volume fractions. On the other hand, for all packings of circular disks examined, the triangular-lattice packings have the smallest values of B[over ¯]_{V} for the possible range of volume fractions. Among all structures studied here, the triangular-lattice packing of circular disks have the minimal values of the 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 with desirable bulk physical properties by inverse design procedures.

Stone-Wales defects preserve hyperuniformity in amorphous two-dimensional networks

Disordered hyperuniformity (DHU) is a recently discovered novel state of many-body systems that possesses vanishing normalized infinite-wavelength density fluctuations similar to a perfect crystal and an amorphous structure like a liquid or glass. Here, we discover a hyperuniformity-preserving topological transformation in two-dimensional (2D) network structures that involves continuous introduction of Stone-Wales (SW) defects. Specifically, the static structure factor [Formula: see text] of the resulting defected networks possesses the scaling [Formula: see text] for small wave number k, where [Formula: see text] monotonically decreases as the SW defect concentration p increases, reaches [Formula: see text] at [Formula: see text], and remains almost flat beyond this p. Our findings have important implications for amorphous 2D materials since the SW defects are well known to capture the salient feature of disorder in these materials. Verified by recently synthesized single-layer amorphous graphene, our network models reveal unique electronic transport mechanisms and mechanical behaviors associated with distinct classes of disorder in 2D materials.

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 Λ[over ¯]. On the other hand, the disappearance and reemergence of long-range order, depending on whether the system is cloaked 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.

Methodology to construct large realizations of perfectly hyperuniform disordered packings

Disordered hyperuniform packings (or dispersions) are unusual amorphous two-phase materials that are endowed with exotic physical properties. Such hyperuniform systems are characterized by an anomalous suppression of volume-fraction fluctuations at infinitely long-wavelengths, compared to ordinary disordered materials. While there has been growing interest in such singular states of amorphous matter, a major obstacle has been an inability to produce large samples that are perfectly hyperuniform due to practical limitations of conventional numerical and experimental methods. To overcome these limitations, we introduce a general theoretical methodology to construct perfectly hyperuniform packings in d-dimensional Euclidean space R^{d}. Specifically, beginning with an initial general tessellation of space by disjoint cells that meets a "bounded-cell" condition, hard particles of general shape are placed inside each cell such that the local-cell particle packing fractions are identical to the global packing fraction. We prove that the constructed packings with a polydispersity in size are perfectly hyperuniform in the infinite-sample-size limit, regardless of particle shapes, positions, and numbers per cell. We use this theoretical formulation to devise an efficient and tunable algorithm to generate extremely large realizations of such packings. We employ two distinct initial tessellations: Voronoi as well as sphere tessellations. Beginning with Voronoi tessellations, we show that our algorithm can remarkably convert extremely large nonhyperuniform packings into hyperuniform ones in R^{2} and R^{3}. Implementing our theoretical methodology on sphere tessellations, we establish the hyperuniformity of the classical Hashin-Shtrikman multiscale coated-spheres structures, which are known to be two-phase media microstructures that possess optimal effective transport and elastic properties. A consequence of our work is a rigorous demonstration that packings that have identical tessellations can either be nonhyperuniform or hyperuniform by simply tuning local characteristics. It is noteworthy that our computationally designed hyperuniform two-phase systems can easily be fabricated via state-of-the-art methods, such as 2D photolithographic and 3D printing technologies. In addition, the tunability of our methodology offers a route for the discovery of novel disordered hyperuniform two-phase materials.

Structural disorder and anomalous diffusion in random packing of spheres

Nowadays Nuclear Magnetic Resonance diffusion (dNMR) measurements of water molecules in heterogeneous systems have broad applications in material science, biophysics and medicine. Up to now, microstructural rearrangement in media has been experimentally investigated by studying the diffusion coefficient (D(t)) behavior in the tortuosity limit. However, this method is not able to describe structural disorder and transitions in complex systems. Here we show that, according to the continuous time random walk framework, the dNMR measurable parameter α, quantifying the anomalous regime of D(t), provides a quantitative characterization of structural disorder and structural transition in heterogeneous systems. To demonstrate this, we compare α measurements obtained in random packed monodisperse micro-spheres with Molecular Dynamics simulations of disordered porous media and 3D Monte Carlo simulation of particles diffusion in these kind of systems. Experimental results agree well with simulations that correlate the most used parameters and functions characterizing the disorder in porous media.

One-dimensional gravity in infinite point distributions

The dynamics of infinite asymptotically uniform distributions of purely self-gravitating particles in one spatial dimension provides a simple and interesting toy model for the analogous three dimensional problem treated in cosmology. In this paper we focus on a limitation of such models as they have been treated so far in the literature: the force, as it has been specified, is well defined in infinite point distributions only if there is a centre of symmetry (i.e., the definition requires explicitly the breaking of statistical translational invariance). The problem arises because naive background subtraction (due to expansion, or by "Jeans swindle" for the static case), applied as in three dimensions, leaves an unregulated contribution to the force due to surface mass fluctuations. Following a discussion by Kiessling of the Jeans swindle in three dimensions, we show that the problem may be resolved by defining the force in infinite point distributions as the limit of an exponentially screened pair interaction. We show explicitly that this prescription gives a well defined (finite) force acting on particles in a class of perturbed infinite lattices, which are the point processes relevant to cosmological N -body simulations. For identical particles the dynamics of the simplest toy model (without expansion) is equivalent to that of an infinite set of points with inverted harmonic oscillator potentials which bounce elastically when they collide. We discuss and compare with previous results in the literature and present new results for the specific case of this simplest (static) model starting from "shuffled lattice" initial conditions. These show qualitative properties of the evolution (notably its "self-similarity") like those in the analogous simulations in three dimensions, which in turn resemble those in the expanding universe.

Tilings of space and superhomogeneous point processes

We consider the construction of point processes from tilings, with equal-volume tiles, of d -dimensional Euclidean space R;{d} . We show that one can generate, with simple algorithms ascribing one or more points to each tile, point processes which are "superhomogeneous" (or "hyperuniform")-i.e., for which the structure factor S(k) vanishes when the wave vector k tends to zero. The exponent gamma characterizing the leading small- k behavior, S(k-->0) proportional, variant k(gamma), depends in a simple manner on the nature of the correlation properties of the specific tiling and on the conservation of the mass moments of the tiles. Assigning one point to the center of mass of each tile gives the exponent gamma=4 for any tiling in which the shapes and orientations of the tiles are short-range correlated. Smaller exponents in the range 4-d<gamma<4 (and thus always superhomogeneous for d< or =4 ) may be obtained in the case that the latter quantities have long-range correlations. Assigning more than one point to each tile in an appropriate way, we show that one can obtain arbitrarily higher exponents in both cases. We illustrate our results with explicit constructions using known deterministic tilings, as well as some simple stochastic tilings for which we can calculate S(k) exactly. Our results provide an explicit analytical construction of point processes with gamma>4 . Applications to condensed matter physics, and also to cosmology, are briefly discussed.

Two-point correlation properties of stochastic splitting processes

We study how the two-point density correlation properties of a point particle distribution are modified when each particle is divided, by a stochastic process, into an equal number of identical "daughter" particles. We consider generically that there may be nontrivial correlations in the displacement fields describing the positions of the different daughters of the same "mother" particle and then treat separately the cases in which there are, or are not, correlations also between the displacements of daughters belonging to different mothers. For both cases exact formulas are derived relating the structure factor (power spectrum) of the daughter distribution to that of the mothers. An application of these results is that they give explicit algorithms for generating, starting from regular lattice arrays, stochastic particle distributions with an arbitrarily high degree of large-scale uniformity. Such distributions are of interest, in particular, in the context of studies of self-gravitating systems in cosmology.

Gravitational dynamics of an infinite shuffled lattice: Particle coarse-graining, nonlinear clustering, and the continuum limit

We study the evolution under their self-gravity of particles evolving from infinite "shuffled lattice" initial conditions. We focus here specifically on the comparison between the evolution of such a system and that of "daughter" coarse-grained particle distributions. These are sparser (i.e., lower density) particle distributions, defined by a simple coarse-graining procedure, which share the same large-scale mass fluctuations. We consider both the case that such coarse-grainings are performed (i) on the initial conditions, and (ii) at a finite time with a specific additional prescription. In numerical simulations we observe that, to a first approximation, these coarse-grainings represent well the evolution of the two-point correlation properties over a significant range of scales. We note, in particular, that the form of the two-point correlation function in the original system, when it is evolving in the asymptotic "self-similar" regime, may be reproduced well in a daughter coarse-grained system in which the dynamics are still dominated by two-body (nearest neighbor) interactions. This provides a simple physical description of the origin of the form of part of the asymptotic nonlinear correlation function. Using analytical results on the early time evolution of these systems, however, we show that small observed differences between the evolved system and its coarse-grainings at the initial time will in fact diverge as the ratio of the coarse-graining scale to the original interparticle distance increases. The second coarse-graining studied, performed at a finite time in a specified manner, circumvents this problem. It also makes it more physically transparent why gravitational dynamics from these initial conditions tends toward a self-similar evolution. We finally discuss the precise definition of a limit in which a continuum (specifically Vlasov-type) description of the observed linear and nonlinear evolution should be applicable. This requires the introduction of an additional intrinsic length scale (e.g., a physical smoothing in the force at small scales), which is kept fixed as the particle density diverges. In this limit the different coarse-grainings are equivalent and leave the evolution of the "mother" system invariant.

Gravitational dynamics of an infinite shuffled lattice of particles

We study, using numerical simulations, the dynamical evolution of self-gravitating point particles in static Euclidean space, starting from a simple class of infinite "shuffled lattice" initial conditions. These are obtained by applying independently to each particle on an infinite perfect lattice a small random displacement, and are characterized by a power spectrum (structure factor) of density fluctuations which is quadratic in the wave number k, at small k. For a specified form of the probability distribution function of the "shuffling" applied to each particle, and zero initial velocities, these initial configurations are characterized by a single relevant parameter: the variance delta(2) of the "shuffling" normalized in units of the lattice spacing l. The clustering, which develops in time starting from scales around l, is qualitatively very similar to that seen in cosmological simulations, which begin from lattices with applied correlated displacements and incorporate an expanding spatial background. From very soon after the formation of the first nonlinear structures, a spatiotemporal scaling relation describes well the evolution of the two-point correlations. At larger times the dynamics of these correlations converges to what is termed "self-similar" evolution in cosmology, in which the time dependence in the scaling relation is specified entirely by that of the linearized fluid theory. Comparing simulations with different delta, different resolution, but identical large scale fluctuations, we are able to identify and study features of the dynamics of the system in the transient phase leading to this behavior. In this phase, the discrete nature of the system explicitly plays an essential role.

Collective coordinate control of density distributions

Real collective density variables C(k) [cf. Eq. 1 3 ] in many-particle systems arise from nonlinear transformations of particle positions, and determine the structure factor S(k) , where k denotes the wave vector. Our objective is to prescribe C(k) and then to find many-particle configurations that correspond to such a target C(k) using a numerical optimization technique. Numerical results reported here extend earlier one- and two-dimensional studies to include three dimensions. In addition, they demonstrate the capacity to control S(k) in the neighborhood of |k|=0. The optimization method employed generates multiparticle configurations for which S(k) proportional, |k|alpha, |k|<or=K, and alpha=1, 2, 4, 6, 8, and 10. The case alpha=1 is relevant for the Harrison-Zeldovich model of the early universe, for superfluid 4He, and for jammed amorphous sphere packings. The analysis also provides specific examples of interaction potentials whose classical ground states are configurationally degenerate and disordered.

Force distribution in a randomly perturbed lattice of identical particles with 1/r2 pair interaction

We study the statistics of the force felt by a particle in the class of a spatially correlated distribution of identical pointlike particles, interacting via a 1/r2 pair force (i.e., gravitational or Coulomb), and obtained by randomly perturbing an infinite perfect lattice. We specify the conditions under which the force on a particle is a well-defined stochastic quantity. We then study the small displacements approximation, giving both the limitations of its validity and, when it is valid, an expression for the force variance. The method introduced by Chandrasekhar to find the force probability density function for the homogeneous Poisson particle distribution is extended to shuffled lattices of particles. In this way, we can derive an approximate expression for the probability distribution of the force over the full range of perturbations of the lattice, i.e., from very small (compared to the lattice spacing) to very large where the Poisson limit is recovered. We show in particular the qualitative change in the large-force tail of the force distribution between these two limits. Excellent accuracy of our analytic results is found on detailed comparison with results from numerical simulations. These results provide basic statistical information about the fluctuations of the interactions (i) of the masses in self-gravitating systems like those encountered in the context of cosmological N -body simulations, and (ii) of the charges in the ordered phase of the one-component plasma, the so-called Coulomb or Wigner crystal.

Scale invariant forces in one-dimensional shuffled lattices

We present a detailed and exact study of the probability density function P(F) of the total force F acting on a point particle belonging to a perturbed lattice of identical point sources of a power-law pair interaction. The main results concern the large-F tail of P(F) for which two cases are mainly distinguished: (i) Gaussian-like fast decreasing P(F) for lattice with perturbations forbidding any pair of particles to be found arbitrarily close to one each other and (ii) Lévy-like power-law decreasing P(F) when this possibility is instead permitted. It is important to note that in the second case the exponent of the power-law tail of P(F) is the same for all perturbations (apart from very singular cases) and is in a one-to-one correspondence with the exponent characterizing the behavior of the pair interaction with the distance between the two particles.

Gravitational evolution of a perturbed lattice and its fluid limit

We apply a simple linearization, well known in solid state physics, to approximate the evolution at early times of cosmological N-body simulations of gravity. In the limit that the initial perturbations, applied to an infinite perfect lattice, are at wavelengths much greater than the lattice spacing l, the evolution is exactly that of a pressureless self-gravitating fluid treated in the analogous (Lagrangian) linearization, with the Zeldovich approximation as a subclass of asymptotic solutions. Our less restricted approximation allows one to trace the evolution of the discrete distribution until the time when particles approach one another (i.e., "shell crossing"). We calculate modifications of the fluid evolution, explicitly dependent on l, i.e., discreteness effects in the N-body simulations. We note that these effects become increasingly important as the initial redshift is increased at fixed l.