N-body decomposition of bipartite networks

In this paper, we present a method to project co-authorship networks, that accounts in detail for the geometrical structure of scientists collaborations. By restricting the scope to 3-body interactions, we focus on the number of triangles in the system, and show the importance of multi-scientists (more than 2) collaborations in the social network. This motivates the introduction of generalized networks, where basic connections are not binary, but involve arbitrary number of components. We focus on the 3-body case, and study numerically the percolation transition.Comment: 5 pages, submitted to PR

Generalized Fokker-Planck equation, Brownian motion, and ergodicity

Microscopic theory of Brownian motion of a particle of mass $M$ in a bath of molecules of mass $m\ll M$ is considered beyond lowest order in the mass ratio $m/M$. The corresponding Langevin equation contains nonlinear corrections to the dissipative force, and the generalized Fokker-Planck equation involves derivatives of order higher than two. These equations are derived from first principles with coefficients expressed in terms of correlation functions of microscopic force on the particle. The coefficients are evaluated explicitly for a generalized Rayleigh model with a finite time of molecule-particle collisions. In the limit of a low-density bath, we recover the results obtained previously for a model with instantaneous binary collisions. In general case, the equations contain additional corrections, quadratic in bath density, originating from a finite collision time. These corrections survive to order $(m/M)^2$ and are found to make the stationary distribution non-Maxwellian. Some relevant numerical simulations are also presented

Strong Shock Waves and Nonequilibrium Response in a One-dimensional Gas: a Boltzmann Equation Approach

We investigate the nonequilibrium behavior of a one-dimensional binary fluid on the basis of Boltzmann equation, using an infinitely strong shock wave as probe. Density, velocity and temperature profiles are obtained as a function of the mixture mass ratio \mu. We show that temperature overshoots near the shock layer, and that heavy particles are denser, slower and cooler than light particles in the strong nonequilibrium region around the shock. The shock width w(\mu), which characterizes the size of this region, decreases as w(\mu) ~ \mu^{1/3} for \mu-->0. In this limit, two very different length scales control the fluid structure, with heavy particles equilibrating much faster than light ones. Hydrodynamic fields relax exponentially toward equilibrium, \phi(x) ~ exp[-x/\lambda]. The scale separation is also apparent here, with two typical scales, \lambda_1 and \lambda_2, such that \lambda_1 ~ \mu^{1/2} as \mu-->0$, while \lambda_2, which is the slow scale controlling the fluid's asymptotic relaxation, increases to a constant value in this limit. These results are discussed at the light of recent numerical studies on the nonequilibrium behavior of similar 1d binary fluids.Comment: 9 pages, 8 figs, published versio

Singular forces and point-like colloids in lattice Boltzmann hydrodynamics

We present a second-order accurate method to include arbitrary distributions of force densities in the lattice Boltzmann formulation of hydrodynamics. Our method may be used to represent singular force densities arising either from momentum-conserving internal forces or from external forces which do not conserve momentum. We validate our method with several examples involving point forces and find excellent agreement with analytical results. A minimal model for dilute sedimenting particles is presented using the method which promises a substantial gain in computational efficiency.Comment: 22 pages, 9 figures. Submitted to Phys. Rev.

Coupling of thermal and mass diffusion in regular binary thermal lattice-gases

We have constructed a regular binary thermal lattice-gas in which the thermal diffusion and mass diffusion are coupled and form two nonpropagating diffusive modes. The power spectrum is shown to be similar in structure as for the one in real fluids, in which the central peak becomes a combination of coupled entropy and concentration contributions. Our theoretical findings for the power spectra are confirmed by computer simulations performed on this model.Comment: 5 pages including 3 figures in RevTex

Field induced stationary state for an accelerated tracer in a bath

Our interest goes to the behavior of a tracer particle, accelerated by a constant and uniform external field, when the energy injected by the field is redistributed through collision to a bath of unaccelerated particles. A non equilibrium steady state is thereby reached. Solutions of a generalized Boltzmann-Lorentz equation are analyzed analytically, in a versatile framework that embeds the majority of tracer-bath interactions discussed in the literature. These results --mostly derived for a one dimensional system-- are successfully confronted to those of three independent numerical simulation methods: a direct iterative solution, Gillespie algorithm, and the Direct Simulation Monte Carlo technique. We work out the diffusion properties as well as the velocity tails: large v, and either large -v, or v in the vicinity of its lower cutoff whenever the velocity distribution is bounded from below. Particular emphasis is put on the cold bath limit, with scatterers at rest, which plays a special role in our model.Comment: 20 pages, 6 figures v3:minor corrections in sec.III and added reference

Generalized dynamical density functional theory for classical fluids and the significance of inertia and hydrodynamic interactions

We study the dynamics of a colloidal fluid including inertia and hydrodynamic interactions, two effects which strongly influence the non-equilibrium properties of the system. We derive a general dynamical density functional theory (DDFT) which shows very good agreement with full Langevin dynamics. In suitable limits, we recover existing DDFTs and a Navier-Stokes-like equation with additional non-local terms.Comment: 5 pages, 4 figures, 4 supplementary movie files, I supplementary pd

The electrical double layer for a fully asymmetric electrolyte around a spherical colloid: an integral equation study

The hypernetted chain/mean spherical approximation (HNC/MSA) integral equation is obtained and solved numerically for a totally asymmetric primitive model electrolyte around a spherical macroparticle. The ensuing radial distribution functions show a very good agreement when compared to our Monte Carlo and molecular dynamics simulations for spherical geometry and with respect to previous anisotropic reference HNC calculations in the planar limit. We report an analysis of the potential vs charge relationship, radial distribution functions, mean electrostatic potential and cumulative reduced charge for representative cases of 1:1 and 2:2 salts with a size asymmetry ratio of 2. Our results are collated with those of the Modified Gouy-Chapman (MGC) and unequal radius Modified Gouy-Chapman (URMGC) theories and with those of HNC/MSA in the restricted primitive model (RPM) to assess the importance of size asymmetry effects. One of the most striking characteristics found is that,\textit{contrary to the general belief}, away from the point of zero charge the properties of an asymmetric electrical double layer (EDL) are not those corresponding to a symmetric electrolyte with the size and charge of the counterion, i.e. \textit{counterions do not always dominate}. This behavior suggests the existence of a new phenomenology in the EDL that genuinely belongs to a more realistic size-asymmetric model where steric correlations are taken into account consistently. Such novel features can not be described by traditional mean field theories like MGC, URMGC or even by enhanced formalisms, like HNC/MSA, if they are based on the RPM.Comment: 29 pages, 13 figure

On the velocity distributions of the one-dimensional inelastic gas

We consider the single-particle velocity distribution of a one-dimensional fluid of inelastic particles. Both the freely evolving (cooling) system and the non-equilibrium stationary state obtained in the presence of random forcing are investigated, and special emphasis is paid to the small inelasticity limit. The results are obtained from analytical arguments applied to the Boltzmann equation along with three complementary numerical techniques (Molecular Dynamics, Direct Monte Carlo Simulation Methods and iterative solutions of integro-differential kinetic equations). For the freely cooling fluid, we investigate in detail the scaling properties of the bimodal velocity distribution emerging close to elasticity and calculate the scaling function associated with the distribution function. In the heated steady state, we find that, depending on the inelasticity, the distribution function may display two different stretched exponential tails at large velocities. The inelasticity dependence of the crossover velocity is determined and it is found that the extremely high velocity tail may not be observable at ``experimentally relevant'' inelasticities.Comment: Latex, 14 pages, 12 eps figure

Kinetics and scaling in ballistic annihilation

We study the simplest irreversible ballistically-controlled reaction, whereby particles having an initial continuous velocity distribution annihilate upon colliding. In the framework of the Boltzmann equation, expressions for the exponents characterizing the density and typical velocity decay are explicitly worked out in arbitrary dimension. These predictions are in excellent agreement with the complementary results of extensive Monte Carlo and Molecular Dynamics simulations. We finally discuss the definition of universality classes indexed by a continuous parameter for this far from equilibrium dynamics with no conservation laws