Effective pair potentials for spherical nanoparticles

An effective description for spherical nanoparticles in a fluid of point particles is presented. The points inside the nanoparticles and the point particles are assumed to interact via spherically symmetric additive pair potentials, while the distribution of points inside the nanoparticles is taken to be spherically symmetric and smooth. The resulting effective pair interactions between a nanoparticle and a point particle, as well as between two nanoparticles, are then given by spherically symmetric potentials. If overlap between particles is allowed, the effective potential generally has non-analytic points, but for each effective potential the expressions for different overlapping cases can be written in terms of one analytic auxiliary potential. Effective potentials for hollow nanoparticles (appropriate e.g. for buckyballs) are also considered, and shown to be related to those for solid nanoparticles. Finally, explicit expressions are given for the effective potentials derived from basic pair potentials of power law and exponential form, as well as from the commonly used London-Van der Waals, Morse, Buckingham, and Lennard-Jones potential. The applicability of the latter is demonstrated by comparison with an atomic description of nanoparticles with an internal face centered cubic structure.Comment: 27 pages, 12 figures. Unified description of overlapping and nonoverlapping particles added, as well as a comparison with an idealized atomic descriptio

Crucial role of sidewalls in velocity distributions in quasi-2D granular gases

Our experiments and three-dimensional molecular dynamics simulations of particles confined to a vertical monolayer by closely spaced frictional walls (sidewalls) yield velocity distributions with non-Gaussian tails and a peak near zero velocity. Simulations with frictionless sidewalls are not peaked. Thus interactions between particles and their container are an important determinant of the shape of the distribution and should be considered when evaluating experiments on a tightly constrained monolayer of particles.Comment: 4 pages, 4 figures, Added reference, model explanation charified, other minor change

Velocity distributions in dissipative granular gases

Motivated by recent experiments reporting non-Gaussian velocity distributions in driven dilute granular materials, we study by numerical simulation the properties of 2D inelastic gases. We find theoretically that the form of the observed velocity distribution is governed primarily by the coefficient of restitution $\eta$ and $q=N_H/N_C$, the ratio between the average number of heatings and the average number of collisions in the gas. The differences in distributions we find between uniform and boundary heating can then be understood as different limits of $q$, for $q \gg 1$ and $q \lesssim 1$ respectively.Comment: 5 figure

Extended Heat-Fluctuation Theorems for a System with Deterministic and Stochastic Forces

Heat fluctuations over a time \tau in a non-equilibrium stationary state and in a transient state are studied for a simple system with deterministic and stochastic components: a Brownian particle dragged through a fluid by a harmonic potential which is moved with constant velocity. Using a Langevin equation, we find the exact Fourier transform of the distribution of these fluctuations for all \tau. By a saddle-point method we obtain analytical results for the inverse Fourier transform, which, for not too small \tau, agree very well with numerical results from a sampling method as well as from the fast Fourier transform algorithm. Due to the interaction of the deterministic part of the motion of the particle in the mechanical potential with the stochastic part of the motion caused by the fluid, the conventional heat fluctuation theorem is, for infinite and for finite \tau, replaced by an extended fluctuation theorem that differs noticeably and measurably from it. In particular, for large fluctuations, the ratio of the probability for absorption of heat (by the particle from the fluid) to the probability to supply heat (by the particle to the fluid) is much larger here than in the conventional fluctuation theorem.Comment: 23 pages, 6 figures. Figures are now in color, Eq. (67) was corrected and a footnote was added on the d-dimensional cas

Mode-coupling theory for multiple-time correlation functions of tagged particle densities and dynamical filters designed for glassy systems

The theoretical framework for higher-order correlation functions involving multiple times and multiple points in a classical, many-body system developed by Van Zon and Schofield [Phys. Rev. E 65, 011106 (2002)] is extended here to include tagged particle densities. Such densities have found an intriguing application as proposed measures of dynamical heterogeneities in structural glasses. The theoretical formalism is based upon projection operator techniques which are used to isolate the slow time evolution of dynamical variables by expanding the slowly-evolving component of arbitrary variables in an infinite basis composed of the products of slow variables of the system. The resulting formally exact mode-coupling expressions for multiple-point and multiple-time correlation functions are made tractable by applying the so-called N-ordering method. This theory is used to derive for moderate densities the leading mode coupling expressions for indicators of relaxation type and domain relaxation, which use dynamical filters that lead to multiple-time correlations of a tagged particle density. The mode coupling expressions for higher order correlation functions are also succesfully tested against simulations of a hard sphere fluid at relatively low density.Comment: 15 pages, 2 figure

Onsager-Machlup theory and work fluctuation theorem for a harmonically driven Brownian particle

We extend Tooru-Cohen analysis for nonequilirium steady state(NSS) of a Brownian particle to nonequilibrium oscillatory state (NOS) of Brownian particle by considering time dependent external drive protocol. We consider an unbounded charged Brownian particle in the presence of an oscillating electric field and prove work fluctuation theorem, which is valid for any initial distribution and at all times. For harmonically bounded and constantly dragged Brownian particle considered by Tooru and Cohen, work fluctuation theorem is valid for any initial condition(also NSS), but only in large time limit. We use Onsager-Machlup Lagrangian with a constraint to obtain frequency dependent work distribution function, and describe entropy production rate and properties of dissipation functions for the present system using Onsager-Machlup functional.Comment: 6 pages, 1 figur

Grand canonical ensemble in generalized thermostatistics

We study the grand-canonical ensemble with a fluctuating number of degrees of freedom in the context of generalized thermostatistics. Several choices of grand-canonical entropy functional are considered. The ideal gas is taken as an example.Comment: 14 pages, no figure

Lyapunov Exponent Pairing for a Thermostatted Hard-Sphere Gas under Shear in the Thermodynamic Limit

We demonstrate why for a sheared gas of hard spheres, described by the SLLOD equations with an iso-kinetic Gaussian thermostat in between collisions, deviations of the conjugate pairing rule for the Lyapunov spectrum are to be expected, employing a previous result that for a large number of particles $N$, the iso-kinetic Gaussian thermostat is equivalent to a constant friction thermostat, up to $1/\sqrt{N}$ fluctuations. We also show that these deviations are at most of the order of the fourth power in the shear rate.Comment: 4 pages, to appear in Rapid Comm., Phys. Rev.

Comment on ``Lyapunov Exponent of a Many Body System and Its Transport Coefficients''

In a recent Letter, Barnett, Tajima, Nishihara, Ueshima and Furukawa obtained a theoretical expression for the maximum Lyapunov exponent $\lambda_1$ of a dilute gas. They conclude that $\lambda_1$ is proportional to the cube root of the self-diffusion coefficient $D$, independent of the range of the interaction potential. They validate their conjecture with numerical data for a dense one-component plasma, a system with long-range forces. We claim that their result is highly non-generic. We show in the following that it does not apply to a gas of hard spheres, neither in the dilute nor in the dense phase.Comment: 1 page, Revtex - 1 PS Figs - Submitted to Physical Review Letter

Largest Lyapunov Exponent for Many Particle Systems at Low Densities

The largest Lyapunov exponent $\lambda^+$ for a dilute gas with short range interactions in equilibrium is studied by a mapping to a clock model, in which every particle carries a watch, with a discrete time that is advanced at collisions. This model has a propagating front solution with a speed that determines $\lambda^+$, for which we find a density dependence as predicted by Krylov, but with a larger prefactor. Simulations for the clock model and for hard sphere and hard disk systems confirm these results and are in excellent mutual agreement. They show a slow convergence of $\lambda^+$ with increasing particle number, in good agreement with a prediction by Brunet and Derrida.Comment: 4 pages, RevTeX, 2 Figures (encapsulated postscript). Submitted to Phys. Rev. Let