Fully coupled simulations of non-colloidal monodisperse sheared suspensions

In this work we investigate numerically the dynamics of sheared suspensions in the limit of vanishingly small fluid and particle inertia. The numerical model we used is able to handle the multi-body hydrodynamic interactions between thousands of particles embedded in a linear shear flow. The presence of the particles is modeled by momentum source terms spread out on a spherical envelop forcing the Stokes equations of the creeping flow. Therefore all the velocity perturbations induced by the moving particles are simultaneously accounted for. The statistical properties of the sheared suspensions are related to the velocity fluctuation of the particles. We formed averages for the resulting velocity fluctuation and rotation rate tensors. We found that the latter are highly anisotropic and that all the velocity fluctuation terms grow linearly with particle volume fraction. Only one off-diagonal term is found to be non zero (clearly related to trajectory symmetry breaking induced by the non-hydrodynamic repulsion force). We also found a strong correlation of positive/negative velocities in the shear plane, on a time scale controlled by the shear rate (direct interaction of two particles). The time scale required to restore uncorrelated velocity fluctuations decreases continuously as the concentration increases. We calculated the shear induced self-diffusion coefficients using two different methods and the resulting diffusion tensor appears to be anisotropic too. The microstructure of the suspension is found to be drastically modified by particle interactions. First the probability density function of velocity fluctuations showed a transition from exponential to Gaussian behavior as particle concentration varies. Second the probability of finding close pairs while the particles move under shear flow is strongly enhanced by hydrodynamic interactions when the concentration increases

Ion pump activity generates fluctuating electrostatic forces in biomembranes

We study the non-equilibrium dynamics of lipid membranes with proteins that actively pump ions across the membrane. We find that the activity leads to a fluctuating force distribution due to electrostatic interactions arising from variation in dielectric constant across the membrane. By applying a multipole expansion we find effects on both the tension and bending rigidity dominated parts of the membranes fluctuation spectrum. We discuss how our model compares with previous studies of force-multipole models.Comment: 6 pages, 2 figures, to appear in EP

Descriptions of membrane mechanics from microscopic and effective two-dimensional perspectives

Mechanics of fluid membranes may be described in terms of the concepts of mechanical deformations and stresses, or in terms of mechanical free-energy functions. In this paper, each of the two descriptions is developed by viewing a membrane from two perspectives: a microscopic perspective, in which the membrane appears as a thin layer of finite thickness and with highly inhomogeneous material and force distributions in its transverse direction, and an effective, two-dimensional perspective, in which the membrane is treated as an infinitely thin surface, with effective material and mechanical properties. A connection between these two perspectives is then established. Moreover, the functional dependence of the variation in the mechanical free energy of the membrane on its mechanical deformations is first studied in the microscopic perspective. The result is then used to examine to what extent different, effective mechanical stresses and forces can be derived from a given, effective functional of the mechanical free energy.Comment: 37 pages, 3 figures, minor change

Subdiffusion and weak ergodicity breaking in the presence of a reactive boundary

We derive the boundary condition for a subdiffusive particle interacting with a reactive boundary with finite reaction rate. Molecular crowding conditions, that are found to cause subdiffusion of larger molecules in biological cells, are shown to effect long-tailed distributions with identical exponent for both the unbinding times from the boundary to the bulk and the rebinding times from the bulk. This causes a weak ergodicity breaking: typically, an individual particle either stays bound or remains in the bulk for very long times. We discuss why this may be beneficial for in vivo gene regulation by DNA-binding proteins, whose typical concentrations are nanomolarComment: 4 pages, 1 figure, REVTeX4, accepted to Phys Rev Lett, some typos correcte

Open String Fluctuations in AdS with and without Torsion

The equations of motion and boundary conditions for the fluctuations around a classical open string, in a curved space-time with torsion, are considered in compact and world-sheet covariant form. The rigidly rotating open strings in Anti de Sitter space with and without torsion are investigated in detail. By carefully analyzing the tangential fluctuations at the boundary, we show explicitly that the physical fluctuations (which at the boundary are combinations of normal and tangential fluctuations) are finite, even though the world-sheet is singular there. The divergent 2-curvature thus seems less dangerous than expected, in these cases. The general formalism can be straightforwardly used also to study the (bosonic part of the) fluctuations around the closed strings, recently considered in connection with the AdS/CFT duality, on AdS_5 \times S^5 and AdS_3 \times S^3 \times T^4.Comment: 19 pages, Late

Mean first-passage time of surface-mediated diffusion in spherical domains

We present an exact calculation of the mean first-passage time to a target on the surface of a 2D or 3D spherical domain, for a molecule alternating phases of surface diffusion on the domain boundary and phases of bulk diffusion. The presented approach is based on an integral equation which can be solved analytically. Numerically validated approximation schemes, which provide more tractable expressions of the mean first-passage time are also proposed. In the framework of this minimal model of surface-mediated reactions, we show analytically that the mean reaction time can be minimized as a function of the desorption rate from the surface.Comment: to appear in J. Stat. Phy

How Landscape Heterogeneity Frames Optimal Diffusivity in Searching Processes

Theoretical and empirical investigations of search strategies typically have failed to distinguish the distinct roles played by density versus patchiness of resources. It is well known that motility and diffusivity of organisms often increase in environments with low density of resources, but thus far there has been little progress in understanding the specific role of landscape heterogeneity and disorder on random, non-oriented motility. Here we address the general question of how the landscape heterogeneity affects the efficiency of encounter interactions under global constant density of scarce resources. We unveil the key mechanism coupling the landscape structure with optimal search diffusivity. In particular, our main result leads to an empirically testable prediction: enhanced diffusivity (including superdiffusive searches), with shift in the diffusion exponent, favors the success of target encounters in heterogeneous landscapes

Weakly non-ergodic Statistical Physics

We find a general formula for the distribution of time averaged observables for weakly non-ergodic systems. Such type of ergodicity breaking is known to describe certain systems which exhibit anomalous fluctuations, e.g. blinking quantum dots and the sub-diffusive continuous time random walk model. When the fluctuations become normal we recover usual ergodic statistical mechanics. Examples of a particle undergoing fractional dynamics in a binding force field are worked out in detail. We briefly discuss possible physical applications in single particle experiments

Lattice Boltzmann simulations of soft matter systems

This article concerns numerical simulations of the dynamics of particles immersed in a continuum solvent. As prototypical systems, we consider colloidal dispersions of spherical particles and solutions of uncharged polymers. After a brief explanation of the concept of hydrodynamic interactions, we give a general overview over the various simulation methods that have been developed to cope with the resulting computational problems. We then focus on the approach we have developed, which couples a system of particles to a lattice Boltzmann model representing the solvent degrees of freedom. The standard D3Q19 lattice Boltzmann model is derived and explained in depth, followed by a detailed discussion of complementary methods for the coupling of solvent and solute. Colloidal dispersions are best described in terms of extended particles with appropriate boundary conditions at the surfaces, while particles with internal degrees of freedom are easier to simulate as an arrangement of mass points with frictional coupling to the solvent. In both cases, particular care has been taken to simulate thermal fluctuations in a consistent way. The usefulness of this methodology is illustrated by studies from our own research, where the dynamics of colloidal and polymeric systems has been investigated in both equilibrium and nonequilibrium situations.Comment: Review article, submitted to Advances in Polymer Science. 16 figures, 76 page