Approach to a Stationary State in an External Field

We study relaxation towards a stationary out-of-equilibrium state by analyzing a one-dimensional stochastic process followed by a particle accelerated by an external field and propagating through a thermal bath. The effect of collisions is described within one-dimensional formulation of Boltzmann’s kinetic theory. We present analytical solutions for the Maxwell gas and for the very hard particle model. The exponentially fast relaxation of the velocity distribution towards the stationary form is demonstrated. In the reference frame moving with constant drift velocity the hydrodynamic diffusive mode is shown to govern the distribution in the position space. We show that the exact value of the diffusion coefficient for any value of the field is correctly predicted by the Green-Kubo autocorrelation formula generalized to the stationary state.     

Deterministic single-file dynamics in collisional representation

We re-examine numerically the diffusion of a deterministic, or ballistic single file with preassigned velocity distribution (Jepsen's gas) from a collisional viewpoint. For a two-modal velocity distribution, where half the particles have velocity +/-c, the collisional statistics is analytically proven to reproduce the continuous time representation. For a three-modal velocity distribution with equal fractions, where less than 12 of the particles have velocity +/-c, with the remaining particles at rest, the collisional process is shown to be inhomogeneous; its stationary properties are discussed here by combining exact and phenomenological arguments. Collisional memory effects are then related to the negative power-law tails in the velocity autocorrelation functions, predicted earlier in the continuous time formalism. Numerical and analytical results for Gaussian and four-modal Jepsen's gases are also reported for the sake of a comparison.     

Influence of wall motion on particle sedimentation using hybrid LB-IBM scheme

We integrate the lattice Boltzmann method (LBM) and immersed boundary method (IBM) to capture the coupling between a rigid boundary surface and the hydrodynamic response of an enclosed particle laden fluid. We focus on a rigid box filled with a Newtonian fluid where the drag force based on the slip velocity at the wall and settling particles induces the interaction. We impose an external harmonic oscillation on the system boundary and found interesting results in the sedimentation behavior. Our results reveal that the sedimentation and particle locations are sensitive to the boundary walls oscillation amplitude and the subsequent changes on the enclosed flow field. Two different particle distribution analyses were performed and showed the presence of an agglomerate structure of particles. Despite the increase in the amplitude of wall motion, the turbulence level of the flow field and distribution of particles are found to be less in quantity compared to the stationary walls. The integrated LBM-IBM methodology promised the prospect of an efficient and accurate dynamic coupling between a non-compliant bounding surface and flow field in a wide-range of systems. Understanding the dynamics of the fluid-filled box can be particularly important in a simulation of particle deposition within biological systems and other engineering applications.     

Effective coupling between two Brownian particles

We use the system-plus-reservoir approach to study the dynamics of a system composed of two independent Brownian particles. We present an extension of the well-known model of a bath of oscillators which is capable of inducing an effective coupling between the two particles depending on the choice made for the spectral function of the bath oscillators. The coupling is nonlinear in the variables of interest, and an exponential dependence on these variables is imposed in order to guarantee the translational invariance of the model if the two particles are not subject to any external potential. The effective equations of motion for the particles are obtained by the Laplace transform method, and, besides recovering all the local dynamical properties for each particle, we end up with an effective interaction potential between them. We explicitly analyze one of its possible forms.     

Relaxation Times for Hamiltonian Systems

Usually, the relaxation times of a gas are estimated in the frame of the Boltzmann equation. In this paper, instead, we deal with the relaxation problem in the frame of the dynamical theory of Hamiltonian systems, in which the definition itself of a relaxation time is an open question. We introduce a lower bound for the relaxation time, and give a general theorem for estimating it. Then we give an application to a concrete model of an interacting gas, in which the lower bound turns out to be of the order of magnitude of the relaxation times observed in dilute gases.     

Phase locking in networks of synaptically coupled McKean relaxation oscillators

We use geometric dynamical systems methods to derive phase equations for networks of weakly connected McKean relaxation oscillators. We derive an explicit formula for the connection function when the oscillators are coupled with chemical synapses modeled as the convolution of some input spike train with an appropriate synaptic kernel. The theory allows the systematic investigation of the way in which a slow recovery variable can interact with synaptic time scales to produce phase-locked solutions in networks of pulse coupled neural relaxation oscillators. The theory is exact in the singular limit that the fast and slow time scales of the neural oscillator become effectively independent. By focusing on a pair of mutually coupled McKean oscillators with alpha function synaptic kernels, we clarify the role that fast and slow synapses of excitatory and inhibitory type can play in producing stable phase-locked rhythms. In particular we show that for fast excitatory synapses there is coexistence of a stable synchronous, a stable anti-synchronous, and one stable asynchronous solution. For slower synapses the anti-synchronous solution can lose stability, whilst for even slower synapses it can regain stability. The case of inhibitory synapses is similar up to a reversal of the stability of solution branches. Using a return-map analysis the case of strong pulsatile coupling is also considered. In this case it is shown that the synchronous solution can co-exist with a continuum of asynchronous states.     

Equilibrium and stationary nonequilibrium states in a chain of colliding harmonic oscillators

Equilibrium and nonequilibrium properties of a chain of colliding harmonic oscillators (ding-dong model) are investigated. Our chain is modeled as harmonically bounded particles that can only interact with neighboring particles by hard-core interaction. Between the collisions, particles are just independent harmonic oscillators. We are especially interested in the stationary nonequilibrium state of the ding-dong model coupled with two stochastic heat reservoirs (not thermostated) at the ends, whose temperature is different. We check the Gallavotti-Cohen fluctuation theorem [G. Gallavoti and E. G. D. Cohen, Phys. Rev. Lett. 74, 2694 (1995)] and also the Evans-Searles identity [D. Evans and D. Searles, Phys. Rev. E. 50, 1994 (1994)] numerically. It is verified that the former theorem is satisfied for this system, although the system is not a thermostated system.     

Noise-induced synchronization and clustering in ensembles of uncoupled limit-cycle oscillators

We study synchronization properties of general uncoupled limit-cycle oscillators driven by common and independent Gaussian white noises. Using phase reduction and averaging methods, we analytically derive the stationary distribution of the phase difference between oscillators for weak noise intensity. We demonstrate that in addition to synchronization, clustering, or more generally coherence, always results from arbitrary initial conditions, irrespective of the details of the oscillators.     

Tagged Particle Diffusion in One-Dimensional Gas with Hamiltonian Dynamics

We consider a one-dimensional gas of hard point particles in a finite box that are in thermal equilibrium and evolving under Hamiltonian dynamics. Tagged particle correlation functions of the middle particle are studied. For the special case where all particles have the same mass, we obtain analytic results for the velocity auto-correlation function in the short time diffusive regime and the long time approach to the saturation value when finite-size effects become relevant. In the case where the masses are unequal, numerical simulations indicate sub-diffusive behaviour with mean square displacement of the tagged particle growing as t/ln(t) with time t. Also various correlation functions, involving the velocity and position of the tagged particle, show damped oscillations at long times that are absent for the equal mass case.     

Particle injection into a chain: decoherence versus relaxation for Hermitian and non‐Hermitian dynamics

A model system for the injection of fermionic particles from filled source sites into an empty chain is investigated. The ensuing dynamics for Hermitian as well as for non‐Hermitian time evolution, where the particles cannot return to the bath sites (quantum ratchet), is studied. A non‐homogeneous hybridization between bath and chain sites permits transient currents in the chain. Non‐interacting particles show decoherence in the thermodynamic limit: the average particle number and the average current density in the chain become stationary for long times, whereas the single‐particle density matrix displays large fluctuations around its mean value. Using the numerical time‐dependent density‐matrix renormalization group (t‐DMRG) method it is demonstrated, on the other hand, that sizable density‐density interactions between the particles introduce relaxation which is by orders of magnitudes faster than the decoherence processes.     

Statistical Theory of Two-Sided Multipactor

We develop a statistical theory of secondary-emission discharge (SED) taking the energy distribution of secondary electrons into account. The theory allows one to describe quantitatively the initial stage of development of a two-sided multipactor. For an arbitrary probability density of normal components of the ejection velocity and an arbitrary distance between the walls enclosing the microwave discharge plasma, we construct an analytical solution for the electron distribution function over transit times. The performed analysis is based on the results of a detailed study of conditions under which an electron reaches the opposite side. With allowance for the spread in thermal velocities, we derive a recurrence relation between the electron distribution functions over emission phases and formulate a general integral equation from which the resulting stationary distribution and the threshold of SED onset are determined.     

Non-Equilibrium Statistical Physics of Currents in Queuing Networks

We consider a stable open queuing network as a steady non-equilibrium system of interacting particles. The network is completely specified by its underlying graphical structure, type of interaction at each node, and the Markovian transition rates between nodes. For such systems, we ask the question “What is the most likely way for large currents to accumulate over time in a network ?”, where time is large compared to the system correlation time scale. We identify two interesting regimes. In the first regime, in which the accumulation of currents over time exceeds the expected value by a small to moderate amount (moderate large deviation), we find that the large-deviation distribution of currents is universal (independent of the interaction details), and there is no long-time and averaged over time accumulation of particles (condensation) at any nodes. In the second regime, in which the accumulation of currents over time exceeds the expected value by a large amount (severe large deviation), we find that the large-deviation current distribution is sensitive to interaction details, and there is a long-time accumulation of particles (condensation) at some nodes. The transition between the two regimes can be described as a dynamical second order phase transition. We illustrate these ideas using the simple, yet non-trivial, example of a single node with feedback.     

Stable regimes for hard disks in a channel with twisting walls

We study a gas of N hard disks in a box with semi-periodic boundary conditions. The unperturbed gas is hyperbolic and ergodic (these facts are proved for N=2 and expected to be true for all N≥2). We study various perturbations by twisting the outgoing velocity at collisions with the walls. We show that the dynamics tends to collapse to various stable regimes, however we define the perturbations, and however small they are.     

Turbulence and coarsening in active and passive binary mixtures

Phase separation between two fluids in two dimensions is investigated by means of direct numerical simulations of coupled Navier-Stokes and Cahn-Hilliard equations. We study the phase ordering process in the presence of an external stirring acting on the velocity field. For both active and passive mixtures we find that, for a sufficiently strong stirring, coarsening is arrested in a stationary dynamical state characterized by a continuous rupture and formation of finite domains. Coarsening arrest is shown to be independent of the chaotic or regular nature of the flow.     

Stationary nonequilibrium states in boundary-driven Hamiltonian systems: Shear flow

We investigate stationary nonequilibrium states of systems of particles moving according to Hamiltonian dynamics with specified potentials. The systems are driven away from equilibrium by Maxwell-demon reflection rules at the walls. These deterministic rules conserve energy but not phase space volume, and the resulting global dynamics may or may not be time reversible (or even invertible). Using rules designed to simulate moving walls, we can obtain a stationary shear flow. Assuming that for macroscopic systems this flow satisfies the Navier-Stokes equations, we compare the hydrodynamic entropy production with the average rate of phase-space volume compression. We find that they are equalwhen the velocity distribution of particles incident on the walls is a local Maxwellian. An argument for a general equality of this kind, based on the assumption of local thermodynamic equilibrium, is given. Molecular dynamic simulations of hard disks in a channel produce a steady shear flow with the predicted behavior.     

Polydispersity and Optimal Relaxation in the Hard Sphere Fluid

We consider the mass heterogeneity in a gas of polydisperse hard particles as a key to optimizing a dynamical property: the kinetic relaxation rate. Using the framework of the Boltzmann equation, we study the long time approach of a perturbed velocity distribution toward the equilibrium Maxwellian solution. We work out the cases of discrete as well as continuous distributions of masses, as found in dilute fluids of mesoscopic particles such as granular matter and colloids. On the basis of analytical and numerical evidence, we formulate a dynamical equipartition principle that leads to the result that no such continuous dispersion in fact minimizes the relaxation time, as the global optimum is characterized by a finite number of species. This optimal mixture is found to depend on the dimension $d$ of space, ranging from five species for $d=1$ to a single one for $d\ge 4$ . The role of the collisional kernel is also discussed, and extensions to dissipative systems are shown to be possible.