Pair correlation function characteristics of nearly jammed disordered and ordered hard-sphere packings

Aleksandar Donev, Salvatore Torquato, Frank H. Stillinger

Research output: Contribution to journalArticle

Abstract

We study the approach to jamming in hard-sphere packings and, in particular, the pair correlation function g 2(r) around contact, both theoretically and computationally. Our computational data unambiguously separate the narrowing δ-function contribution to g 2 due to emerging interparticle contacts from the background contribution due to near contacts. The data also show with unprecedented accuracy that disordered hard-sphere packings are strictly isostatic: i.e., the number of exact contacts in the jamming limit is exactly equal to the number of degrees of freedom, once rattlers are removed. For such isostatic packings, we derive a theoretical connection between the probability distribution of interparticle forces P f(f), which we measure computationally, and the contact contribution to g 2. We verify this relation for computationally generated isostatic packings that are representative of the maximally random jammed state. We clearly observe a maximum in P f and a nonzero probability of zero force, shedding light on long-standing questions in the granular-media literature. We computationally observe an unusual power-law divergence in the near-contact contribution to g 2, persistent even in the jamming limit, with exponent -0.4 clearly distinguishable from previously proposed inverse-square-root divergence. Additionally, we present high-quality numerical data on the two discontinuities in the split-second peak of g 2 and use a shared-neighbor analysis of the graph representing the contact network to study the local particle clusters responsible for the peculiar features. Finally, we present the computational data on the contact contribution to g 2 for vacancy-diluted fcc crystal packings and also investigate partially crystallized packings along the transition from maximally disordered to fully ordered packings. We find that the contact network remains isostatic even when ordering is present. Unlike previous studies, we find that ordering has a significant impact on the shape of P f for small forces.

Original languageEnglish (US)
Article number011105
JournalPhysical Review E
Volume71
Issue number1
DOIs
StatePublished - Jan 2005

Fingerprint

Sphere packing
Pair Correlation Function
Hard Spheres
Contact
Packing
jamming
Jamming
divergence
Divergence
Granular Media
Vacancy
emerging
discontinuity
degrees of freedom
Square root
exponents
Discontinuity
Power Law
Crystal
Probability Distribution

ASJC Scopus subject areas

  • Physics and Astronomy(all)
  • Condensed Matter Physics
  • Statistical and Nonlinear Physics
  • Mathematical Physics

Cite this

Pair correlation function characteristics of nearly jammed disordered and ordered hard-sphere packings. / Donev, Aleksandar; Torquato, Salvatore; Stillinger, Frank H.

In: Physical Review E, Vol. 71, No. 1, 011105, 01.2005.

Research output: Contribution to journalArticle

@article{daded2995b4947e5a86f96e420c8a12e,
title = "Pair correlation function characteristics of nearly jammed disordered and ordered hard-sphere packings",
abstract = "We study the approach to jamming in hard-sphere packings and, in particular, the pair correlation function g 2(r) around contact, both theoretically and computationally. Our computational data unambiguously separate the narrowing δ-function contribution to g 2 due to emerging interparticle contacts from the background contribution due to near contacts. The data also show with unprecedented accuracy that disordered hard-sphere packings are strictly isostatic: i.e., the number of exact contacts in the jamming limit is exactly equal to the number of degrees of freedom, once rattlers are removed. For such isostatic packings, we derive a theoretical connection between the probability distribution of interparticle forces P f(f), which we measure computationally, and the contact contribution to g 2. We verify this relation for computationally generated isostatic packings that are representative of the maximally random jammed state. We clearly observe a maximum in P f and a nonzero probability of zero force, shedding light on long-standing questions in the granular-media literature. We computationally observe an unusual power-law divergence in the near-contact contribution to g 2, persistent even in the jamming limit, with exponent -0.4 clearly distinguishable from previously proposed inverse-square-root divergence. Additionally, we present high-quality numerical data on the two discontinuities in the split-second peak of g 2 and use a shared-neighbor analysis of the graph representing the contact network to study the local particle clusters responsible for the peculiar features. Finally, we present the computational data on the contact contribution to g 2 for vacancy-diluted fcc crystal packings and also investigate partially crystallized packings along the transition from maximally disordered to fully ordered packings. We find that the contact network remains isostatic even when ordering is present. Unlike previous studies, we find that ordering has a significant impact on the shape of P f for small forces.",
author = "Aleksandar Donev and Salvatore Torquato and Stillinger, {Frank H.}",
year = "2005",
month = "1",
doi = "10.1103/PhysRevE.71.011105",
language = "English (US)",
volume = "71",
journal = "Physical Review E - Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics",
issn = "1063-651X",
publisher = "American Physical Society",
number = "1",

}

TY - JOUR

T1 - Pair correlation function characteristics of nearly jammed disordered and ordered hard-sphere packings

AU - Donev, Aleksandar

AU - Torquato, Salvatore

AU - Stillinger, Frank H.

PY - 2005/1

Y1 - 2005/1

N2 - We study the approach to jamming in hard-sphere packings and, in particular, the pair correlation function g 2(r) around contact, both theoretically and computationally. Our computational data unambiguously separate the narrowing δ-function contribution to g 2 due to emerging interparticle contacts from the background contribution due to near contacts. The data also show with unprecedented accuracy that disordered hard-sphere packings are strictly isostatic: i.e., the number of exact contacts in the jamming limit is exactly equal to the number of degrees of freedom, once rattlers are removed. For such isostatic packings, we derive a theoretical connection between the probability distribution of interparticle forces P f(f), which we measure computationally, and the contact contribution to g 2. We verify this relation for computationally generated isostatic packings that are representative of the maximally random jammed state. We clearly observe a maximum in P f and a nonzero probability of zero force, shedding light on long-standing questions in the granular-media literature. We computationally observe an unusual power-law divergence in the near-contact contribution to g 2, persistent even in the jamming limit, with exponent -0.4 clearly distinguishable from previously proposed inverse-square-root divergence. Additionally, we present high-quality numerical data on the two discontinuities in the split-second peak of g 2 and use a shared-neighbor analysis of the graph representing the contact network to study the local particle clusters responsible for the peculiar features. Finally, we present the computational data on the contact contribution to g 2 for vacancy-diluted fcc crystal packings and also investigate partially crystallized packings along the transition from maximally disordered to fully ordered packings. We find that the contact network remains isostatic even when ordering is present. Unlike previous studies, we find that ordering has a significant impact on the shape of P f for small forces.

AB - We study the approach to jamming in hard-sphere packings and, in particular, the pair correlation function g 2(r) around contact, both theoretically and computationally. Our computational data unambiguously separate the narrowing δ-function contribution to g 2 due to emerging interparticle contacts from the background contribution due to near contacts. The data also show with unprecedented accuracy that disordered hard-sphere packings are strictly isostatic: i.e., the number of exact contacts in the jamming limit is exactly equal to the number of degrees of freedom, once rattlers are removed. For such isostatic packings, we derive a theoretical connection between the probability distribution of interparticle forces P f(f), which we measure computationally, and the contact contribution to g 2. We verify this relation for computationally generated isostatic packings that are representative of the maximally random jammed state. We clearly observe a maximum in P f and a nonzero probability of zero force, shedding light on long-standing questions in the granular-media literature. We computationally observe an unusual power-law divergence in the near-contact contribution to g 2, persistent even in the jamming limit, with exponent -0.4 clearly distinguishable from previously proposed inverse-square-root divergence. Additionally, we present high-quality numerical data on the two discontinuities in the split-second peak of g 2 and use a shared-neighbor analysis of the graph representing the contact network to study the local particle clusters responsible for the peculiar features. Finally, we present the computational data on the contact contribution to g 2 for vacancy-diluted fcc crystal packings and also investigate partially crystallized packings along the transition from maximally disordered to fully ordered packings. We find that the contact network remains isostatic even when ordering is present. Unlike previous studies, we find that ordering has a significant impact on the shape of P f for small forces.

UR - http://www.scopus.com/inward/record.url?scp=37649027929&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=37649027929&partnerID=8YFLogxK

U2 - 10.1103/PhysRevE.71.011105

DO - 10.1103/PhysRevE.71.011105

M3 - Article

VL - 71

JO - Physical Review E - Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics

JF - Physical Review E - Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics

SN - 1063-651X

IS - 1

M1 - 011105

ER -