Reconstruction of the three-dimensional beat pattern underlying swimming behaviors of sperm.

 The eukaryotic flagellum propels sperm cells and simultaneously detects physical and chemical cues that modulate the waveform of the flagellar beat. Most previous studies have characterized the flagellar beat and swimming trajectories in two space dimensions (2D) at a water/glass interface. Here, using refined holographic imaging methods, we report high-quality recordings of three-dimensional (3D) flagellar bending waves. As predicted by theory, we observed that an asymmetric and planar flagellar beat results in a circular swimming path, whereas a symmetric and non-planar flagellar beat results in a twisted-ribbon swimming path. During swimming in 3D, human sperm flagella exhibit torsion waves characterized by maxima at the low curvature regions of the flagellar wave. We suggest that these torsion waves are common in nature and that they are an intrinsic property of beating axonemes. We discuss how 3D beat patterns result in twisted-ribbon swimming paths. This study provides new insight into the axoneme dynamics, the 3D flagellar beat, and the resulting swimming behavior.

Physics of microswimmers--single particle motion and collective behavior: a review.

Locomotion and transport of microorganisms in fluids is an essential aspect of life. Search for food, orientation toward light, spreading of off-spring, and the formation of colonies are only possible due to locomotion. Swimming at the microscale occurs at low Reynolds numbers, where fluid friction and viscosity dominates over inertia. Here, evolution achieved propulsion mechanisms, which overcome and even exploit drag. Prominent propulsion mechanisms are rotating helical flagella, exploited by many bacteria, and snake-like or whip-like motion of eukaryotic flagella, utilized by sperm and algae. For artificial microswimmers, alternative concepts to convert chemical energy or heat into directed motion can be employed, which are potentially more efficient. The dynamics of microswimmers comprises many facets, which are all required to achieve locomotion. In this article, we review the physics of locomotion of biological and synthetic microswimmers, and the collective behavior of their assemblies. Starting from individual microswimmers, we describe the various propulsion mechanism of biological and synthetic systems and address the hydrodynamic aspects of swimming. This comprises synchronization and the concerted beating of flagella and cilia. In addition, the swimming behavior next to surfaces is examined. Finally, collective and cooperate phenomena of various types of isotropic and anisotropic swimmers with and without hydrodynamic interactions are discussed.

Dynamics of a polymer chain confined in a membrane.

We present a Brownian dynamics theory with full hydrodynamics (Stokesian dynamics) for a Gaussian polymer chain embedded in a liquid membrane which is surrounded by bulk solvent and walls. The mobility tensors are derived in Fourier space for the two geometries, namely, a free membrane embedded in a bulk fluid, and a membrane sandwiched by the two walls. Within the preaveraging approximation, a new expression for the diffusion coefficient of the polymer is obtained for the free-membrane geometry. We also carry out a Rouse normal mode analysis to obtain the relaxation time and the dynamical structure factor. For large polymer size, both quantities show Zimm-like behavior in the free-membrane case, whereas they are Rouse-like for the sandwiched membrane geometry. We use the scaling argument to discuss the effect of excluded-volume interactions on the polymer relaxation time.

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'.

The conformation of fluid membranes: Monte Carlo simulations.

The conformation and scaling properties of self-avoiding fluid membranes with an extrinsic bending rigidity kappa were studied with the use of Monte Carlo methods. For kappa = 0, the results are consistent with branched polymer behavior at large length scales. There is a smooth crossover from a crumpled to an extended state with increasing kappa, with a peak in the specific heat when the persistence length reaches the system size. The scale-dependent effective bending rigidity is a decreasing function of system size for all bare rigidities. These results indicate that fluid membranes are always crumpled at sufficiently long length scales.

Nano- and microparticles at fluid and biological interfaces.

Systems with interfaces are abundant in both technological applications and biology. While a fluid interface separates two fluids, membranes separate the inside of vesicles from the outside, the interior of biological cells from the environment, and compartmentalize cells into organelles. The physical properties of interfaces are characterized by interface tension, those of membranes are characterized by bending and stretching elasticity. Amphiphilic molecules like surfactants that are added to a system with two immiscible fluids decrease the interface tension and induce a bending rigidity. Lipid bilayer membranes of vesicles can be stretched or compressed by osmotic pressure; in biological cells, also the presence of a cytoskeleton can induce membrane tension. If the thickness of the interface or the membrane is small compared with its lateral extension, both can be described using two-dimensional mathematical surfaces embedded in three-dimensional space. We review recent work on the interaction of particles with interfaces and membranes. This can be micrometer-sized particles at interfaces that stabilise emulsions or form colloidosomes, as well as typically nanometer-sized particles at membranes, such as viruses, parasites, and engineered drug delivery systems. In both cases, we first discuss the interaction of single particles with interfaces and membranes, e.g. particles in external fields, non-spherical particles, and particles at curved interfaces, followed by interface-mediated interaction between two particles, many-particle interactions, interface and membrane curvature-induced phenomena, and applications.

Dynamic regimes of fluids simulated by multiparticle-collision dynamics.

We investigate the hydrodynamic properties of a fluid simulated with a mesoscopic solvent model. Two distinct regimes are identified, the "particle regime" in which the dynamics is gaslike and the "collective regime" where the dynamics is fluidlike. This behavior can be characterized by the Schmidt number, which measures the ratio between viscous and diffusive transport. Analytical expressions for the tracer diffusion coefficient, which have been derived on the basis of a molecular-chaos assumption, are found to describe the simulation data very well in the particle regime, but important deviations are found in the collective regime. These deviations are due to hydrodynamic correlations. The model is then extended in order to investigate self-diffusion in colloidal dispersions. We study first the transport properties of heavy pointlike particles in the mesoscopic solvent, as a function of their mass and number density. Second, we introduce excluded-volume interactions among the colloidal particles and determine the dependence of the diffusion coefficient on the colloidal volume fraction for different solvent mean-free paths. In the collective regime, the results are found to be in good agreement with previous theoretical predictions based on Stokes hydrodynamics and the Smoluchowski equation.

Stability of inverse bicontinuous cubic phases in lipid-water mixtures

We investigate the stability of seven inverse bicontinuous cubic phases [ G, D, P, C(P), S, I-WP, F-RD] in lipid-water mixtures based on a curvature model of membranes. Lipid monolayers are described by parallel surfaces to triply periodic minimal surfaces. The phase behavior is determined by the distribution of the Gaussian curvature on the minimal surface and the porosity of each structure. Only G, D, and P are found to be stable, and to coexist along a triple line. The calculated phase diagram agrees very well with experimental results for 2:1 lauric acid/DLPC.

Elastic properties of polymer interfaces: aggregation of pure diblock, mixed diblock, and triblock copolymers.

Block copolymers adsorbing to an interface between two immiscible homopolymers modify the elastic constants of this interface. Within self-consistent field calculations for Gaussian chains, we determine how the bending constants vary in dependence on the block copolymer concentration and architecture. Four phenomena are discussed. (i) When a tricritical or isotropic Lifshitz critical point is approached in a ternary mixture by varying the concentration of diblock copolymers or changing temperature, the elastic constants vanish. We determine the corresponding power laws, and show that the de Gennes-Taupin criterium for the stability of lamellar phases against undulations and the Ginzburg-Landau criterium for bulk fluctuations yield identical predictions for the validity of the mean-field approximation. (ii) Addition of a small amount of diblock copolymers modifies the bending rigidity. If the diblock copolymers are comparable in length to the homopolymers, adsorption of the diblocks reduces the bending rigidity. If the diblocks are much longer, they increase the bending rigidity. Only for an extreme ratio of chain lengths (>100), the predictions for polymers tethered to an infinitely thin, impenetrable sheet become accurate. (iii) Mixtures of short and long symmetric diblock copolymers are studied, as well as mixtures of two asymmetric diblocks, which are obtained by exchanging the long and short ends. The saddle-splay modulus is found to be the same in both mixtures, while the bending rigidity is significantly smaller in the latter case. (iv) The role of the block copolymer architecture is studied by comparing the effect of triblock copolymers with the effect of diblocks with the same overall length and composition. We propose that triblock copolymers are a very efficient way to control the spontaneous curvature of an interface.

Mesoscopic solvent simulations: multiparticle-collision dynamics of three-dimensional flows.

A recently developed mesoscopic solvent model with multiparticle-collision dynamics is applied to three-dimensional solvent flows in a channel with and without a spherical obstacle. The advantage of a gravitationally driven flow of the solvent over the flow induced by a pressure gradient in the calculation of the solvent viscosity is demonstrated. Three different algorithms for stochastic collision steps are investigated and compared. In particular, we have examined an alternative algorithm with relative velocities drawn from a Maxwell-Boltzmann distribution at each collision step. This algorithm increases the numerical efficiency of the mesoscopic model for solvent flows with low and intermediate Reynolds numbers. Our simulation results for the recirculation length of stationary vortices behind a spherical obstacle are in good agreement with the previous experimental measurements.

Adsorption of monovalent and multivalent cations and anions on DNA molecules.

Adsorption of monovalent and multivalent cations and anions on a deoxyribose nucleic acid (DNA) molecule from a salt solution is investigated by computer simulation. The ions are modeled as charged hard spheres, the DNA molecule as a point charge pattern following the double-helical phosphate strands. The geometrical shape of the DNA molecules is modeled on different levels ranging from a simple cylindrical shape to structured models which include the major and minor grooves between the phosphate strands. The densities of the ions adsorbed on the phosphate strands in the major and in the minor grooves are calculated. First, we find that the adsorption pattern on the DNA surface depends strongly on its geometrical shape: counterions adsorb preferentially along the phosphate strands for a cylindrical model shape, but in the minor groove for a geometrically structured model. Second, we find that an addition of monovalent salt ions results in an increase of the charge density in the minor groove while the total charge density of ions adsorbed in the major groove stays unchanged. The adsorbed ion densities are highly structured along the minor groove while they are almost smeared along the major groove. Furthermore, for a fixed amount of added salt, the major-groove cationic charge is independent of the counterion valency. For increasing salt concentration the major groove is neutralized while the total charge adsorbed in the minor groove is constant. DNA overcharging is detected for multivalent salts. Simulations for larger ion radii, which mimic the effect of ion hydration, indicate an increased adsorbtion of cations in the major groove.

Attraction between DNA molecules mediated by multivalent ions.

The effective force between two parallel DNA molecules is calculated as a function of their mutual separation for different valencies of counterion and salt ions and different salt concentrations. Computer simulations of the primitive model are used and the shape of the DNA molecules is accurately modeled using different geometrical shapes. We find that multivalent ions induce a significant attraction between the DNA molecules whose strength can be tuned by the averaged valency of the ions. The physical origin of the attraction is traced back either to electrostatics or to entropic contributions. For multivalent counterions and monovalent salt ions, we find a salt-enhanced repulsion effect: the force is first attractive but gets repulsive with increasing salt concentration. Furthermore, we show that the multivalent-ion-induced attraction does not necessarily correlate with DNA overcharging.

Budding of crystalline domains in fluid membranes.

Crystalline domains embedded in fluid membrane vesicles are studied by Monte Carlo simulations of dynamically triangulated surfaces and by scaling arguments. A budding transition from a caplike state to a budded shape is observed for increasing spontaneous curvature C0 of the crystalline domain as well as increasing line tension lambda. The location of the budding transition is determined as a function of C0, lambda, and the radius R(A) of the crystalline domain. In contrast to previous theoretical predictions, it is found that budding occurs at a value of the spontaneous curvature C0, that is always a decreasing function of the domain size R(A). Several characteristic scaling regimes are predicted. The distribution of five- and sevenfold disclinations as the budding transition is approached is determined, and the dynamics of the generation of defects is studied.

Helical polymer in cylindrical confining geometries.

Using an algorithm for simulating equilibrium configurations, we study a fluctuating helical polymer either (i) contained in a cylindrical pore or (ii) wound around a cylindrical rod. We work in the regime where both the contour length and the persistence length of the helical polymer are much larger than the diameter of the cylinder. In case (i) we calculate the free energy of confinement and interpret it in terms of a wormlike chain in a pore with an effective diameter that depends on the parameters of the helix. In case (ii) we consider the possibility that one end of the helical polymer escapes from the rod and wanders away. The average numbers of turns at which the helix escapes or intersects the rod are measured in the simulations, as a function of the pitch p(0) . The behavior for large and small p(0) is explained with simple scaling arguments.

Semiflexible polymer in a uniform force field in two dimensions.

The conformational properties of a semiflexible polymer chain, anchored at one end in a uniform force field, are studied in a simple two-dimensional model. Recursion relations are derived for the partition function and then iterated numerically. We calculate the angular fluctuations of the polymer about the direction of the force field and the average polymer configuration as functions of the bending rigidity, chain length, chain orientation at the anchoring point, and field strength.

Near-surface structure of a bicontinuous microemulsion with a transition region.

The lamellar ordering of bicontinuous microemulsions adjacent to a planar hydrophilic wall is investigated experimentally by grazing-incidence small-angle neutron scattering and theoretically by computer simulations. It is shown that precise depth information in neutron scattering can be obtained by tuning the scattering length density of the overall microemulsion. Neutron reflectometry completes the characterization. The nucleation of a lamellar phase at the wall is observed, and a perforated lamellar transition region is identified at the lamellar-microemulsion interface. The thickness of the lamellar region is about 400 Å, which corresponds to two bilayers.