
Networks of channels conveying particles are often subject to blockages due
to the limited carrying capacity of the individual channels. If the channels
are coupled, blockage of one causes an increase in the flux entering the
remaining open channels leading to a cascade of failures. Once all channels are
blocked no additional particle can enter the system. If the blockages are of
finite duration, however, the system reaches a steady state with an exiting
flux that is reduced compared to the incoming one. We propose a stochastic
model consisting of $N_c$ channels each with a blocking threshold of $N$
particles. Particles enter the system according to a Poisson process with the
entering flux of intensity $\Lambda$ equally distributed over the open
channels. Any particle in an open channel exits at a rate $\mu$ and a blocked
channel unblocks at a rate $\mu^*$. We present a method to obtain the exiting
flux in the steady state, and other properties, for arbitrary $N_c$ and $N$ and
we present explicit solutions for $N_c=2,3$. We apply these results to compare
the efficiency of conveying a particulate stream of intensity $\Lambda$ using
different channel configurations. We compare a single "robust" channel with a
large capacity with multiple "fragile" channels with a proportionately reduced
capacity. The "robust" channel is more efficient at low intensity, while
multiple, "fragile" channels have a higher throughput at large intensity. We
also compare $N_c$ coupled channels with $N_c$ independent channels, both with
threshold $N=2$. For $N_c=2$ if $\mu^*/\mu>1/4$, the coupled channels are
always more efficient. Otherwise the independent channels are more efficient
for sufficiently large $\Lambda$.

We study by Molecular Dynamics simulation a dense onecomponent system of
particles confined on a spherical substrate. We more specifically investigate
the evolution of the structural and dynamical properties of the system when
changing the control parameters, the temperature and the curvature of the
substrate. We find that the dynamics becomes glassy at low temperature, with a
strong slowdown of the relaxation and the emergence of dynamical heterogeneity.
The prevalent local $6$fold order is frustrated by curvature and we analyze in
detail the role of the topological defects in the statics and the dynamics of
the particle assembly.

We consider the two dimensional motion of a particle into a confining
potential, subjected to Brownian forces, associated with two different
temperatures on the orthogonal directions. Exact solutions are obtained for an
asymmetric harmonic potential in the overdamped and underdamped regimes,
whereas perturbative approaches are used for more general potentials. The
resulting non equilibrium stationary state is characterized with a nonzero
orthoradial mean current, corresponding to a global rotation of the particle
around the center. The rotation is due to two symmetry breaking: two different
temperatures and a mismatch between the principal axes of the confining
asymmetric potential and the temperature axes. We confirm our predictions by
performing Brownian dynamics simulation. Finally, we propose to observe this
effect on a laser cooled atomic system.

Longrange interacting Hamiltonian systems are believed to relax generically
towards nonequilibrium states called "quasistationary" because they evolve
towards thermodynamic equilibrium very slowly, on a timescale diverging with
particle number. We show here that, by applying a suitable perturbation
operator for a finite time interval, we obtain, in a family of longrange
systems, nonequilibrium states which appear to be strictly stationary. They
exist even in the case of a harmonic potential, and are characterised by an
ordered microscopic phase space structure. We give some simple heuristic
arguments which predict reasonably well some properties of these states.

Isolated longrange interacting particle systems appear generically to relax
to nonequilibrium states ("quasistationary states" or QSS) which are
stationary in the thermodynamic limit. A fundamental open question concerns the
"robustness" of these states when the system is not isolated. In this paper we
explore, using both analytical and numerical approaches to a paradigmatic one
dimensional model, the effect of a simple class of perturbations. We call them
"internal local perturbations" in that the particle energies are perturbed at
collisions in a way which depends only on the local properties. Our central
finding is that the effect of the perturbations is to drive all the very
different QSS we consider towards a unique QSS. The latter is thus independent
of the initial conditions of the system, but determined instead by both the
longrange forces and the details of the perturbations applied. Thus in the
presence of such a perturbation the longrange system evolves to a unique
nonequilibrium stationary state, completely different to its state in absence
of the perturbation, and it remains in this state when the perturbation is
removed. We argue that this result may be generic for longrange interacting
systems subject to perturbations which are dependent on the local properties
(e.g. spatial density or velocity distribution) of the system itself.

We study the dynamics of a onecomponent liquid constrained on a spherical
substrate, a 2sphere, and investigate how the modecoupling theory (MCT) can
describe the new features brought by the presence of curvature. To this end we
have derived the MCT equations in a spherical geometry. We find that, as seen
from the MCT, the slow dynamics of liquids in curved space at low temperature
does not qualitatively differ from that of glassforming liquids in Euclidean
space. The MCT predicts the right trend for the evolution of the relaxation
slowdown with curvature but is dramatically off at a quantitative level.

We investigate by Molecular Dynamics simulation a system of $N$ particles
moving on the surface of a twodimensional sphere and interacting by a
LennardJones potential. We detail the way to account for the changes brought
by a nonzero curvature, both at a methodological and at a physical level. When
compared to a twodimensional LennardJones liquid on the Euclidean plane,
where a phase transition to an ordered hexagonal phase takes place, we find
that the presence of excess defects imposed by the topology of the sphere
frustrates the hexagonal order. We observe at high density a rapid increase of
the relaxation time when the temperature is decreased, whereas in the same
range of temperature the pair correlation function of the system evolves only
moderately.

Hamiltonian systems with longrange interactions give rise to long lived out
of equilibrium macroscopic states, socalled quasistationary states. We show
here that, in a suitably generalized form, this result remains valid for many
such systems in the presence of dissipation. Using an appropriate meanfield
kinetic description, we show that models with dissipation due to a viscous
damping or due to inelastic collisions admit "scaling quasistationary states",
i.e., states which are quasistationary in rescaled variables. A numerical
study of one dimensional selfgravitating systems confirms both the relevance
of these solutions, and gives indications of their regime of validity in line
with theoretical predictions. We underline that the velocity distributions
never show any tendency to evolve towards a MaxwellBoltzmann form.

We investigate models in which blocking can interrupt a particulate flow
process at any time. Filtration, and flow in micro/nanochannels and traffic
flow are examples of such processes. We first consider concurrent flow models
where particles enter a channel randomly. If at any time two particles are
simultaneously present in the channel, failure occurs. The key quantities are
the survival probability and the distribution of the number of particles that
pass before failure. We then consider a counterflow model with two opposing
Poisson streams. There is no restriction on the number of particles passing in
the same direction, but blockage occurs if, at any time, two opposing particles
are simultaneously present in the passage.

We investigate the influence of dry friction on an asymmetric, granular
piston of mass $M$ composed of two materials undergoing inelastic collisions
with bath particles of mass $m$. Numerical simulations of the BoltzmannLorentz
equation reveal the existence of two scaling regimes depending on the strength
of friction. In the large friction limit, we introduce an exact model giving
the asymptotic behavior of the BoltzmannLorentz equation. For small friction
and for large mass ratio $M/m$, we derive a FokkerPlanck equation for which
the exact solution is also obtained. Static friction attenuates the motor
effect and results in a discontinuous velocity distribution.

Within the framework of a BoltzmannLorentz equation, we analyze the dynamics
of a granular rotor immersed in a bath of thermalized particles in the presence
of a frictional torque on the axis. In numerical simulations of the equation,
we observe two scaling regimes at low and high bath temperatures. In the large
friction limit, we obtain the exact solution of a model corresponding to
asymptotic behavior of the BoltzmannLorentz equation. In the limit of large
rotor mass and small friction, we derive a FokkerPlanck equation for which the
exact solution is also obtained.

We study the properties of a heterogeneous, chiral granular rotor that is
capable of performing useful work when immersed in a bath of thermalized
particles. The dynamics can be obtained in general from a numerical solution of
the BoltzmannLorentz equation. We show that a mechanical approach gives the
exact mean angular velocity in the limit of an infinitely massive rotor. We
examine the dependence of the mean angular velocity on the coefficients of
restitution of the two materials composing the motor. We compute the power and
efficiency and compare with numerical simulations. We also perform a realistic
numerical simulation of a granular rotor which shows that the presence of non
uniformity of the bath density within the region where the motor rotates, and
that the ratchet effect is slightly weakened, but qualitatively sustained.
Finally we discuss the results in connection with recent experiments.

In this article, we review the progress made on the statistical mechanics of
liquids and fluids embedded in curved space. Our main focus will be on
twodimensional manifolds of constant nonzero curvature and on the influence of
the latter on the phase behavior, thermodynamics and structure of simple
liquids. Reference will also be made to existing work on threedimensional
curved space and twodimensional manifolds with varying curvature.

We provide a consistent statisticalmechanical treatment for describing the
thermodynamics and the structure of fluids embedded in the hyperbolic plane. In
particular, we derive a generalization of the virial equation relating the bulk
thermodynamic pressure to the pair correlation function and we develop the
appropriate setting for extending the integralequation approach of
liquidstate theory in order to describe the fluid structure. We apply the
formalism and study the influence of negative space curvature on two types of
systems that have been recently considered: Coulombic systems, such as the one
and twocomponent plasma models, and fluids interacting through shortrange
pair potentials, such as the harddisk and the LennardJones models.

We show by using a Discrete Element Method that a bilayer of vibrated
granular bidisperse spheres exhibits the striking feature that the horizontal
velocity distribution of the top layer particles has a quasiGaussian shape,
whereas that of the bottom layer is far from Gaussian. We examine in detail the
relevance of all physical parameters (acceleration of the bottom plate, mass
ratio, layer coverage). Moreover, a microscopic analysis of the trajectories
and the collision statistics reveal how the mechanism of randomization.

We investigate the influence of space curvature, and of the associated
"frustration", on the dynamics of a model glassformer: a monatomic liquid on
the hyperbolic plane. We find that the system's fragility, i.e. the sensitivity
of the relaxation time to temperature changes, increases as one decreases the
frustration. As a result, curving space provides a way to tune fragility and
make it as large as wanted. We also show that the nature of the emerging
"dynamic heterogeneities", another distinctive feature of slowly relaxing
systems, is directly connected to the presence of frustrationinduced
topological defects.

We investigate the velocitycorrelation distributions after $n$ collisions of
a tagged particle undergoing binary collisions. Analytical expressions are
obtained in any dimension for the velocitycorrelation distribution after the
firstcollision as well as for the velocitycorrelation function after an
infinite number of collisions, in the limit of Gaussian velocity distributions.
It appears that the decay of the firstcollision velocitycorrelation
distribution for negative argument is exponential in any dimension with a
coefficient that depends on the mass and on the coefficient of restitution. We
also obtained the velocitycorrelation distribution when the velocity
distributions are not Gaussian: by inserting Sonine corrections of the velocity
distributions, we derive the corrections to the velocitycorrelation
distribution which agree perfectly with a DSMC (Direct Simulation Monte Carlo)
simulation. We emphasize that these new quantities can be easily obtained in
simulations and likely in experiments: they could be an efficient probe of the
local environment and of the degree of inelasticity of the collisions.

We introduce a minimal model describing the physics of classical
twodimensional (2D) frustrated Heisenberg systems, where spins order in a
nonplanar way at T=0. This model, consisting of coupled trihedra (or
Ising$\mathbb{R}P^3$ model), encompasses Ising (chiral) degrees of freedom,
spinwave excitations and $\Z_2$ vortices. Extensive Monte Carlo simulations
show that the T=0 chiral order disappears at finite temperature in a continuous
phase transition in the 2D Ising universality class, despite misleading
intermediatesize effects observed at the transition. The analysis of
configurations reveals that shortrange spin fluctuations and $\Z_2$ vortices
proliferate near the chiral domain walls explaining the strong renormalization
of the transition temperature. Chiral domain walls can themselves carry an
unlocalized $\Z_2$ topological charge, and vortices are then preferentially
paired with charged walls. Further, we conjecture that the anomalous
sizeeffects suggest the proximity of the present model to a tricritical point.
A body of results is presented, that all support this claim: (i) Firstorder
transitions obtained by Monte Carlo simulations on several related models (ii)
Approximate mapping between the Ising$\mathbb{R}P^3$ model and a dilute Ising
model (exhibiting a tricritical point) and, finally, (iii) Meanfield results
obtained for Isingmultispin Hamiltonians, derived from the hightemperature
expansion for the vector spins of the Ising$\mathbb{R}P^3$ model.

We propose that the salient feature to be explained about the glass
transition of supercooled liquids is the temperaturecontrolled superArrhenius
activated nature of the viscous slowing down, more strikingly seen in
weaklybonded, fragile systems. In the light of this observation, the relevance
of simple models of spherically interacting particles and that of models based
on freevolume congested dynamics are questioned. Finally, we discuss how the
main aspects of the phenomenology of supercooled liquids, including the
crossover from Arrhenius to superArrhenius activated behavior and the
heterogeneous character of the $\alpha$ relaxation, can be described by an
approach based on frustrationlimited domains.

Many experimental studies of protein deposition on solid surfaces involve
alternating adsorption/desorption steps. In this paper, we investigate the
effect of a desorption step (separating two adsorption steps) on the kinetics,
the adsorbedlayer structure, and the saturation density. Our theoretical
approach involves a density expansion of the pair distribution function and an
application of an interpolation formula to estimate the saturation density as a
function of the density at which the desorption process commences, $\rho_1$,
and the density of the depleted configuration, $\rho_2$. The theory predicts an
enhancement of the saturation density compared with that of a simple,
uninterrupted RSA process and a maximum in the saturation density when
$\rho_2={2/3}\rho_1$. The theoretical results are in qualitative and in
semiquantitative agreement with the results of numerical simulations.