Arxiv:1804.04927V2 [Cond-Mat.Dis-Nn] 20 Aug 2018
Total Page:16
File Type:pdf, Size:1020Kb
Absence of hyperuniformity in amorphous hard-sphere packings of nonvanishing complexity M. J. Godfrey and M. A. Moore School of Physics and Astronomy, University of Manchester, Manchester M13 9PL, UK (Dated: August 21, 2018) We relate the structure factor S(k ! 0) in a system of jammed hard spheres of number density ρ to its complexity per particle Σ(ρ) by the formula S(k ! 0) = −1=[ρ2Σ00(ρ) + 2ρΣ0(ρ)]. We have verified this formula for the case of jammed disks in a narrow channel, for which it is possible to find Σ(ρ) and S(k) analytically. Hyperuniformity, which is the vanishing of S(k ! 0), will therefore not occur if the complexity is nonzero. An example is given of a jammed state of hard disks in a narrow channel which is hyperuniform when generated by dynamical rules that produce a non-extensive complexity. Systems that are hyperuniform have small long-wavelength in the first place”. In this paper we argue for the correctness density fluctuations. In other words, their structure factor of this view, by showing that at least in the regime where the S(k) vanishes as the wavevector k ! 0. The structure logarithm of the number of jammed states (a quantity usually factor of N spheres whose centers are at rj is defined as referred to as the complexity ΣN;V ) is extensive and propor- 2 S(k) = jρkj =N, where the Fourier transform of the num- tional to N one can derive a simple formula for the finite value PN of S(k ! 0) in terms of derivatives of Σ . Our formula ber density ρk = j=1 exp(−ik · rj). The whole subject N;V of hyperuniformity has been recently reviewed at length by for S(k ! 0) is the analog of that for an equilibrium system Torquato [1]. In this paper we focus on one of the most ex- at temperature T : tensively studied questions connected with hyperuniformity – hM 2i − hMi2 the extent to which jammed systems of hard spheres or disks S(k ! 0) = = ρkBT κT ; (1) are hyperuniform. Torquato and Stillinger [2] have made the hMi following conjecture: where ρ is the number density N=V , and κ is the isother- Any strictly jammed saturated infinite packing of identical T mal compressibility [1]; M is the number of particles within spheres in d is hyperuniform. R a large sub-volume (window) of the system of volume Ω so The restrictions on the packing are important. A “strictly” hMi = ρΩ. The name “hyperuniformity” means that the vari- jammed packing is one such that each particle is firstly “lo- ance of the number of particles in Ω is vanishingly small com- cally” trapped by its neighbors so that it cannot be trans- pared to the natural expectation that it should be proportional lated while their positions are held fixed and secondly is also to M. Hard sphere fluids at a temperature T are not hyperuni- “collectively jammed” so that no subset of the particles can form as their compressibility κ is finite. be simultaneously moved so that its members move out of T We shall first re-write Eq. (1) in a more useful form for our contact with one another with the remainder set fixed. Fi- discussion of jammed states. Let p be the pressure in the sys- nally, the condition of strict jamming precludes any uniform tem. The compressibility κ is related to the second derivative volume-decreasing strains of the boundary, which means that T of the Helmholtz free energy F of the system, F = E − TS, their bulk and shear moduli are infinitely large [1]. The word where E is the energy of the system and S is its entropy: “saturated” in the conjecture eliminates packings that contain voids large enough to accommodate a “rattler” [1]. Such voids 1 @V 1 1 κ ≡ − = : (2) would be a source of diffuse scattering, leading to a nonzero T V @p V (@2F=@V 2) limiting value for the structure factor, S(k ! 0) 6= 0. T;N T;N 1 Many investigators have reported that the limit S(k ! 0) For a system of hard spheres in d dimensions, E = 2 NdkBT , 2 2 2 2 seems to be non-zero for jammed packings of hard spheres so that we can write @ F=@V T;N = −T @ S=@V T;N . in three dimensions; Refs. [3,4] are recent examples. In re- Then Eq. (1) becomes sponse, Torquato has asserted that the “difficulty of ensur- ρ k ing jamming as N becomes sufficiently large to access the S(k ! 0) = − B : (3) V (@2S=@V 2) arXiv:1804.04927v2 [cond-mat.dis-nn] 20 Aug 2018 small-wavenumber regime in the structure factor as well as T;N the presence of rattlers that degrade hyperuniformity makes it We turn now to generalizing Eq. (3) to jammed states. We virtually impossible to test the Torquato-Stillinger jamming- imagine that an averaging process is used. The simplest is the hyperuniformity conjecture via current numerical packing Edwards procedure of giving equal weight to all the jammed protocols” [1]. states of a given number density ρ. Let us call the complex- Ikeda and Berthier [3], who found that S(k ! 0) 6= 0 in ity ΣN;V , the logarithm of the number of jammed states of their work, observed that it is “difficult to provide a physical N spheres in a volume V . Numerical methods have been explanation for the existence of the observed deviations from recently devised for determining ΣN;V [5]. When it is ex- strict hyperuniformity, mainly because there is no deep physi- tensive, it is expected to be of the form ΣN;V = NΣ(ρ) as cal reason to expect perfect hyperuniformity in these systems N ! 1 at fixed density ρ ≡ N=V . It is then our claim that 2 the average of the structure factor S(k ! 0) over the jammed The structure factor is given by the usual expression states can be obtained from Eq. (3) by making the replacement S(k) = hρ~kρ~−ki=N, where ρ~k is the Fourier transform of ρ~(x) = ρ(x) − hρi, the deviation of ρ(x) away from its aver- S=N ! kBΣ(ρ): (4) age value. The subtraction removes from S(k) a δ-function at k = 0 and makes it possible to show that This kind of idea is clearly related to those that have been suc- cessfully used in the theory of granular systems [6]. In terms var M S(k ! 0) = ; (9) of derivatives of the function Σ(ρ), Eqs. (3) and (4) yield hMi 1 where M is the number of particles in a window of volume S(k ! 0) = − : (5) ρ2Σ00(ρ) + 2ρΣ0(ρ) Ω V , placed at random in V , in the limit where Ω and hMi are large [1]. Hyperuniformity can arise if the function Σ(ρ) has divergent To find the distribution of M we adopt an approach that is first or second derivatives, which is only likely to happen at similar to the elementary one used to derive the Gibbs distri- exceptional values of ρ; we shall see later that for a system bution in statistical physics. In a typical state drawn from the of disks in a narrow channel this happens at the two packing Edwards ensemble, the number of spheres M whose centers fractions, φmin and φmax, for which the system has crystalline lie within a randomly-placed window of volume Ω follows a order. In general, the right hand side of Eq. (5) will not be certain distribution P (M) which we suppose to be indepen- zero and there will be no hyperuniformity. However, if the dent of the particular state chosen, in the limit V ! 1. We complexity is not extensive, Eq. (5) cannot be used and hy- expect P (M) to be unchanged by further averaging over the peruniformity may then occur. Later we shall give an explicit Edwards ensemble, and to be identical to the probability that example of such behavior for a system of hard disks in a nar- a fixed window of volume Ω contains M spheres. Provided row channel for which the configuration has been generated the window is much larger than any correlation length for the by a set of dynamical rules that result in a non-extensive com- system, we can use the extensivity of the complexity to ap- plexity. proximate the probabilities by The result in Eq. (5), obtained by analogy with the equi- Σ +Σ librium theory, can be derived directly for jammed systems. P (M) / e N−M;V −Ω M;Ω ; (10) Denote by V [C ] the volume of a jammed configuration C N N where the terms in the exponent are given by expressions sim- of N spheres. Then the complexity ΣN;V is given by ilar to (6). Expanding ΣN−M;V −Ω to first order in M then ΣN;V X gives e = δ(V − V [CN ]); (6) Σ +µM CN P (M) / e M;Ω : (11) where we follow the Edwards prescription [6] of giving equal Thus we are led to consider the partition function weight to each of the jammed states. (We discuss other aver- X ΣM;Ω+µM aging possibilities later.) Often it is easier to work with the ZM;Ω = e : (12) Legendre transform of Eq. (6), and discuss M Z The steepest-descents condition for this sum recovers the ex- X −fV [CN ] ΣN;V −fV ZN;f = e = e dV: (7) pected µ = −@ΣM;Ω=@M, so that by writing ΣM;Ω = CN MΣ(ρ), with ρ = M=Ω, we obtain The reciprocal of f is what Edwards and Oakeshott have re- µ = −d[ρΣ(ρ)]/dρ.