BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Umberto Marini Bettolo Marconi

Umberto Marini Bettolo Marconi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechcond-mat.softcond-mat.mes-hallBiological PhysicsBiomoleculesphysics.chem-phphysics.comp-phphysics.flu-dyn

Catalog footprint

What is connected

34works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

34 published item(s)

preprint2022arXiv

Active particles driven by competing spatially dependent self-propulsion and external force

We investigate how the competing presence of a nonuniform motility landscape and an external confining field affects the properties of active particles. We employ the active Ornstein-Uhlenbeck particle (AOUP) model with a periodic swim velocity profile to derive analytical approximations for the steady-state probability distribution of position and velocity, encompassing both the Unified Colored Noise Approximation and the theory of potential-free active particles with spatially dependent swim velocity recently developed. We test the theory by confining an active particle in a harmonic trap, which gives rise to interesting properties, such as a transition from a unimodal to a bimodal (and, eventually multimodal) spatial density, induced by decreasing the spatial period of the self propulsion. Correspondingly, the velocity distribution shows pronounced deviations from the Gaussian shape, even displaying a bimodal profile in the high-motility regions. Our results can be confirmed by real-space experiments on active colloidal Janus particles in external fields.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Collective effects in confined Active Brownian Particles

We investigate a two-dimensional system of active particles confined to a narrow annular domain. Despite the absence of explicit interactions among the velocities or the active forces of different particles, the system displays a transition from a disordered and stuck state to an ordered state of global collective motion where the particles rotate persistently clockwise or anticlockwise. We describe this behavior by introducing a suitable order parameter, the velocity polarization, measuring the global alignment of the particles' velocities along the tangential direction of the ring. We also measure the spatial velocity correlation function and its correlation length to characterize the two states. In the rotating phase, the velocity correlation displays an algebraic decay that is analytically predicted together with its correlation length while in the stuck regime the velocity correlation decays exponentially with a correlation length that increases with the persistence time. In the first case, the correlation (and, in particular, its correlation length) does not depend on the active force but the system size only. The global collective motion, an effect caused by the interplay between finite-size, periodicity, and persistent active forces, disappears as the size of the ring becomes infinite, suggesting that this phenomenon does not correspond to a phase transition in the usual thermodynamic sense.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Correlated escape of active particles across a potential barrier

We study the dynamics of one-dimensional active particles confined in a double-well potential, focusing on the escape properties of the system, such as the mean escape time from a well. We first consider a single-particle both in near and far-from-equilibrium regimes by varying the persistence time of the active force and the swim velocity. A non-monotonic behavior of the mean escape time is observed with the persistence time of the activity, revealing the existence of an optimal choice of the parameters favoring the escape process. For small persistence times, a Kramers-like formula with an effective potential obtained within the Unified Colored Noise Approximation is shown to hold. Instead, for large persistence times, we developed a simple theoretical argument based on the first passage theory which explains the linear dependence on the escape time with the persistence of the active force. In the second part of the work, we consider the escape of two active particles mutually repelling. Interestingly, the subtle interplay of active and repulsive forces may lead to a correlation between particles favoring the simultaneous jump across the barrier. This mechanism cannot be observed in the escape process of two passive particles. Finally, we find that, in the small-persistence regime, the repulsion favors the escape, like in passive systems, in agreement with our theoretical predictions, while for large persistence times, the repulsive and active forces produce an effective attraction which hinders the barrier crossing.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Dynamics of active particles with space-dependent swim velocity

We study the dynamical properties of an active particle subject to a swimming speed explicitly depending on the particle position. The oscillating spatial profile of the swim velocity considered in this paper takes inspiration from experimental studies based on Janus particles whose speed can be modulated by an external source of light. We suggest and apply an appropriate model of an active Ornstein Uhlenbeck particle to the present case. This allows us to predict the stationary properties, by finding the exact solution of the steady-state probability distribution of particle position and velocity. From this, we obtain the spatial density profile and show that its form is consistent with the one found in the framework of other popular models. The reduced velocity distribution highlights the emergence of non-Gaussianity in our generalized AOUP model which becomes more evident as the spatial dependence of the velocity profile becomes more pronounced. Then, we focus on the time-dependent properties of the system. Velocity autocorrelation functions are studied in the steady-state combining numerical and analytical methods derived under suitable approximations. We observe a non-monotonic decay in the temporal shape of the velocity autocorrelation function which depends on the ratio between the persistence length and the spatial period of the swim velocity. Finally, we numerically and analytically study the mean square displacement and the long-time diffusion coefficient. The ballistic regime, observed in the small-time region, is deeply affected by the properties of the swim velocity landscape which induces also a crossover to a sub-ballistic but superdiffusive regime for intermediate times. Finally, the long-time diffusion coefficient decreases as the amplitude of the swim velocity oscillations increases because the diffusion is determined by those regions where the particles are slow.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Hydrodynamics of simple active liquids: the emergence of velocity correlations

We derive the Hydrodynamics for a system of N active, spherical, underdamped particles, interacting through conservative forces. At the microscopic level, we represent the evolution of the particles in terms of the Kramers equation for the probability density distribution of their positions, velocities, and orientations, while at a mesoscopic level we switch to a coarse-grained description introducing an appropriate set of hydrodynamic fields given by the lower-order moments of the distribution. In addition to the usual density and polarization fields, the Hydrodynamics developed in this paper takes into account the velocity and kinetic temperature fields, which are crucial to understanding new aspects of the behavior of active liquids. By imposing a suitable closure of the hydrodynamic moment equations and truncation of the Born-Bogolubov-Green-Kirkwood-Yvon hierarchy, we obtain a closed set of mesoscopic balance equations. At this stage, we focus our interest on the small deviations of the hydrodynamic fields from their averages and apply the methods of the theory of linear hydrodynamic fluctuations. Our treatment sheds light on the peculiar properties of isotropic active liquids and their emergent dynamical collective phenomena, such as the spontaneous alignment of the particle velocities. We predict the existence within the liquid phase of spatial equal-time Ornstein-Zernike-like velocity correlations both for the longitudinal and the transverse modes. At variance with active solids, in active liquids, the correlation length of the transverse velocity fluctuations is sensibly shorter than the length of the longitudinal fluctuations. In particular, the latter depends on the sound speed and increases with the persistence time, while the former displays a weaker dependence on these parameters. Finally, within the same framework, we derive the dynamical structure factors ...

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Active Matter at high density: velocity distribution and kinetic temperature

We consider the solid or hexatic non-equilibrium phases of an interacting two-dimensional system of Active Brownian Particles at high density and investigate numerically and theoretically the properties of the velocity distribution function and the associated kinetic temperature. We obtain approximate analytical predictions for the shape of the velocity distribution and find a transition from a Mexican-hat-like to a Gaussian-like distribution as the persistence time of the active force changes from the small to the large persistence regime. Through a detailed numerical and theoretical analysis of the single-particle velocity variance, we report an exact analytical expression for the kinetic temperature of dense spherical self-propelled particles that holds also in the non-equilibrium regimes with large persistence times and discuss its range of validity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Activity-controlled clogging and unclogging of micro-channels

We propose a mechanism to control the formation of stable obstructions in two-dimensional micro-channels of variable sections taking advantage of the peculiar clustering property of active systems. Under the activation of the self-propulsion by external stimuli, the system behaves as a switch according to the following principle: by turning-on the self-propulsion the particles become active and even at very low densities stick to the walls and form growing layers eventually blocking the channel bottleneck, while the obstruction dissolves when the self-propulsion is turned off. We construct the phase diagram distinguishing clogged and open states in terms of density and bottleneck width. The study of the average clogging time, as a function of density and bottleneck width, reveals the marked efficiency of the active clogging that swiftly responds to the self-propulsion turning on. The resulting picture shows a profound difference with respect to the clogging obtained through the slow diffusive dynamics of attractive passive Brownian disks. This numerical work suggests a novel method to use particles with externally tunable self-propulsion to create or destroy plugs in micro-channels.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Hidden velocity ordering in dense suspensions of self-propelled disks

Recent studies of the phase diagram for spherical, purely repulsive, active particles established the existence of a transition from a liquid-like to a solid-like phase analogous to the one observed in colloidal systems at thermal equilibrium, in particular in two dimensions an intermediate hexatic phase is observed. Here, we present evidence that the active dense phases (solid, hexatic and liquid) exhibit interesting dynamical anomalies. First, we unveil the growth - with density and activity - of ordered domains where the particles' velocities align in parallel or vortex-like patches. Second, when activity is strong, the spatial distribution of kinetic energy becomes heterogeneous with high energy regions correlated to defects of the crystalline structure. This spatial heterogeneity is accompanied by temporal intermittency, with sudden peaks in the time-series of kinetic energy. The observed dynamical anomalies are not present in a dense equilibrium system and cannot be detected by considering only the structural properties of the system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Inertial self-propelled particles

We study a self-propelled particle moving in a solvent with the active Ornstein Uhlenbeck dynamics in the underdamped regime to evaluate the influence of the inertia. We focus on the properties of potential-free and harmonically confined underdamped active particles, studying how the single-particle trajectories modify for different values of the drag coefficient. In both cases, we solve the dynamics in terms of correlation matrices and steady-state probability distribution functions revealing the explicit correlations between velocity and active force. We also evaluate the influence of the inertia on the time-dependent properties of the system, discussing the mean square displacement and the time-correlations of particle positions and velocities. Particular attention is devoted to the study of the Virial active pressure unveiling the role of the inertia on this observable.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Spontaneous velocity alignment in Motility-induced Phase Separation

We study a system of purely repulsive spherical self-propelled particles in the minimal set-up inducing Motility-Induced Phase Separation (MIPS). We show that, even if explicit alignment interactions are absent, a growing order in the velocities of the clustered particles accompanies MIPS. Particles arrange into aligned or vortex-like domains. Their sizes increase as the persistence of the self-propulsion grows, an effect that is quantified studying the spatial correlation function of the velocities. We explain the velocity-alignment by unveiling a hidden alignment interaction of the Vicsek-like form, induced by the interplay between steric interactions and self-propulsion. As a consequence, we argue that the MIPS transition cannot be fully understood in terms of a scalar field, the density, since the collective orientation of the velocities should be included in effective coarse-grained descriptions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Critical Phenomena in Active Matter

We investigate the effect of self-propulsion on a mean-field order-disorder transition. Starting from a $φ^4$ scalar field theory subject to an exponentially correlated noise, we exploit the Unified Colored Noise Approximation to map the non-equilibrium active dynamics onto an effective equilibrium one. This allows us to follow the evolution of the second-order critical point as a function of the noise parameters: the correlation time $τ$ and the noise strength $D$. Our results suggest that $τ$ is a crucial ingredient that changes the location of the critical point but, remarkably, not the universality class of the model. We also estimate the effect of Gaussian fluctuations on the mean-field approximation finding an Ornstein-Zernike like expression for the static structure factor at long wave lengths. Finally, to assess the validity of our predictions, we compare the mean-field theoretical results with numerical simulations of active Lennard-Jones particles in two and three dimensions, finding a good qualitative agreement at small $τ$ values.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Effective potential method for active particles

We investigate the steady state properties of an active fluid modeled as an assembly of soft repulsive spheres subjected to Gaussian colored noise. Such a noise captures one of the salient aspects of active particles, namely the persistence of their motion and determines a variety of novel features with respect to familiar passive fluids. We show that within the so-called multidimensional unified colored noise approximation, recently introduced in the field of active matter, the model can be treated by methods similar to those employed in the study of standard molecular fluids. The system shows a tendency of the particles to aggregate even in the presence of purely repulsive forces because the combined action of colored noise and interactions enhances the the effective friction between nearby particles. We also discuss whether an effective two-body potential approach, which would allow to employ methods similar to those of density functional theory, is appropriate. The limits of such an approximation are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Pressure and surface tension of an active simple liquid: a comparison between kinetic, mechanical and free-energy based approaches

We discuss different definitions of pressure for a system of active spherical particles driven by a non-thermal coloured noise. We show that mechanical, kinetic and free-energy based approaches lead to the same result up to first order in the non-equilibrium expansion parameter. The first prescription is based on a generalisation of the kinetic mesoscopic virial equation and expresses the pressure exerted on the walls in terms of the average of the virial of the inter-particle forces. In the second approach, the pressure and the surface tension are identified with the volume and area derivatives, respectively, of the partition function associated with the known stationary non-equilibrium distribution of the model. The third method is a mechanical approach and is related to the work necessary to deform the system. The pressure is obtained by comparing the expression of the work in terms of local stress and strain with the corresponding expression in terms of microscopic distribution. This is determined from the force balance encoded in the Born-Green-Yvon equation. Such a method has the advantage of giving a formula for the local pressure tensor and the surface tension even in inhomogeneous situations. By direct inspection, we show that the three procedures lead to the same values of the pressure, and give support to the idea that the partition function, obtained via the unified coloured noise approximation, is more than a formal property of the system, but determines the stationary non-equilibrium thermodynamics of the model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Electro-osmotic flow in coated nanocapillaries: a theoretical investigation

Motivated by recent experiments, we present a theoretical investigation of how the electro-osmotic flow occurring in a capillary is modified when its charged surfaces are coated by charged polymers. The theoretical treatment is based on a three dimensional model consisting of a ternary fluid-mixture, representing the solvent and two species for the ions, confined between two parallel charged plates decorated by a fixed array of scatterers representing the polymer coating. The electro-osmotic flow, generated by a constant electric field applied in a direction parallel to the plates, is studied numerically by means of Lattice Boltzmann simulations. In order to gain further understanding we performed a simple theoretical analysis by extending the Stokes-Smoluchowski equation to take into account the porosity induced by the polymers in the region adjacent the walls. We discuss the nature of the velocity profiles by focusing on the competing effects of the polymer charges and the frictional forces they exert. We show evidence of the flow reduction and of the flow inversion phenomenon when the polymer charge is opposite to the surface charge. By using the density of polymers and the surface charge as control variables, we propose a phase diagram that discriminates the direct and the reversed flow regimes and determine its dependence on the ionic concentration.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Kinetic Density Functional Theory: A microscopic approach to fluid mechanics

In the present paper we give a brief summary of some recent theoretical advances in the treatment of inhomogeneous fluids and methods which have applications in the study of dynamical properties of liquids in situations of extreme confinement, such as nanopores, nanodevices, etc. The approach obtained by combining kinetic and density functional methods is microscopic, fully self-consistent and allows to determine both configurational and flow properties of dense fluids. The theory predicts the correct hydrodynamic behavior and provides a practical and numerical tool to determine how the transport properties are modified when the length scales of the confining channels are comparable with the size of the molecules. The applications range from the dynamics of simple fluids under confinement, to that of neutral binary mixtures and electrolytes where the theory in the limit of slow gradients reproduces the known phenomenological equations such as the Planck-Nernst-Poisson and the Smoluchowski equations. The approach here illustrated allows for fast numerical solution of the evolution equations for the one-particle phase-space distributions by means of the weighted density lattice Boltzmann method and is particularly useful when one considers flows in complex geometries.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Lattice Boltzmann Method for mixtures at variable Schmidt number

When simulating multicomponent mixtures via the Lattice Boltzmann Method, it is desirable to control the mutual diffusivity between species while maintaining the viscosity of the solution fixed. This goal is herein achieved by a modification of the multicomponent Bhatnagar-Gross-Krook (BGK) evolution equations by introducing two different timescales for mass and momentum diffusion. Diffusivity is thus controlled by an effective drag force acting between species. Numerical simulations confirm the accuracy of the method for neutral binary and charged ternary mixtures in bulk conditions. The simulation of a charged mixture in a charged slit channel show that the conductivity and electro-osmotic mobility exhibit a departure from the Helmholtz-Smoluchowski prediction at high diffusivity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Modulation of current through a nanopore induced by a charged globule: implications for DNA-docking

The passage of DNA through a nanopore can be effectively decomposed into two distinct phases, docking and actual translocation. In experiments each phase is characterized by a distinct current signature which allows the discrimination of the two events. However, at low voltages a clear distinction of the two phases is lost. By using numerical simulations we clarify how the current signature associated to the docking events depends on the applied voltage. The simulations show that at small voltage the DNA globule enhances the pore conductance due to an enrichment of charge carriers. At high voltage, the globule drains substantial charge carriers from the pore region, thereby reducing the overall conductance. The results provide a new interpretation to the experimental data on conductance and show how docking interferes with the translocation signal, of potential interest for sequencing applications.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multidimensional Stationary Probability Distribution for Interacting Active Particles

We derive the stationary probability distribution for a non-equilibrium system composed by an arbitrary number of degrees of freedom that are subject to Gaussian colored noise and a conservative potential. This is based on a multidimensional version of the Unified Colored Noise Approximation. By comparing theory with numerical simulations we demonstrate that the theoretical probability density quantitatively describes the accumulation of active particles around repulsive obstacles. In particular, for two particles with repulsive interactions, the probability of close contact decreases when one of the two particle is pinned. Moreover, in the case of isotropic confining potentials, the radial density profile shows a non trivial scaling with radius. Finally we show that the theory well approximates the "pressure" generated by the active particles allowing to derive an equation of state for a system of non-interacting colored noise-driven particles.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Steric modulation of ionic currents in DNA translocation through nanopores

Ionic currents accompanying DNA translocation strongly depend on molarity of the electrolyte solution and the shape and surface charge of the nanopore. By means of the Poisson-Nernst-Planck equations it is shown how conductance is modulated by the presence of the DNA intruder and as a result of competing electrostatic and confinement factors. The theoretical results reproduce quantitatively the experimental ones and are summarized in a conductance diagram that allows distinguishing the region of reduced conductivity from the region of enhanced conductivity as a function of molarity and the pore dimension.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Towards a statistical mechanical theory of active fluids

We present a stochastic description of a model of N mutually repelling active spheres in the presence of external fields and characterize its steady state behavior. To reproduce the effects of the experimentally observed persistence of the trajectories of the active particles we consider a Gaussian forcing having a non vanishing correlation time $τ$, whose finiteness is a measure of the activity of the system. With these ingredients we show that it is possible to develop a statistical mechanical approach similar to the one employed in the study of equilibrium liquids and to obtain the explicit form of the many-particle distribution function by means of the multidimensional unified colored noise approximation. Such a distribution plays a role analogous to the Gibbs distribution in equilibrium statistical mechanics and provides a complete information about the microscopic state of the system. From here we develop a method to determine the one and two-particle distribution functions in the spirit of the Born-Green-Yvon (BGY) equations of equilibrium statistical mechanics. The resulting equations which contain extra-correlations induced by the activity allow to determine the stationary density profiles in the presence of external fields, the pair correlations and the pressure of active fluids. In the low density regime we obtain the effective pair potential $ϕ(r)$ acting between two isolated particles separated by a distance, $r$, showing the existence an effective attraction between them induced by activity. Based on these results, in the second half of the paper we propose a mean field theory as an approach simpler than the BGY hierarchy and use it to derive a van der Waals expression of the equation of state, which can serve as the basis to understand the phase behavior of active fluids.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Tracer diffusion of hard-sphere binary mixtures under nano-confinement

The physics of diffusion phenomena in nano and micro channels has attracted a lot of attention in recent years, due to its close connection with many technological, medical and industrial applications. In the present paper we employ a kinetic approach to investigate how the confinement in nanostructured geometries affects the diffusive properties of fluid mixtures and leads to the appearance of properties different from those of bulk systems. In particular, we derive an expression for the friction tensor in the case of a bulk fluid mixture confined to a narrow slit having undulated walls. The boundary roughness leads to a new mechanism for transverse diffusion, and can even lead to an effective diffusion along the channel larger than the one corresponding to a planar channel of equivalent section. Finally we discuss a reduction of the previous equation to a one dimensional effective diffusion equation in which an entropic term encapsulates the geometrical information on the channel shape.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Velocity distribution in active particles systems

We derive an analytic expression for the distribution of velocities of multiple interacting active particles which we test by numerical simulations. In clear contrast with equilibrium we find that the velocities are coupled to positions. Our model shows that, even for two particles only, the individual velocities display a variance depending on the interparticle separation and the emergence of correlations between the velocities of the particles. When considering systems composed of many particles we find an analytic expression connecting the overall velocity variance to density, at the mean-field level, and to the pair distribution function valid in the limit of small noise correlation times. Finally we discuss the intriguing analogies and main differences between our effective free energy functional and the theoretical scenario proposed so far for phase-separating active particles.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

About an H-theorem for systems with non-conservative interactions

We exhibit some arguments in favour of an H-theorem for a generalization of the Boltzmann equation including non-conservative interactions and a linear Fokker-Planck-like thermostatting term. Such a non-linear equation describing the evolution of the single particle probability $P_i(t)$ of being in state $i$ at time $t$, is a suitable model for granular gases and is indicated here as Boltzmann-Fokker-Planck (BFP) equation. The conjectured H-functional, which appears to be non-increasing, is $H_C(t)=\sum_i P_i(t) \ln P_i(t)/Π_i$ with $Π_i = \lim_{t \to \infty} P_i(t)$, in analogy with the H-functional of Markov processes. The extension to continuous states is straightforward. A simple proof can be given for the elastic BFP equation. A semi-analytical proof is also offered for the BFP equation for so-called inelastic Maxwell molecules. Other evidence is obtained by solving particular BFP cases through numerical integration or through "particle schemes" such as the Direct Simulation Monte Carlo.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Weighted density Lattice Boltzmann approach to fluids under confinement

The Enskog-like kinetic approach, recently introduced by us to study strongly inhomogeneous flu- ids, is reconsidered in order to improve the description of the transport coefficients. The approach is based on a separation of the interaction between hydrodynamic and non-hydrodynamic parts. The latter is treated within a simple relaxation approximation. We show that, by considering the non-hydrodynamic part via a weighted density approximation, we obtain a better prediction of the transport coefficients. By virtue of the simplicity of the kinetic equation we are able to solve numer- ically the phase space distribution in the presence of inhomogeneities, such as confining surfaces, via a Lattice Boltzmann method. Analytical estimates of the importance of these corrections to the transport coefficients in bulk conditions is provided. Poiseuille flow of the hard-sphere fluid confined between two parallel smooth walls is studied and their pore-averaged properties are determined.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Charge transport in nanochannels: a molecular theory

We introduce a theoretical and numerical method to investigate the flow of charged fluid mixtures under extreme confinement. We model the electrolyte solution as a ternary mixture, comprising two ionic species of opposite charge and a third uncharged component. The microscopic approach is based on kinetic theory and is fully self-consistent. It allows to determine configurational prop- erties, such as layering near the confining walls, and the flow properties. We show that, under appropriate assumptions, the approach reproduces the phenomenological equations used to describe electrokinetic phenomena, without requiring the introduction of constitutive equations to determine the fluxes. Moreover, we model channels of arbitrary shape and nanometric roughness, features that have important repercussions on the transport properties of these systems. Numerical simulations are obtained by solving the evolution dynamics of the one-particle phase- space distributions of each species by means of a Lattice Boltzmann method for flows in straight and wedged channels. Results are presented for the microscopic density, the velocity profiles and for the volumetric and charge flow-rates. Strong departures from electroneutrality are shown to appear at molecular level.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Stabilized Lattice Boltzmann-Enskog method for compressible flows and its application to one and two-component fluids in nanochannels

A numerically stable method to solve the discretized Boltzmann-Enskog equation describing the behavior of non ideal fluids under inhomogeneous conditions is presented. The algorithm employed uses a Lagrangian finite-difference scheme for the treatment of the convective term and a forcing term to account for the molecular repulsion together with a Bhatnagar-Gross-Krook relaxation term. In order to eliminate the spurious currents induced by the numerical discretization procedure, we use a trapezoidal rule for the time integration together with a version of the two-distribution method of He et al. (J. Comp. Phys 152, 642 (1999)). Numerical tests show that, in the case of one component fluid in the presence of a spherical potential well, the proposed method reduces the numerical error by several orders of magnitude. We conduct another test by considering the flow of a two component fluid in a channel with a bottleneck and provide information about the density and velocity field in this structured geometry.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Thermally induced directed currents in hard rod systems

We study the non equilibrium statistical properties of a one dimensional hard-rod fluid undergoing collisions and subject to a spatially non uniform Gaussian heat-bath and periodic potential. The system is able to sustain finite currents when the spatially inhomogeneous heat-bath and the periodic potential profile display an appropriate relative phase shift, $ϕ$. By comparison with the collisionless limit, we determine the conditions for the most efficient transport among inelastic, elastic and non interacting rods. We show that the situation is complex as, depending on shape of the temperature profile, the current of one system may outperform the others.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamics of Fluid Mixtures in Nanospaces

A multicomponent extension of our recent theory of simple fluids [ U.M.B. Marconi and S. Melchionna, Journal of Chemical Physics, 131, 014105 (2009) ] is proposed to describe miscible and immiscible liquid mixtures under inhomogeneous, non steady conditions typical of confined fluid flows. We first derive from a microscopic level the evolution equations of the phase space distribution function of each component in terms of a set of self consistent fields, representing both body forces and viscous forces (forces dependent on the density distributions in the fluid and on the velocity distributions). Secondly, we solve numerically the resulting governing equations by means of the Lattice Boltzmann method whose implementation contains novel features with respect to existing approaches. Our model incorporates hydrodynamic flow, diffusion, surface tension, and the possibility for global and local viscosity variations. We validate our model by studying the bulk viscosity dependence of the mixture on concentration, packing fraction and size ratio. Finally we consider inhomogeneous systems and study the dynamics of mixtures in slits of molecular thickness and relate structural and flow properties.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Electro-osmotic flows under nanoconfinement: a self-consistent approach

We introduce a theoretical and numerical method to investigate the properties of electro-osmotic flows under conditions of extreme confinement. The present approach, aiming to provide a simple modeling of electrolyte solutions described as ternary mixtures, which comprises two ionic species and a third uncharged component, is an extension of our recent work on binary neutral mixtures. The approach, which combines elements of kinetic theory, density functional theory with Lattice-Boltzmann algorithms, is microscopic and self-consistent and does not require the us e of constitutive equations to determine the fluxes. Numerical solutions are obtained by solving the resulting coupled equations for the one-particle phase-space distributions of the species by means of a Lattice Boltzmann discretization procedure. Results are given for the microscopic density and velocity profiles and for the volumetric and charge flow.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Multicomponent Diffusion in Nanosystems

We present the detailed analysis of the diffusive transport of spatially inhomogeneous fluid mixtures and the interplay between structural and dynamical properties varying on the atomic scale. The present treatment is based on different areas of liquid state theory, namely kinetic and density functional theory and their implementation as an effective numerical method via the Lattice Boltzmann approach. By combining the first two methods it is possible to obtain a closed set of kinetic equations for the singlet phase space distribution functions of each species. The interactions among particles are considered within a self-consistent approximation and the resulting effective molecular fields are analyzed. We focus on multispecies diffusion in systems with short-range hard-core repulsion between particles of unequal sizes and weak attractive long-range interactions. As a result, the attractive part of the potential does not contribute explicitly to viscosity but to diffusivity and the thermodynamic properties. Finally, we obtain a practical scheme to solve the kinetic equations by employing a discretization procedure derived from the Lattice Boltzmann approach. Within this framework, we present numerical data concerning the mutual diffusion properties both in the case of a quiescent bulk fluid and shear flow inducing Taylor dispersion.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dynamic density functional theory versus Kinetic theory of simple fluids

By combining methods of kinetic and density functional theory, we present a description of molecular fluids which accounts for their microscopic structure and thermodynamic properties as well as for the hydrodynamic behavior. We focus on the evolution of the one particle phase space distribution, rather than on the evolution of the average particle density, which features in dynamic density functional theory. The resulting equation can be studied in two different physical limits: diffusive dynamics, typical of colloidal fluids without hydrodynamic interaction, where particles are subject to overdamped motion resulting from the coupling with a solvent at rest, and inertial dynamics, typical of molecular fluids. Finally, we propose an algorithm to solve numerically and efficiently the resulting kinetic equation by employing a discretization procedure analogous to the one used in the Lattice Boltzmann method.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Thermal convection in mono-disperse and bi-disperse granular gases: A simulation study

We present results of a simulation study of inelastic hard-disks vibrated in a vertical container. An Event-Driven Molecular Dynamics method is developed for studying the onset of convection. Varying the relevant parameters (inelasticity, number of layers at rest, intensity of the gravity) we are able to obtain a qualitative agreement of our results with recent hydrodynamical predictions. Increasing the inelasticity, a first continuous transition from the absence of convection to one convective roll is observed, followed by a discontinuous transition to two convective rolls, with hysteretic behavior. At fixed inelasticity and increasing gravity, a transition from no convection to one roll can be evidenced. If the gravity is further increased, the roll is eventually suppressed. Increasing the number of monolayers the system eventually localizes mostly at the bottom of the box: in this case multiple convective rolls as well as surface waves appear. We analyze the density and temperature fields and study the existence of symmetry breaking in these fields in the direction perpendicular to the injection of energy. We also study a binary mixture of grains with different properties (inelasticity or diameters). The effect of changing the properties of one of the components is analyzed, together with density, temperature and temperature ratio fields. Finally, the presence of a low-fraction of quasi-elastic impurities is shown to determine a sharp transition between convective and non-convective steady states.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Velocity fluctuations in cooling granular gases

We study the formation and the dynamics of correlations in the velocity field for 1D and 2D cooling granular gases with the assumption of negligible density fluctuations (``Homogeneous Velocity-correlated Cooling State'', HVCS). It is shown that the predictions of mean field models fail when velocity fluctuations become important. The study of correlations is done by means of molecular dynamics and introducing an Inelastic Lattice Maxwell Models. This lattice model is able to reproduce all the properties of the Homogeneous Cooling State and several features of the HVCS. Moreover it allows very precise measurements of structure functions and other crucial statistical indicators. The study suggests that both the 1D and the 2D dynamics of the velocity field are compatible with a diffusive dynamics at large scale with a more complex behavior at small scale. In 2D the issue of scale separation, which is of interest in the context of kinetic theories, is addressed.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Dynamic Density Functional Theory of Fluids

We present a new time-dependent Density Functional approach to study the relaxational dynamics of an assembly of interacting particles subject to thermal noise. Starting from the Langevin stochastic equations of motion for the velocities of the particles we are able by means of an approximated closure to derive a self-consistent deterministic equation for the temporal evolution of the average particle density. The closure is equivalent to assuming that the equal-time two-point correlation function out of equilibrium has the same properties as its equilibrium version. The changes in time of the density depend on the functional derivatives of the grand canonical free energy functional $F[ρ]$ of the system. In particular the static solutions of the equation for the density correspond to the exact equilibrium profiles provided one is able to determine the exact form of $F[ρ]$. In order to assess the validity of our approach we performed a comparison between the Langevin dynamics and the dynamic density functional method for a one-dimensional hard-rod system in three relevant cases and found remarkable agreement, with some interesting exceptions, which are discussed and explained. In addition, we consider the case where one is forced to use an approximate form of $F[ρ]$. Finally we compare the present method with the stochastic equation for the density proposed by other authors [Kawasaki,Kirkpatrick etc.] and discuss the role of the thermal fluctuations.

0 citations0 reviewsNo code linkedOpen access