Source author record

Rodrigo Soto

Rodrigo Soto 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

28works
9topics
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

28 published item(s)

preprint2021arXiv

Navier--Stokes transport coefficients for a model of a confined quasi-two-dimensional granular binary mixture

The Navier--Stokes transport coefficients for a model of a confined quasi-two-dimensional granular binary mixture of inelastic hard spheres are determined from the Boltzmann kinetic equation. A normal or hydrodynamic solution to the Boltzmann equation is obtained via the Chapman--Enskog method for states near the local version of the homogeneous time-dependent state. The mass, momentum, and heat fluxes are determined to first order in the spatial gradients of the hydrodynamic fields, and the associated transport coefficients are identified. As expected, they are given in terms of the solutions of a set of coupled linear integral equations. In addition, in contrast to previous results obtained for low-density granular mixtures, there are also nonzero contributions to the first-order approximations to the partial temperatures $T_i^{(1)}$ and the cooling rate $ζ^{(1)}$. Explicit forms for the diffusion transport coefficients, the shear viscosity coefficient, and the quantities $T_i^{(1)}$ and $ζ^{(1)}$ are obtained by assuming the steady-state conditions and by considering the leading terms in a Sonine polynomial expansion. The above transport coefficients are given in terms of the coefficients of restitution, concentration, and the masses and diameters of the components of the mixture. The results apply in principle for arbitrary degree of inelasticity and are not limited to specific values of concentration, mass and/or size ratios. As a simple application of these results, the violation of the Onsager reciprocal relations for a confined granular mixture is quantified in terms of the parameter space of the problem.

preprint2020arXiv

Run-and-tumble bacteria slowly approaching the diffusive regime

The run-and-tumble (RT) dynamics followed by bacterial swimmers gives rise first to a ballistic motion due to their persistence, and later, through consecutive tumbles, to a diffusive process. Here we investigate how long it takes for a dilute swimmer suspension to reach the diffusive regime as well as what is the amplitude of the deviations from the diffusive dynamics. A linear time dependence of the mean-squared displacement (MSD) is insufficient to characterize diffusion and thus we also focus on the excess kurtosis of the displacement distribution. Four swimming strategies are considered: (i) the conventional RT model with complete reorientation after tumbling, (ii) the case of partial reorientation, characterized by a distribution of tumbling angles, (iii) a run-and-reverse model with rotational diffusion, and (iv) a RT particle where the tumbling rate depends on the stochastic concentration of an internal protein. By analyzing the associated kinetic equations for the probability density function and simulating the models, we find that for models (ii), (iii), and (iv) the relaxation to diffusion can take much longer than the mean time between tumble events, evidencing the existence of large tails in the particle displacements. Moreover, the excess kurtosis can assume large positive values. In model (ii) it is possible for some distributions of tumbling angles that the MSD reaches a linear time dependence but, still, the dynamics remains non-Gaussian for long times. This is also the case in model (iii) for small rotational diffusivity. For all models, the long-time diffusion coefficients are also obtained. The theoretical approach, which relies on eigenvalue and angular Fourier expansions of the van Hove function, is in excellent agreement with the simulations.

preprint2019arXiv

3D spatial exploration by E. coli echoes motor temporal variability

Unraveling bacterial strategies for spatial exploration is crucial for understanding the complexity in the organization of life. Bacterial motility determines the spatio-temporal structure of microbial communities, controls infection spreading and the microbiota organization in guts or in soils. Most theoretical approaches for modeling bacterial transport rely on their run-and-tumble motion. For Escherichia coli, the run time distribution was reported to follow a Poisson process with a single characteristic time related to the rotational switching of the flagellar motors. However, direct measurements on flagellar motors show heavy-tailed distributions of rotation times stemming from the intrinsic noise in the chemotactic mechanism. Currently, there is no direct experimental evidence that the stochasticity in the chemotactic machinery affect the macroscopic motility of bacteria. In stark contrast with the accepted vision of run-and-tumble, here we report a large behavioral variability of wild-type \emph{E. coli}, revealed in their three-dimensional trajectories. At short observation times, a large distribution of run times is measured on a population and attributed to the slow fluctuations of a signaling protein triggering the flagellar motor reversal. Over long times, individual bacteria undergo significant changes in motility. We demonstrate that such a large distribution of run times introduces measurement biases in most practical situations. Our results reconcile the notorious conundrum between run time observations and motor switching statistics. We finally propose that statistical modeling of transport properties currently undertaken in the emerging framework of active matter studies, should be reconsidered under the scope of this large variability of motility features.

preprint2019arXiv

E. coli "super-contaminates" narrow ducts fostered by broad run-time distribution

One striking feature of bacterial motion is their ability to swim upstream along corners and crevices, by leveraging hydrodynamic interactions. This motion through anatomic ducts or medical devices might be at the origin of serious infections. However, it remains unclear how bacteria can maintain persistent upstream motion while exhibiting run-and-tumble dynamics. Here we demonstrate that E. coli can travel upstream in microfluidic devices over distances of 15 millimeters in times as short as 15 minutes. Using a stochastic model relating the run times to the time bacteria spend on surfaces, we quantitatively reproduce the evolution of the contamination profiles when considering a broad distribution of run times. Interestingly, the experimental data cannot be reproduced using the usually accepted exponential distribution of run times. Our study demonstrates that the run-and-tumble statistics determine macroscopic bacterial transport properties. This effect, that we name "super-contamination", could explain the fast onset of some life-threatening medical emergencies.

preprint2016arXiv

Coarsening and clustering in run-and-tumble dynamics with short-range exclusion

The emergence of clustering and coarsening in crowded ensembles of self-propelled agents is studied using a lattice model in one-dimension. The persistent exclusion process, where particles move at directions that change randomly at a low tumble rate $α$, is extended allowing sites to be occupied by more than one particle, with a maximum $n_\text{max}$ per site. Three phases are distinguished. For $n_\text{max}=1$ a gas of clusters form, with sizes distributed exponentially and no coarsening takes place. For $n_\text{max}\geq 3$ and small values of $α$, coarsening takes place and few large clusters appear, with a large fraction of the total number of particles in them. In the same range of $n_\text{max}$ but for larger values of $α$, a gas phase where a negligible fraction of particles takes part of clusters. Finally, $n_\text{max}=2$ corresponds to a crossover phase. The character of the transitions between phases is studied extending the model to allow $n_\text{max}$ to take real values and jumps to an occupied site are probabilistic. The transition from the gas of clusters to the coarsening phase is continuous and the mass of the large clusters grows continuously when varying the maximum occupancy, and the crossover found corresponds to values close to the transition. The second transition, from the coarsening to the gaseous phase can be either continuous or discontinuous depending on the parameters, with a critical point separating both cases.

preprint2016arXiv

Dynamics of noncohesive confined granular media

Despite the ubiquitousness and technological and scientific importance of granular matter, our understanding is still very poor compared to molecular fluids and solids. Until today, there is no unified description, which indeed seems unreachable. However, it has been proposed that important advances could be attained for noncohesive, hard-sphere like systems, by combining fluid dynamics with phase-field modeling through an appropriate order parameter (Aranson and Tsimring, 2006). Here, we present a review of the dynamics of confined granular matter, for which this systematic approach has proven its value. Motivated by the pioneering work of Olafsen and Urbach (1998), many experimental, theoretical and numerical studies of model confined granular systems have been realized, which have unveiled a very large variety of fundamental phenomena. In this review, we focus on few of these fundamental aspects, namely phase coexistence, effective surface tension, and a detailed description of the liquid state.

preprint2016arXiv

Effect of the vibration profile on shallow granular systems

We describe the collective behavior of a system of many inelastic spherical particles inside a box which is being periodically vibrated. The box is shallow, with large horizontal dimensions, while the height is less than two particle diameters. The vibrations are not symmetric: the time the box is moving up is, in general, different to the time it is moving down. The limit cycles of isolated grains are largely affected by the asymmetry of the vibration mode, increasing the size in phase space of the chaotic regions. When many grains are placed in the box, the phase separation between dense, solid-like regions, coexisting with fluid-like regions takes place at smaller global densities for asymmetric vibration profiles. Besides, the order parameter of the transition takes larger values when asymmetric forcing is used.

preprint2015arXiv

Self-assembly of Active Colloidal Molecules with Dynamic Function

Catalytically active colloids maintain non-equilibrium conditions in which they produce and deplete chemicals and hence effectively act as sources and sinks of molecules. While individual colloids that are symmetrically coated do not exhibit any form of dynamical activity, the concentration fields resulting from their chemical activity decay as $1/r$ and produce gradients that attract or repel other colloids depending on their surface chemistry and ambient variables. This results in a non-equilibrium analogue of ionic systems, but with the remarkable novel feature of action-reaction symmetry breaking. We study solutions of such chemically active colloids in dilute conditions when they join up to form molecules via generalized ionic bonds, and discuss how we can achieve structures with time dependent functionality. In particular, we study a molecule that adopts a spontaneous oscillatory pattern of conformations, and another that exhibits a run-and-tumble dynamics similar to bacteria. Our study shows that catalytically active colloids could be used for designing self-assembled structures that posses dynamical functionalities that are determined by their prescribed 3D structures, a strategy that follows the design principle of proteins.

preprint2015arXiv

Universality and criticality of a second-order granular solid-liquid-like phase transition

We experimentally study the critical properties of the non-equilibrium solid-liquid-like transition that takes place in vibrated granular matter. The critical dynamics is characterized by the coupling of the density field with the bond-orientational order parameter $Q_4$, which measures the degree of local crystallization. Two setups are compared, which present the transition at different critical accelerations as a a result of modifying the energy dissipation parameters. In both setups five independent critical exponents are measured, associated to different properties of $Q_4$: the correlation length, relaxation time, vanishing wavenumber limit (static susceptibility), the hydrodynamic regime of the pair correlation function, and the amplitude of the order parameter. The respective critical exponents agree in both setups and are given by $ν_{\perp} = 1$, $ν_{\parallel} = 2$, $γ= 1$, $η\approx 0.6 - 0.67$, and $β=1/2$, whereas the dynamical critical exponent is $z = ν_{\parallel}/ν_{\perp} = 2$. The agreement on five exponents is an exigent test for the universality of the transition. Thus, while dissipation is strictly necessary to form the crystal, the path the system undergoes towards the phase separation is part of a well defined universality class. In fact, the local order shows critical properties while density does not. Being the later conserved, the appropriate model that couples both is model C in the Hohenberg and Halperin classification. The measured exponents are in accord with the non-equilibrium extension to model C if we assume that $α$, the exponent associated in equilibrium to the specific heat divergence but with no counterpart in this non-equilibrium experiment, vanishes.

preprint2014arXiv

Dynamics of a first order transition to an absorbing state

A granular system confined in a quasi two-dimensional box that is vertically vibrated can transit to an absorbing state in which all particles bounce vertically in phase with the box, with no horizontal motion. In principle, this state can be reached for any density lower than the one corresponding to one complete monolayer, which is then the critical density. Below this critical value, the transition to the absorbing state is of first order, with long metastable periods, followed by rapid transitions driven by homogeneous nucleation. Molecular dynamics simulations and experiments show that there is a dramatic increase on the metastable times far below the critical density; in practice, it is impossible to observe spontaneous transitions close to the critical density. This peculiar feature is a consequence of the non-equilibrium nature of this first order transition to the absorbing state. A Ginzburg-Landau model, with multiplicative noise, describes qualitatively the observed phenomena and explains the macroscopic size of the critical nuclei. The nuclei become of small size only close to a second critical point where the active phase becomes unstable via a saddle node bifurcation. It is only close to this second critical point that experiments and simulations can evidence spontaneous transitions to the absorbing state while the metastable times grow dramatically moving away from it.

preprint2014arXiv

Self-assembly of catalytically active colloidal molecules: Tailoring activity through surface chemistry

A heterogeneous and dilute suspension of catalytically active colloids is studied as a non-equilibrium analogue of ionic systems, which has the remarkable feature of action-reaction symmetry breaking. Symmetrically coated colloids are found to join up to form self-assembled molecules that could be inert or have spontaneous activity in the form of net translational velocity and spin depending on their symmetry properties and their constituents. The type of activity can be adjusted by changing the surface chemistry and ambient variables that control the surface reactions and the phoretic drift.

preprint2014arXiv

Shear viscosity of a model for confined granular media

The shear viscosity in the dilute regime of a model for confined granular matter is studied by simulations and kinetic theory. The model consists on projecting into two dimensions the motion of vibrofluidized granular matter in shallow boxes by modifying the collision rule: besides the restitution coefficient that accounts for the energy dissipation, there is a separation velocity that is added in each collision in the normal direction. The two mechanisms balance on average, producing stationary homogeneous states. Molecular dynamics simulations show that in the steady state the distribution function departs from a Maxwellian, with cumulants that remain small in the whole range of inelasticities. The shear viscosity normalized with stationary temperature presents a clear dependence with the inelasticity, taking smaller values compared to the elastic case. A Boltzmann-like equation is built and analyzed using linear response theory. It is found that the predictions show an excellent agreement with the simulations when the correct stationary distribution is used but a Maxwellian approximation fails in predicting the inelasticity dependence of the viscosity. These results confirm that transport coefficients depend strongly on the mechanisms that drive them to stationary states.

preprint2013arXiv

Capillary-like Fluctuations of a Solid-Liquid Interface in a Non-Cohesive Granular System

One of the most noticeable collective motion of non-cohesive granular matter is clustering under certain conditions. In particular, when a quasi-two-dimensional monolayer of mono-disperse non-cohesive particles is vertically vibrated, a solid-liquid-like transition occurs when the driving amplitude exceeds a critical value. Here, the physical mechanism underlying particle clustering relies on the strong interactions mediated by grain collisions, rather than on grain-grain cohesive forces. In average, the solid cluster resembles a drop, with a striking circular shape. We experimentally investigate the coarse-grained solid-liquid interface fluctuations, which are characterized through the static and dynamic correlation functions in the Fourier space. These fluctuations turn out to be well described by the capillary wave theory, which allows us to measure the solid-liquid interface surface tension and mobility once the granular "thermal" kinetic energy is determined. Despite the system is strongly out of equilibrium and that the granular temperature is not uniform, there is energy equipartition at the solid-liquid interface, for a relatively large range of angular wave-numbers. Furthermore, both surface tension and mobility are consistent with a simple order of magnitude estimation considering the characteristic energy, length and time scales, which is very similar to what can be done for atomic systems. Please find the supplementary material of this article here: arXiv:1304.2646

preprint2013arXiv

Hydrodynamic modes in a confined granular fluid

Confined granular fluids, placed in a shallow box that is vibrated vertically, can achieve homogeneous stationary states thanks to energy injection mechanisms that take place throughout the system. These states can be stable even at high densities and inelasticities allowing for a detailed analysis of the hydrodynamic modes that govern the dynamics of granular fluids. Analyzing the decay of the time correlation functions it is shown that there is a crossover between a quasielastic regime in which energy evolves as a slow mode, to a inelastic regime, with energy slaved to the other conserved fields. The two regimes have well differentiated transport properties and, in the inelastic regime, the dynamics can be described by a reduced hydrodynamics with modified longitudinal viscosity and sound speed. The crossover between the two regimes takes place at a wavevector that is proportional to the inelasticity. A two dimensional granular model, with collisions that mimic the energy transfers that take place in a confined system is studied by means of microscopic simulations. The results show excellent agreement with the theoretical framework and allows the validation of hydrodynamic-like models.

preprint2013arXiv

Run-and-tumble in a crowded environment: persistent exclusion process for swimmers

The effect of crowding on the run-and-tumble dynamics of swimmers such as bacteria is studied using a discrete lattice model of mutually excluding particles that move with constant velocity along a direction that is randomized at a rate $α$. In stationary state, the system is found to break into dense clusters in which particles are trapped or stopped from moving. The characteristic size of these clusters predominantly scales as $α^{-0.5}$ both in 1D and 2D. For a range of densities, due to cooperative effects, the stopping time scales as ${\cal T}_{1d}^{0.85}$ and as ${\cal T}_{2d}^{0.8}$, where ${\cal T}_d$ is the diffusive time associated with the motion of cluster boundaries. Our findings might be helpful in understanding the early stages of biofilm formation.

preprint2012arXiv

Fluctuations and Criticality of a Granular Solid-Liquid-like Phase Transition

We present an experimental study of density and order fluctuations in the vicinity of the solid-liquid-like transition that occurs in a vibrated quasi-two-dimensional granular system. The two-dimensional projected static and dynamic correlation functions are studied. We show that density fluctuations, characterized through the structure factor, increase in size and intensity as the transition is approached, but they do not change significantly at the transition itself. The dense, metastable clusters, which present square symmetry, also increase their local order in the vicinity of the transition. This is characterized through the bond-orientational order parameter $Q_4$, which in Fourier space obeys an Ornstein-Zernike behavior. Depending on filling density and vertical height, the transition can be of first or second order type. In the latter case, the associated correlation length $ξ_4$, relaxation time $τ_4$, zero $k$ limit of $Q_4$ fluctuations (static susceptibility), the pair correlation function of $Q_4$, and the amplitude of the order parameter obey critical power laws, with saturations due to finite size effects. Their respective critical exponents are $ν_{\bot} = 1$, $ν_{||} = 2$, $γ= 1$, $η=0.67$, and $β=1/2$, whereas the dynamical critical exponent $z = ν_{||}/ν_{\bot} = 2$. These results are consistent with model C of dynamical critical phenomena, valid for a non-conserved critical order parameter (bond-orientation order) coupled to a conserved field (density).

preprint2012arXiv

Stochastic Quantization and Casimir Forces: Pistons of Arbitrary Cross Section

Recently, a method based on stochastic quantization has been proposed to compute the Casimir force and its fluctuations in arbitrary geometries. It relies on the spectral decomposition of the Laplacian operator in the given geometry. Both quantum and thermal fluctuations are considered. Here we use such method to compute the Casimir force on the plates of a finite piston of arbitrary cross section. Asymptotic expressions valid at low and high temperatures and short and long distances are obtained. The case of a piston with triangular cross section is analysed in detail. The regularization of the divergent stress tensor is described.

preprint2012arXiv

Stochastic resonance on the transverse displacement of swimmers in an oscillatory shear flow

Self-propelled microorganisms, such as unicellular algae or bacteria, swim along their director relative to the fluid velocity. Under a steady shear flow the director rotates in close orbit, a periodic structure that is preserved under an oscillatory shear flow. If the shear flow is subjected to small fluctuations produced by small irregularities in the microchannel or by other swimmers nearby, the director dynamics becomes stochastic. Numerical integration of the swimmer motion shows that there is stochastic resonance: The displacement in the vorticity direction is maximized for a finite noise intensity. This transverse displacement resonance is observed when the displacement is coarse grained over several periods, although the director is preferentially oriented along the flow. The resonant noise intensity is proportional to the oscillation frequency and independent of the shear rate. The enhanced displacement can have effects on the transverse diffusion of swimmers and the rheology of the suspension.

preprint2011arXiv

A two-sphere model for bacteria swimming near solid surfaces

We present a simple model for bacteria like \emph{Escherichia coli} swimming near solid surfaces. It consists of two spheres of different radii connected by a dragless rod. The effect of the flagella is taken into account by imposing a force on the tail sphere and opposite torques exerted by the rod over the spheres. The hydrodynamic forces and torques on the spheres are computed by considering separately the interaction of a single sphere with the surface and with the flow produced by the other sphere. Numerically, we solve the linear system which contains the geometrical constraints and the force-free and torque-free conditions. The dynamics of this swimmer near a solid boundary is very rich, showing three different behaviors depending on the initial conditions: (1) swimming in circles in contact with the wall, (2) swimming in circles at a finite distance from the wall, and (3) swimming away from it. Furthermore, the order of magnitude of the radius of curvature for the circular motion is in the range $8-50\,μ$m, close to values observed experimentally.

preprint2011arXiv

Characterization of the Melting Transition in Two Dimensions at Vanishing External Pressure Using Molecular Dynamics Simulations

A molecular dynamics study of a two dimensional system of particles interacting through a Lennard-Jones pairwise potential is performed at fixed temperature and vanishing external pressure. As the temperature is increased, a solid-to-liquid transition occurs. When the melting temperature $T_c$ is approached from below, there is a proliferation of dislocation pairs and the elastic constant approaches the value predicted by the KTHNY theory. In addition, as $T_c$ is approached from above, the relaxation time increases, consistent with an approach to criticality. However, simulations fail to produce a stable hexatic phase using systems with up to 90,000 particles. A significant jump in enthalpy at $T_c$ is observed, consistent with either a first order or a continuous transition. The role of external pressure is discussed.

preprint2011arXiv

Stochastic Quantization and Casimir Forces

In this paper we show how the stochastic quantization method developed by Parisi and Wu can be used to obtain Casimir forces. Both quantum and thermal fluctuations are taken into account by a Langevin equation for the field. The method allows the Casimir force to be obtained directly, derived from the stress tensor instead of the free energy. It only requires the spectral decomposition of the Laplacian operator in the given geometry. The formalism provides also an expression for the fluctuations of the force. As an application we compute the Casimir force on the plates of a finite piston of arbitrary cross section. Fluctuations of the force are also directly obtained, and it is shown that, in the piston case, the variance of the force is twice the force squared.

preprint2011arXiv

Sudden chain energy transfer events in vibrated granular media

In a mixture of two species of grains of equal size but different mass, placed in a vertically vibrated shallow box, there is spontaneous segregation. Once the system is at least partly segregated and clusters of the heavy particles have formed, there are sudden peaks of the horizontal kinetic energy of the heavy particles, that is otherwise small. Together with the energy peaks the clusters rapidly expand and the segregation is partially lost. The process repeats once segregation has taken place again. Depending on the experimental or numerical parameters, the energy bursts can occur either randomly or with some regularity in time. An explanation for these events is provided based on the existence of a fixed point for an isolated particle bouncing with only vertical motion. The horizontal energy peaks occur when the energy stored in the vertical motion is partly transferred into horizontal energy through a chain reaction of collisions between heavy particles. A necessary condition for the observed regularity of the events is that chain reactions involve most of the heavy particles.

preprint2010arXiv

Competition of Brazil nut effect, buoyancy, and inelasticity induced segregation in a granular mixture

It has been recently reported that a granular mixture in which grains differ in their restitution coefficients presents segregation: the more inelastic particles sink to the bottom. When other segregation mechanisms as buoyancy and the Brazil nut effect are present, the inelasticity induced segregation can compete with them. First, a detailed analysis, based on numerical simulations of two dimensional systems, of the competition between buoyancy and the inelasticity induced segregation is presented, finding that there is a transition line in the parameter space that determines which mechanism is dominant. In the case of neutrally buoyant particles having different sizes the inelasticity induced segregation can compete with the Brazil nut effect (BNE). Reverse Brazil nut effect (RBNE) could be obtained at large inelasticities of the intruder. At intermediate values, BNE and RBNE coexist and large inelastic particles are found both near the bottom and at the top of the system.

preprint2010arXiv

Dynamical approach to the Casimir effect

Casimir forces can appear between intrusions placed in different media driven by several fluctuation mechanisms, either in equilibrium or out of it. Herein, we develop a general formalism to obtain such forces from the dynamical equations of the fluctuating medium, the statistical properties of the driving noise, and the boundary conditions of the intrusions (which simulate the interaction between the intrusions and the medium). As a result, an explicit formula for the Casimir force over the intrusions is derived. This formalism contains the thermal Casimir effect as a particular limit and generalizes the study of the Casimir effect to such systems through their dynamical equations, with no appeal to their Hamiltonian, if any exists. In particular, we study the Casimir force between two infinite parallel plates with Dirichlet or Neumann boundary conditions, immersed in several media with finite correlation lengths (reaction--diffusion system, liquid crystals, and two coupled fields with non-Hermitian evolution equations). The driving Gaussian noises have vanishing or finite spatial or temporal correlation lengths; in the first case, equilibrium is reobtained and finite correlations produce nonequilibrium dynamics. The results obtained show that, generally, nonequilibrium dynamics leads to Casimir forces, whereas Casimir forces are obtained in equilibrium dynamics if the stress tensor is anisotropic.

preprint2010arXiv

Enhanced diffusion due to active swimmers at a solid surface

We consider two systems of active swimmers moving close to a solid surface, one being a living population of wild-type \textit{E. coli} and the other being an assembly of self-propelled Au-Pt rods. In both situations, we have identified two different types of motion at the surface and evaluated the fraction of the population that displayed ballistic trajectories (active swimmers) with respect to those showing random-like behavior. We studied the effect of this complex swimming activity on the diffusivity of passive tracers also present at the surface. We found that the tracer diffusivity is enhanced with respect to standard Brownian motion and increases linearly with the activity of the fluid, defined as the product of the fraction of active swimmers and their mean velocity. This result can be understood in terms of series of elementary encounters between the active swimmers and the tracers.

preprint2010arXiv

Extended event driven molecular dynamics for simulating dense granular matter

A new numerical method is presented to efficiently simulate the inelastic hard sphere (IHS) model for granular media, when fluid and frozen regions coexist in the presence of gravity. The IHS model is extended by allowing particles to change their dynamics into either a frozen state or back to the normal collisional state, while computing the dynamics only for the particles in the normal state. Careful criteria, local in time and space, are designed such that particles become frozen only at mechanically stable positions. The homogeneous deposition over a static surface and the dynamics of a rotating drum are studied as test cases. The simulations agree with previous experimental results. The model is much more efficient than the usual event driven method and allows to overcome some of the difficulties of the standard IHS model, such as the existence of a static limit.

preprint2002arXiv

Temperature inversion in granular fluids under gravity

We study, via hydrodynamic equations, the granular temperature profile of a granular fluid under gravity and subjected to energy injection from a base. It is found that there exists a turn-up in the granular temperature and that, far from the base, it increases linearly with height. We show that this phenomenon, observed previously in experiments and computer simulations, is a direct consequence of the heat flux law, different form Fourier's, in granular fluids. The positive granular temperature gradient is proportional to gravity and a transport coefficient $μ_0$, relating the heat flux to the density gradients, that is characteristic of granular systems. Our results provide a method to compute the value $μ_0$ for different restitution coefficients. The theoretical predictions are verified by means of molecular dynamics simulations, and the value of $μ_0$ is computed for the two dimensional inelastic hard sphere model. We provide, also, a boundary condition for the temperature field that is consistent with the modified Fourier's law.