Arxiv:2103.14989V2 [Cond-Mat.Soft]
Total Page:16
File Type:pdf, Size:1020Kb
Structural Characterization of Many-Particle Systems on Approach to Hyperuniform States Salvatore Torquato∗ Department of Chemistry, Department of Physics, Princeton Institute for the Science and Technology of Materials, and Program in Applied and Computational Mathematics, Princeton University, Princeton, New Jersey 08544, USA (Dated: April 1, 2021) 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 Rd 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 infor- mation 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 hyperunifor- mity, 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 sys- tems 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. I. INTRODUCTION An equivalent definition of hyperuniformity is based on the local number variance σ2(R) N(R)2 N(R) 2 ≡ h i − h i A hyperuniform point configuration in d-dimensional associated with the number N(R) of points within a d- Euclidean space Rd is characterized by an anomalous dimensional spherical observation window of radius R, suppression of large-scale density fluctuations relative where angular brackets denote an ensemble average. A Rd to those in typical disordered systems, such as liquids point pattern in is hyperuniform if its variance grows d and structural glasses [1, 2]. More precisely, a hyperuni- in the large-R limit slower than R . This behavior is to form point pattern is one in which the structure factor be contrasted with those of typical disordered systems, S(k) 1+ρh˜(k) tends to zero as the wavenumber k k such as Poisson point patterns, gases and liquids, where arXiv:2103.14989v2 [cond-mat.soft] 30 Mar 2021 tends≡ to zero [1, 2], i.e., ≡ | | the number variance scales like the volume v1(R) of the observation window, which is given by lim S(k)=0, (1) |k|→0 πd/2Rd ˜ k v1(R)= . (2) where h( ) is the Fourier transform of the total corre- Γ(1 + d/2) lation function h(r) g (r) 1 and g (r) is the pair- ≡ 2 − 2 correlation function [3]. The hyperuniformity concept Consider systems that are characterized by a structure generalizes the traditional notion of long-range order in factor with a radial power-law form in the vicinity of the many-particle systems to include all perfect crystals, per- origin, i.e., fect quasicrystals, and exotic amorphous states of matter. Disordered hyperuniform materials can have advantages S(k) k α for k 0. (3) over crystalline ones, such as unique or nearly optimal, ∼ | | | | → direction-independent physical properties and robustness For hyperuniform systems, the exponent α is positive against defects [4–16]. (α > 0) and its value determines three different large- 2 R scaling behaviors of the number variance [1, 2, 17]: phenomenon, a hyperuniformity length scale ξc and a hy- peruniformity index H. The latter two quantities provide Rd−1, α> 1 (class I) measures of nearness to hyperuniformity. σ2(R) Rd−1 ln R, α = 1 (class II) . (4) In the remainder of the paper, we obtain a variety of ∼ Rd−α, α< 1 (class III) theoretical results to study the problem at hand. First, 2 for a general system, we establish quantitative criteria These scalings of σ (R) define three classes of hype- to ascertain the extent of hyperuniform and nonhyper- runiformity [2], with classes I and III describing the uniform distance-scaling regimes as well as the crossover strongest and weakest forms of hyperuniformity, respec- point between them in terms of the “volume” coeffi- tively. States of matter that belong to class I include all cient A and “surface-area” coefficient B associated with perfect crystals [1, 17], many perfect quasicrystals [17– the variance σ2(R) (Sec. III). Specifically, the ratio 19], and “randomly” perturbed crystal structures [20– B/A determines the crossover length scale Rc. The 23], classical disordered ground states of matter [1, 24, 25] larger the ratio B/A, the larger the hyperuniform scal- as well as systems out of equilibrium [26, 27]. Class II hy- ing regime, which becomes of infinite extent in the limit peruniform systems include some quasicrystals [19], the B/A . This capacity to determine hyperuniform positions of the prime numbers [28], and many disordered scaling→ regimes ∞ is expected to be particularly useful in classical [26, 29–32] and quantum [33–35] states of mat- analyzing experimentally- or computationally-generated ter. Examples of class III hyperuniform systems include samples that are necessarily of finite size. classical disordered ground states [36], random organiza- tion models [37] and perfect glasses [26]. Second, to complement the known direct-space repre- By contrast, for any nonhyperuniform system, it is sentation of the coefficient B in terms of the total corre- shown in Appendix A that the local variance has the lation function h(r)[1], we derive here its corresponding following large-R scaling behaviors: Fourier representation in terms of the structure factor S(k) (Sec. IV). The latter representation is particularly d + useful when the scattering intensity is available experi- 2 R , α =0 (typical nonhyperuniform) σ (R) d−α mentally or if the structure factor is known analytically. ∼ (R , α< 0 (anti-hyperuniform). (5) Third, we show that the free-volume theory of the pres- sure of equilibrium packings of identical hard spheres that For a “typical” nonhyperuniform system, S(0) is approach either a strictly jammed crystal or disordered bounded [2]. In anti-hyperuniform systems, S(0) is un- state dictates that such jammed states be perfectly hype- bounded, i.e., runiform (Sec. V). We describe why this outcome implies that such jammed states must be defect-free. lim S(k)=+ , (6) |k|→0 ∞ Fourth, motivated by a desire to rely on analytical rather than numerical methods, we structurally charac- and hence are diametrically opposite to hyperuniform terize three different disordered-system models as a func- systems. Anti-hyperuniform systems include fractals, tion of the relevant control parameter, either density or systems at thermal critical points (e.g., liquid-vapor and temperature, with end states that are perfectly hyperuni- magnetic critical points) [38–42] as well as certain sub- form. These models are distinguished from most other stitution tilings [43]. models with hyperuniform states in that their pair corre- Our main concern in this paper is the exploration lation functions and structure factors are known exactly of quantitative descriptors that herald when a many- for all values of the control parameter in the thermody- particle system is nearly hyperuniform or approaching namic limit. We purposely avoid the use of simulations a hyperuniform state, whether ordered or not. Eluci- of many-particle systems in finite boxes to draw con- dating such questions not only is expected to lead to a clusions, since hyperuniformity in an infinite-wavelength deeper fundamental understanding of the nature and for- property. The first model that we characterize is an equi- mation of hyperuniform systems but has great practical librium system of hard rods (Sec. VI). Here the control value. For example, since hyperuniformity can endow parameter is the number density ρ (or packing fraction) a system with novel or optimal physical properties, it and its terminal value corresponds to the jammed state is essential to know how close the system must be to that is the integer lattice. The second model studied is perfect hyperuniformity without significantly degrading a certain “sticky” hard-sphere system in arbitrary space its ideal performance. In practice, perfect hyperunifor- dimension d as a function of the number density (or pack- mity is never achieved due to defects that are inevitably ing fraction) (Sec. VII). The third model that we char- present in any real finite-sized system [44], whether crys- acterize are low-temperature excited states of so-called talline, quasicrystalline or disordered [23].