Local structure of percolating gels at very low volume fractions

The formation of colloidal gels is strongly dependent on the volume fraction of the system and the strength of the interactions between the colloids. Here we explore very dilute solutions by the means of numerical simulations, and show that, in the absence of hydrodynamic interactions and for sufficiently strong interactions, percolating colloidal gels can be realised at very low values of the volume fraction. Characterising the structure of the network of the arrested material we find that, when reducing the volume fraction, the gels are dominated by low-energy local structures, analogous to the isolated clusters of the interaction potential. Changing the strength of the interaction allows us to tune the compactness of the gel as characterised by the fractal dimension, with low interaction strength favouring more chain-like structures

Steady State of microemulsions in shear flow

Steady-state properties of microemulsions in shear flow are studied in the context of a Ginzburg-Landau free-energy approach. Explicit expressions are given for the structure factor and the time correlation function at the one loop level of approximation. Our results predict a four-peak pattern for the structure factor, implying the simultaneous presence of interfaces aligned with two different orientations. Due to the peculiar interface structure a non-monotonous relaxation of the time correlator is also found.Comment: 5 pages, 3 figure

Are stress-free membranes really 'tensionless'?

In recent years it has been argued that the tension parameter driving the fluctuations of fluid membranes, differs from the imposed lateral stress, the 'frame tension'. In particular, stress-free membranes were predicted to have a residual fluctuation tension. In the present paper, this argument is reconsidered and shown to be inherently inconsistent -- in the sense that a linearized theory, the Monge model, is used to predict a nonlinear effect. Furthermore, numerical simulations of one-dimensional stiff membranes are presented which clearly demonstrate, first, that the internal 'intrinsic' stress in membranes indeed differs from the frame tension as conjectured, but second, that the fluctuations are nevertheless driven by the frame tension. With this assumption, the predictions of the Monge model agree excellently with the simulation data for stiffness and tension values spanning several orders of magnitude

Budding and vesiculation induced by conical membrane inclusions

Conical inclusions in a lipid bilayer generate an overall spontaneous curvature of the membrane that depends on concentration and geometry of the inclusions. Examples are integral and attached membrane proteins, viruses, and lipid domains. We propose an analytical model to study budding and vesiculation of the lipid bilayer membrane, which is based on the membrane bending energy and the translational entropy of the inclusions. If the inclusions are placed on a membrane with similar curvature radius, their repulsive membrane-mediated interaction is screened. Therefore, for high inclusion density the inclusions aggregate, induce bud formation, and finally vesiculation. Already with the bending energy alone our model allows the prediction of bud radii. However, in case the inclusions induce a single large vesicle to split into two smaller vesicles, bending energy alone predicts that the smaller vesicles have different sizes whereas the translational entropy favors the formation of equal-sized vesicles. Our results agree well with those of recent computer simulations.Comment: 11 pages, 12 figure

Micellar Aggregates of Gemini Surfactants: Monte Carlo Simulation of a Microscopic Model

We propose a "microscopic" model of gemini surfactants in aqueous solution. Carrying out extensive Monte Carlo simulations, we study the variation of the critical micellar concentration (CMC) of these model gemini surfactants with the variation of the (a) length of the spacer connecting the two hydrophilic heads, (b) length of the hydrophobic tail and (c) the bending rigidity of the hydrocarbon chains forming the spacer and the tail; some of the trends of variation are counter-intuitive but are in excellent agreement with the available experimental results. Our simulations also elucidate the dependence of the shapes of the micellar aggregates and the magnitude of the CMC on the geometrical shape and size of the surfactant molecules and the electrical charge on the hydrophilic heads

Traction force microscopy with optimized regularization and automated Bayesian parameter selection for comparing cells

Adherent cells exert traction forces on to their environment, which allows them to migrate, to maintain tissue integrity, and to form complex multicellular structures. This traction can be measured in a perturbation-free manner with traction force microscopy (TFM). In TFM, traction is usually calculated via the solution of a linear system, which is complicated by undersampled input data, acquisition noise, and large condition numbers for some methods. Therefore, standard TFM algorithms either employ data filtering or regularization. However, these approaches require a manual selection of filter- or regularization parameters and consequently exhibit a substantial degree of subjectiveness. This shortcoming is particularly serious when cells in different conditions are to be compared because optimal noise suppression needs to be adapted for every situation, which invariably results in systematic errors. Here, we systematically test the performance of new methods from computer vision and Bayesian inference for solving the inverse problem in TFM. We compare two classical schemes, L1- and L2-regularization, with three previously untested schemes, namely Elastic Net regularization, Proximal Gradient Lasso, and Proximal Gradient Elastic Net. Overall, we find that Elastic Net regularization, which combines L1 and L2 regularization, outperforms all other methods with regard to accuracy of traction reconstruction. Next, we develop two methods, Bayesian L2 regularization and Advanced Bayesian L2 regularization, for automatic, optimal L2 regularization. Using artificial data and experimental data, we show that these methods enable robust reconstruction of traction without requiring a difficult selection of regularization parameters specifically for each data set. Thus, Bayesian methods can mitigate the considerable uncertainty inherent in comparing cellular traction forces

Stress Tensors of Multiparticle Collision Dynamics Fluids

Stress tensors are derived for the multiparticle collision dynamics algorithm, a particle-based mesoscale simulation method for fluctuating fluids, resembling those of atomistic or molecular systems. Systems with periodic boundary conditions as well as fluids confined in a slit are considered. For every case, two equivalent expressions for the tensor are provided, the internal stress tensor, which involves all degrees of freedom of a system, and the external stress, which only includes the interactions with the confining surfaces. In addition, stress tensors for a system with embedded particles are determined. Based on the derived stress tensors, analytical expressions are calculated for the shear viscosity. Simulations illustrate the difference in fluctuations between the various derived expressions and yield very good agreement between the numerical results and the analytically derived expression for the viscosity

The steering gaits of sperm

Sperm are highly specialized cells, which have been subject to substantial evolutionary pressure. Whereas some sperm features are highly conserved, others have undergone major modifications. Some of these variations are driven by adaptation to mating behaviours or fitness at the organismic level. Others represent alternative solutions to the same task. Sperm must find the egg for fertilization. During this task, sperm rely on long slender appendages termed flagella that serve as sensory antennas, propellers and steering rudders. The beat of the flagellum is periodic. The resulting travelling wave generates the necessary thrust for propulsion in the fluid. Recent studies reveal that, for steering, different species rely on different fundamental features of the beat wave. Here, we discuss some examples of unity and diversity across sperm from different species with a particular emphasis on the steering mechanisms. This article is part of the Theo Murphy meeting issue ‘Unity and diversity of cilia in locomotion and transport’

Implicit and explicit solvent models for the simulation of a single polymer chain in solution: Lattice Boltzmann vs Brownian dynamics

We present a comparative study of two computer simulation methods to obtain static and dynamic properties of dilute polymer solutions. The first approach is a recently established hybrid algorithm based upon dissipative coupling between Molecular Dynamics and lattice Boltzmann (LB), while the second is standard Brownian Dynamics (BD) with fluctuating hydrodynamic interactions. Applying these methods to the same physical system (a single polymer chain in a good solvent in thermal equilibrium) allows us to draw a detailed and quantitative comparison in terms of both accuracy and efficiency. It is found that the static conformations of the LB model are distorted when the box length L is too small compared to the chain size. Furthermore, some dynamic properties of the LB model are subject to an $L^{-1}$ finite size effect, while the BD model directly reproduces the asymptotic $L \to \infty$ behavior. Apart from these finite size effects, it is also found that in order to obtain the correct dynamic properties for the LB simulations, it is crucial to properly thermalize all the kinetic modes. Only in this case, the results are in excellent agreement with each other, as expected. Moreover, Brownian Dynamics is found to be much more efficient than lattice Boltzmann as long as the degree of polymerization is not excessively large.Comment: 11 figures, submitted to J. Chem. Phy

Phase Behaviour of Amphiphilic Monolayers: Theory and Simulation

Coarse grained models of monolayers of amphiphiles (Langmuir monolayers) have been studied theoretically and by computer simulations. We discuss some of the insights obtained with this approach, and present new simulation results which show that idealised models can successfully reproduce essential aspects of the generic phase behaviour of Langmuir monolayers.Comment: To appear in J. Phys.: Cond. Matte