Analogies between granular jamming and the liquid-glass transition

Environmentally Sensitive Luminescence Reveals Spatial Confinement, Dynamics, and Their Molecular Weight Dependence in a Polymer Glass

Polymer glasses have an irregular structure. Among the causes for such complexity are the chemically distinct chain end groups that are the most abundant irregularities in any linear polymer. In this work, we demonstrate that chain end induced defects allow polymer glasses to create confined environments capable of hosting small emissive molecules. Using environmentally sensitive luminescent complexes, we show that the size of these confinements depends on molecular weight and can dramatically affect the photoluminescence of free or covalently bound emissive complexes. We confirm the impact of chain end confinement on the bulk glass transition in poly(methyl acrylate) (pMA) and show that commonly observed T g changes induced by the chain ends should have a structural origin. Finally, we demonstrate that the size and placement of luminescent molecular probes in pMA can dramatically affect the probe luminescence and its temperature dependence, suggesting that polymer glass is a highly irregular and complex environment, marking its difference with conventional small molecule solvents. Considering the ubiquity of luminescent glassy materials, our work lays down a blueprint for designing them with structural considerations in mind, ones where packing density and chain end size are key factors.

Micromechanical description of the compaction of soft pentagon assemblies

We analyze the isotropic compaction of assemblies composed of soft pentagons interacting through classical Coulomb friction via numerical simulations. The effect of the initial particle shape is discussed by comparing packings of pentagons with packings of soft circular particles. We characterize the evolution of the packing fraction, the elastic modulus, and the microstructure (particle rearrangement, connectivity, contact force, and particle stress distributions) as a function of the applied stresses. Both systems behave similarly: the packing fraction increases and tends asymptotically to a maximum value ϕ_{max}, where the bulk modulus diverges. At the microscopic scale we show that particle rearrangements occur even beyond the jammed state, the mean coordination increases as a square root of the packing fraction, and the force and stress distributions become more homogeneous as the packing fraction increases. Soft pentagons experience larger particle rearrangements than circular particles, and such behavior decreases proportionally to the friction. Interestingly, the friction between particles also contributes to a better homogenization of the contact force network in both systems. From the expression of the granular stress tensor we develop a model that describes the compaction behavior as a function of the applied pressure, the Young modulus, and the initial shape of the particles. This model, settled on the joint evolution of the particle connectivity and the contact stress, provides outstanding predictions from the jamming point up to very high densities.

Biomechanics of Collective Cell Migration in Cancer Progression: Experimental and Computational Methods

Cell migration is essential for regulating many biological processes in physiological or pathological conditions, including embryonic development and cancer invasion. In vitro and in silico studies suggest that collective cell migration is associated with some biomechanical particularities such as restructuring of extracellular matrix (ECM), stress and force distribution profiles, and reorganization of the cytoskeleton. Therefore, the phenomenon could be understood by an in-depth study of cells' behavior determinants, including but not limited to mechanical cues from the environment and from fellow "travelers". This review article aims to cover the recent development of experimental and computational methods for studying the biomechanics of collective cell migration during cancer progression and invasion. We also summarized the tested hypotheses regarding the mechanism underlying collective cell migration enabled by these methods. Together, the paper enables a broad overview on the methods and tools currently available to unravel the biophysical mechanisms pertinent to cell collective migration as well as providing perspectives on future development toward eventually deciphering the key mechanisms behind the most lethal feature of cancer.

Application of low-order potential solutions to higher-order vertical traction boundary problems in an elastic half-space

New solutions of potential functions for the bilinear vertical traction boundary condition are derived and presented. The discretization and interpolation of higher-order tractions and the superposition of the bilinear solutions provide a method of forming approximate and continuous solutions for the equilibrium state of a homogeneous and isotropic elastic half-space subjected to arbitrary normal surface tractions. Past experimental measurements of contact pressure distributions in granular media are reviewed in conjunction with the application of the proposed solution method to analysis of elastic settlement in shallow foundations. A numerical example is presented for an empirical 'saddle-shaped' traction distribution at the contact interface between a rigid square footing and a supporting soil medium. Non-dimensional soil resistance is computed as the reciprocal of normalized surface displacements under this empirical traction boundary condition, and the resulting internal stresses are compared to classical solutions to uniform traction boundary conditions.

Force fluctuations on a wall in interaction with a granular lid-driven cavity flow

The force fluctuations experienced by a boundary wall subjected to a lid-driven cavity flow are investigated by means of numerical simulations based on the discrete-element method. The time-averaged dynamics inside the cavity volume and the resulting steady force on the wall are governed by the boundary macroscopic inertial number, the latter being derived from the shearing velocity and the confinement pressure imposed at the top. The force fluctuations are quantified through measuring both the autocorrelation of force time series and the distributions of grain-wall forces, at distinct spatial scales from particle scale to wall scale. A key result is that the grain-wall force distributions are entirely driven by the boundary macroscopic inertial number, whatever the spatial scale considered. In particular, when the wall scale is considered, the distributions are found to evolve from nearly exponential to nearly Gaussian distributions by decreasing the macroscopic inertial number. The transition from quasistatic to dense inertial flow is well identified through remarkable changes in the shapes of the distributions of grain-wall forces, accompanied by a loss of system memory in terms of the mesoscale force transmitted toward the wall.

Modifying self-assembly and species separation in three-dimensional systems of shape-anisotropic particles

The behaviors of large, dynamic assemblies of macroscopic particles are of direct relevance to geophysical and industrial processes and may also be used as easily studied analogs to micro- or nano-scale systems, or model systems for microbiological, zoological, and even anthropological phenomena. We study vibrated mixtures of elongated particles, demonstrating that the inclusion of differing particle "species" may profoundly alter a system's dynamics and physical structure in various diverse manners. The phase behavior observed suggests that our system, despite its athermal nature, obeys a minimum free energy principle analogous to that observed for thermodynamic systems. We demonstrate that systems of exclusively spherical objects, which form the basis of numerous theoretical frameworks in many scientific disciplines, represent only a narrow region of a wide, multidimensional phase space. Thus, our results raise significant questions as to whether such models can accurately describe the behaviors of systems outside this highly specialized case.

Physics of active jamming during collective cellular motion in a monolayer

Although collective cell motion plays an important role, for example during wound healing, embryogenesis, or cancer progression, the fundamental rules governing this motion are still not well understood, in particular at high cell density. We study here the motion of human bronchial epithelial cells within a monolayer, over long times. We observe that, as the monolayer ages, the cells slow down monotonously, while the velocity correlation length first increases as the cells slow down but eventually decreases at the slowest motions. By comparing experiments, analytic model, and detailed particle-based simulations, we shed light on this biological amorphous solidification process, demonstrating that the observed dynamics can be explained as a consequence of the combined maturation and strengthening of cell-cell and cell-substrate adhesions. Surprisingly, the increase of cell surface density due to proliferation is only secondary in this process. This analysis is confirmed with two other cell types. The very general relations between the mean cell velocity and velocity correlation lengths, which apply for aggregates of self-propelled particles, as well as motile cells, can possibly be used to discriminate between various parameter changes in vivo, from noninvasive microscopy data.

Regime transitions of granular flow in a shear cell: a micromechanical study

The regime transitions of granular flow in a model shear cell are investigated numerically with a stress-controlled boundary condition. The correlations between the elastically and kinetically scaled stresses and the packing fraction are examined, and two packing fractions (0.58 and 0.50) are identified for the quasistatic to intermediate and intermediate to inertial regime transitions. The profiles and structures of contact networks and force chains among particles in different flow regimes are investigated. It is shown that the connectivity (coordination number) among particles and the homogeneity in the shear flow increase as the system goes through the inertial, intermediate, and then quasistatic regimes, and there is only little variation in the internal structure after the system has entered the quasistatic regime. Short-range force chains start to appear in the inertial regime, which also depend on the magnitude of the shear rate. The percolation of system-spanning force chains through the whole system is a characteristic of the onset of the quasistatic regime, which happens at a packing fraction that is close to the glass transition, i.e., about random loose packing (0.58) but far below the isotropic quasistatic (athermal) jamming packing fraction of random close packing (0.64). The tails of the probability density distribution P(f) of the scaled normal contact forces for the flows in different regimes are quantified by a stretched exponential P(f)=exp(-cf^{n}) with a remarkable finding that n ∼ 1.1 may be a potential demarcation point separating the quasistatic regime and the inertial or intermediate regimes.

Entangled granular media

We study the geometrically induced cohesion of ensembles of granular "u particles" that mechanically entangle through particle interpenetration. We vary the length-to-width ratio l/w of the u particles and form them into freestanding vertical columns. In a laboratory experiment, we monitor the response of the columns to sinusoidal vibration (with peak acceleration Γ). Column collapse occurs in a characteristic time τ which follows the relation τ∝exp(Γ/Δ). Δ resembles an activation energy and is maximal at intermediate l/w. A simulation reveals that optimal strength results from competition between packing and entanglement.

Power-law flow statistics in anisometric (wedge) hoppers

We find the probability for N particles to exit an anisometric (having unequal dimensions) hopper before jamming to have a broad power-law decay with exponent α = -2, in marked contrast to the exponential decay seen in hoppers with symmetric apertures. The transition from exponential to power law is explained by amodel that assumes particle motion is correlated over a distinct length scale. Hoppers with lengths larger than this length are modeled as a series of adjacent, statistically independent "cells." Experiments with apertures 27-37 particle diameters D long are well fit by a three-cell model, implying that the correlation length is ≈ 9-12D.

Steady flow of smooth, inelastic particles on a bumpy inclined plane: hard and soft particle simulations

We study smooth, slightly inelastic particles flowing under gravity on a bumpy inclined plane using event-driven and discrete-element simulations. Shallow layers (ten particle diameters) are used to enable simulation using the event-driven method within reasonable computational times. Steady flows are obtained in a narrow range of angles (13 degrees-14.5 degrees); lower angles result in stopping of the flow and higher angles in continuous acceleration. The flow is relatively dense with the solid volume fraction, nu approximately 0.5 , and significant layering of particles is observed. We derive expressions for the stress, heat flux, and dissipation for the hard and soft particle models from first principles. The computed mean velocity, temperature, stress, dissipation, and heat flux profiles of hard particles are compared to soft particle results for different values of stiffness constant (k). The value of stiffness constant for which results for hard and soft particles are identical is found to be k>or=2x10(6) mg/d, where m is the mass of a particle, g is the acceleration due to gravity, and d is the particle diameter. We compare the simulation results to constitutive relations obtained from the kinetic theory of Jenkins and Richman [J. T. Jenkins and M. W. Richman, Arch. Ration. Mech. Anal. 87, 355 (1985)] for pressure, dissipation, viscosity, and thermal conductivity. We find that all the quantities are very well predicted by kinetic theory for volume fractions nu<0.5. At higher densities, obtained for thicker layers (H=15d and H=20d), the kinetic theory does not give accurate prediction. Deviations of the kinetic theory predictions from simulation results are relatively small for dissipation and heat flux and most significant deviations are observed for shear viscosity and pressure. The results indicate the range of applicability of soft particle simulations and kinetic theory for dense flows.

Rheology of confined granular flows: scale invariance, glass transition, and friction weakening

We study fully developed, steady granular flows confined between parallel flat frictional sidewalls using numerical simulations and experiments. Above a critical rate, sidewall friction stabilizes the underlying heap at an inclination larger than the angle of repose. The shear rate is constant and independent of inclination over much of the flowing layer. In the direction normal to the free surface, the solid volume fraction increases on a scale equal to half the flowing layer depth. Beneath a critical depth at which internal friction is invariant, grains exhibit creeping and intermittent cage motion similar to that in glasses, causing gradual weakening of friction at the walls.

Spatial force correlations in granular shear flow. I. Numerical evidence

We investigate the emergence of correlations in granular shear flow. By increasing the density of a simulated granular flow, we observe a transition from a dilute regime, where interactions are dominated by binary collisions, to a dense regime characterized by large force networks and collective motions. With increasing density, interacting grains tend to form networks of simultaneous contacts due to the dissipative nature of collisions. We quantify the size of these networks by measuring two-point force correlations and find dramatic changes in the statistics of contact forces as the size of the networks increases.

Velocity correlations in dense granular shear flows: effects on energy dissipation and normal stress

We study the effect of precollisional velocity correlations on granular shear flow by molecular dynamics simulations of an inelastic hard sphere system. Comparison of the simulations with kinetic theory reveals that the theory overestimates both the energy dissipation rate and the normal stress in the dense flow region. We find that the relative normal velocity of colliding particles is smaller than that expected from random collisions, and the discrepancies in the dissipation and the normal stress can be adjusted by introducing the idea of the collisional temperature, from which we conclude that the velocity correlation neglected in the kinetic theory is responsible for the discrepancies. Our analysis of the distributions of the precollisional velocity suggests that the correlation grows through multiple inelastic collisions during the time scale of the inverse of the shear rate. As for the shear stress, the discrepancy is also found in the dense region, but it depends strongly on the particle inelasticity.

Towards a theoretical picture of dense granular flows down inclines

Unlike most fluids, granular materials include coexisting solid, liquid or gaseous regions, which produce a rich variety of complex flows. Dense flows down inclines preserve this complexity but remain simple enough for detailed analysis. In this review we survey recent advances in this rapidly evolving area of granular flow, with the aim of providing an organized, synthetic review of phenomena and a characterization of the state of understanding. The perspective that we adopt is influenced by the hope of obtaining a theory for dense, inclined flows that is based on assumptions that can be tested in physical experiments and numerical simulations, and that uses input parameters that can be independently measured. We focus on dense granular flows over three kinds of inclined surfaces: flat-frictional, bumpy-frictional and erodible. The wealth of information generated by experiments and numerical simulations for these flows has led to meaningful tests of relatively simple existing theories.

Thermodynamics and statistical mechanics of dense granular media

By detailed molecular dynamics and Monte Carlo simulations of a model system we show that granular materials at rest can be described as thermodynamics systems. First, we show that granular packs can be characterized by few parameters, as much as fluids or solids. Then, in a second independent step, we demonstrate that these states can be described in terms of equilibrium distributions which coincide with the statistical mechanics of powders first proposed by Edwards. We also derive the system equation of state as a function of the "configurational temperature," its new intensive thermodynamic parameter.

Structural signatures of the unjamming transition at zero temperature

We study the pair correlation function g(r) for zero-temperature, disordered, soft-sphere packings just above the onset of jamming. We find distinct signatures of the transition in both the first and split second peaks of this function. As the transition is approached from the jammed side (at higher packing fraction) the first peak diverges and narrows on the small-r side to a delta function. On the high-r side of this peak, g(r) decays as a power law. In the split second peak, the two subpeaks are both singular at the transition, with power-law behavior on their low-r sides and step-function drop-offs on their high-r sides. These singularities at the transition are reminiscent of empirical criteria that have previously been used to distinguish glassy structures from liquid ones.

Sheared force networks: anisotropies, yielding, and geometry

A scenario for the yielding of granular matter is presented by considering the ensemble of force networks for a given contact network and applied shear stress tau. As tau is increased, the probability distribution of contact forces becomes highly anisotropic, the difference between average contact forces along minor and major axes grows, and the allowed networks span a shrinking subspace of all force networks. Eventually, contacts start to break, and at the maximal shear stress the packing becomes effectively isostatic. The size of the allowed subspace exhibits simple scaling properties, which lead to a prediction for the yield stress for packings of an arbitrary contact number.

A Lattice Model of Vitrification and Gelation

In this work we introduce a simple lattice model with T-shaped molecules in two dimensions that exhibits a rich range of morphological behaviors. Depending on the volume fraction and quench path, this system can adopt uniform liquid, solution, and phase-separated states, as well as inhomogeneous glass or gel-like states, as revealed by dynamic mean-field simulations. An important characteristic of this system is the existence of a large number of degenerate low-energy states with small barriers that leads to a broad, kinetically explored landscape. The mean-field stability and phase diagram of this model is constructed and provides a useful guide for understanding the complex behaviors of the system. One striking feature is that there is a cascade of instabilities that converge to mark the onset of what we identify as the glass transition. Both dynamic mean-field and Monte Carlo simulations reveal glass-like relaxation dynamics. Our results lead to a picture of gelation as a continuation of the glass transition into the two-phase region, or equivalently, as an incomplete phase separation arrested by the onset of the glass transition.

Shear-induced overaging in a polymer glass

A phenomenon recently coined as overaging implies a slowdown in the collective (slow) relaxation modes of a glass when a transient shear strain is imposed. We are able to reproduce this behavior in simulations of a supercooled polymer melt by imposing instantaneous shear deformations. The increase in relaxation times Delta(tau(1/2)) rises rapidly with deformation, becoming exponential in the plastic regime, and is accompanied by significant changes in the distribution of these relaxation times throughout the system. This overaging is distinct from standard aging. We find increases in pressure, bond-orientational order, and in the average energy of the inherent structures (<e(IS)>) of the system, all dependent on the size of the deformation. The observed change in behavior from elastic to plastic deformation suggests a link to the physics of the "jammed state."

Isoenergetic jamming transition in particle-filled systems

A full understanding of the jamming transition remains elusive, but recent advances which draw upon the common features of frustrated systems are encouraging. Herein, we show that, for mixtures of oil and silica particles, the dependence of the dejamming stress on filler volume fraction, phi is consistent with the shape of a reported jamming phase diagram [Trappe, Nature (London) 411, 772 (2001)]. We discover for the first time, however, that the role of phi disappears when mechanical energy input, defined as stress multiplied by strain, is used instead of stress as the critical parameter. We also examine literature results for aqueous suspensions of boehmite alumina powders, latex dispersions of polystyrene particles, and carbon black-filled elastomers in order to illustrate the universality of our finding. This study provides evidence for a thermodynamic interpretation of the jamming transition.

Structural signature of jamming in granular media

Glasses are rigid, but flow when the temperature is increased. Similarly, granular materials are rigid, but become unjammed and flow if sufficient shear stress is applied. The rigid and flowing phases are strikingly different, yet measurements reveal that the structures of glass and liquid are virtually indistinguishable. It is therefore natural to ask whether there is a structural signature of the jammed granular state that distinguishes it from its flowing counterpart. Here we find evidence for such a signature, by measuring the contact-force distribution between particles during shearing. Because the forces are sensitive to minute variations in particle position, the distribution of forces can serve as a microscope with which to observe correlations in the positions of nearest neighbours. We find a qualitative change in the force distribution at the onset of jamming. If, as has been proposed, the jamming and glass transitions are related, our observation of a structural signature associated with jamming hints at the existence of a similar structural difference at the glass transition--presumably too subtle for conventional scattering techniques to uncover. Our measurements also provide a determination of a granular temperature that is the counterpart in granular systems to the glass-transition temperature in liquids.

Jamming transition in emulsions and granular materials

We investigate the jamming transition in packings of emulsions and granular materials via molecular dynamics simulations. The emulsion model is composed of frictionless droplets interacting via nonlinear normal forces obtained using experimental data acquired by confocal microscopy of compressed emulsions systems. Granular materials are modeled by Hertz-Mindlin deformable spherical grains with Coulomb friction. In both cases, we find power-law scaling for the vanishing of pressure and excess number of contacts as the system approaches the jamming transition from high volume fractions. We find that the construction history parametrized by the compression rate during the preparation protocol has a strong effect on the micromechanical properties of granular materials but not on emulsions. This leads the granular system to jam at different volume fractions depending on the histories. Isostaticity is found in the packings close to the jamming transition in emulsions and in granular materials at slow compression rates and infinite friction. Heterogeneity of interparticle forces increases as the packings approach the jamming transition which is demonstrated by the exponential tail in force distributions and the small values of the participation number measuring spatial localization of the forces. However, no signatures of the jamming transition are observed in structural properties, like the radial distribution functions and the distributions of contacts.

Temporally heterogeneous dynamics in granular flows

Granular simulations are used to probe the particle scale dynamics at short, intermediate, and long time scales for gravity-driven, dense granular flows down an inclined plane. On approach to the angle of repose, where motion ceases, the dynamics become intermittent over intermediate times, with strong temporal correlations between particle motions-temporally heterogeneous dynamics. This intermittency is characterized through large-scale structural events whereby the contact network periodically spans the system. A characteristic time scale associated with these processes increases as the stopped state is approached. These features are discussed in the context of the dynamics of supercooled liquids near the glass transition.

Ensemble theory for force networks in hyperstatic granular matter

An ensemble approach for force networks in static granular packings is developed. The framework is based on the separation of packing and force scales, together with an a priori flat measure in the force phase space under the constraints that the contact forces are repulsive and balance on every particle. In this paper we will give a general formulation of this force network ensemble, and derive the general expression for the force distribution P(f). For small regular packings these probability densities are obtained in closed form, while for larger packings we present a systematic numerical analysis. Since technically the problem can be written as a noninvertible matrix problem (where the matrix is determined by the contact geometry), we study what happens if we perturb the packing matrix or replace it by a random matrix. The resulting P(f) 's differ significantly from those of normal packings, which touches upon the deep question of how network statistics is related to the underlying network structure. Overall, the ensemble formulation opens up a different perspective on force networks that is analytically accessible, and which may find applications beyond granular matter.

Packing geometry and statistics of force networks in granular media

The relation between packing geometry and force network statistics is studied for granular media. Based on simulations of two-dimensional packings of Hertzian spheres, we develop a geometrical framework relating the distribution of interparticle forces P(f) to the weight distribution P(w), which is measured in experiments. We apply this framework to reinterpret recent experimental data on strongly deformed packings and suggest that the observed changes of P(w) are dominated by changes in contact network while P(f) remains relatively unaltered. We furthermore investigate the role of packing disorder in the context of the q model and address the question of how force fluctuations build up as a function of the distance beneath the top surface.

NMR experiments on a three-dimensional vibrofluidized granular medium

A three-dimensional granular system fluidized by vertical container vibrations was studied using pulsed field gradient NMR coupled with one-dimensional magnetic resonance imaging. The system consisted of mustard seeds vibrated vertically at 50 Hz, and the number of layers N(l)<or=4 was sufficiently low to achieve a nearly time-independent granular fluid. Using NMR, the vertical profiles of density and granular temperature were directly measured, along with the distributions of vertical and horizontal grain velocities. The velocity distributions showed modest deviations from Maxwell-Boltzmann statistics, except for the vertical velocity distribution near the sample bottom, which was highly skewed and non-Gaussian. Data taken for three values of N(l) and two dimensionless accelerations Gamma=15,18 were fitted to a hydrodynamic theory, which successfully models the density and temperature profiles away from the vibrating container bottom. A temperature inversion near the free upper surface is observed, in agreement with predictions based on the hydrodynamic parameter micro which is nonzero only in inelastic systems.

Force network ensemble: a new approach to static granular matter

An ensemble approach for force distributions in static granular packings is developed. This framework is based on the separation of packing and force scales, together with an a priori flat measure in the force phase space under the constraints that the contact forces are repulsive and balance on every particle. We show how the formalism yields realistic results, both for disordered and regular triangular "snooker ball" configurations, and obtain a shear-induced unjamming transition of the type proposed recently for athermal media.

Partially fluidized shear granular flows: continuum theory and molecular dynamics simulations

The continuum theory of partially fluidized shear granular flows is tested and calibrated using two-dimensional soft particle molecular dynamics simulations. The theory is based on the relaxational dynamics of the order parameter that describes the transition between static and flowing regimes of granular material. We define the order parameter as a fraction of static contacts among all contacts between particles. We also propose and verify by direct simulations the constitutive relation based on the splitting of the shear stress tensor into a"fluid part" proportional to the strain rate tensor, and a remaining "solid part." The ratio of these two parts is a function of the order parameter. The rheology of the fluid component agrees well with the kinetic theory of granular fluids even in the dense regime. Based on the hysteretic bifurcation diagram for a thin shear granular layer obtained in simulations, we construct the "free energy" for the order parameter. The theory calibrated using numerical experiments with the thin granular layer is applied to the surface-driven stationary two-dimensional granular flows in a thick granular layer under gravity.

Jamming at zero temperature and zero applied stress: the epitome of disorder

We have studied how two- and three-dimensional systems made up of particles interacting with finite range, repulsive potentials jam (i.e., develop a yield stress in a disordered state) at zero temperature and zero applied stress. At low packing fractions phi, the system is not jammed and each particle can move without impediment from its neighbors. For each configuration, there is a unique jamming threshold phi(c) at which particles can no longer avoid each other, and the bulk and shear moduli simultaneously become nonzero. The distribution of phi(c) values becomes narrower as the system size increases, so that essentially all configurations jam at the same packing fraction in the thermodynamic limit. This packing fraction corresponds to the previously measured value for random close packing. In fact, our results provide a well-defined meaning for "random close packing" in terms of the fraction of all phase space with inherent structures that jam. The jamming threshold, point J, occurring at zero temperature and applied stress and at the random-close-packing density, has properties reminiscent of an ordinary critical point. As point J is approached from higher packing fractions, power-law scaling is found for the divergence of the first peak in the pair correlation function and in the vanishing of the pressure, shear modulus, and excess number of overlapping neighbors. Moreover, near point J, certain quantities no longer self-average, suggesting the existence of a length scale that diverges at J. However, point J also differs from an ordinary critical point: the scaling exponents do not depend on dimension but do depend on the interparticle potential. Finally, as point J is approached from high packing fractions, the density of vibrational states develops a large excess of low-frequency modes. Indeed, at point J, the density of states is a constant all the way down to zero frequency. All of these results suggest that point J is a point of maximal disorder and may control behavior in its vicinity-perhaps even at the glass transition.

Order parameter description of stationary partially fluidized shear granular flows

We carry out a detailed comparison of soft-particle molecular dynamics simulations with the theory of partially fluidized shear granular flows. We verify by direct simulations a constitutive relation based on the separation of the shear stress tensor into a fluid part proportional to the strain rate tensor, and a remaining solid part. The ratio of these two components is determined by the order parameter. Based on results of the simulations we construct the "free energy" function for the order parameter. We also present the simulations of the stationary deep 2D granular flows driven by an upper wall and compare it with the continuum theory.

Role of rigidity in the fluid-solid transition

We examine the fluid-solid transition for a hard-disk system. By counting the near neighbors in the average configurations of a grand-canonical Monte Carlo simulation, this enables us to relate the thermodynamic transition with the rigidity theory, since we find that the coordination number in the fluid-solid transition is close to the coordination number predicted by a mean field rigidity theory, due to dynamical jamming of particles, where the contact region between disks is the radial ring outside a disk with a maximum allowed coordination number that is not bigger than six. Using these ideas, we were able to produce a continuous glass-like transition when nucleation of rigidity is suppressed.

Confined granular packings: structure, stress, and forces

The structure and stresses of static granular packs in cylindrical containers are studied by using large-scale discrete element molecular dynamics simulations in three dimensions. We generate packings by both pouring and sedimentation and examine how the final state depends on the method of construction. The vertical stress becomes depth independent for deep piles and we compare these stress depth profiles to the classical Janssen theory. The majority of the tangential forces for particle-wall contacts are found to be close to the Coulomb failure criterion, in agreement with the theory of Janssen, while particle-particle contacts in the bulk are far from the Coulomb criterion. In addition, we show that a linear hydrostaticlike region at the top of the packings unexplained by the Janssen theory arises because most of the particle-wall tangential forces in this region are far from the Coulomb yield criterion. The distributions of particle-particle and particle-wall contact forces P(f) exhibit exponential-like decay at large forces in agreement with previous studies.