Shear-induced self-diffusion in non-colloidal suspensions

Self-diffusion in a monodisperse suspension of non-Brownian particles in simple shear flow is studied using accelerated Stokesian dynamics (ASD) simulation. The availability of a much faster computational algorithm allows the study of large systems (typically of 1000 particles) and the extraction of accurate results for the complete shear-induced self-diffusivity tensor. The finite, and often large, autocorrelation time requires the mean-square displacements to be followed for very long times, which is now possible with ASD. The self-diffusivities compare favourably with the available experimental measurements when allowance is made for the finite strains sampled in the experiments. The relationship between the mean-square displacements and the diffusivities appearing in a Fokker–Planck equation when advection couples to diffusion is discussed

Anisotropic anomalous diffusion modulated by log-periodic oscillations

We introduce finite ramified self-affine substrates in two dimensions with a set of appropriate hopping rates between nearest-neighbor sites, where the diffusion of a single random walk presents an anomalous {\it anisotropic} behavior modulated by log-periodic oscillations. The anisotropy is revealed by two different random walk exponents, $\nu_x$ and $\nu_y$, in the {\it x} and {\it y} direction, respectively. The values of these exponents, as well as the period of the oscillation, are analytically obtained and confirmed by Monte Carlo simulations.Comment: 7 pages, 7 figure

Microscopic elasticity of complex systems

Lecture Notes for the Erice Summer School 2005 Computer Simulations in Condensed Matter: from Materials to Chemical Biology. Perspectives in celebration of the 65th Birthday of Mike Klein organized by Kurt Binder, Giovanni Ciccotti and Mauro Ferrari

Amorphous Systems in Athermal, Quasistatic Shear

We present results on a series of 2D atomistic computer simulations of amorphous systems subjected to simple shear in the athermal, quasistatic limit. The athermal quasistatic trajectories are shown to separate into smooth, reversible elastic branches which are intermittently broken by discrete catastrophic plastic events. The onset of a typical plastic event is studied with precision, and it is shown that the mode of the system which is responsible for the loss of stability has structure in real space which is consistent with a quadrupolar source acting on an elastic matrix. The plastic events themselves are shown to be composed of localized shear transformations which organize into lines of slip which span the length of the simulation cell, and a mechanism for the organization is discussed. Although within a single event there are strong spatial correlations in the deformation, we find little correlation from one event to the next, and these transient lines of slip are not to be confounded with the persistent regions of localized shear -- so-called "shear bands" -- found in related studies. The slip lines gives rise to particular scalings with system length of various measures of event size. Strikingly, data obtained using three differing interaction potentials can be brought into quantitative agreement after a simple rescaling, emphasizing the insensitivity of the emergent plastic behavior in these disordered systems to the precise details of the underlying interactions. The results should be relevant to understanding plastic deformation in systems such as metallic glasses well below their glass temperature, soft glassy systems (such as dense emulsions), or compressed granular materials.Comment: 21 pages, 18 figure

Lumpy Structures in Self-Gravitating Disks

Following Toomre & Kalnajs (1991), local models of slightly dissipative self-gravitating disks show how inhomogeneous structures can be maintained over several galaxy rotations. Their basic physical ingredients are self-gravity, dissipation and differential rotation. In order to explore the structures resulting from these processes on the kpc scale, local simulation of self-gravitating disks are performed in this paper in 2D as well as in 3D. The third dimension becomes a priori important as soon as matter clumping causes a tight coupling of the 3D equations of motion. The physically simple and general framework of the model permits to make conclusions beyond the here considered scales. A time dependent affine coordinate system is used, allowing to calculate the gravitational forces via a particle-mesh FFT-method, increasing the performance with respect to previous direct force calculations. Persistent patterns, formed by transient structures, whose intensity and morphological characteristic depend on the dissipation rate are obtained and described. Some of our simulations reveal first signs of mass-size and velocity dispersion-size power-law relations, but a clear scale invariant behavior will require more powerful computer techniques.Comment: 28 pages, 32 figures. Accepted for publication in A&A. Full resolution paper available at http://obswww.unige.ch/Preprints/dyn_art.htm

Dense Packings of Superdisks and the Role of Symmetry

We construct the densest known two-dimensional packings of superdisks in the plane whose shapes are defined by |x^(2p) + y^(2p)| <= 1, which contains both convex-shaped particles (p > 0.5, with the circular-disk case p = 1) and concave-shaped particles (0 < p < 0.5). The packings of the convex cases with p 1 generated by a recently developed event-driven molecular dynamics (MD) simulation algorithm [Donev, Torquato and Stillinger, J. Comput. Phys. 202 (2005) 737] suggest exact constructions of the densest known packings. We find that the packing density (covering fraction of the particles) increases dramatically as the particle shape moves away from the "circular-disk" point (p = 1). In particular, we find that the maximal packing densities of superdisks for certain p 6 = 1 are achieved by one of the two families of Bravais lattice packings, which provides additional numerical evidence for Minkowski's conjecture concerning the critical determinant of the region occupied by a superdisk. Moreover, our analysis on the generated packings reveals that the broken rotational symmetry of superdisks influences the packing characteristics in a non-trivial way. We also propose an analytical method to construct dense packings of concave superdisks based on our observations of the structural properties of packings of convex superdisks.Comment: 15 pages, 8 figure