Athermal jamming of soft frictionless Platonic solids

Random close packing of semi-flexible polymers in two dimensions: Emergence of local and global order

Through extensive Monte Carlo simulations, we systematically study the effect of chain stiffness on the packing ability of linear polymers composed of hard spheres in extremely confined monolayers, corresponding effectively to 2D films. First, we explore the limit of random close packing as a function of the equilibrium bending angle and then quantify the local and global order by the degree of crystallinity and the nematic or tetratic orientational order parameter, respectively. A multi-scale wealth of structural behavior is observed, which is inherently absent in the case of athermal individual monomers and is surprisingly richer than its 3D counterpart under bulk conditions. As a general trend, an isotropic to nematic transition is observed at sufficiently high surface coverages, which is followed by the establishment of the tetratic state, which in turn marks the onset of the random close packing. For chains with right-angle bonds, the incompatibility of the imposed bending angle with the neighbor geometry of the triangular crystal leads to a singular intra- and inter-polymer tiling pattern made of squares and triangles with optimal local filling at high surface concentrations. The present study could serve as a first step toward the design of hard colloidal polymers with a tunable structural behavior for 2D applications.

Three-dimensional compaction of soft granular packings

This paper analyzes the compaction behavior of assemblies composed of soft (elastic) spherical particles beyond the jammed state, using three-dimensional non-smooth contact dynamic simulations. The assemblies of particles are characterized using the evolution of the packing fraction, the coordination number, and the von Misses stress distribution within the particles as the confining stress increases. The packing fraction increases and tends toward a maximum value close to 1, and the mean coordination number increases as a square root of the packing fraction. As the confining stress increases, a transition is observed from a granular-like material with exponential tails of the shear stress distributions to a continuous-like material characterized by Gaussian-like distributions of the shear stresses. We develop an equation that describes the evolution of the packing fraction as a function of the applied pressure. This equation, based on the micromechanical expression of the granular stress tensor, the limit of the Hertz contact law for small deformation, and the power-law relation between the packing fraction and the coordination of the particles, provides good predictions from the jamming point up to very high densities without the need for tuning any parameters.

Structural characterization and statistical properties of jammed soft ellipsoid packing

The jamming transition and jammed packing structures of hydrogel soft ellipsoids are studied using magnetic resonance imaging techniques. As the packing fraction increases, the fluctuation of local free volume decreases and the fluctuation of particle deformation increases. Effective thermodynamic quantities are obtained by characterizing these fluctuations using k-gamma distributions based on an underlying statistical model for granular materials. Surprisingly, the two granular temperatures measuring the relative fluctuations of both free volume and particle deformation remain basically unchanged as the packing fraction increases. The total configurational entropy is also approximately constant for packing with different packing fractions. The significantly different behaviors of these effective thermodynamic quantities compared with hard sphere systems are further attributed to a statistically affine structural transformation of the packing structures along with particle deformations when the packing fraction changes.

Predicting maximally random jammed packing density of non-spherical hard particles via analytical continuation of fluid equation of state

We developed a formalism for accurately predicting the density of MRJ packing state of a wide spectrum of congruent non-spherical hard particles in 3D via analytical fluid EOS.

Jammed packings of 3D superellipsoids with tunable packing fraction, coordination number, and ordering

We carry out numerical studies of static packings of frictionless superellipsoidal particles in three spatial dimensions. We consider more than 200 different particle shapes by varying the three shape parameters that define superellipsoids. We characterize the structural and mechanical properties of both disordered and ordered packings using two packing-generation protocols. We perform athermal quasi-static compression simulations starting from either random, dilute configurations (Protocol 1) or thermalized, dense configurations (Protocol 2), which allows us to tune the orientational order of the packings. In general, we find that superellipsoid packings are hypostatic, with coordination number zJ < ziso, where ziso = 2df and df = 5 or 6 depending on whether the particles are axi-symmetric or not. Over the full range of orientational order, we find that the number of quartic modes of the dynamical matrix for the packings always matches the number of missing contacts relative to the isostatic value. This result suggests that there are no mechanically redundant contacts for ordered, yet hypostatic packings of superellipsoidal particles. Additionally, we find that the packing fraction at jamming onset for disordered packings of superellipsoidal depends on at least two particle shape parameters, e.g. the asphericity A and reduced aspect ratio β of the particles.

Inertial shear flow of assemblies of frictionless polygons: Rheology and microstructure

Motivated by the understanding of shape effects in granular materials, we numerically investigate the macroscopic and microstructural properties of anisotropic dense assemblies of frictionless polydisperse rigid pentagons in shear flow, and compare them with similar systems of disks. Once subjected to large cumulative shear strains their rheology and microstructure are investigated in uniform steady states, depending on inertial number I, which ranges from the quasistatic limit ([Formula: see text]) to 0.2. In the quasistatic limit both systems are devoid of Reynolds dilatancy, i.e., flow at their random close packing density. Both macroscopic friction angle [Formula: see text], an increasing function of I , and solid fraction [Formula: see text], a decreasing function of I, are larger with pentagons than with disks at small I, but the differences decline for larger I and, remarkably, nearly vanish for [Formula: see text]. Under growing I , the depletion of contact networks is considerably slower with pentagons, in which increasingly anisotropic, but still well-connected force-transmitting structures are maintained throughout the studied range. Whereas contact anisotropy and force anisotropy contribute nearly equally to the shear strength in disk assemblies, the latter effect dominates with pentagons at small I, while the former takes over for I of the order of 10^-2. The size of clusters of grains in side-to-side contact, typically comprising more than 10 pentagons in the quasistatic limit, very gradually decreases for growing I.

Geometrical properties of rigid frictionless granular packings as a function of particle size and shape

Three-dimensional discrete numerical simulation is used to investigate the properties of close-packed frictionless granular assemblies as a function of particle polydispersity and shape. Unlike some experimental results, simulations show that disordered packings of pinacoids (eight-face convex polyhedra) achieve higher solid fraction values than amorphous packings of spherical or rounded particles, thus fulfilling the analog of Ulam's conjecture stated by Jiao and co-workers for random packings [Y. Jiao and S. Torquato, Phys. Rev. E 84, 041309 (2011)PLEEE81539-375510.1103/PhysRevE.84.041309]. This seeming discrepancy between experimental and numerical results is believed to result from difficulties in overcoming inter particle friction through experimental densification processes. Moreover, solid fraction is shown to increase further with bidispersity and peak when the volume proportion of small particles reaches 30%. Contrarily, substituting up to 50% of flat pinacoids for isometric ones yields solid fraction decrease, especially when flat particles are also elongated. Nevertheless, particle shape seems to play a minor role in packing solid fraction compared to polydispersity. Additional investigations focused on the packing microstructure confirm that pinacoid packings fulfill the isostatic conjecture and that they are free of order except beyond 30% to 50% of flat or flat-elongated polyhedra in the packing. This order increase progressively takes the form of a nematic phase caused by the reorientation of flat or flat-elongated particles to minimize the packing potential energy. Simultaneously, this reorientation seems to increase the solid fraction value slightly above the maximum achieved by monodisperse isometric pinacoids, as well as the coordination number. Finally, partial substitution of elongated pinacoids for isometric ones has limited effect on packing solid fraction or order.

Maximally dense random packings of cubes and cuboids via a novel inverse packing method

The packings of cubes and cuboids (i.e., "elongated" or "compressed" cubes) are ubiquitous in nature. The high symmetry and space-tiling nature of such particles make them easily packable in dense configurations with a high degree of orientational and translational order. In this paper, we devise a novel inverse packing method that enables one to generate dense hard-particle packings with a controllable degree of disorder/order quantified by predefined order metrics via stochastic Monte Carlo optimizations. We employ the inverse packing method to generate and investigate the maximally dense random packings (MDRPs) of hard cubes and cuboids with aspect ratio α, in which a series of newly introduced normalized local cubatic order parameters sensitive to the onset of any spatial order in packings of cubes and cuboids is minimized. The density of the MDRP of cubes is φ ≈ 0.637, which increases as the shape deviates from the cube limit (α = 1) and reaches the maximal values for cuboids with aspect ratios α = 0.7 or 1.5. These special α values associated with local density extrema are almost identical for those associated with the random packings of spherocylinders, spheroids and superellipsoids, suggesting a universal influence of shape elongation on random packing density. Our inverse packing method can be readily utilized to study the MDRPs of other hard particles and the normalized local cubatic order parameter introduced here is applicable to other shaped particles characterized by three principal axes.

Evolution of the dense packings of spherotetrahedral particles: from ideal tetrahedra to spheres

Particle shape plays a crucial role in determining packing characteristics. Real particles in nature usually have rounded corners. In this work, we systematically investigate the rounded corner effect on the dense packings of spherotetrahedral particles. The evolution of dense packing structure as the particle shape continuously deforms from a regular tetrahedron to a sphere is investigated, starting both from the regular tetrahedron and the sphere packings. The dimer crystal and the quasicrystal approximant are used as initial configurations, as well as the two densest sphere packing structures. We characterize the evolution of spherotetrahedron packings from the ideal tetrahedron (s = 0) to the sphere (s = 1) via a single roundness parameter s. The evolution can be partitioned into seven regions according to the shape variation of the packing unit cell. Interestingly, a peak of the packing density Φ is first observed at s ≈ 0.16 in the Φ-s curves where the tetrahedra have small rounded corners. The maximum density of the deformed quasicrystal approximant family (Φ ≈ 0.8763) is slightly larger than that of the deformed dimer crystal family (Φ ≈ 0.8704), and both of them exceed the densest known packing of ideal tetrahedra (Φ ≈ 0.8563).

Cluster and constraint analysis in tetrahedron packings

The disordered packings of tetrahedra often show no obvious macroscopic orientational or positional order for a wide range of packing densities, and it has been found that the local order in particle clusters is the main order form of tetrahedron packings. Therefore, a cluster analysis is carried out to investigate the local structures and properties of tetrahedron packings in this work. We obtain a cluster distribution of differently sized clusters, and peaks are observed at two special clusters, i.e., dimer and wagon wheel. We then calculate the amounts of dimers and wagon wheels, which are observed to have linear or approximate linear correlations with packing density. Following our previous work, the amount of particles participating in dimers is used as an order metric to evaluate the order degree of the hierarchical packing structure of tetrahedra, and an order map is consequently depicted. Furthermore, a constraint analysis is performed to determine the isostatic or hyperstatic region in the order map. We employ a Monte Carlo algorithm to test jamming and then suggest a new maximally random jammed packing of hard tetrahedra from the order map with a packing density of 0.6337.

Internal friction and absence of dilatancy of packings of frictionless polygons

By means of numerical simulations, we show that assemblies of frictionless rigid pentagons in slow shear flow possess an internal friction coefficient (equal to 0.183±0.008 with our choice of moderately polydisperse grains) but no macroscopic dilatancy. In other words, despite side-side contacts tending to hinder relative particle rotations, the solid fraction under quasistatic shear coincides with that of isotropic random close packings of pentagonal particles. Properties of polygonal grains are thus similar to those of disks in that respect. We argue that continuous reshuffling of the force-bearing network leads to frequent collapsing events at the microscale, thereby causing the macroscopic dilatancy to vanish. Despite such rearrangements, the shear flow favors an anisotropic structure that is at the origin of the ability of the system to sustain shear stress.

Effect of size polydispersity versus particle shape in dense granular media

We present a detailed analysis of the morphology of granular systems composed of frictionless pentagonal particles by varying systematically both the size span and particle shape irregularity, which represent two polydispersity parameters of the system. The microstructure is characterized in terms of various statistical descriptors such as global and local packing fractions, radial distribution functions, coordination number, and fraction of floating particles. We find that the packing fraction increases with the two parameters of polydispersity, but the effect of shape polydispersity for all the investigated structural properties is significant only at low size polydispersity where the positional and/or orientational ordering of the particles prevail. We focus in more detail on the class of side/side contacts, which is the interesting feature of our system as compared to a packing of disks. We show that the proportion of such contacts has weak dependence on the polydispersity parameters. The side- side contacts do not percolate but they define clusters of increasing size as a function of size polydispersity and decreasing size as a function of shape polydispersity. The clusters have anisotropic shapes but with a decreasing aspect ratio as polydispersity increases. This feature is argued to be a consequence of strong force chains (forces above the mean), which are mainly captured by side-side contacts. Finally, the force transmission is intrinsically multiscale, with a mean force increasing linearly with particle size.

Variable-cell method for stress-controlled jamming of athermal, frictionless grains

A method is introduced to simulate jamming of polyhedral grains under controlled stress that incorporates global degrees of freedom through the metric tensor of a periodic cell containing grains. Jamming under hydrostatic (isotropic) stress and athermal conditions leads to a precise definition of the ideal jamming point at zero shear stress. The structures of tetrahedra jammed hydrostatically exhibit less translational order and lower jamming-point density than previously described maximally random jammed hard tetrahedra. Under the same conditions, cubes jam with negligible nematic order. Grains with octahedral symmetry having s>0.5 (where s interpolates from octahedra [s=0] to cubes [s=1]) jam with an abundance of face-face contacts in the absence of nematic order. For sufficiently large face-face contact number, percolating clusters form that span the entire simulation box. The response of hydrostatically jammed tetrahedra and cubes to shear-stress perturbation is also demonstrated with the variable-cell method.

Jammed frictional tetrahedra are hyperstatic

We prepare packings of frictional tetrahedra with volume fractions ϕ ranging from 0.469 to 0.622 using three different experimental protocols under isobaric conditions. Analysis via x-ray microtomography reveals that the contact number Z grows with ϕ, but does depend on the preparation protocol. While there exist four different types of contacts in tetrahedra packings, our analysis shows that the edge-to-face contacts contribute about 50% of the total increase in Z. The number of constraints per particle C increases also with ϕ and even the loosest packings are strongly hyperstatic, i.e., mechanically overdetermined with C approximately twice the degrees of freedom each particle possesses.

Packings of irregular polyhedral particles: strength, structure, and effects of angularity

We present a systematic numerical investigation of the shear strength and structure of granular packings composed of irregular polyhedral particles. The angularity of the particles is varied by increasing the number of faces from 8 (octahedronlike shape) to 596. We find that the shear strength increases with angularity up to a maximum value and saturates as the particles become more angular (below 46 faces). At the same time, the packing fraction increases to a peak value but declines for more angular particles. We analyze the connectivity and anisotropy of the microstructure by considering both the contacts and branch vectors joining particle centers. The increase of the shear strength with angularity is shown to be due to a net increase of the fabric and force anisotropies but at higher particle angularity a rapid falloff of the fabric anisotropy is compensated by an increase of force anisotropy, leading thus to the saturation of shear strength.

Spontaneous crystallization in athermal polymer packings

We review recent results from extensive simulations of the crystallization of athermal polymer packings. It is shown that above a certain packing density, and for sufficiently long simulations, all random assemblies of freely-jointed chains of tangent hard spheres of uniform size show a spontaneous transition into a crystalline phase. These polymer crystals adopt predominantly random hexagonal close packed morphologies. An analysis of the local environment around monomers based on the shape and size of the Voronoi polyhedra clearly shows that Voronoi cells become more spherical and more symmetric as the system transits to the ordered state. The change in the local environment leads to an increase in the monomer translational contribution to the entropy of the system, which acts as the driving force for the phase transition. A comparison of the crystallization of hard-sphere polymers and monomers highlights similarities and differences resulting from the constraints imposed by chain connectivity.

Effect of confinement on dense packings of rigid frictionless spheres and polyhedra

We study numerically the influence of confinement on the solid fraction and on the structure of three-dimensional random close-packed granular materials subject to gravity. The effects of grain shape (spherical or polyhedral), material polydispersity, and confining wall friction on this dependence are investigated. In agreement with a simple geometrical model, the solid fraction is found to decrease linearly for increasing confinement no matter the grain shape. Furthermore, this decrease remains valid for bidisperse sphere packings, although the gradient seems to reduce significantly when the proportion of small particles reaches 40% by volume. The confinement effect on the coordination number is also captured by an extension of the aforementioned model.

Shape effects on the random-packing density of tetrahedral particles

Regular tetrahedra have been demonstrated recently giving high packing density in random configurations. However, it is unknown whether the random-packing density of tetrahedral particles with other shapes can reach an even higher value. A numerical investigation on the random packing of regular and irregular tetrahedral particles is carried out. Shape effects of rounded corner, eccentricity, and height on the packing density of tetrahedral particles are studied. Results show that altering the shape of tetrahedral particles by rounding corners and edges, by altering the height of one vertex, or by lateral displacement of one vertex above its opposite face, all individually have the effect of reducing the random-packing density. In general, the random-packing densities of irregular tetrahedral particles are lower than that of regular tetrahedra. The ideal regular tetrahedron should be the shape which has the highest random-packing density in the family of tetrahedra, or even among convex bodies. An empirical formula is proposed to describe the rounded corner effect on the packing density, and well explains the density deviation of tetrahedral particles with different roundness ratios. The particles in the simulations are verified to be randomly packed by studying the pair correlation functions, which are consistent with previous results. The spherotetrahedral particle model with the relaxation algorithm is effectively applied in the simulations.

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. All of our organizing principles are applied to and tested against the most comprehensive set of both convex and concave particle shapes examined 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.

Rheology of granular materials composed of nonconvex particles

By means of contact dynamics simulations, we investigate the shear strength and internal structure of granular materials composed of two-dimensional nonconvex aggregates. We find that the packing fraction first grows as the nonconvexity is increased but declines at higher nonconvexity. This unmonotonic dependence reflects the competing effects of pore size reduction between convex borders of aggregates and gain in porosity at the nonconvex borders that are captured in a simple model fitting nicely the simulation data both in the isotropic and sheared packings. On the other hand, the internal angle of friction increases linearly with nonconvexity and saturates to a value independent of nonconvexity. We show that fabric anisotropy, force anisotropy, and friction mobilization, all enhanced by multiple contacts between aggregates, govern the observed increase of shear strength and its saturation with increasing nonconvexity. The main effect of interlocking is to dislocate frictional dissipation from the locked double and triple contacts between aggregates to the simple contacts between clusters of aggregates. This self-organization of particle motions allows the packing to keep a constant shear strength at high nonconvexity.

Maximally random jammed packings of Platonic solids: hyperuniform long-range correlations and isostaticity

We generate maximally random jammed (MRJ) packings of the four nontiling Platonic solids (tetrahedra, octahedra, dodecahedra, and icosahedra) using the adaptive-shrinking-cell method [S. Torquato and Y. Jiao, Phys. Rev. E 80, 041104 (2009)]. Such packings can be viewed as prototypical glasses in that they are maximally disordered while simultaneously being mechanically rigid. The MRJ packing fractions for tetrahedra, octahedra, dodecahedra, and icosahedra are, respectively, 0.763±0.005, 0.697±0.005, 0.716±0.002, and 0.707±0.002. We find that as the number of facets of the particles increases, the translational order in the packings increases while the orientational order decreases. Moreover, we show that the MRJ packings are hyperuniform (i.e., their infinite-wavelength local-number-density fluctuations vanish) and possess quasi-long-range pair correlations that decay asymptotically with scaling r(-4). This provides further evidence that hyperuniform quasi-long-range correlations are a universal feature of MRJ packings of frictionless particles of general shape. However, unlike MRJ packings of ellipsoids, superballs, and superellipsoids, which are hypostatic, MRJ packings of the nontiling Platonic solids are isostatic. We provide a rationale for the organizing principle that the MRJ packing fractions for nonspherical particles with sufficiently small asphericities exceed the corresponding value for spheres (∼0.64). We also discuss how the shape and symmetry of a polyhedron particle affects its MRJ packing fraction.

Isostaticity of constraints in amorphous jammed systems of soft frictionless Platonic solids

The average number of constraints per particle <C(total)> in mechanically stable amorphous systems of Platonic solids approaches the isostatic limit at the jamming point (<C(total)>→12), though average number of contacts are hypostatic. By introducing angular alignment metrics to classify the degree of constraint imposed by each contact, constraints are shown to arise as a direct result of local orientational order reflected in edge-face and face-face alignment angle distributions. With approximately one face-face contact per particle at jamming, chainlike face-face clusters form with finite extent--a signature of amorphous jammed systems.