Vibrations and diverging length scales near the unjamming transition

We numerically study the vibrations of jammed packings of particles interacting with finite-range, repulsive potentials at zero temperature. As the packing fraction $\phi$ is lowered towards the onset of unjamming at $\phi_{c}$, the density of vibrational states approaches a non-zero value in the limit of zero frequency. For $\phi>\phi_{c}$, there is a crossover frequency, $\omega^{*}$ below which the density of states drops towards zero. This crossover frequency obeys power-law scaling with $\phi-\phi_{c}$. Characteristic length scales, determined from the dominant wavevector contributing to the eigenmode at $\omega^{*}$, diverge as power-laws at the unjamming transition.Comment: Submitted to PRL, 4 pages + 7 .eps figure

Granular flow down a rough inclined plane: transition between thin and thick piles

The rheology of granular particles in an inclined plane geometry is studied using molecular dynamics simulations. The flow--no-flow boundary is determined for piles of varying heights over a range of inclination angles $\theta$. Three angles determine the phase diagram: $\theta_{r}$, the angle of repose, is the angle at which a flowing system comes to rest; $\theta_{m}$, the maximum angle of stability, is the inclination required to induce flow in a static system; and $\theta_{max}$ is the maximum angle for which stable, steady state flow is observed. In the stable flow region $\theta_{r}<\theta<\theta_{max}$, three flow regimes can be distinguished that depend on how close $\theta$ is to $\theta_{r}$: i) $\theta>>\theta_{r}$: Bagnold rheology, characterized by a mean particle velocity $v_{x}$ in the direction of flow that scales as $v_{x}\propto h^{3/2}$, for a pile of height $h$, ii) $\theta\gtrsim\theta_{r}$: the slow flow regime, characterized by a linear velocity profile with depth, and iii) $\theta\approx\theta_{r}$: avalanche flow characterized by a slow underlying creep motion combined with occasional free surface events and large energy fluctuations. We also probe the physics of the initiation and cessation of flow. The results are compared to several recent experimental studies on chute flows and suggest that differences between measured velocity profiles in these experiments may simply be a consequence of how far the system is from jamming.Comment: 19 pages, 14 figs, submitted to Physics of Fluid

Normal Modes in Model Jammed Systems in Three Dimensions

Vibrational spectra and normal modes of mechanically stable particle packings in three dimensions are analyzed over a range of compressions, from near the jamming transition, where the packings lose their rigidity, to far above it. At high frequency, the normal modes are localized at all compressions. At low frequency, the nature of the modes depends somewhat on compression. At large compressions, far from the transition, the lowest-frequency normal modes have some plane-wave character, though less than one would expect for a crystalline or isotropic solid. At low compressions near the jamming transition, the lowest-frequency modes are neither plane-wave-like nor localized. We characterize these differences, highlighting the unusual dispersion behavior that emerges for marginally jammed solids.Comment: Under review at Phys. Rev. E. Lower resolution figures her

Geometric origin of excess low-frequency vibrational modes in amorphous solids

Glasses have a large excess of low-frequency vibrational modes in comparison with crystalline solids. We show that such a feature is a necessary consequence of the geometry generic to weakly connected solids. In particular, we analyze the density of states of a recently simulated system, comprised of weakly compressed spheres at zero temperature. We account for the observed a) constancy of the density of modes with frequency, b) appearance of a low-frequency cutoff, and c) power-law increase of this cutoff with compression. We predict a length scale below which vibrations are very different from those of a continuous elastic body.Comment: 4 pages, 2 figures. Argument rewritten, identical result

Density of states in random lattices with translational invariance

We propose a random matrix approach to describe vibrational excitations in disordered systems. The dynamical matrix M is taken in the form M=AA^T where A is some real (not generally symmetric) random matrix. It guaranties that M is a positive definite matrix which is necessary for mechanical stability of the system. We built matrix A on a simple cubic lattice with translational invariance and interaction between nearest neighbors. We found that for certain type of disorder phonons cannot propagate through the lattice and the density of states g(w) is a constant at small w. The reason is a breakdown of affine assumptions and inapplicability of the elasticity theory. Young modulus goes to zero in the thermodynamic limit. It strongly reminds of the properties of a granular matter at the jamming transition point. Most of the vibrations are delocalized and similar to diffusons introduced by Allen, Feldman et al., Phil. Mag. B v.79, 1715 (1999).Comment: 4 pages, 5 figure

Stress response inside perturbed particle assemblies

The effect of structural disorder on the stress response inside three dimensional particle assemblies is studied using computer simulations of frictionless sphere packings. Upon applying a localised, perturbative force within the packings, the resulting {\it Green's} function response is mapped inside the different assemblies, thus providing an explicit view as to how the imposed perturbation is transmitted through the packing. In weakly disordered arrays, the resulting transmission of forces is of the double-peak variety, but with peak widths scaling linearly with distance from the source of the perturbation. This behaviour is consistent with an anisotropic elasticity response profile. Increasing the disorder distorts the response function until a single-peak response is obtained for fully disordered packings consistent with an isotropic description.Comment: 8 pages, 7 figure captions To appear in Granular Matte

Sparse random matrices and vibrational spectra of amorphous solids

A random matrix approach is used to analyze the vibrational properties of amorphous solids. We investigated a dynamical matrix M=AA^T with non-negative eigenvalues. The matrix A is an arbitrary real NxN sparse random matrix with n independent non-zero elements in each row. The average values =0 and dispersion =V^2 for all non-zero elements. The density of vibrational states g(w) of the matrix M for N,n >> 1 is given by the Wigner quarter circle law with radius independent of N. We argue that for n^2 << N this model can be used to describe the interaction of atoms in amorphous solids. The level statistics of matrix M is well described by the Wigner surmise and corresponds to repulsion of eigenfrequencies. The participation ratio for the major part of vibrational modes in three dimensional system is about 0.2 - 0.3 and independent of N. Together with term repulsion it indicates clearly to the delocalization of vibrational excitations. We show that these vibrations spread in space by means of diffusion. In this respect they are similar to diffusons introduced by Allen, Feldman, et al., Phil. Mag. B 79, 1715 (1999) in amorphous silicon. Our results are in a qualitative and sometimes in a quantitative agreement with molecular dynamic simulations of real and model glasses.Comment: 24 pages, 7 figure

Patterns and Collective Behavior in Granular Media: Theoretical Concepts

Granular materials are ubiquitous in our daily lives. While they have been a subject of intensive engineering research for centuries, in the last decade granular matter attracted significant attention of physicists. Yet despite a major efforts by many groups, the theoretical description of granular systems remains largely a plethora of different, often contradicting concepts and approaches. Authors give an overview of various theoretical models emerged in the physics of granular matter, with the focus on the onset of collective behavior and pattern formation. Their aim is two-fold: to identify general principles common for granular systems and other complex non-equilibrium systems, and to elucidate important distinctions between collective behavior in granular and continuum pattern-forming systems.Comment: Submitted to Reviews of Modern Physics. Full text with figures (2Mb pdf) avaliable at http://mti.msd.anl.gov/AransonTsimringReview/aranson_tsimring.pdf Community responce is appreciated. Comments/suggestions send to [email protected]

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 (PFG) NMR coupled with one-dimensional magnetic resonance imaging (MRI). The system consisted of mustard seeds vibrated vertically at 50 Hz, and the number of layers N_ell <= 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_ell and two dimensionless accelerations Gamma=15,18 were fit to a hydrodynamic theory, which successfully models the density and temperature profiles including a temperature inversion near the free upper surface.Comment: 14 pages, 15 figure

Planetary Dynamics and Habitable Planet Formation In Binary Star Systems

Whether binaries can harbor potentially habitable planets depends on several factors including the physical properties and the orbital characteristics of the binary system. While the former determines the location of the habitable zone (HZ), the latter affects the dynamics of the material from which terrestrial planets are formed (i.e., planetesimals and planetary embryos), and drives the final architecture of the planets assembly. In order for a habitable planet to form in a binary star system, these two factors have to work in harmony. That is, the orbital dynamics of the two stars and their interactions with the planet-forming material have to allow terrestrial planet formation in the habitable zone, and ensure that the orbit of a potentially habitable planet will be stable for long times. We have organized this chapter with the same order in mind. We begin by presenting a general discussion on the motion of planets in binary stars and their stability. We then discuss the stability of terrestrial planets, and the formation of potentially habitable planets in a binary-planetary system.Comment: 56 pages, 29 figures, chapter to appear in the book: Planets in Binary Star Systems (Ed. N. Haghighipour, Springer publishing company