Source author record

Emmanuel Trizac

Emmanuel Trizac appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

31works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

31 published item(s)

preprint2023arXiv

Control of Active Brownian Particles: An exact solution

Control of stochastic systems is a challenging open problem in statistical physics, with potential applications in a wealth of systems from biology to granulates. Unlike most cases investigated so far, we aim here at controlling a genuinely out-of-equilibrium system, the two dimensional Active Brownian Particles model in a harmonic potential, a paradigm for the study of self-propelled bacteria. We search for protocols for the driving parameters (stiffness of the potential and activity of the particles) bringing the system from an initial passive-like stationary state to a final active-like one, within a chosen time interval. The exact analytical results found for this prototypical system of self-propelled particles brings control techniques to a wider class of out-of-equilibrium systems.

preprint2023arXiv

Metropolis Monte Carlo sampling: convergence, localization transition and optimality

Among random sampling methods, Markov Chain Monte Carlo algorithms are foremost. Using a combination of analytical and numerical approaches, we study their convergence properties towards the steady state, within a random walk Metropolis scheme. Analysing the relaxation properties of some model algorithms sufficiently simple to enable analytic progress, we show that the deviations from the target steady-state distribution can feature a localization transition as a function of the characteristic length of the attempted jumps defining the random walk. While the iteration of the Monte Carlo algorithm converges to equilibrium for all choices of jump parameters, the localization transition changes drastically the asymptotic shape of the difference between the probability distribution reached after a finite number of steps of the algorithm and the target equilibrium distribution. We argue that the relaxation before and after the localisation transition is respectively limited by diffusion and rejection rates.

preprint2022arXiv

Shortcuts To Adiabaticity for Lévy processes in harmonic traps

Lévy stochastic processes, with noise distributed according to a Lévy stable distribution, are ubiquitous in science. Focusing on the case of a particle trapped in an external harmonic potential, we address the problem of finding "shortcuts to adiabaticity": after the system is prepared in a given initial stationary state, we search for time-dependent protocols for the driving external potential, such that a given final state is reached in a given, finite time. These techniques, usually used for stochastic processes with additive Gaussian noise, are typically based on a inverse-engineering approach. We generalise the approach to the wider class of Lévy stochastic processes, both in the overdamped and in the underdamped regime, by finding exact equations for the relevant characteristic functions in Fourier space.

preprint2020arXiv

Electroosmosis as a probe for electrostatic correlations

We study the role of ionic correlations on the electroosmotic flow in planar double-slit channels, without salt. We propose an analytical theory, based on recent advances in the understanding of correlated systems. We compare the theory with mean-field results and validate it by means of dissipative particle dynamics simulations. Interestingly, for some surface separations, correlated systems exhibit a larger flow than predicted by mean-field. We conclude that the electroosmotic properties of a charged system can be used, in general, to infer and weight the importance of electrostatic correlations therein.

preprint2019arXiv

Lattice Boltzmann Electrokinetics simulation of nanocapacitors

We propose a method to model metallic surfaces in Lattice Boltzmann Electrokinetics simulations (LBE), a lattice-based algorithm rooted in kinetic theory which captures the coupled solvent and ion dynamics in electrolyte solutions. This is achieved by a simple rule to impose electrostatic boundary conditions, in a consistent way with the location of the hydrodynamic interface for stick boundary conditions. The proposed method also provides the local charge induced on the electrode by the instantaneous distribution of ions under voltage. We validate it in the low voltage regime by comparison with analytical results in two model nanocapacitors: parallel plate and coaxial electrodes. We examine the steady-state ionic concentrations and electric potential profiles (and corresponding capacitance), the time-dependent response of the charge on the electrodes, as well as the steady-state electro-osmotic profiles in the presence of an additional, tangential electric field. The LBE method further provides the time-dependence of these quantities, as illustrated on the electro-osmotic response. While we do not consider this case in the present work, which focuses on the validation of the method, the latter readily applies to large voltages between the electrodes, as well as to time-dependent voltages. This work opens the way to the LBE simulation of more complex systems involving electrodes and metallic surfaces, such as sensing devices based on nanofluidic channels and nanotubes, or porous electrodes.

preprint2019arXiv

Optimizing Brownian escape rates by potential shaping

Brownian escape is key to a wealth of physico-chemical processes, including polymer folding, and information storage. The frequency of thermally activated energy barrier crossings is assumed to generally decrease exponentially with increasing barrier height. Here, we show experimentally that higher, fine-tuned barrier profiles result in significantly enhanced escape rates in breach of the intuition relying on the above scaling law, and address in theory the corresponding conditions for maximum speed-up. Importantly, our barriers end on the same energy on which they start. For overdamped dynamics, the achievable boost of escape rates is, in principle, unbounded so that the barrier optimization has to be regularized. We derive optimal profiles under two different regularizations, and uncover the efficiency of N-shaped barriers. We then demonstrate the viability of such a potential in automated microfluidic Brownian dynamics experiments using holographic optical tweezers and achieve a doubling of escape rates compared to unhindered Brownian motion. Finally, we show that this escape rate boost extends into the low-friction inertial regime.

preprint2016arXiv

Engineered swift equilibration of a Brownian particle

A fundamental and intrinsic property of any device or natural system is its relaxation time relax, which is the time it takes to return to equilibrium after the sudden change of a control parameter [1]. Reducing $tau$ relax , is frequently necessary, and is often obtained by a complex feedback process. To overcome the limitations of such an approach, alternative methods based on driving have been recently demonstrated [2, 3], for isolated quantum and classical systems [4--9]. Their extension to open systems in contact with a thermostat is a stumbling block for applications. Here, we design a protocol,named Engineered Swift Equilibration (ESE), that shortcuts time-consuming relaxations, and we apply it to a Brownian particle trapped in an optical potential whose properties can be controlled in time. We implement the process experimentally, showing that it allows the system to reach equilibrium times faster than the natural equilibration rate. We also estimate the increase of the dissipated energy needed to get such a time reduction. The method paves the way for applications in micro and nano devices, where the reduction of operation time represents as substantial a challenge as miniaturization [10]. The concepts of equilibrium and of transformations from an equilibrium state to another, are cornerstones of thermodynamics. A textbook illustration is provided by the expansion of a gas, starting at equilibrium and expanding to reach a new equilibrium in a larger vessel. This operation can be performed either very slowly by a piston, without dissipating energy into the environment, or alternatively quickly, letting the piston freely move to reach the new volume.

preprint2016arXiv

Fast equilibrium switch of a micro mechanical oscillator

We demonstrate an accurate method to control the motion of a micromechanical oscillator in contact with a thermal bath. The experiment is carried out on the cantilever tip of an Atomic Force Microscope (AFM). Applying an appropriate time dependent external force, we decrease the time necessary to reach equilibrium by two orders of magnitude compared to the intrinsic equilibration time. Finally, we analyze the energetic cost of such a fast equilibration, by measuring with kBT accuracy the energy exchanges along the process.

preprint2016arXiv

Microscopic origin of self-similarity in granular blast waves

The self-similar expansion of a blast wave, well-studied in air, has peculiar counterparts in dense and dissipative media such as granular gases. Recent results have shown that, while the traditional Taylor-von Neumann-Sedov (TvNS) derivation is not applicable to such granular blasts, they can nevertheless be well understood via a combination of microscopic and hydrodynamic insights. In this article, we provide a detailed analysis of these methods associating Molecular Dynamics simulations and continuum equations, which successfully predict hydrodynamic profiles, scaling properties and the instability of the self-similar solution. We also present new results for the energy conserving case, including the particle-level analysis of the classic TvNS solution and its breakdown at higher densities.

preprint2016arXiv

Renormalized jellium model for colloidal mixtures

In an attempt to quantify the role of polydispersity in colloidal suspensions, we present an efficient implementation of the renormalized jellium model for a mixture of spherical charged colloids. The different species may have different size, charge and density. Advantage is taken from the fact that the electric potential pertaining to a given species obeys a Poisson's equation that is species independent; only boundary conditions do change from a species to the next. All species are coupled through the renormalized background (jellium) density, that is determined self-consistently. The corresponding predictions are compared to the results of Monte Carlo simulations of binary mixtures, where Coulombic interactions are accounted for exactly, at the primitive model level (structureless solvent with fixed dielectric permittivity). An excellent agreement is found.

preprint2016arXiv

Tracer diffusion coefficients in a sheared inelastic Maxwell gas

We study the transport properties of an impurity in a sheared granular gas, in the framework of the Boltzmann equation for inelastic Maxwell models. We investigate here the impact of a nonequilibrium phase transition found in such systems, where the tracer species carries a finite fraction of the total kinetic energy (ordered phase). To this end, the diffusion coefficients are first obtained for a granular binary mixture in spatially inhomogeneous states close to the simple shear flow. In this situation, the set of coupled Boltzmann equations are solved by means of a Chapman-Enskog-like expansion around the (local) shear flow distributions for each species, thereby retaining all the hydrodynamic orders in the shear rate $a$. Due to the anisotropy induced by the shear flow, three tensorial quantities $D_{ij}$, $D_{p,ij}$, and $D_{T,ij}$ are required to describe the mass transport process instead of the conventional scalar coefficients. These tensors are given in terms of the solutions of a set of coupled algebraic equations, which can be \emph{exactly} solved as functions of the shear rate $a$, the coefficients of restitution $α_{sr}$ and the parameters of the mixture (masses and composition). Once the forms of $D_{ij}$, $D_{p,ij}$, and $D_{T,ij}$ are obtained for arbitrary mole fraction $x_1=n_1/(n_1+n_2)$ (where $n_r$ is the number density of species $r$), the tracer limit ($x_1\to 0$) is carefully considered for the above three diffusion tensors. Explicit forms for these coefficients are derived showing that their shear rate dependence is significantly affected by the order-disorder transition.

preprint2015arXiv

Anomalous transport of impurities in inelastic Maxwell gases

A mixture of dissipative hard grains generically exhibits a breakdown of kinetic energy equipartition. The undriven and thus freely cooling binary problem, in the tracer limit where the density of one species becomes minute, may exhibit an extreme form of this breakdown, with the minority species carrying a finite fraction of the total kinetic energy of the system. We investigate the fingerprint of this non-equilibrium phase transition, akin to an ordering process, on transport properties. The analysis, performed by solving the Boltzmann kinetic equation from a combination of analytical and Monte Carlo techniques, hints at the possible failure of hydrodynamics in the ordered region. As a relevant byproduct of the study, the behaviour of the second and fourth-degree velocity moments is also worked out.

preprint2015arXiv

Blast dynamics in a dissipative gas

The blast caused by an intense explosion has been extensively studied in conservative fluids, where the Taylor-von Neumann-Sedov hydrodynamic solution is a prototypical example of self-similarity driven by conservation laws. In dissipative media however, energy conservation is violated, yet a distinctive self-similar solution appears. It hinges on the decoupling of random and coherent motion permitted by a broad class of dissipative mechanisms. This enforces a peculiar layered structure in the shock, for which we derive the full hydrodynamic solution, validated by a microscopic approach based on Molecular Dynamics simulations. We predict and evidence a succession of temporal regimes, as well as a long-time corrugation instability, also self-similar, which disrupts the blast boundary. These generic results may apply from astrophysical systems to granular gases, and invite further cross-fertilization between microscopic and hydrodynamic approaches of shockwaves.

preprint2015arXiv

Generalized transport coefficients for inelastic Maxwell mixtures under shear flow

The Boltzmann equation framework for inelastic Maxwell models is considered to determine the transport coefficients associated with the mass, momentum and heat fluxes of a granular binary mixture in spatially inhomogeneous states close to the simple shear flow. The Boltzmann equation is solved by means of a Chapman-Enskog-like expansion around the (local) shear flow distributions $f_r^{(0)}$ for each species that retain all the hydrodynamic orders in the shear rate. Due to the anisotropy induced by the shear flow, tensorial quantities are required to describe the transport processes instead of the conventional scalar coefficients. These tensors are given in terms of the solutions of a set of coupled equations, which can be analytically solved as functions of the shear rate $a$, the coefficients of restitution $α_{rs}$ and the parameters of the mixture (masses, diameters and composition). Since the reference distribution functions $f_r^{(0)}$ apply for arbitrary values of the shear rate and are not restricted to weak dissipation, the corresponding generalized coefficients turn out to be nonlinear functions of both $a$ and $α_{rs}$. The dependence of the relevant elements of the three diffusion tensors on both the shear rate and dissipation is illustrated in the tracer limit case, the results showing that the deviation of the generalized transport coefficients from their forms for vanishing shear rates is in general significant. A comparison with the previous results obtained analytically for inelastic hard spheres by using Grad's moment method is carried out showing a good agreement over a wide range of values for the coefficients of restitution. Finally, as an application of the theoretical expressions derived here for the transport coefficients, thermal diffusion segregation of an intruder immersed in a granular gas is also studied.

preprint2015arXiv

Screening like-charges in one-dimensional Coulomb systems: Exact results

The possibility that like-charges can attract each other under the mediation of mobile counterions is by now well documented experimentally, numerically, and analytically. Yet, obtaining exact results is in general impossible, or restricted to some limiting cases. We work out here in detail a one dimensional model that retains the essence of the phenomena present in higher dimensional systems. The partition function is obtained explicitly, from which a wealth of relevant quantities follow, such as the effective force between the charges or the counterion profile in their vicinity. Isobaric and canonical ensembles are distinguished. The case of two equal charges screened by an arbitrary number $N$ of counterions is first studied, before the more general asymmetric situation is addressed. It is shown that the parity of $N$ plays a key role in the long range physics.

preprint2014arXiv

Interplay of anisotropy in shape and interactions in charged platelet suspensions

Motivated by the intriguing phase behavior of charged colloidal platelets, we investigate the structure and dynamics of charged repulsive disks by means of Monte-Carlo simulations. The electrostatic interactions are taken into account through an effective two-body potential, obtained within the non-linear Poisson-Boltzmann formalism, which has the form of anisotropic screened Coulomb potential. Recently, we showed that the original intrinsic anisotropy of the electrostatic potential in competition with excluded volume effects leads to a rich phase behavior that not only includes various liquid-crsytalline phases but also predicts the existence of novel structures composed of alternating nematic-antinematic sheets. Here, we examine the structural and dynamical signatures of each of the observed structures for both translational and rotational degrees of freedom. Finally, we discuss the influence of effective charge value and our results in relation to experimental findings on charged platelet suspensions.

preprint2014arXiv

Taking one charge off a two-dimensional Wigner crystal

A planar array of identical charges at vanishing temperature forms a Wigner crystal with hexagonal symmetry. We take off one (reference) charge in a perpendicular direction, hold it fixed, and search for the ground state of the whole system. The planar projection of the reference charge should then evolve from a six-fold coordination (center of a hexagon) for small distances to a three-fold arrangement (center of a triangle), at large distances $d$ from the plane. The aim of this paper is to describe the corresponding non-trivial lattice transformation. For that purpose, two numerical methods (direct energy minimization and Monte Carlo simulations), together with an analytical treatment, are presented. Our results indicate that the $d=0$ and $d\to\infty$ limiting cases extend for finite values of $d$ from the respective starting points into two sequences of stable states, with intersecting energies at some value $d_t$; beyond this value the branches continue as metastable states.

preprint2014arXiv

The contact theorem for charged fluids: from planar to curved geometries

When a Coulombic fluid is confined between two parallel charged plates, an exact relation links the difference of ionic densities at contact with the plates, to the surface charges of these boundaries. It no longer applies when the boundaries are curved, and we work out how it generalizes when the fluid is confined between two concentric spheres (or cylinders), in two and in three space dimensions. The analysis is thus performed within the cell model picture. The generalized contact relation opens the possibility to derive new exact expressions, of particular interest in the regime of strong coulombic couplings. Some emphasis is put on cylindrical geometry, for which we discuss in depth the phenomenon of counter-ion evaporation/condensation, and obtain novel results. Good agreement is found with Monte Carlo simulation data.

preprint2014arXiv

Thinking outside the box: fluctuations and finite size effects

The isothermal compressibility of an interacting or non interacting system may be extracted from the fluctuations of the number of particles in a well chosen control volume. Finite size effects are prevalent and should then be accounted for to obtain a meaningful, thermodynamic compressibility. In the traditional computational setup where a given simulation box is replicated with periodic boundary conditions, we study particle number fluctuations outside the box (i.e. when the control volume exceeds the box itself), which bear relevant thermodynamic information. We also investigate the related problem of extracting the compressibility from the structure factor in the small wave-vector limit ($k\to 0$). The calculation should be restricted to the discrete set of wave-vectors $k$ that are compatible with the periodicity of the system, and we assess the consequences of considering other $k$ values, a widespread error among beginners.

preprint2013arXiv

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>=4. The role of the collisional kernel is also discussed, and extensions to dissipative systems are shown to be possible.

preprint2012arXiv

Counter-ion density profile around charged cylinders: the strong-coupling needle limit

Charged rod-like polymers are not able to bind all their neutralizing counter-ions: a fraction of them evaporates while the others are said to be condensed. We study here counter-ion condensation and its ramifications, both numerically by means of Monte Carlo simulations employing a previously introduced powerful logarithmic sampling of radial coordinates, and analytically, with special emphasis on the strong-coupling regime. We focus on the thin rod, or needle limit, that is naturally reached under strong coulombic couplings, where the typical inter-particle spacing $a'$ along the rod is much larger than its radius R. This regime is complementary and opposite to the simpler thick rod case where $a'\ll R$. We show that due account of counter-ion evaporation, a universal phenomenon in the sense that it occurs in the same clothing for both weakly and strongly coupled systems, allows to obtain excellent agreement between the numerical simulations and the strong-coupling calculations.

preprint2012arXiv

Dynamic Monte Carlo Simulations of Anisotropic Colloids

We put forward a simple procedure for extracting dynamical information from Monte Carlo simulations, by appropriate matching of the short-time diffusion tensor with its infinite-dilution limit counterpart, which is supposed to be known. This approach --discarding hydrodynamics interactions-- first allows us to improve the efficiency of previous Dynamic Monte Carlo algorithms for spherical Brownian particles. In a second step, we address the case of anisotropic colloids with orientational degrees of freedom. As an illustration, we present a detailed study of the dynamics of thin platelets, with emphasis on long-time diffusion and orientational correlations.

preprint2012arXiv

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.

preprint2012arXiv

Phase behaviour of colloidal assemblies on 2D corrugated substrates

We investigate - with Monte Carlo computer simulations - the phase behaviour of dimeric colloidal molecules on periodic substrates with square symmetry. The molecules are formed in a two-dimensional suspension of like charged colloids subject to periodic external confinement, which can be experimentally realized by optical methods. We study the evolution of positional and orientational order by varying the temperature across the melting transition. We propose and evaluate appropriate order parameters as well as the specific heat capacity and show that the decay of positional correlations belongs to a class of crossover transitions while the orientational melting is a second-order phase transition.

preprint2012arXiv

Self-assembly of spherical interpolyelectrolyte complexes from oppositely charged polymers

The formation of inter-polyelectrolyte complexes from the association of oppositely charged polymers in an electrolyte is studied. The charged polymers are linear oppositely charged polyelectrolytes, with possibly a neutral block. This leads to complexes with a charged core, and a more dilute corona of dangling chains, or of loops (flower-like structure). The equilibrium aggregation number of the complexes (number of polycations m+ and polyanions m-) is determined by minimizing the relevant free energy functional, the Coulombic contribution of which is worked out within Poisson-Boltzmann theory. The complexes can be viewed as colloids that are permeable to micro-ionic species, including salt. We find that the complexation process can be highly specific, giving rise to very localized size distribution in composition space (m+,m-).

preprint2011arXiv

Colloidal ionic complexes on periodic substrates: ground state configurations and pattern switching

We theoretically and numerically studied ordering of "colloidal ionic clusters" on periodic substrate potentials as those generated by optical trapping. Each cluster consists of three charged spherical colloids: two negatively and one positively charged. The substrate is a square or rectangular array of traps, each confining one such cluster. By varying the lattice constant from large to small, the observed clusters are first rod-like and form ferro- and antiferro-like phases, then they bend into a banana-like shape and finally condense into a percolated structure. Remarkably, in a broad parameter range between single-cluster and percolated structures, we have found stable supercomplexes composed of six colloids forming grape-like or rocket-like structures. We investigated the possibility of macroscopic pattern switching by applying external electrical fields.

preprint2011arXiv

Interaction regimes for oppositely charged plates with multivalent counterions

Within a mean field treatment of the interaction between two oppositely charged plates in a salt free solution, the distance at which a transition from an attractive to a repulsive regime appears can be computed analytically. The mean field description however breaks down under strong coulombic couplings, that can be achieved at room temperature with multivalent counter-ions and highly charged surfaces. Making use of the contact theorem and simple physical arguments, we propose explicit expressions for the equation of state in several situations at short distances. The possibility of Bjerrum pair formation is addressed and is shown to have profound consequences on the interactions. To complete the picture, we finally consider the large distance limit, from which schematic phase diagrams discriminating attractive from repulsive regions can be proposed.

preprint2009arXiv

Rheological properties for inelastic Maxwell mixtures under shear flow

The Boltzmann equation for inelastic Maxwell models is considered to determine the rheological properties in a granular binary mixture in the simple shear flow state. The transport coefficients (shear viscosity and viscometric functions) are {\em exactly} evaluated in terms of the coefficients of restitution, the (reduced) shear rate and the parameters of the mixture (particle masses, diameters and concentration). The results show that in general, for a given value of the coefficients of restitution, the above transport properties decrease with increasing shear rate.

preprint2005arXiv

Dynamics of a tracer granular particle as a non-equilibrium Markov process

The dynamics of a tracer particle in a stationary driven granular gas is investigated. We show how to transform the linear Boltzmann equation describing the dynamics of the tracer into a master equation for a continuous Markov process. The transition rates depend upon the stationary velocity distribution of the gas. When the gas has a Gaussian velocity probability distribution function (pdf), the stationary velocity pdf of the tracer is Gaussian with a lower temperature and satisfies detailed balance for any value of the restitution coefficient $α$. As soon as the velocity pdf of the gas departs from the Gaussian form, detailed balance is violated. This non-equilibrium state can be characterized in terms of a Lebowitz-Spohn action functional $W(τ)$ defined over trajectories of time duration $τ$. We discuss the properties of this functional and of a similar functional $\bar{W}(τ)$ which differs from the first for a term which is non-extensive in time. On the one hand we show that in numerical experiments, i.e. at finite times $τ$, the two functionals have different fluctuations and $\bar{W}$ always satisfies an Evans-Searles-like symmetry. On the other hand we cannot observe the verification of the Lebowitz-Spohn-Gallavotti-Cohen (LS-GC) relation, which is expected for $W(τ)$ at very large times $τ$. We give an argument for the possible failure of the LS-GC relation in this situation. We also suggest practical recipes for measuring $W(τ)$ and $\bar{W}(τ)$ in experiments.

preprint2004arXiv

Hydrodynamics within the Electric Double Layer on slipping surfaces

We show, using extensive Molecular Dynamics simulations, that the dynamics of the electric double layer (EDL) is very much dependent on the wettability of the charged surface on which the EDL develops. For a wetting surface, the dynamics, characterized by the so-called Zeta potential, is mainly controlled by the electric properties of the surface, and our work provides a clear interpretation for the traditionally introduced immobile Stern layer. In contrast, the immobile layer disappears for non-wetting surfaces and the Zeta potential deduced from electrokinetic effects is considerably amplified by the existence of a slippage at the solid substrate.

preprint2003arXiv

Kinetics of Fragmenting Freely Evolving Granular Gases

We investigate the effect of fragmentation on the homogeneous free cooling of inelastic hard spheres, using Boltzmann kinetic theory and Direct Monte Carlo simulations. We analyze in detail a model where dissipative collisions may subsequently lead to a break-up of the grains. With a given probability, two off-springs are then created from one of the two colliding partners, with conservation of mass, momentum and kinetic energy. We observe a scaling regime characterized by a single collisional average, that quantifies the deviations from Gaussian behaviour for the joint size and velocity distribution function. We also discuss the possibility of a catastrophe whereby the number of particles diverges in a finite time. This phenomenon appears correlated to a ``shattering'' transition marked by a delta singularity at vanishingly small grains for the rescaled size distribution.