Jamming transition in a two-dimensional open granular pile with rolling resistance

We present a molecular dynamics study of the jamming/unjamming transition in two-dimensional granular piles with open boundaries. The grains are modeled by viscoelastic forces, Coulomb friction and resistance to rolling. Two models for the rolling resistance interaction were assessed: one considers a constant rolling friction coefficient, and the other one a strain dependent coefficient. The piles are grown on a finite size substrate and subsequently discharged through an orifice opened at the center of the substrate. Varying the orifice width and taking the final height of the pile after the discharge as the order parameter, one can devise a transition from a jammed regime (when the grain flux is always clogged by an arch) to a catastrophic regime, in which the pile is completely destroyed by an avalanche as large as the system size. A finite size analysis shows that there is a finite orifice width associated with the threshold for the unjamming transition, no matter the model used for the microscopic interactions. As expected, the value of this threshold width increases when rolling resistance is considered, and it depends on the model used for the rolling friction.Comment: 9 pages, 6 figure

Scale separation in granular packings: stress plateaus and fluctuations

It is demonstrated, by numerical simulations of a 2D assembly of polydisperse disks, that there exists a range (plateau) of coarse graining scales for which the stress tensor field in a granular solid is nearly resolution independent, thereby enabling an `objective' definition of this field. Expectedly, it is not the mere size of the the system but the (related) magnitudes of the gradients that determine the widths of the plateaus. Ensemble averaging (even over `small' ensembles) extends the widths of the plateaus to sub-particle scales. The fluctuations within the ensemble are studied as well. Both the response to homogeneous forcing and to an external compressive localized load (and gravity) are studied. Implications to small solid systems and constitutive relations are briefly discussed.Comment: 4 pages, 4 figures, RevTeX 4, Minor corrections to match the published versio

A knowledge-based view of the extending enterprise for enhancing a collaborative innovation advantage

In animal societies as well as in human crowds, many observed collective behaviours result from self-organized processes based on local interactions among individuals. However, models of crowd dynamics are still lacking a systematic individual-level experimental verification, and the local mechanisms underlying the formation of collective patterns are not yet known in detail. We have conducted a set of well-controlled experiments with pedestrians performing simple avoidance tasks in order to determine the laws ruling their behaviour during interactions. The analysis of the large trajectory dataset was used to compute a behavioural map that describes the average change of the direction and speed of a pedestrian for various interaction distances and angles. The experimental results reveal features of the decision process when pedestrians choose the side on which they evade, and show a side preference that is amplified by mutual interactions. The predictions of a binary interaction model based on the above findings were then compared to bidirectional flows of people recorded in a crowded street. Simulations generate two asymmetric lanes with opposite directions of motion, in quantitative agreement with our empirical observations. The knowledge of pedestrian behavioural laws is an important step ahead in the understanding of the underlying dynamics of crowd behaviour and allows for reliable predictions of collective pedestrian movements under natural conditions

Relaxation kinetics in two-dimensional structures

We have studied the approach to equilibrium of islands and pores in two dimensions. The two-regime scenario observed when islands evolve according to a set of particular rules, namely relaxation by steps at low temperature and smooth at high temperature, is generalized to a wide class of kinetic models and the two kinds of structures. Scaling laws for equilibration times are analytically derived and confirmed by kinetic Monte Carlo simulations.Comment: 6 pages, 7 figures, 1 tabl

The contamination of the surface of Vesta by impacts and the delivery of the dark material

The Dawn spacecraft observed the presence of dark material, which in turn proved to be associated with OH and H-rich material, on the surface of Vesta. The source of this dark material has been identified with the low albedo asteroids, but it is still a matter of debate whether the delivery of the dark material is associated with a few large impact events, to micrometeorites or to the continuous, secular flux of impactors on Vesta. The continuous flux scenario predicts that a significant fraction of the exogenous material accreted by Vesta should be due to non-dark impactors likely analogous to ordinary chondrites, which instead represent only a minor contaminant in the HED meteorites. We explored the continuous flux scenario and its implications for the composition of the vestan regolith, taking advantage of the data from the Dawn mission and the HED meteorites. We used our model to show that the stochastic events scenario and the micrometeoritic flux scenario are natural consequences of the continuous flux scenario. We then used the model to estimate the amounts of dark and hydroxylate materials delivered on Vesta since the LHB and we showed how our results match well with the values estimated by the Dawn mission. We used our model to assess the amount of Fe and siderophile elements that the continuous flux of impactors would mix in the vestan regolith: concerning the siderophile elements, we focused our attention on the role of Ni. The results are in agreement with the data available on the Fe and Ni content of the HED meteorites and can be used as a reference frame in future studies of the data from the Dawn mission and of the HED meteorites. Our model cannot yet provide an answer to the fate of the missing non-carbonaceous contaminants, but we discuss possible reasons for this discrepancy.Comment: 31 pages, 7 figures, 4 tables. Accepted for publication on the journal ICARUS, "Dark and Bright Materials on Vesta" special issu

Internal states of model isotropic granular packings. I. Assembling process, geometry and contact networks

This is the first paper of a series of three, reporting on numerical simulation studies of geometric and mechanical properties of static assemblies of spherical beads under an isotropic pressure. Frictionless systems assemble in the unique random close packing (RCP) state in the low pressure limit if the compression process is fast enough, slower processes inducing traces of crystallization, and exhibit specific properties directly related to isostaticity of the force-carrying structure. The different structures of frictional packings assembled by various methods cannot be classified by the sole density. While lubricated systems approach RCP densities and coordination number z^*~=6 on the backbone in the rigid limit, an idealized "vibration" procedure results in equally dense configurations with z^*~=4.5. Near neighbor correlations on various scales are computed and compared to available laboratory data, although z^* values remain experimentally inaccessible. Low coordination packings have many rattlers (more than 10% of the grains carry no force), which should be accounted for on studying position correlations, and a small proportion of harmless "floppy modes" associated with divalent grains. Frictional packings, however slowly assembled under low pressure, retain a finite level of force indeterminacy, except in the limit of infinite friction.Comment: 29 pages. Published in Physical Review

Devil's staircase in kinetically limited growth of Ising model

The devil's staircase is a term used to describe surface or an equilibrium phase diagram in which various ordered facets or phases are infinitely closely packed as a function of some model parameter. A classic example is a 1-D Ising model [bak] wherein long-range and short range forces compete, and the periodicity of the gaps between minority species covers all rational values. In many physical cases, crystal growth proceeds by adding surface layers which have the lowest energy, but are then frozen in place. The emerging layered structure is not the thermodynamic ground state, but is uniquely defined by the growth kinetics. It is shown that for such a system, the grown structure tends to the equilibrium ground state via a devil's staircase traversing an infinity of intermediate phases. It would be extremely difficult to deduce the simple growth law based on measurement made on such an grown structure.Comment: 4 pages, PRL submitte

Group Irregularity Strength of Connected Graphs

We investigate the group irregularity strength ($s_g(G)$) of graphs, i.e. the smallest value of $s$ such that taking any Abelian group \gr of order $s$, there exists a function f:E(G)\rightarrow \gr such that the sums of edge labels at every vertex are distinct. We prove that for any connected graph $G$ of order at least 3, $s_g(G)=n$ if $n\neq 4k+2$ and $s_g(G)\leq n+1$ otherwise, except the case of some infinite family of stars

Diffusion of gold nanoclusters on graphite

We present a detailed molecular-dynamics study of the diffusion and coalescence of large (249-atom) gold clusters on graphite surfaces. The diffusivity of monoclusters is found to be comparable to that for single adatoms. Likewise, and even more important, cluster dimers are also found to diffuse at a rate which is comparable to that for adatoms and monoclusters. As a consequence, large islands formed by cluster aggregation are also expected to be mobile. Using kinetic Monte Carlo simulations, and assuming a proper scaling law for the dependence on size of the diffusivity of large clusters, we find that islands consisting of as many as 100 monoclusters should exhibit significant mobility. This result has profound implications for the morphology of cluster-assembled materials