Search

This is a comment to a letter by D. Mandal, K. Klymko and M. R. DeWeese published as Phys. Rev. Lett. 119, 258001 (2017)., Comment: 2 pages without figures, in press as a Comment on Physical Review Letters

Self-propulsion (SP) is a main feature of active particles (AP), such as bacteria or biological micromotors, distinguishing them from passive colloids. A renowned consequence of SP is accumulation at static interfaces, even in the absence of hydrodynamic interactions. Here we address the role of SP in the interaction between AP and a moving semipermeable membrane. In particular, we implement a model of noninteracting AP in a channel crossed by a partially penetrable wall, moving at a constant velocity $c$. With respect to both the cases of passive colloids with $c>0$ and AP with $c=0$, the AP with finite $c$ show enhancement of accumulation in front of the obstacle and experience a largely increased drag force. This effect is understood in terms of an effective potential localised at the interface between particles and membrane, of height proportional to $c\tau/\xi$, where $\tau$ is the AP's re-orientation time and $\xi$ the width characterising the surface's smoothness ($\xi\to 0$ for hard core obstacles). An approximate analytical scheme is able to reproduce the observed density profiles and the measured drag force, in very good agreement with numerical simulations. The effects discussed here can be exploited for automatic selection and filtering of AP with desired parameters., Comment: 13 pages, 3 figures

We study the stationary state of a one-dimensional kinetic model where a probe particle is driven by an external field E and collides, elastically or inelastically, with a bath of particles at temperature T. We focus on the stationary distribution of the velocity of the particle, and of two estimates of the total entropy production \Delta s_tot. One is the entropy production of the medium \Delta s_m, which is equal to the energy exchanged with the scatterers, divided by a parameter \theta, coinciding with the particle temperature at E=0. The other is the work W done by the external field, again rescaled by \theta. At small E, a good collapse of the two distributions is found: in this case the two quantities also verify the Fluctuation Relation (FR), indicating that both are good approximations of \Delta s_tot. Differently, for large values of E, the fluctuations of W violate the FR, while \Delta s_m still verifies it., Comment: 6 pages, 4 figures

We study a general model of granular Brownian ratchet consisting of an asymmetric object moving on a line and surrounded by a two-dimensional granular gas, which in turn is coupled to an external random driving force. We discuss the two resulting Boltzmann equations describing the gas and the object in the dilute limit and obtain a closed system for the first few moments of the system velocity distributions. Predictions for the net ratchet drift, the variance of its velocity fluctuations and the transition rates in the Markovian limit, are compared to numerical simulations and a fair agreement is observed., Comment: 15 pages, 4 figures, to be published on Journal of Statistical Mechanics: Theory and Experiments

We consider a massive inelastic piston, whose opposite faces have different coefficients of restitution, moving under the action of an infinitely dilute gas of hard disks maintained at a fixed temperature. The dynamics of the piston is Markovian and obeys a continuous Master Equation: however, the asymmetry of restitution coefficients induces a violation of detailed balance and a net drift of the piston, as in a Brownian ratchet. Numerical investigations of such non-equilibrium stationary state show that the velocity fluctuations of the piston are symmetric around the mean value only in the limit of large piston mass, while they are strongly asymmetric in the opposite limit. Only taking into account such an asymmetry, i.e. including a third parameter in addition to the mean and the variance of the velocity distribution, it is possible to obtain a satisfactory analytical prediction for the ratchet drift velocity., Comment: 6 pages, 5 figures, to be published on Europhysics Letters; some references have been added

We present a lattice-based numerical method to describe the non equilibrium behavior of a simple fluid under non-uniform spatial conditions. The evolution equation for the one-particle phase-space distribution function is derived starting from a microscopic description of the system. It involves a series of approximations which are similar to those employed in theories of inhomogeneous fluids, such as Density Functional theory. Among the merits of the present approach: the possibility to determine the equation of state of the model, the transport coefficients and to provide an efficient method of numerical solution under non-uniform conditions. The algorithm is tested in a particular non equilibrium situation, namely the steady flow of a hard-sphere fluid across a narrow slit. Pronounced non-hydrodynamic oscillations in the density and velocity profiles are found.

We consider the velocity fluctuations of a system of particles described by the Inelastic Maxwell Model. The present work extends the methods, previously employed to obtain the one-particle velocity distribution function, to the study of the two particle correlations. Results regarding both the homogeneous cooling process and the steady state driven regime are presented. In particular we obtain the form of the pair correlation function in the scaling region of the homogeneous cooling process and show that some of its moments diverge. This fact has repercussions on the behavior of the energy fluctuations of the model., Comment: 16 pages, 1 figure, to be published on Journal of Statistical Mechanics: Theory and Experiments

In the present work we propose a method to determine fluctuation induced forces in non equilibrium systems. These forces are the analogue of the well known Casimir forces, which were originally introduced in Quantum Field theory and later extended to the area of Critical Phenomena. The procedure starts from the observation that many non equilibrium systems exhibit long-range correlations and the associated structure factors diverge in the long wavelength limit. The introduction of external bodies into such systems in general modifies the spectrum of these fluctuations and leads to the appearance of a net force between these bodies. The mechanism is illustrated by means of a simple example: a reaction diffusion equation with random noises., Comment: Submitted to Europhysics Letters. 7 pages, 2 figures

We have numerically investigated the behavior of driven non-cohesive granular media and found that two fixed large intruder particles, immersed in a sea of small particles, experience, in addition to a short range depletion force, a long range repulsive force. The observed long range interaction is fluctuation-induced and we propose a mechanism similar to the Casimir effect that generates it: the hydrodynamic fluctuations are geometrically confined between the intruders, producing an unbalanced renormalized pressure. An estimation based on computing the possible Fourier modes explains the repulsive force and is in qualitative agreement with the simulations., Comment: 4 pages, 3 figures. Accepted in Phys. Rev. Letters

The behavior of a driven granular gas in a container consisting of $M$ connected compartments is studied employing a microscopic kinetic model. After obtaining the governing equations for the occupation numbers and the granular temperatures of each compartment we consider the various dynamical regimes. The system displays interesting analogies with the ordering processes of phase separating mixtures quenched below the their critical point. In particular, we show that below a certain value of the driving intensity the populations of the various compartments become unequal and the system clusterizes. Such a phenomenon is not instantaneous, but is characterized by a time scale, $\tau$, which follows a Vogel-Vulcher exponential behavior. On the other hand, the reverse phenomenon which involves the ``evaporation'' of a cluster due to the driving force is also characterized by a second time scale which diverges at the limit of stability of the cluster., Comment: 11 pages, 17 figures

We study the behavior of an assembly of $N$ granular particles contained in two compartments within a simple kinetic approach. The particles belonging to each compartment collide inelastically with each other and are driven by a stochastic heat bath. In addition, the fastest particles can change compartment at a rate which depends on their kinetic energy. Via a Boltzmann velocity distribution approach we first study the dynamics of the model in terms of a coupled set of equations for the populations in the containers and their granular temperatures and find a crossover from a symmetric high-temperature phase to an asymmetric low temperature phase. Finally, in order to include statistical fluctuations, we solve the model within the Direct Simulation Monte Carlo approach. Comparison with previous studies are presented., Comment: 13 pages, 12 figures

We investigate the phase diagram of a two-component associating fluid mixture in the presence of selectively adsorbing substrates. The mixture is characterized by a bulk phase diagram which displays peculiar features such as closed loops of immiscibility. The presence of the substrates may interfere the physical mechanism involved in the appearance of these phase diagrams, leading to an enhanced tendency to phase separate below the lower critical solution point. Three different cases are considered: a planar solid surface in contact with a bulk fluid, while the other two represent two models of porous systems, namely a slit and an array on infinitely long parallel cylinders. We confirm that surface transitions, as well as capillary transitions for a large area/volume ratio, are stabilized in the one-phase region. Applicability of our results to experiments reported in the literature is discussed., Comment: 12 two-column pages, 12 figures, accepted for publication in Physical Review E; corrected version

Motivated by recent experiments we have carried out an Event Driven computer simulation of a diluted binary mixture of granular particles vertically vibrated in the presence of gravity. The simulations not only confirm that the kinetic energies of the two species are not equally distributed, as predicted by various theoretical models, but also seem to reproduce rather well the density and temperature profiles measured experimentally. Rotational degrees of freedom do not seem to play any important qualitative role. Instead, simulation shows the onset of a clustering instability along the horizontal direction. At the interior of the cluster we observe a secondary instability with respect to the perfect mixing situation, so that segregation of species is observed within the cluster., Comment: 12 pages, 11 figures, submitted for publication

We study the steady state properties of a 2D granular mixture in the presence of energy driving by employing simple analytical estimates and Direct Simulation Monte Carlo. We adopt two different driving mechanisms: a) a homogeneous heat bath with friction and b) a vibrating boundary (thermal or harmonic) in the presence of gravity. The main findings are: the appearance of two different granular temperatures, one for each species; the existence of overpopulated tails in the velocity distribution functions and of non trivial spatial correlations indicating the spontaneous formation of cluster aggregates. In the case of a fluid subject to gravity and to a vibrating boundary, both densities and temperatures display non uniform profiles along the direction normal to the wall, in particular the temperature profiles are different for the two species while the temperature ratio is almost constant with the height. Finally, we obtained the velocity distributions at different heights and verified the non gaussianity of the resulting distributions., Comment: 19 pages, 12 figures, submitted for publication

The free evolution of inelastic particles in one dimension is studied by means of Molecular Dynamics (MD), of an inelastic pseudo-Maxwell model and of a lattice model, with emphasis on the role of spatial correlations. We present an exact solution of the 1d granular pseudo-Maxwell model for the scaling distribution of velocities and discuss how this model fails to describe correctly the homogeneous cooling stage of the 1d granular gas. Embedding the pseudo-Maxwell gas on a lattice (hence allowing for the onset of spatial correlations), we find a much better agreement with the MD simulations even in the inhomogeneous regime. This is seen by comparing the velocity distributions, the velocity profiles and the structure factors of the velocity field., Comment: Latex file: 6 pages, 5 figures (.eps). See also http://axtnt3.phys.uniroma1.it/Maxwell

We present a new microscopic model of granular medium to study the role of dynamical correlations and the onset of spatial order induced by the inelasticity of the interactions. In spite of its simplicity, it features several different aspects of the rich phenomenology observed in granular materials and allows to make contact with other topics of statistical mechanics such as diffusion processes, domain growth, persistence, aging phenomena. Interestingly, while local observables being controlled by the largest wavelength fluctuations seem to suggest a purely diffusive behavior, the formation of spatially extended structures and topological defects, such as vortices and shocks, reveals a more complex scenario., Comment: 4 pages, 4 figures

We study fluidized granular gases in a stationary state determined by the balance between an external driving and the bulk dissipation. The two considered situations are inspired by recent experiments, where the gravity plays a major role as a driving mechanism: in the first case gravity acts only in one direction and the bottom wall is vibrated, in the second case gravity acts in both directions and no vibrating walls are present. Simulations performed under the molecular chaos assumption show averaged profiles of density, velocity and granular temperature which are in good agreement with the experiments. Moreover we measure the velocity distributions which show strong non-Gaussian behavior, as experiments pointed out, but also density correlations accounting for clustering, at odds with the experimental results. The hydrodynamics of the first model is discussed and an exact solution is found for the density and granular temperature as functions of the distance from the vibrating wall. The limitations of such a solution, in particular in a broad layer near the wall injecting energy, are discussed., Comment: Revised version accepted for publication. New results added and discussions considering tangential forces. 27 pages (19 figures included), to appear in Phys.Rev.E

We present a novel picture of a non isothermal solidification process starting from a molecular level, where the microscopic origin of the basic mechanisms and of the instabilities characterizing the approach to equilibrium is rendered more apparent than in existing approaches based on coarse grained free energy functionals \`a la Landau. The system is composed by a lattice of Potts spins, which change their state according to the stochastic dynamics proposed some time ago by Creutz. Such a method is extended to include the presence of latent heat and thermal conduction. Not only the model agrees with previous continuum treatments, but it allows to introduce in a consistent fashion the microscopic stochastic fluctuations. These play an important role in nucleating the growing solid phase in the melt. The approach is also very satisfactory from the quantitative point of view since the relevant growth regimes are fully characterized in terms of scaling exponents., Comment: 7 pages Latex +3 figures.ps

We address the problem of the so-called ``granular gases'', i.e. gases of massive particles in rapid movement undergoing inelastic collisions. We introduce a class of models of driven granular gases for which the stationary state is the result of the balance between the dissipation and the random forces which inject energies. These models exhibit a genuine thermodynamic limit, i.e. at fixed density the mean values of kinetic energy and dissipated energy per particle are independent of the number $N$ of particles, for large values of $N$. One has two regimes: when the typical relaxation time $\tau$ of the driving Brownian process is small compared with the mean collision time $\tau_c$ the spatial density is nearly homogeneous and the velocity probability distribution is gaussian. In the opposite limit $\tau \gg \tau_c$ one has strong spatial clustering, with a fractal distribution of particles, and the velocity probability distribution strongly deviates from the gaussian one. Simulations performed in one and two dimensions under the {\it Stosszahlansatz} Boltzmann approximation confirm the scenario. Furthermore we analyze the instabilities bringing to the spatial and the velocity clusterization. Firstly, in the framework of a mean-field model, we explain how the existence of the inelasticity can bring to a spatial clusterization; on the other side we discuss, in the framework of a Langevin dynamics treating the collisions in a mean-field way, how a non-gaussian distribution of velocity can arise. The comparison between the numerical and the analytical results exhibits an excellent agreement., Comment: 33 pages including 25 figures, submitted to Phys.Rev.E

We study the problem of phase separation in systems with a positive definite order parameter, and in particular, in systems with absorbing states. Owing to the presence of a single minimum in the free energy driving the relaxation kinetics, there are some basic properties differing from standard phase separation. We study analytically and numerically this class of systems; in particular we determine the phase diagram, the growth laws in one and two dimensions and the presence of scale invariance. Some applications are also discussed., Comment: Submitted to Europhysics Lett

The present paper focuses on the order-disorder transition of an Ising model on a self-similar lattice. We present a detailed numerical study, based on the Monte Carlo method in conjunction with the finite size scaling method, of the critical properties of the Ising model on some two dimensional deterministic fractal lattices with different Hausdorff dimensions. Those with finite ramification order do not display ordered phases at any finite temperature, whereas the lattices with infinite connectivity show genuine critical behavior. In particular we considered two Sierpinski carpets constructed using different generators and characterized by Hausdorff dimensions d_H=log 8/log 3 = 1.8927.. and d_H=log 12/log 4 = 1.7924.., respectively. The data show in a clear way the existence of an order-disorder transition at finite temperature in both Sierpinski carpets. By performing several Monte Carlo simulations at different temperatures and on lattices of increasing size in conjunction with a finite size scaling analysis, we were able to determine numerically the critical exponents in each case and to provide an estimate of their errors. Finally we considered the hyperscaling relation and found indications that it holds, if one assumes that the relevant dimension in this case is the Hausdorff dimension of the lattice., Comment: 21 pages, 7 figures; a new section has been added with results for a second fractal; there are other minor changes

We investigate the properties of a model of granular matter consisting of $N$ Brownian particles on a line subject to inelastic mutual collisions. This model displays a genuine thermodynamic limit for the mean values of the energy and the energy dissipation. When the typical relaxation time $\tau$ associated with the Brownian process is small compared with the mean collision time $\tau_c$ the spatial density is nearly homogeneous and the velocity probability distribution is gaussian. In the opposite limit $\tau \gg \tau_c$ one has strong spatial clustering, with a fractal distribution of particles, and the velocity probability distribution strongly deviates from the gaussian one., Comment: 4 pages including 3 eps figures, LaTex, added references, corrected typos, minimally changed contents and abstract, to published in Phys.Rev.Lett. (tentatively on 28th of October, 1998)

We present a theoretical study for the intermediate stages of the growth of membranes and vesicles in supersaturated solutions of amphiphilic molecules. The problem presents important differences with the growth of droplets in the classical theory of Lifshitz-Slyozov-Wagner, because the aggregates are extensive only in two dimensions, but still grow in a three dimensional bath. The balance between curvature and edge energy favours the nucleation of small planar membranes, but as they grow beyond a critical size they close themselves to form vesicles. We obtain a system of coupled equations describing the growth of planar membranes and vesicles, which is solved numerically for different initial conditions. Finally, the range of parameters relevant in experimental situations is discussed., Comment: 13 pages and 5 postscript figures. To appear in Phys. Rev E

We discuss the static and kinetic properties of a Ginzburg-Landau spherically symmetric $O(N)$ model recently introduced (Phys. Rev. Lett. {\bf 75}, 2176, (1995)) in order to generalize the so called Phase field model of Langer. The Hamiltonian contains two $O(N)$ invariant fields $\phi$ and $U$ bilinearly coupled. The order parameter field $\phi$ evolves according to a non conserved dynamics, whereas the diffusive field $U$ follows a conserved dynamics. In the limit $N \to \infty$ we obtain an exact solution, which displays an interesting kinetic behavior characterized by three different growth regimes. In the early regime the system displays normal scaling and the average domain size grows as $t^{1/2}$, in the intermediate regime one observes a finite wavevector instability, which is related to the Mullins-Sekerka instability; finally, in the late stage the structure function has a multiscaling behavior, while the domain size grows as $t^{1/4}$., Comment: 9 pages RevTeX, 9 figures included, files packed with uufiles to appear on Phy. Rev. E

To study the effect of slow heat conduction during phase separation, we discuss the relaxation properties of an $O(N)$ symmetric model with Phase field type dynamics, where a non conserved order parameter field couples bilinearly to a diffusive field. In the limit $N\to\infty$ we obtain an exact solution. The analysis reveals three different types of growth regimes and a very rich dynamical behavior. Finally the connection with the Mullins-Sekerka instability is expounded., Comment: 4 pages, 4 figures included, Postscript file packed with uufiles to appear on Phy. Rev. Lett

We extend the discussion of the growth kinetics of the large-N time-dependent Ginzburg-Landau model with an order parameter described by a $\Phi^6$ free energy functional, to the conserved case. Quenches from a high temperature initial state to a zero temperature state are studied for different selections of parameters entering the effective potential. In all cases we obtain the asymptotic form of the structure factor. As expected in analogy with the well studied $\Phi^4$ model, we find multiscaling behavior whenever stable equilibrium is reached in the ordering region. On the other hand the present model also displays a novel feature, namely the occurrence of metastable relaxation., Comment: 20 pages,Plain Latex

In this paper we present a formulation of the nonlinear stochastic differential equation which allows for systematic approximations. The method is not restricted to the asymptotic, i.e., stationary, regime but can be applied to derive effective equations describing the relaxation of the system from arbitrary initial conditions. The basic idea is to reduce the nonlinear Langevin equation to an equivalent equilibrium problem, which can then be studied with the methods of conventional equilibrium statistical field theory. A particular well suited perturbative scheme is that developed in quantum field theory by Cornwall, Jackiw and Tomboulis. We apply this method to the study of $N$ component Ginzburg-Landau equation in zero spatial dimension. In the limit of $N\to\infty$ we can solve the effective equations and obtain closed forms for the time evolutions of the average field and of the two-time connected correlation functions., Comment: 19 Pages + 4 Figures, RevTeX 3.0, sub. Phys. Rev. E

We compute the corrections of next to leading order in the ${1 \over N}$ expansion to the effective potential of a system described by a Ginzburg-Landau model with $N$ components and quartic interaction, in the case of spontaneous symmetry breaking. The method we apply allows to generalize in a simple way the so-called Self-Consistent Screened Approximation (SCSA)., Comment: p. 8, LATEX, DFF 193/9/1993