Search

We characterize a system of hard spheres with a simple collision rule that breaks time reversal symmetry, but conserves energy. The collisions lead to an a-chiral, isotropic, and homogeneous stationary state, whose properties are determined in simulations and compared to an approximate theory originally developed for elastic hard spheres. In the nonequilibrium fluid state, velocities are correlated, a phenomenon known from other nonequilibrium stationary states. The correlations are long-ranged decaying like $1/r^d$ in $d$ dimensions. Such correlations are expected on general grounds far from equilibrium and had previously been observed in driven or non-stationary systems., 14 pages, 17 Figures

It has recently been reported that bacteria, such as E.coli and P. putida, perform distinct modes of motion when placed in porous media as compared to dilute regions or free space. This has led us to suggest an efficient strategy for active particles in a disordered environment: reorientations are suppressed in locally dilute regions and intensified in locally dense ones. Thereby the local geometry determines the optimal path of the active agent and substantially accelerates the dynamics for up to two orders of magnitude. We observe a non-monotonic behavior of the diffusion coefficient in dependence on the tumbling rate and identify a localisation transition, either by increasing the density of obstacles or by decreasing the reorientation rate., Comment: 5 pages, 5 figure. Movies can be accessed via: phi_obs0.4-sigmoid.webm: see https://nxcloud.omid.land/s/iDkDYGHwjczpZg9 phi_obs0.6-sigmoid.webm: see https://nxcloud.omid.land/s/eJQo8btyWQZ9nrA phi_obs0.4-uniform.webm: see https://nxcloud.omid.land/s/Fnp7cpGnm8JYqxX phi_obs0.6-uniform.webm: see https://nxcloud.omid.land/s/cBMgtL5QfYiXWZD

A liquid droplet, immersed into a Newtonian fluid, can be propelled solely by internal flow. In a simple model, this flow is generated by a collection of point forces, which represent externally actuated devices or model autonomous swimmers. We work out the general framework to compute the self-propulsion of the droplet as a function of the actuating forces and their positions within the droplet. A single point force, F with general orientation and position, r_0, gives rise to both, translational and rotational motion of the droplet. We show that the translational mobility is anisotropic and the rotational mobility can be nonmonotonic as a function of | r_0|, depending on the viscosity contrast. Due to the linearity of the Stokes equation, superposition can be used to discuss more complex arrays of point forces. We analyse force dipoles, such as a stresslet, a simple model of a biflagellate swimmer and a rotlet, representing a helical swimmer, driven by an external magnetic field. For a general force distribution with arbitrary high multipole moments the propulsion properties of the droplet depend only on a few low order multipoles: up to the quadrupole for translational and up to a special octopole for rotational motion. The coupled motion of droplet and device is discussed for a few exemplary cases. We show in particular that a biflagellate swimmer, modeled as a stresslet, achieves a steady comoving state, where the position of the device relative to the droplet remains fixed. In fact there are two fixpoints, symmetric with respect to the center of the droplet. A tiny external force selects one of them and allows to switch between forward and backward motion., Comment: arXiv admin note: text overlap with arXiv:1902.01159

Abstract We work out the propulsion of a viscous drop which is driven by two mechanisms: the active velocity of an encapsulated squirmer and an externally applied force acting on the squirmer. Of particular interest is the existence of a stable comoving state of drop and squirmer, allowing for controlled manipulation of the viscous drop by external forcing. The velocities of droplet and squirmer, as well as the conditions for a stable comoving state are worked out analytically for the axisymmetric configuration with a general displacement of the squirmer from the center of the droplet Graphic abstract

Dense granular media can be prepared in a stationary state by suitable driving. Such driving can be given by a random, momentum-conserving external force acting upon, say, a fluid comprised of inelastic hard spheres. While this out-of-equilibrium stationary state violates time reversal symmetry, it can still be investigated by means similar to ordinary fluids. For high enough density, the driven granular fluid undergoes a glass transition, and for this transition an extension to the mode-coupling theory can be derived. In addition to the quiescent stationary state, a kinetic theory as well as experiments in 2D for the active microrheology can be devised, where a selected intruder is pulled through the system as a probe for either constant velocity or force.

We consider self-propelled droplets which are driven by internal flow. Tracer particles, which are advected by the flow, in general follow chaotic trajectories, even though the motion of the autonomous swimmer is completely regular. The flow is mixing, and for P\'{e}clet and Batchelor numbers, which are realized e.g. in eucaryotic cells, advective mixing can substantially accelerate and even dominate transport by diffusion.

published

We analyze the flow curves of a two-dimensional assembly of granular particles which are interacting via frictional contact forces. For packing fractions slightly below jamming, the fluid undergoes a large scale instability, implying a range of stress and strain rates where no stationary flow can exist. Whereas small systems were shown previously to exhibit hysteretic jumps between the low and high stress branches, large systems exhibit continuous shear thickening arising from averaging unsteady, spatially heterogeneous flows. The observed large scale patterns as well as their dynamics are found to depend on strain rate: At the lower end of the unstable region, force chains merge to form giant bands that span the system in the compressional direction and propagate in the dilational direction. At the upper end, we observe large scale clusters which extend along the dilational direction and propagate along the compressional direction. Both patterns, bands and clusters, come in with infinite correlation length similar to the sudden onset of system-spanning plugs in impact experiments.

We suggest a simple model for reversible cross-links, binding and unbinding to/from a network of semiflexible polymers. The resulting frequency dependent response of the network to an applied shear is calculated via Brownian dynamics simulations. It is shown to be rather complex with the timescale of the linkers competing with the excitations of the network. If the lifetime of the linkers is the longest timescale, as is indeed the case in most biological networks, then a distinct low frequency peak of the loss modulus develops. The storage modulus shows a corresponding decay from its plateau value, which for irreversible cross-linkers extends all the way to the static limit. This additional relaxation mechanism can be controlled by the relative weight of reversible and irreversible linkers., Comment: accepted J. Chem. Phys

We discuss the flow field and propulsion velocity of active droplets, which are driven by body forces residing on a rigid gel. The latter is modelled as a porous medium which gives rise to permeation forces. In the simplest model, the Brinkman equation, the porous medium is characterised by a single lengthscale [Formula: see text] --the square root of the permeability. We compute the flow fields inside and outside of the droplet as well as the energy dissipation as a function of [Formula: see text]. We furthermore show that there are optimal gel fractions, giving rise to maximal linear and rotational velocities. In the limit [Formula: see text], corresponding to a very dilute gel, we recover Stokes flow. The opposite limit, [Formula: see text], corresponding to a space filling gel, is singular and not equivalent to Darcy's equation, which cannot account for self-propulsion.

A two-dimensional bidisperse granular fluid is shown to exhibit pronounced long-ranged dynamical heterogeneities as dynamical arrest is approached. Here we focus on the most direct approach to study these heterogeneities: we identify clusters of slow particles and determine their size, $N_c$, and their radius of gyration, $R_G$. We show that $N_c\propto R_G^{d_f}$, providing direct evidence that the most immobile particles arrange in fractal objects with a fractal dimension, $d_f$, that is observed to increase with packing fraction $\phi$. The cluster size distribution obeys scaling, approaching an algebraic decay in the limit of structural arrest, i.e., $\phi\to\phi_c$. Alternatively, dynamical heterogeneities are analyzed via the four-point structure factor $S_4(q,t)$ and the dynamical susceptibility $\chi_4(t)$. $S_4(q,t)$ is shown to obey scaling in the full range of packing fractions, $0.6\leq\phi\leq 0.805$, and to become increasingly long-ranged as $\phi\to\phi_c$. Finite size scaling of $\chi_4(t)$ provides a consistency check for the previously analyzed divergences of $\chi_4(t)\propto (\phi-\phi_c)^{-\gamma_{\chi}}$ and the correlation length $\xi\propto (\phi-\phi_c)^{-\gamma_{\xi}}$. We check the robustness of our results with respect to our definition of mobility. The divergences and the scaling for $\phi\to\phi_c$ suggest a non-equilibrium glass transition which seems qualitatively independent of the coefficient of restitution., Comment: 14 pages, 25 figures

The actin and microtubule networks form the dynamic cytoskeleton. Network dynamics is driven by molecular motors applying force onto the networks and the interactions between the networks. Here we assay the dynamics of centrosomes in the scale of seconds as a proxy for the movement of microtubule asters. With this assay we want to detect the role of specific motors and of network interaction. During interphase of syncytial embryos of Drosophila, cortical actin and the microtubule network depend on each other. Centrosomes induce cortical actin to form caps, whereas F-actin anchors microtubules to the cortex. In addition, lateral interactions between microtubule asters are assumed to be important for regular spatial organization of the syncytial embryo. The functional interaction between the microtubule asters and cortical actin has been largely analyzed in a static manner, so far. We recorded the movement of centrosomes at 1 Hz and analyzed their fluctuations for two processes—pair separation and individual movement. We found that F-actin is required for directional movements during initial centrosome pair separation, because separation proceeds in a diffusive manner in latrunculin-injected embryos. For assaying individual movement, we established a fluctuation parameter as the deviation from temporally and spatially slowly varying drift movements. By analysis of mutant and drug-injected embryos, we found that the fluctuations were suppressed by both cortical actin and microtubules. Surprisingly, the microtubule motor Kinesin-1 also suppressed fluctuations to a similar degree as F-actin. Kinesin-1 may mediate linkage of the microtubule (+)-ends to the actin cortex. Consistent with this model is our finding that Kinesin-1-GFP accumulates at the cortical actin caps.

A theory for the non-local shear stress correlations in supercooled liquids is derived from first principles. It captures the crossover from viscous to elastic dynamics at an idealized liquid to glass transition and explains the emergence of long-ranged stress correlations in glass, as expected from classical continuum elasticity. The long-ranged stress correlations can be traced to the coupling of shear stress to transverse momentum, which is ignored in the classic Maxwell model. To rescue this widely used model, we suggest a generalization in terms of a single relaxation time $\tau$ for the fast degrees of freedom only. This generalized Maxwell model implies a divergent correlation length $\xi\propto\tau$ as well as dynamic critical scaling and correctly accounts for the far-field stress correlations. It can be rephrased in terms of generalized hydrodynamic equations, which naturally couple stress and momentum and furthermore allow to connect to fluidity and elasto-plastic models., Comment: 5 pages plus 1 page SI

We compute the rheological properties of inelastic hard spheres in steady shear flow for general shear rates and densities. Starting from the microscopic dynamics we generalise the Integration Through Transients (\textsc{itt}) formalism to a fluid of dissipative, randomly driven granular particles. The stress relaxation function is computed approximately within a mode-coupling theory---based on the physical picture, that relaxation of shear is dominated by slow structural relaxation, as the glass transition is approached. The transient build-up of stress in steady shear is thus traced back to transient density correlations which are computed self-consistently within mode-coupling theory. The glass transition is signalled by the appearance of a yield stress and a divergence of the Newtonian viscosity, characterizing linear response. For shear rates comparable to the structural relaxation time, the stress becomes independent of shear rate and we observe shear thinning, while for the largest shear rates Bagnold scaling, i.e., a quadratic increase of shear stress with shear rate, is recovered. The rheological properties are qualitatively similar for all values of $\varepsilon$, the coefficient of restitution; however, the magnitude of the stress as well as the range of shear thinning and thickening show significant dependence on the inelasticity., Comment: 37 pages, 7 figures

Considering a granular fluid of inelastic smooth hard spheres we discuss the conditions delineating the rheological regimes comprising Newtonian, Bagnoldian, shear thinning, and shear thickening behavior. Developing a kinetic theory, valid at finite shear rates and densities around the glass transition density, we predict the viscosity and Bagnold coefficient at practically relevant values of the control parameters. The determination of full flow curves relating the shear stress $\sigma$ to the shear rate $\dot\gamma$, and predictions of the yield stress complete our discussion of granular rheology derived from first principles., Comment: 7 pages, minor fixes and extended appendix

We study the self-propulsion of spherical droplets as simplified hydrodynamic models of swimming micro-organisms or artificial micro-swimmers. In contrast to approaches that start from active velocity fields produced by the system, we consider active interface tractions, body force densities and active stresses as the origin of autonomous swimming. For negligible Reynolds number and given activity, we compute the external and internal flow fields as well as the centre of mass velocity and angular velocity of the droplet at fixed time. To construct trajectories from single time snapshots, the evolution of active forces or stresses must be determined in the laboratory frame. Here, we consider the case of active matter, which is carried by a continuously distributed rigid but sparse (cyto)-skeleton that is immersed in the droplet interior. We calculate examples of trajectories of a droplet and its skeleton from force densities or stresses, which may be explicitly time-dependent in a frame fixed within the skeleton.

We develop a generalized hydrodynamic theory, which can account for the build-up of long-ranged and long-lived shear stress correlations in supercooled liquids as the glass transition is approached. Our theory is based on the decomposition of tensorial stress relaxation into fast microscopic processes and slow dynamics due to conservation laws. In the fluid, anisotropic shear stress correlations arise from the tensorial nature of stress. By approximating the fast microscopic processes by a single relaxation time in the spirit of Maxwell, we find viscoelastic precursors of the Eshelby-type correlations familiar in an elastic medium. The spatial extent of shear stress fluctuations is characterized by a correlation length ξ which grows like the viscosity η or time scale τ ∼ η, whose divergence signals the glass transition. In the solid, the correlation length is infinite and stress correlations decay algebraically as r−d in d dimensions. published

In a granular gas of rough particles the spin of a grain is correlated with its linear velocity. We develop an analytical theory to account for these correlations and compare its predictions to numerical simulations, using Direct Simulation Monte Carlo as well as Molecular Dynamics. The system is shown to relax from an arbitrary initial state to a steady-state, which is characterized by time-independent, finite correlations of spin and linear velocity. The latter are analyzed systematically for a wide range of system parameters, including the coefficients of tangential and normal restitution as well as the moment of inertia of the particles. For most parameter values the axis of rotation and the direction of linear momentum are perpendicular like in a sliced tennis ball, while parallel orientation, like in a rifled bullet, occurs only for a small range of parameters. The limit of smooth spheres is singular: any arbitrarily small roughness unavoidably causes significant translation-rotation correlations, whereas for perfectly smooth spheres the rotational degrees of freedom are completely decoupled from the dynamic evolution of the gas. peerReviewed

We analyze the stretching elasticity of a wormlike chain with a tension discontinuity resulting from a Hookean spring connecting its backbone to a fixed point. The elasticity of isolated semiflexible filaments has been the subject in a significant body of literature, primarily because of its relevance to the mechanics of biological matter. In real systems, however, these filaments are usually part of supramolecular structures involving cross-linkers or molecular motors which cause tension discontinuities. Our model is intended as a minimal structural element incorporating such a discontinuity. We obtain analytical results in the weakly bending limit of the filament, concerning its force-extension relation and the response of the two parts in which the filament is divided by the spring. For a small tension discontinuity, the linear response of the filament extension to this discontinuity strongly depends on the external tension. For large external tension $f$, the spring force contributes a subdominant correction $\sim 1/f^{3/2}$ to the well known $\sim 1/\sqrt{f}$ dependence of the end-to-end extension.

Spider silk is well-known for its outstanding mechanical properties. However, there is a significant variation of these properties in literature and studies analyzing large numbers of silk samples to explain these variations are still lacking. To fill this gap, the following work examines the mechanical properties of major ampullate silk based on a large ensemble of threads from Nephila clavipes and Nephila senegalensis. In addition, the effect of relative humidity (RH) on the mechanical properties was quantified. The large effect of RH on the mechanical properties makes it plausible that the variation in the literature values can to a large extent be attributed to changes in RH. Spider silk’s most remarkable property—its high tenacity—remains unchanged. In addition, this work also includes hysteresis as well as relaxation measurements. It is found that the relaxation process is well described by a stretched exponential decay.

We study the kinetic theory of driven and undriven granular gases, taking into account both translational and rotational degrees of freedom. We obtain the high-energy tail of the stationary bivariate energy distribution, depending on the total energy E and the ratio \(x=\sqrt{E_{w}/E}\) of rotational energy Ew to total energy. Extremely energetic particles have a unique and well-defined distribution f(x) which has several remarkable features: x is not uniformly distributed as in molecular gases; f(x) is not smooth but has multiple singularities. The latter behavior is sensitive to material properties such as the collision parameters, the moment of inertia and the collision rate. Interestingly, there are preferred ratios of rotational-to-total energy. In general, f(x) is strongly correlated with energy and the deviations from a uniform distribution grow with energy. We also solve for the energy distribution of freely cooling Maxwell Molecules and find qualitatively similar behavior.

Cells moving on a two dimensional substrate generate motion by polymerizing actin filament networks inside a flat membrane protrusion. New filaments are generated by branching off existing ones, giving rise to branched network structures. We investigate the force-extension relation of branched filaments, grafted on an elastic structure at one end and pushing with the free ends against the leading edge cell membrane. Single filaments are modeled as worm-like chains, whose thermal bending fluctuations are restricted by the leading edge cell membrane, resulting in an effective force. Branching can increase the stiffness considerably; however the effect depends on branch point position and filament orientation, being most pronounced for intermediate tilt angles and intermediate branch point positions. We describe filament networks without cross-linkers to focus on the effect of branching. We use randomly positioned branch points, as generated in the process of treadmilling, and orientation distributions as measured in lamellipodia. These networks reproduce both the weak and strong force response of lamellipodia as measured in force-velocity experiments. We compare properties of branched and unbranched networks. The ratio of the network average of the force per branched filament to the average force per unbranched filament depends on the orientation distribution of the filaments. The ratio exhibits compression dependence and may go up to about 4.5 in networks with a narrow orientation distribution. With orientation distributions measured in lamellipodia, it is about two and essentially independent from network compression, graft elasticity and filament persistence length.

The dynamics of randomly crosslinked liquids is addressed via a Rouse- and a Zimm-type model with crosslink statistics taken either from bond percolation or Erdoes-Renyi random graphs. While the Rouse-type model isolates the effects of the random connectivity on the dynamics of molecular clusters, the Zimm-type model also accounts for hydrodynamic interactions on a preaveraged level. The incoherent intermediate scattering function is computed in thermal equilibrium, its critical behaviour near the sol-gel transition is analysed and related to the scaling of cluster diffusion constants at the critical point. Second, non-equilibrium dynamics is studied by looking at stress relaxation in a simple shear flow. Anomalous stress relaxation and critical rheological properties are derived. Some of the results contradict long-standing scaling arguments, which are shown to be flawed by inconsistencies., Comment: 21 pages, 3 figures; Dedicated to Lothar Schaefer on the occasion of his 60th birthday; Changes: added comments on the gel phase and some references

Starting from a Zimm model we study selfdiffusion in a solution of crosslinked monomers. We focus on the effects of the hydrodynamic interaction on the dynamics and the critical behaviour at the sol-gel-point. Hydrodynamic interactions cause the clusters' diffusion constant to depend not only on the cluster's size but also on the cluster's shape -- in contrast to the Rouse model. This gives rise to a nontrivial scaling of the Kirkwood diffusion constant averaged over all clusters of fixed size $n$, ${\hat D}_n\sim n^{-{\hat b}}$ with ${\hat b}=1/d_s-1/2$ given in terms of the spectral dimension $d_s$ of critical percolation clusters. The long-time decay of the incoherent scattering function is determined by the diffusive motion of the largest clusters. This implies the critical vanishing $D_{\rm eff}\sim \epsilon^a$ of the cluster-averaged effective diffusion constant at the gel point with exponent $a = (3/2 -\tau +1/d_{s})/\sigma$., Comment: 7 pages, 4 figures; slightly expanded version, typos corrected; to appear in Eur. Phys. J. E

We discuss a model of presynaptic vesicle dynamics, which allows for heterogeneity in release probability among vesicles. Specifically, we explore the possibility that synaptic activity is carried by two types of vesicles; first, a readily releasable pool and, second, a reluctantly releasable pool. The pools differ regarding their probability of release and time scales on which released vesicles are replaced by new ones. Vesicles of both pools increase their release probability during repetitive stimulation according to the buildup of Ca2+ concentration in the terminal. These properties are modeled to fit data from the calyx of Held, a giant synapse in the auditory pathway. We demonstrate that this arrangement of two pools of releasable vesicles can account for a variety of experimentally observed patterns of synaptic depression and facilitation at this synapse. We conclude that synaptic transmission cannot be accurately described unless heterogeneity of synaptic release probability is taken into account.

A two-dimensional dense fluid of frictional grains is shown to exhibit time-chaotic, spatially heterogeneous flow in a range of stress values, $\sigma$, chosen in the unstable region of s-shaped flow curves. Stress controlled simulations reveal a phase diagram with reentrant stationary flow for small and large stress $\sigma$. In between no steady flow state can be reached, instead the system either jams or displays time dependent heterogeneous strain rates $\dot\gamma({\bf r},t)$. The results of simulations are in agreement with the stability analysis of a simple hydrodynamic model, coupling stress and microstructure which we tentatively associate with the frictional contact network.

We consider active particles in a heterogeneous medium, modeled by static, random obstacles. In accordance with the known tendency of active particles to cluster, we observe accumulation and crystallization of active particles around the obstacles which serve as nucleation sites. In the limit of high activity, the crystals start to rotate spontaneously, resembling a rotating rigid body. We trace the occurrence of these oscillations to the enhanced attraction of particles whose orientation points along the rotational velocity as compared to those whose orientation points in the opposite direction.

Molecular dynamics simulations are used to generate an ensemble of saddles of the potential energy of a Lennard-Jones liquid. Classifying all extrema by their potential energy u and number of unstable directions k, a well defined relation k(u) is revealed. The degree of instability of typical stationary points vanishes at a threshold potential energy, which lies above the energy of the lowest glassy minima of the system. The energies of the inherent states, as obtained by the Stillinger-Weber method, approach the threshold energy at a temperature close to the mode-coupling transition temperature Tc., Comment: 4 RevTeX pages, 6 eps figures. Revised version

Recently, we have modelled inelastic collisions of one-dimensional rods [1,2] by the absorption of translational energy E tr through internal degrees of freedom, in particular elastic vibrations. We arrived at a stochastic description of collision processes, characterised by a stochastic coefficient of restitution e . In this paper, we construct an analytic approximation for the transition probability E tr → E tr ′=(1− e 2 ) E tr . This allows us to perform much longer simulations of large, strongly inelastic granular systems and study relaxation to the true equilibrium state. If the internal vibrations are undamped, equilibrium is characterised by propagating sound waves. In the case of damping, the system develops towards a final state which consists of one big cluster, containing all particles at rest.

We investigate a binary Lennard-Jones mixture with molecular dynamics simulations. We consider first a system cooled linearly in time with the cooling rate gamma. By varying gamma over almost four decades we study the influence of the cooling rate on the glass transition and on the resulting glass. We find for all investigated quantities a cooling rate dependence; with decreasing cooling rate the system falls out of equilibrium at decreasing temperatures, reaches lower enthalpies and obtains increasing local order. Next we study the dynamics of the melting process by investigating the most immobile and most mobile particles in the glass. We find that their spatial distribution is heterogeneous and that the immobile/mobile particles are surrounded by denser/less dense cages than an average particle., 9 pages, latex, 10 figures

We investigate the static shear viscosity on the sol side of the vulcanization transition within a minimal mesoscopic model for the Rouse-dynamics of a randomly crosslinked melt of phantom polymers. We derive an exact relation between the viscosity and the resistances measured in a corresponding random resistor network. This enables us to calculate the viscosity exactly for an ensemble of crosslinks without correlations. The viscosity diverges logarithmically as the critical point is approached. For a more realistic ensemble of crosslinks amenable to the scaling description of percolation, we prove the scaling relation $k=\phi-\beta$ between the critical exponent $k$ of the viscosity, the thermal exponent $\beta$ associated with the gel fraction and the crossover exponent $\phi$ of a random resistor network., Comment: 8 pages, uses Europhysics Letters style; Revisions: results extended

A detailed mathematical analysis of the diffusion process of neurotransmitter inside the synaptic cleft is presented and the spatio-temporal concentration profile is calculated. Using information about the experimentally observed time course of glutamate in the cleft the effective diffusion coefficient Dnet is estimated as Dnet approximately 20-50 nm(2) microseconds(-1), implying a strong reduction compared with free diffusion in aqueous solution. The tortuosity of the cleft and interactions with transporter molecules are assumed to affect the transmitter motion. We estimate the transporter density to be 5170 to 8900 micrometer(-2) in the synaptic cleft and its vicinity, using the experimentally observed time constant of glutamate. Furthermore a theoretical model of synaptic transmission is presented, taking the spatial distribution of post-synaptic (AMPA-) receptors into account. The transmitter diffusion and receptor dynamics are modeled by Monte Carlo simulations preserving the typically observed noisy character of post-synaptic responses. Distributions of amplitudes, rise and decay times are calculated and shown to agree well with experiments. Average open probabilities are computed from a novel kinetic model and are shown to agree with averages over many Monte Carlo runs. Our results suggest that post-synaptic currents are only weakly potentiated by clustering of post-synaptic receptors, but increase linearly with the total number of receptors. Distributions of amplitudes and rise times are used to discriminate between different morphologies, e.g. simple and perforated synapses. A skew in the miniature amplitude distribution can be caused by multiple release of pre-synaptic vesicles at perforated synapses.

We calculate the coefficient of restitution, $\epsilon$, starting from a microscopic model of elastic disks. The theory is shown to agree with the approach of Hertz in the quasistatic limit, but predicts inelastic collisions for finite relative velocities of two approaching disks. The velocity dependence of $\epsilon$ is calculated numerically for a wide range of velocities. The coefficient of restitution furthermore depends on the elastic constants of the material via Poisson's number. The elastic vibrations absorb kinetic energy more effectively for materials with low values of the shear modulus., Comment: 25 pages, 12 Postscript figures, LaTex2e

We investigate freely cooling systems of rough spheres in two and three dimensions. Simulations using an event driven algorithm are compared with results of an approximate kinetic theory, based on the assumption of a generalized homogeneous cooling state. For short times $t$, translational and rotational energy are found to change linearly with $t$. For large times both energies decay like $t^{-2}$ with a ratio independent of time, but not corresponding to equipartition. Good agreement is found between theory and simulations, as long as no clustering instability is observed. System parameters, i.e. density, particle size, and particle mass can be absorbed in a rescaled time, so that the decay of translational and rotational energy is solely determined by normal restitution and surface roughness., Comment: 10 pages, 10 eps-figures

We present a minimal kinetic model for excitatory synaptic transmission to cerebellar Purkinje cells. The main components are a kinetic model for a single glutamate receptor, which is calibrated with the help of patch clamp data, and a mean field approximation for the dynamics of a population of channels, which generate an EPSC. The resulting minimal model of the parallel fiber–Purkinje cell synapse is used to estimate the dynamics of glutamate in the synaptic cleft and to clarify the role of receptor desensitization in synaptic transmission. We also apply the model to different aspects of synaptic modulation, like long-term depression and potentiation by pharmacological application of ampakines. In the framework of the minimal model these effects can be understood as the result of modified receptor kinetics.

Amorphous alloys can be modelled as random networks of atoms or molecules of differing valencies. A statistical-mechanical theory of such networks is presented for the case of binary amorphous alloys. The existence of a continuous equilibrium phase transition, from the liquid state to the amorphous solid state driven by increasing the density of permanent random covalent bonds, is demonstrated. The structural and elastic properties of the amorphous solid state are calculated through the use of a variational Ansatz, and a discussion is given of the dependence of these properties on the relative concentration of the constituent atoms or molecules and their valencies. Implications of the present results in the context of gelation are addressed, especially with regard to the (primarily entropic) elastic properties of gels.

We have formulated a model for the coupled dynamics of the local polarization and the elastic deformations of dipolar glasses. Randomly interacting dipoles undergo a phase transition to a glassy state, which is signalled above Tc by a divergence of the relaxation time of the dipolar autocorrelation. Dipolar dynamics is directly observable in the dynamic dielectric susceptibility, for which we predict a power law dependence on frequency at and close to Tc. The slowing down of dipolar degrees of freedom also gives rise to dynamic critical anomalies in the elastic and structural properties, like the appearance of a central peak and divergence of the sound damping. All these phenomena are a consequence of dynamic scaling of the dipolar autocorrelation and hence interrelated. We discuss two scaling functions in detail: One is calculated within mean field theory and the other one is taken from numerical simulations of short range Ising spin glasses. The predictions of the model are in accordance with most experimental findings for dipolar glasses.

We generalize the statistical mechanical theory of vulcanization to the case of D-dimensional polymerized manifolds. Starting from a continuum model of self-avoiding manifolds, we study the effects of introducing random crosslinks between monomers on (in general) different manifolds. As for the case of linear polymers, one observes a continuous phase transition from a fluid to an amorphous solid state, characterized by a finite fraction of localized monomers. We compute this fraction, as well as the typical localization length near the transition.

When pulling a particle in a driven granular fluid with constant force $F_{ex}$, the probe particle approaches a steady-state average velocity $v$. This velocity and the corresponding friction coefficient of the probe $\zeta=F_{ex}/v$ are obtained within a schematic model of mode-coupling theory and compared to results from event-driven simulations. For small and moderate drag forces, the model describes the simulation results successfully for both the linear as well as the nonlinear region: The linear response regime (constant friction) for small drag forces is followed by shear thinning (decreasing friction) for moderate forces. For large forces, the model demonstrates a subsequent increasing friction in qualitative agreement with the data. The square-root increase of the friction with force found in [Fiege et al., Granular Matter $\boldsymbol{14}$, 247 (2012)] is explained by a simple kinetic theory., Comment: 5 pages, 4 figures

Large scale simulations of two-dimensional bidisperse granular fluids allow us to determine spatial correlations of slow particles via the four-point structure factor $S_4(q,t)$. Both cases, elastic ($\varepsilon=1$) as well as inelastic ($\varepsilon < 1$) collisions, are studied. As the fluid approaches structural arrest, i.e. for packing fractions in the range $0.6 \le \phi \le 0.805$, scaling is shown to hold: $S_4(q,t)/\chi_4(t)=s(q\xi(t))$. Both the dynamic susceptibility, $\chi_4(\tau_{\alpha})$, as well as the dynamic correlation length, $\xi(\tau_{\alpha})$, evaluated at the $\alpha$ relaxation time, $\tau_{\alpha}$, can be fitted to a power law divergence at a critical packing fraction. The measured $\xi(\tau_{\alpha})$ widely exceeds the largest one previously observed for hard sphere 3d fluids. The number of particles in a slow cluster and the correlation length are related by a robust power law, $\chi_4(\tau_{\alpha}) \approx\xi^{d-p}(\tau_{\alpha})$, with an exponent $d-p\approx 1.6$. This scaling is remarkably independent of $\varepsilon$, even though the strength of the dynamical heterogeneity increases dramatically as $\varepsilon$ grows., Comment: 5 pages, 6 figures

We explore the effect of an attractive interaction between parallel-aligned polymers, which are perpendicularly grafted on a substrate. Such an attractive interaction could be due to, e.g., reversible cross-links. The competition between permanent grafting favoring a homogeneous state of the polymer brush and the attraction, which tends to induce in-plane collapse of the aligned polymers, gives rise to an instability of the homogeneous phase to a bundled state. In this latter state the in-plane translational symmetry is spontaneously broken and the density is modulated with a finite wavelength, which is set by the length scale of transverse fluctuations of the grafted polymers. We analyze the instability for two models of aligned polymers: directed polymers with a line tension and weakly bending chains with a bending stiffness., 7 pages, 5 figures, final version as published in PRE

We investigate the long time dynamics of a strong glass former, SiO2, below the glass transition temperature by averaging single particle trajectories over time windows which comprise roughly 100 particle oscillations. The structure on this coarse-grained time scale is very well defined in terms of coordination numbers, allowing us to identify ill-coordinated atoms, called defects in the following. The most numerous defects are OO neighbors, whose lifetimes are comparable to the equilibration time at low temperature. On the other hand SiO and OSi defects are very rare and short lived. The lifetime of defects is found to be strongly temperature dependent, consistent with activated processes. Single-particle jumps give rise to local structural rearrangements. We show that in SiO2 these structural rearrangements are coupled to the creation or annihilation of defects, giving rise to very strong correlations of jumping atoms and defects., 10 pages, 17 figures

Aiming at the mechanical properties of cross-linked biopolymers, we set up and analyze a model of two weakly bending wormlike chains subjected to a tensile force, with regularly spaced inter-chain bonds (cross-links) represented by harmonic springs. Within this model, we compute the force-extension curve and the differential stiffness exactly and discuss several limiting cases. Cross-links effectively stiffen the chain pair by reducing thermal fluctuations transverse to the force and alignment direction. The extra alignment due to cross-links increases both with growing number and with growing strength of the cross-links, and is most prominent for small force f. For large f, the additional, cross-link-induced extension is subdominant except for the case of linking the chains rigidly and continuously along their contour. In this combined limit, we recover asymptotically the elasticity of a weakly bending wormlike chain without constraints, stiffened by a factor four. The increase in differential stiffness can be as large as 100% for small f or large numbers of cross-links., Comment: 11 pages, 6 figures, submitted to PRE

Large scale simulations and analytical theory have been combined to obtain the non-equilibrium velocity distribution, $f(v)$, of randomly accelerated particles in suspension. The simulations are based on an event-driven algorithm, generalised to include friction. They reveal strongly anomalous but largely universal distributions which are independent of volume fraction and collision processes, which suggests a one-particle model should capture all the essential features. We have formulated this one-particle model and solved it analytically in the limit of strong damping, where we find that $f(v)$ decays as $1/v$ for multiple decades, eventually crossing over to a Gaussian decay for the largest velocities. Many particle simulations and numerical solution of the one-particle model agree for all values of the damping., 6 pages, 5 figures

Starting from a microscopic model of interacting elastic rods, we derive an effective dissipative dynamics for the translational motion, where the internal degrees of freedom, i.e., the vibrations of the rods, are represented by a thermalized bath. Focusing on two-particle collisions, we calculate the coefficient of restitution \ensuremath{\epsilon} as a function of the relative length of the colliding rods, the center of mass velocity, and the degree of excitation of the internal vibrations. In general, \ensuremath{\epsilon} is a stochastic quantity; its distribution is interpreted as the transition probability for the Markovian jump process, which describes the loss of translational energy in successive two-particle collisions. \textcopyright{} 1996 The American Physical Society.