Search

More than 30 years ago Edwards and co-authors proposed a model to describe the statistics of granular packings by an ensemble of equiprobable jammed states. Experimental tests of this model remained scarce so far. We introduce a simple system to analyze statistical properties of jammed granular ensembles to test Edwards theory. Identical spheres packed in a nearly two-dimensional geometrical confinement were studied in experiments and numerical simulations. When tapped, the system evolves towards a ground state, but due to incompatible domain structures it gets trapped. Analytical calculations reproduce relatively well our simulation results, which allows us to test Edwards theory on a coupled system of two subsystems with different properties. We find that the joint system can only be described by the Edwards theory if considered as a single system due to the constraints in the stresses. The results show counterintuitive effects as in the coupled system the change in the order parameter is opposite to what is expected from the change in the compactivity.

A trapping mechanism for propelled colloidal particles based on an inhomogeneous drive is presented and studied by means of computer simulations. In experiments this method can be realized using photophoretic Janus particles driven by a light source, which shines through a shading mask and leads to an accumulation of the particles in the passive part. An equation for an accumulation parameter is derived using the effective inhomogeneous diffusion constant generated by the inhomogeneous drive. The impact of particle interaction on the trapping mechanism is studied, as well as the interplay between passivity-induced trapping and the emergent self-clustering of systems containing a high density of active particles. The combination of both effects makes the clusters more controllable for applications., Comment: 7 pages, 4 figures

While the self-learning kinetic Monte Carlo (SLKMC) method enables the calculation of transition rates from a realistic potential, implementations of it were usually limited to one specific surface orientation. An example is the fcc (111) surface in Latz et al. 2012, J. Phys.: Condens. Matter 24, 485005. This work provides an extension by means of detecting the local orientation, and thus allows for the accurate simulation of arbitrarily shaped surfaces. We applied the model to the diffusion of Ag monolayer islands and voids on a Ag(111) and Ag(001) surface, as well as the relaxation of a three-dimensional spherical particle.

A fully parallel version of the contact dynamics (CD) method is presented in this paper. For large enough systems, 100% efficiency has been demonstrated for up to 256 processors using a hierarchical domain decomposition with dynamic load balancing. The iterative scheme to calculate the contact forces is left domain-wise sequential, with data exchange after each iteration step, which ensures its stability. The number of additional iterations required for convergence by the partially parallel updates at the domain boundaries becomes negligible with increasing number of particles, which allows for an effective parallelization. Compared to the sequential implementation, we found no influence of the parallelization on simulation results., Comment: 19 pages, 15 figures, published in Journal of Computational Physics (2011)

Energy dissipation is studied for a hard magnetic tip that scans a soft magnetic substrate. The dynamics of the atomic moments are simulated by solving the Landau-Lifshitz-Gilbert (LLG) equation numerically. The local energy currents are analysed for the case of a Heisenberg spin chain taken as substrate. This leads to an explanation for the velocity dependence of the friction force: The non-linear contribution for high velocities can be attributed to a spin wave front pushed by the tip along the substrate., Comment: 5 pages, 9 figures

True biaxial tests of granular materials are investigated by applying the principle of minimal dissipation and comparing to two dimensional contact dynamics simulations. It is shown that the macroscopic steady state manifested by constant stress ratio and constant volume is the result of the ever changing microscopic structure which minimizes the dissipation rate. The shear band angle in the varying shear band structures is found to be constant. We also show that introducing friction on the walls reduces the degeneracy of the optimal shear band structures to one for a wide range of parameters which gives a non-constant stress ratio curve with varying aspect ratio that can be calculated., Comment: 8 pages, 8 figures

Cohesive powders form agglomerates that can be very porous. Hence they are also very fragile. Consider a process of complete fragmentation on a characteristic length scale $\ell$, where the fragments are subsequently allowed to settle under gravity. If this fragmentation-reagglomeration cycle is repeated sufficiently often, the powder develops a fractal substructure with robust statistical properties. The structural evolution is discussed for two different models: The first one is an off-lattice model, in which a fragment does not stick to the surface of other fragments that have already settled, but rolls down until it finds a locally stable position. The second one is a simpler lattice model, in which a fragment sticks at first contact with the agglomerate of fragments that have already settled. Results for the fragment size distribution are shown as well. One can distinguish scale invariant dust and fragments of a characteristic size. Their role in the process of structure formation will be addressed., Comment: 4 pages, 8 figures

Energy dissipation via spin excitations is investigated for a hard ferromagnetic tip scanning a soft magnetic monolayer. We use the classical Heisenberg model with Landau-Lifshitz-Gilbert (LLG)-dynamics including a stochastic field representing finite temperatures. The friction force depends linearly on the velocity (provided it is small enough) for all temperatures. For low temperatures, the corresponding friction coefficient is proportional to the phenomenological damping constant of the LLG equation. This dependence is lost at high temperatures, where the friction coefficient decreases exponentially. These findings can be explained by properties of the spin polarization cloud dragged along with the tip., Comment: 6 pages, 5 figures

Cohesive powders tend to form porous aggregates which can be compacted by applying an external pressure. This process is modelled using the Contact Dynamics method supplemented with a cohesion law and rolling friction. Starting with ballistic deposits of varying density, we investigate how the porosity of the compacted sample depends on the cohesion strength and the friction coefficients. This allows to explain different pore stabilization mechanisms. The final porosity depends on the cohesion force scaled by the external pressure and on the lateral distance between branches of the ballistic deposit r_capt. Even if cohesion is switched off, pores can be stabilized by Coulomb friction alone. This effect is weak for round particles, as long as the friction coefficient is smaller than 1. However, for nonspherical particles the effect is much stronger., Comment: 10 pages, 15 figures

The systematic errors due to the practical implementation of the Contact Dynamics method for simulation of dense granular media are examined. It is shown that, using the usual iterative solver to simulate a chain of rigid particles, effective elasticity and sound propagation with a finite velocity occur. The characteristics of these phenomena are investigated analytically and numerically in order to assess the limits of applicability of this simulation method and to compare it with soft particle molecular dynamics., Comment: submitted to PRE, 7 pages, 6 figures