ترغب بنشر مسار تعليمي؟ اضغط هنا

Event-driven Molecular Dynamics of Soft Particles

163   0   0.0 ( 0 )
 نشر من قبل Thorsten Poeschel
 تاريخ النشر 2013
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

The dynamics of dissipative soft-sphere gases obeys Newtons equation of motion which are commonly solved numerically by (force-based) Molecular Dynamics schemes. With the assumption of instantaneous, pairwise collisions, the simulation can be accelerated considerably using event-driven Molecular Dynamics, where the coefficient of restitution is derived from the interaction force between particles. Recently it was shown, however, that this approach may fail dramatically, that is, the obtained trajectories deviate significantly from the ones predicted by Newtons equations. In this paper, we generalize the concept of the coefficient of restitution and derive a numerical scheme which, in the case of dilute systems and frictionless interaction, allows us to perform highly efficient event-driven Molecular Dynamics simulations even for non-instantaneous collisions. We show that the particle trajectories predicted by the new scheme agree perfectly with the corresponding (force-based) Molecular Dynamics, except for a short transient period whose duration corresponds to the duration of the contact. Thus, the new algorithm solves Newtons equations of motion like force-based MD while preserving the advantages of event-driven simulations.



قيم البحث

اقرأ أيضاً

Single-file transport in pore-like structures constitute an important topic for both theory and experiment. For hardcore interacting particles, a good understanding of the collective dynamics has been achieved recently. Here we study how softness in the particle interaction affects the emergent transport behavior. To this end, we investigate driven Brownian motion of particles in a periodic potential. The particles interact via a repulsive softcore potential with a shape corresponding to a smoothed rectangular barrier. This shape allows us to elucidate effects of mutual particle penetration and particle crossing in a controlled manner. We find that even weak deviations from the hardcore case can have a strong impact on the particle current. Despite of this fact, the knowledge about the transport in a corresponding hardcore system is shown to be useful to describe and interpret our findings for the softcore case. This is achieved by assigning a thermodynamic effective size to the particles based on the equilibrium density functional of hard spheres.
We study, from first principles, the pressure exerted by an active fluid of spherical particles on general boundaries in two dimensions. We show that, despite the non-uniform pressure along curved walls, an equation of state is recovered upon a prope r spatial averaging. This holds even in the presence of pairwise interactions between particles or when asymmetric walls induce ratchet currents, which are accompanied by spontaneous shear stresses on the walls. For flexible obstacles, the pressure inhomogeneities lead to a modulational instability as well as to the spontaneous motion of short semi-flexible filaments. Finally, we relate the force exerted on objects immersed in active baths to the particle flux they generate around them.
116 - Dominik Lips , Artem Ryabov , 2020
Driven diffusive systems constitute paradigmatic models of nonequilibrium physics. Among them, a driven lattice gas known as the asymmetric simple exclusion process (ASEP) is the most prominent example for which many intriguing exact results have bee n obtained. After summarizing key findings, including the mapping of the ASEP to quantum spin chains, we discuss the recently introduced Brownian asymmetric simple exclusion process (BASEP) as a related class of driven diffusive system with continuous space dynamics. In the BASEP, driven Brownian motion of hardcore-interacting particles through one-dimensional periodic potentials is considered. We study whether current-density relations of the BASEP can be considered as generic for arbitrary periodic potentials and whether repulsive particle interactions other than hardcore lead to similar results. Our findings suggest that shapes of current-density relations are generic for single-well periodic potentials and can always be attributed to the interplay of a barrier reduction, blocking and exchange symmetry effect. This implies that in general up to five different phases of nonequilibrium steady states are possible for such potentials. The phases can occur in systems coupled to particle reservoirs, where the bulk density is the order parameter. For multiple-well periodic potentials, more complex current-density relations are possible and more phases can appear. Taking a repulsive Yukawa potential as an example, we show that the effects of barrier reduction and blocking on the current are also present. The exchange symmetry effect requires hardcore interactions and we demonstrate that it can still be identified when hardcore interactions are combined with weak Yukawa interactions.
Run-and-tumble dynamics is a wide-spread mechanism of swimming bacteria. The accumulation of run-and-tumble microswimmers near impermeable surfaces is studied theoretically and numerically in the low-density limit in two and three spatial dimensions. Both uni-modal and exponential distributions of the run lengths are considered. Constant run lengths lead to {peaks and depletions regions} in the density distribution of particles near the surface, in contrast to {exponentially-distributed run lengths}. Finally, we present a universal accumulation law for large channel widths, which applies not only to run-and-tumble swimmers, but also to many other kinds of self-propelled particles.
We obtain a kinetic description of spatially averaged dynamics of particle systems. Spatial averaging is one of the three types of averaging relevant within the Irwing-Kirkwood procedure (IKP), a general method for deriving macroscopic equations from molecular models. The other two types, ensemble averaging and time averaging, have been extensively studied, while spatial averaging is relatively less understood. We show that the average density, linear momentum, and kinetic energy used in IKP can be obtained from a single average quantity, called the generating function. A kinetic equation for the generating function is obtained and tested numerically on Lennard-Jones oscillator chains.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا