Brownian motion and the hydrodynamic friction tensor for colloidal particles of complex shape

Effective action for dissipative and nonholonomic systems

We show that the action of a dynamical system can be supplemented by an effective action for its environment to reproduce arbitrary coordinate dependent ohmic dissipation and gyroscopic forces. The action is a generalization of the harmonic bath model and describes a set of massless interacting scalar fields in an auxiliary space coupled to the original system at the boundary. A certain limit of the model implements nonholonomic constraints. In the case of dynamics with nonlinearly realized symmetries the effective action takes the form of a two-dimensional nonlinear σ model. It provides a basis for application of path integral methods to general dissipative and nonholonomic systems.

Orientation of graphene nanosheets in suspension under an electric field: theoretical model and molecular dynamic simulations

Orientation regulation of nanoparticles in a suspension by an electric field is a powerful tool to tune its mechanical, thermal, optical, electrical properties etc. However, how molecular modification can affect the orientation of two-dimensional nanoparticles is still unclear. In this paper, the influence of molecular modification on the orientation of graphene nanosheets (GNS) in water was investigated through theoretical analyses and molecular dynamics (MD) simulations. Firstly, a new orientation angle model was proposed, which considers hydration effects, dipole moments and resistance torque. Then, MD simulations were conducted to investigate the effects of position, direction, type, and number of functional groups on the orientation of GNS. The trend observed in MD simulations is consistent with the proposed theoretical model. The results reveal that, under the combined influence of the dipole moment and hydration effects, the modification with hydrophilic functional groups can reduce the orientation angle from 21.31° to 8.34°, while the modification with hydrophobic functional groups increases it to 26.43°. Among the hydrophilic functional groups, orientation of hydroxylated GNS is the best. With an increase in the number of hydroxyl groups, orientation angle is decreased from 12.61° to 8.34°. This work can provide valuable guidance for the design of high-performance suspensions and composites, such as thermal smart materials with adjustable thermal conductivity and intelligent devices with tailored capabilities.

Anisotropic Diffusion of Elongated Particles in Active Coherent Flows

The study of particle diffusion, a classical conundrum in scientific inquiry, holds manifold implications for various real-world applications. Particularly within the domain of active flows, where the motion of self-propelled particles instigates fluid movement, extensive research has been dedicated to unraveling the dynamics of passive spherical particles. This scrutiny has unearthed intriguing phenomena, such as superdiffusion at brief temporal scales and conventional diffusion at longer intervals. In contrast to the spherical counterparts, anisotropic particles, which manifest directional variations, are prevalent in nature. Although anisotropic behavior in passive fluids has been subject to exploration, enigmatic aspects persist in comprehending the interplay of anisotropic particles within active flows. This research delves into the intricacies of anisotropic passive particle diffusion, exposing a notable escalation in translational and rotational diffusion coefficients, as well as the superdiffusion index, contingent upon bacterial concentration. Through a detailed examination of particle coordinates, the directional preference of particle diffusion is not solely dependent on the particle length, but rather determined by the ratio of the particle length to the associated length scale of the background flow field. These revelations accentuate the paramount importance of unraveling the nuances of anisotropic particle diffusion within the context of active flows. Such insights not only contribute to the fundamental understanding of particle dynamics, but also have potential implications for a spectrum of applications.

Decoding the hydrodynamic properties of microscale helical propellers from Brownian fluctuations

The complex motility of bacteria, ranging from single-swimmer behaviors such as chemotaxis to collective dynamics, including biofilm formation and active matter phenomena, is driven by their microscale propellers. Despite extensive study of swimming flagellated bacteria, the hydrodynamic properties of their helical-shaped propellers have never been directly measured. The primary challenges to directly studying microscale propellers are 1) their small size and fast, correlated motion, 2) the necessity of controlling fluid flow at the microscale, and 3) isolating the influence of a single propeller from a propeller bundle. To solve the outstanding problem of characterizing the hydrodynamic properties of these propellers, we adopt a dual statistical viewpoint that connects to the hydrodynamics through the fluctuation-dissipation theorem (FDT). We regard the propellers as colloidal particles and characterize their Brownian fluctuations, described by 21 diffusion coefficients for translation, rotation, and correlated translation-rotation in a static fluid. To perform this measurement, we applied recent advances in high-resolution oblique plane microscopy to generate high-speed volumetric movies of fluorophore-labeled, freely diffusing Escherichia coli flagella. Analyzing these movies with a bespoke helical single-particle tracking algorithm, we extracted trajectories, calculated the full set of diffusion coefficients, and inferred the average propulsion matrix using a generalized Einstein relation. Our results provide a direct measurement of a microhelix's propulsion matrix and validate proposals that the flagella are highly inefficient propellers, with a maximum propulsion efficiency of less than 3%. Our approach opens broad avenues for studying the motility of particles in complex environments where direct hydrodynamic approaches are not feasible.

Dynamics and interactions of Quincke roller clusters: From orbits and flips to excited states

Active matter systems may be characterized by the conversion of energy into active motion, e.g., the self-propulsion of microorganisms. Artificial active colloids form models that exhibit essential properties of more complex biological systems but are amenable to laboratory experiments. While most experimental models consist of spheres, active particles of different shapes are less understood. Furthermore, interactions between these anisotropic active colloids are even less explored. Here, we investigate the motion of active colloidal clusters and the interactions between them. We focus on self-assembled dumbbells and trimers powered by an external dc electric field. For dumbbells, we observe an activity-dependent behavior of spinning, circular, and orbital motions. Moreover, collisions between dumbbells lead to the hierarchical self-assembly of tetramers and hexamers, both of which form rotational excited states. On the other hand, trimers exhibit flipping motion that leads to trajectories reminiscent of a honeycomb lattice.

3-D rotation tracking from 2-D images of spherical colloids with textured surfaces

Tracking the three-dimensional rotation of colloidal particles is essential to elucidate many open questions, e.g. concerning the contact interactions between particles under flow, or the way in which obstacles and neighboring particles affect self-propulsion in active suspensions. In order to achieve rotational tracking, optically anisotropic particles are required. We synthesise here rough spherical colloids that present randomly distributed fluorescent asperities and track their motion under different experimental conditions. Specifically, we propose a new algorithm based on a 3-D rotation registration, which enables us to track the 3-D rotation of our rough colloids at short time-scales, using time series of 2-D images acquired at high frame rates with a conventional wide-field microscope. The method is based on the image correlation between a reference image and rotated 3-D prospective images to identify the most likely angular displacements between frames. We first validate our approach against simulated data and then apply it to the cases of: particles flowing through a capillary, freely diffusing at solid-liquid and liquid-liquid interfaces, and self-propelling above a substrate. By demonstrating the applicability of our algorithm and sharing the code, we hope to encourage further investigations in the rotational dynamics of colloidal systems.

Effect of curvature on the diffusion of colloidal bananas

Anisotropic colloidal particles exhibit complex dynamics which play a crucial role in their functionality, transport, and phase behavior. In this Letter, we investigate the two-dimensional diffusion of smoothly curved colloidal rods-also known as colloidal bananas-as a function of their opening angle α. We measure the translational and rotational diffusion coefficients of the particles with opening angles ranging from 0^{∘} (straight rods) to nearly 360^{∘}(closed rings). In particular, we find that the anisotropic diffusion of the particles varies nonmonotonically with their opening angle and that the axis of fastest diffusion switches from the long to the short axis of the particles when α>180^{∘}. We also find that the rotational diffusion coefficient of nearly closed rings is approximately an order of magnitude higher than that of straight rods of the same length. Finally, we show that the experimental results are consistent with slender body theory, indicating that the dynamical behavior of the particles arises primarily from their local drag anisotropy. These results highlight the impact of curvature on the Brownian motion of elongated colloidal particles, which must be taken into account when seeking to understand the behavior of curved colloidal particles.

Inertial dynamics of an active Brownian particle

Active Brownian motion commonly assumes spherical overdamped particles. However, self-propelled particles are often neither symmetric nor overdamped yet underlie random fluctuations from their surroundings. Active Brownian motion has already been generalized to include asymmetric particles. Separately, recent findings have shown the importance of inertial effects for particles of macroscopic size or in low-friction environments. We aim to consolidate the previous findings into the general description of a self-propelled asymmetric particle with inertia. We derive the Langevin equation of such a particle as well as the corresponding Fokker-Planck equation. Furthermore, a formula is presented that allows reconstructing the hydrodynamic resistance matrix of the particle by measuring its trajectory. Numerical solutions of the Langevin equation show that, independently of the particle's shape, the noise-free trajectory at zero temperature starts with an inertial transition phase and converges to a circular helix. We discuss this universal convergence with respect to the helical motion that many microorganisms exhibit.

Spontaneous chiralization of polar active particles

Polar active particles constitute a wide class of active matter that is able to propel along a preferential direction, given by their polar axis. Here, we demonstrate a generic active mechanism that leads to their spontaneous chiralization through a symmetry-breaking instability. We find that the transition of an active particle from a polar to a chiral symmetry is characterized by the emergence of active rotation and of circular trajectories. The instability is driven by the advection of a solute that interacts differently with the two portions of the particle surface and it occurs through a supercritical pitchfork bifurcation.

Time-dependent inertia of self-propelled particles: The Langevin rocket

Many self-propelled objects are large enough to exhibit inertial effects but still suffer from environmental fluctuations. The corresponding basic equations of motion are governed by active Langevin dynamics, which involve inertia, friction, and stochastic noise for both the translational and orientational degrees of freedom coupled via the self-propulsion along the particle orientation. In this paper, we generalize the active Langevin model to time-dependent parameters and explicitly discuss the effect of time-dependent inertia for achiral and chiral particles. Realizations of this situation are manifold, ranging from minirockets (which are self-propelled by burning their own mass), to dust particles in plasma (which lose mass by evaporating material), to walkers with expiring activity. Here we present analytical solutions for several dynamical correlation functions, such as mean-square displacement and orientational and velocity autocorrelation functions. If the parameters exhibit a slow power law in time, we obtain anomalous superdiffusion with a nontrivial dynamical exponent. Finally, we constitute the "Langevin rocket" model by including orientational fluctuations in the traditional Tsiolkovsky rocket equation. We calculate the mean reach of the Langevin rocket and discuss different mass ejection strategies to maximize it. Our results can be tested in experiments on macroscopic robotic or living particles or in self-propelled mesoscopic objects moving in media of low viscosity, such as complex plasma.

Height distribution and orientation of colloidal dumbbells near a wall

Geometric confinement strongly influences the behavior of microparticles in liquid environments. However, to date, nonspherical particle behaviors close to confining boundaries, even as simple as planar walls, remain largely unexplored. Here, we measure the height distribution and orientation of colloidal dumbbells above walls by means of digital in-line holographic microscopy. We find that while larger dumbbells are oriented almost parallel to the wall, smaller dumbbells of the same material are surprisingly oriented at preferred angles. We determine the total height-dependent force acting on the dumbbells by considering gravitational effects and electrostatic particle-wall interactions. Our modeling reveals that at specific heights both net forces and torques on the dumbbells are simultaneously below the thermal force and energy, respectively, which makes the observed orientations possible. Our results highlight the rich near-wall dynamics of nonspherical particles and can further contribute to the development of quantitative frameworks for arbitrarily shaped microparticle dynamics in confinement.

Collective motion of chiral Brownian particles controlled by a circularly-polarized laser beam

We demonstrate the emergence of circular collective motion in a system of spherical light-propelled Brownian particles. Light-propulsion occurs as consequence of the coupling between the chirality of polymeric particles - left (L)- or right (R)-type - and the circularly-polarized light that irradiates them. Irradiation with light that has the same helicity as the particle material leads to a circular cooperative vortical motion between the chiral Brownian particles. In contrast, opposite circular-polarization does not induce such coupling among the particles but only affects their Brownian motion. The mean angular momentum of each particle has a value and sign that distinguishes between chiral activity dynamics and typical Brownian motion. These outcomes have relevant implications for chiral separation technologies and provide new strategies for optical torque tunability in mesoscopic optical array systems, micro- and nanofabrication of light-activated engines with selective control and collective motion.

Brownian fluctuations and hydrodynamics of a microhelix near a solid wall

We combine two-photon lithography and optical tweezers to investigate the Brownian fluctuations and propeller characteristics of a microfabricated helix. From the analysis of mean squared displacements and time correlation functions we recover the components of the full mobility tensor. We find that Brownian motion displays correlations between angular and translational fluctuations from which we can directly measure the hydrodynamic coupling coefficient that is responsible for thrust generation. By varying the distance of the microhelices from a no-slip boundary we can systematically measure the effects of a nearby wall on the resistance matrix. Our results indicate that a rotated helix moves faster when a nearby no-slip boundary is present, providing a quantitative insight on thrust enhancement in confined geometries for both synthetic and biological microswimmers.

Confinement effects on the dynamics of a rigid particle in a nanochannel

The transport of nanoparticles in confined geometries plays a crucial role in several technological applications ranging from nanopore sensors to filtration membranes. Here we describe a Brownian approach to simulate the motion of a rigid-body nanoparticle of an arbitrary shape under confinement. A quaternion formulation is used for the nanoparticle orientation, and the corresponding overdamped Langevin equation, completed by the proper fluctuation-dissipation relation, is derived. The hydrodynamic mobility matrix is obtained via dissipative particle dynamics simulation equipped with a new method for enforcing the no-slip boundary condition for curved moving solid-liquid interfaces. As an application, we analyzed the motion of a nanoparticle in a cylindrical channel under the action of external fields. We show that both axial effective diffusion and rotational diffusion decrease with confinement.

Diffusion Tensors of Arbitrary-Shaped Nanoparticles in Fluid by Molecular Dynamics Simulation

The anisotropic diffusive behavior of nanoparticles with complex shapes attracts great interest due to its potential applications in many fields ranging from bionics to aeronautic industry. Although molecular dynamics (MD) simulations are used widely to investigate nanoparticle diffusion properties, universal methods to describe the diffusion process comprehensively are still lacking. Here, we address this problem by introducing diffusion tensor as it can describe translational and rotational diffusion in three dimensions both individually and their coupling. We take carbon triple sphere suspended in argon fluid as our model system. The consistency of our results and velocity autocorrelation function(VAF) method validates our simulations. The coupling between translational and rotational diffusion is observed directly from analyzing diffusion tensor, and quantified by coupling diffusion coefficient. Our simulation reveals non-trivial effect of some factors in diffusion at nanoscale, which was not considered in previous theories. In addition to introducing an effective method to calculate the diffusion tensor in MD simulations, our work also provides insights for understanding the diffusion process of arbitrary-shaped particles in nanoengineering.

Mesoscopic simulations of anisotropic chemically powered nanomotors

Chemically powered self-propelled colloids generate a motor force by converting locally a source of energy into directed motion, a process that has been explored both in experiments and in computational models. The use of active colloids as building blocks for nanotechnology opens the doors to interesting applications, provided we understand the behavior of these elementary constituents. We build a consistent mesoscopic simulation model for self-propelled colloids of complex shape with the aim of resolving the coupling between their translational and rotational motion. Considering a passive L-shaped colloidal particle, we study its Brownian dynamics and locate its center of hydrodynamics, the tracking point at which translation and rotation decouple. The active L particle displays the same circling trajectories that have been found experimentally, a result which we compare with the Brownian dynamics model. We put forward the role of hydrodynamics by comparing our results with a fluid model in which the particles' velocities are reset randomly. There, the trajectories only display random orientations. We obtain these original simulation results without any parametrization of the algorithm, which makes it a useful method for the preliminary study of active colloids, prior to experimental work.

Moiré and honeycomb lattices through self-assembly of hard-core/soft-shell microgels: experiment and simulation

Moiré and honeycomb lattices result from the sequential double deposition of monolayers of core/shell microgels in dependence of the drying conditions.

Wet-Chemical Synthesis of Chiral Colloids

We disclose a method for the synthesis of chiral colloids from spontaneously formed hollow sugar-surfactant microtubes with internally confined mobile colloidal spheres. Key feature of our approach is the grafting of colloid surfaces with photoresponsive coumarin moieties, which allow for UV-induced, covalent clicking of colloids into permanent chains, with morphologies set by the colloid-to-tube diameter ratio. Subsequent dissolution of tube confinement yields aqueous suspensions that comprise bulk quantities of a variety of linear chains, including single helical chains of polystyrene colloids. These colloidal equivalents of chiral (DNA) molecules are intended for microscopic study of chiral dynamics on a single-particle level.

Experimental Assessment and Model Validation on How Shape Determines Sedimentation and Diffusion of Colloidal Particles

Understanding the hydrodynamics of colloids with complex shapes is of equal importance to widespread practical applications and fundamental scientific problems, such as gelation, crystallization, and phase behavior. Building on previous work, we present a comprehensive study of sedimentation, diffusion, intrinsic viscosities, and other shape-dependent quantities of clusters built from spherical nanoparticles. Cluster preparation is accomplished by assembling surface-modified polystyrene particles on evaporating emulsion droplets. This results in supracolloids that exhibit well-defined configurations, which are governed by the number of constituent particles. Sorting into uniform cluster fractions is achieved through centrifugation of the cluster mixture in a density gradient. Sedimentation coefficients are elucidated by differential centrifugal sedimentation. Rotational and translational diffusion of the clusters are investigated by polarized and depolarized dynamic light scattering. The experimental results are compared to data obtained via a bead-shell model suitable for predicting hydrodynamic quantities of particles with arbitrary shapes. The experimental data are in excellent agreement with the predictions from hydrodynamic modeling. The variety of investigated shapes shows the robustness of our approach and provides a complete picture of the hydrodynamic behavior of complex particles.

Different bending models predict different dynamics of sedimenting elastic trumbbells

The main goal of this paper is to examine theoretically and numerically the impact of a chosen bending model on the dynamics of elastic filaments settling in a viscous fluid under gravity at low-Reynolds-number. We use the bead-spring approximation of a filament and the Rotne-Prager mobility matrix to describe hydrodynamic interactions between the beads. We analyze the dynamics of trumbbells, for which bending angles are typically larger than for thin and long filaments. Each trumbbell is made of three beads connected by springs and it exhibits a bending resistance, described by the harmonic or - alternatively - by the 'cosine' (also called the Kratky-Porod) bending models, both often used in the literature. Using the harmonic bending potential, and coupling it to the spring potential by the Young's modulus, we find simple benchmark solutions: stable stationary configurations of a single elastic trumbbell and attraction of two elastic trumbbells towards a periodic long-lasting orbit. As the most significant result of this paper, we show that for very elastic trumbbells at the same initial conditions, the Kratky-Porod bending potential can lead to qualitatively and quantitatively different spurious dynamics, with artificially large bending angles and unrealistic shapes. We point out that for the bead models of an elastic filament, the range of applicability of the Kratky-Porod model might not go beyond bending angles smaller than π/2 for touching beads and beyond an even much lower value for beads well-separated from each other. The existence of stable stationary configurations of elastic trumbbells and a family of periodic oscillations of two elastic trumbbells are very important findings on their own.

Re-entrant bimodality in spheroidal chiral swimmers in shear flow

We use a continuum model to report on the behavior of a dilute suspension of chiral swimmers subject to externally imposed shear in a planar channel. Swimmer orientation in response to the imposed shear can be characterized by two distinct phases of behavior, corresponding to unimodal or bimodal distribution functions for swimmer orientation along the channel. These phases indicate the occurrence (or not) of a population splitting phenomenon changing the swimming direction of a macroscopic fraction of active particles to the exact opposite of that dictated by the imposed flow. We present a detailed quantitative analysis elucidating the complexities added to the population splitting behavior of swimmers when they are chiral. In particular, the transition from unimodal to bimodal and vice versa are shown to display a re-entrant behavior across the parameter space spanned by varying the chiral angular speed. We also present the notable effects of particle aspect ratio and self-propulsion speed on system phase behavior and discuss potential implications of our results in applications such as swimmer separation/sorting.

Physical foundation of the fluid particle dynamics method for colloid dynamics simulation

Colloid dynamics is significantly influenced by many-body hydrodynamic interactions mediated by a suspending fluid. However, theoretical and numerical treatments of such interactions are extremely difficult. To overcome this situation, we developed a fluid particle dynamics (FPD) method [H. Tanaka and T. Araki, Phys. Rev. Lett., 2000, 35, 3523], which is based on two key approximations: (i) a colloidal particle is treated as a highly viscous particle and (ii) the viscosity profile is described by a smooth interfacial profile function. Approximation (i) makes our method free from the solid-fluid boundary condition, significantly simplifying the treatment of many-body hydrodynamic interactions while satisfying the incompressible condition without the Stokes approximation. Approximation (ii) allows us to incorporate an extra degree of freedom in a fluid, e.g., orientational order and concentration, as an additional field variable. Here, we consider two fundamental problems associated with these approximations. One is the introduction of thermal noise and the other is the incorporation of coupling of the colloid surface with an order parameter introduced into a fluid component, which is crucial when considering colloidal particles suspended in a complex fluid. Here, we show that our FPD method makes it possible to simulate colloid dynamics properly while including full hydrodynamic interactions, inertia effects, incompressibility, thermal noise, and additional degrees of freedom of a fluid, which may be relevant for wide applications in colloidal and soft matter science.

Dynamics of a mechanical system with multiple degrees of freedom out of thermal equilibrium

Out of thermal equilibrium, an environment imposes effective mechanical forces on nanofabricated devices as well as on microscopic chemical or biological systems. Here we address the question of how to calculate these forces together with the response of the system from first principles. We show that an ideal gaslike environment, even near thermal equilibrium, can enforce a specific steady state on the system by creating effective potentials in otherwise homogeneous space. An example of stable and unstable rectifications of thermal fluctuations is presented using a modified Feynman-Smoluchowski ratchet with two degrees of freedom. Moreover, the stability of a steady configuration depends on its chiral symmetry. The transition rate probabilities and the corresponding kinetic equations are derived for a complex mechanical system with arbitrary degrees of freedom. This work, therefore, extends the applicability of mechanical systems as a toy model playground of statistical physics for active and living matter with multiple degrees of freedom.

The non-Gaussian tops and tails of diffusing boomerangs

Experiments involving the two-dimensional passive diffusion of colloidal boomerangs tracked off their centre of mobility have shown striking non-Gaussian tails in their probability distribution function [Chakrabarty et al., Soft Matter, 2016, 12, 4318]. This in turn can lead to anomalous diffusion characteristics, including mean drift. In this paper, we develop a general theoretical explanation for these measurements. The idea relies on calculating the two-dimensional probability densities at the centre of mobility of the particle, where all distributions are Gaussian, and then transforming them to a different reference point. Our model clearly captures the experimental results, without any fitting parameters, and demonstrates that the one-dimensional probability distributions may also exhibit strongly non-Gaussian tops. These results indicate that the choice of tracking point can cause a considerable departure from Gaussian statistics, potentially causing some common modelling techniques to fail.

Brownian motion of a circle swimmer in a harmonic trap

We study the dynamics of a Brownian circle swimmer with a time-dependent self-propulsion velocity in an external temporally varying harmonic potential. For several situations, the noise-free swimming paths, the noise-averaged mean trajectories, and the mean-square displacements are calculated analytically or by computer simulation. Based on our results, we discuss optimal swimming strategies in order to explore a maximum spatial range around the trap center. In particular, we find a resonance situation for the maximum escape distance as a function of the various frequencies in the system. Moreover, the influence of the Brownian noise is analyzed by comparing noise-free trajectories at zero temperature with the corresponding noise-averaged trajectories at finite temperature. The latter reveal various complex self-similar spiral or rosette-like patterns. Our predictions can be tested in experiments on artificial and biological microswimmers under dynamical external confinement.

Automated tracking of colloidal clusters with sub-pixel accuracy and precision

Quantitative tracking of features from video images is a basic technique employed in many areas of science. Here, we present a method for the tracking of features that partially overlap, in order to be able to track so-called colloidal molecules. Our approach implements two improvements into existing particle tracking algorithms. Firstly, we use the history of previously identified feature locations to successfully find their positions in consecutive frames. Secondly, we present a framework for non-linear least-squares fitting to summed radial model functions and analyze the accuracy (bias) and precision (random error) of the method on artificial data. We find that our tracking algorithm correctly identifies overlapping features with an accuracy below 0.2% of the feature radius and a precision of 0.1 to 0.01 pixels for a typical image of a colloidal cluster. Finally, we use our method to extract the three-dimensional diffusion tensor from the Brownian motion of colloidal dimers.

Rotational friction of dipolar colloids measured by driven torsional oscillations

Despite its prominent role in the dynamics of soft materials, rotational friction remains a quantity that is difficult to determine for many micron-sized objects. Here, we demonstrate how the Stokes coefficient of rotational friction can be obtained from the driven torsional oscillations of single particles in a highly viscous environment. The idea is that the oscillation amplitude of a dipolar particle under combined static and oscillating fields provides a measure for the Stokes friction. From numerical studies we derive a semi-empirical analytic expression for the amplitude of the oscillation, which cannot be calculated analytically from the equation of motion. We additionally demonstrate that this expression can be used to experimentally determine the rotational friction coefficient of single particles. Here, we record the amplitudes of a field-driven dipolar Janus microsphere with optical microscopy. The presented method distinguishes itself in its experimental and conceptual simplicity. The magnetic torque leaves the local environment unchanged, which contrasts with other approaches where, for example, additional mechanical (frictional) or thermal contributions have to be regarded.

Short-time dynamics of monomers and dimers in quasi-two-dimensional colloidal mixtures

We report on the short-time dynamics in colloidal mixtures made up of monomers and dimers highly confined between two glass plates. At low concentrations, the experimental measurements of colloidal motion agree well with the solution of the Navier-Stokes equation at low Reynolds numbers; the latter takes into account the increase in the drag force on a colloidal particle due to wall-particle hydrodynamic forces. More importantly, we find that the ratio of the short-time diffusion coefficient of the monomer and that of the center of mass of the dimmer is almost independent of both the dimer molar fraction, x_{d}, and the total packing fraction, ϕ, up to ϕ≈0.5. At higher concentrations, this ratio displays a small but systematic increase. A similar physical scenario is observed for the ratio between the parallel and the perpendicular components of the short-time diffusion coefficients of the dimer. This dynamical behavior is corroborated by means of molecular dynamics computer simulations that include explicitly the particle-particle hydrodynamic forces induced by the solvent. Our results suggest that the effects of colloid-colloid hydrodynamic interactions on the short-time diffusion coefficients are almost identical and factorable in both species.

Role of Tensorial Electronic Friction in Energy Transfer at Metal Surfaces

An accurate description of nonadiabatic energy relaxation is crucial for modeling atomistic dynamics at metal surfaces. Interfacial energy transfer due to electron-hole pair excitations coupled to motion of molecular adsorbates is often simulated by Langevin molecular dynamics with electronic friction. Here, we present calculations of the full electronic friction tensor by using first order time-dependent perturbation theory at the density functional theory level. We show that the friction tensor is generally anisotropic and nondiagonal, as found for hydrogen atom on Pd(100) and CO on Cu(100) surfaces. This implies that electron-hole pair induced nonadiabatic coupling at metal surfaces leads to friction-induced mode coupling, therefore, opening an additional channel for energy redistribution. We demonstrate the robustness and accuracy of our results by direct comparison to established methods and experimental data.

Effects of translation-rotation coupling on the displacement probability distribution functions of boomerang colloidal particles

Prior studies have shown that low symmetry particles such as micro-boomerangs exhibit behaviour of Brownian motion rather different from that of high symmetry particles because convenient tracking points (TPs) are usually inconsistent with their center of hydrodynamic stress (CoH) where the translational and rotational motions are decoupled. In this paper we study the effects of the translation-rotation coupling on the displacement probability distribution functions (PDFs) of the boomerang colloid particles with symmetric arm length. By tracking the motions of different points on the particle symmetry axis, we show that as the distance between the TP and the CoH is increased, the effects of translation-rotation coupling becomes pronounced, making the short-time 2D PDF for fixed initial orientation to change from elliptical, to bean and then to crescent shape, and the angle averaged PDFs change from ellipsoidal-particle-like PDF to a shape with a Gaussian top and long displacement tails. We also observed that at long times the PDFs revert to Gaussian. These 2D PDF shapes provide a clear physical picture of the non-zero mean displacements observed in boomerangs particles.

Equilibrium fluid-crystal interfacial free energy of bcc-crystallizing aqueous suspensions of polydisperse charged spheres

The interfacial free energy is a central quantity in crystallization from the metastable melt. In suspensions of charged colloidal spheres, nucleation and growth kinetics can be accurately measured from optical experiments. In previous work, from these data effective nonequilibrium values for the interfacial free energy between the emerging bcc nuclei and the adjacent melt in dependence on the chemical potential difference between melt phase and crystal phase were derived using classical nucleation theory (CNT). A strictly linear increase of the interfacial free energy was observed as a function of increased metastability. Here, we further analyze these data for five aqueous suspensions of charged spheres and one binary mixture. We utilize a simple extrapolation scheme and interpret our findings in view of Turnbull's empirical rule. This enables us to present the first systematic experimental estimates for a reduced interfacial free energy, σ(0,bcc), between the bcc-crystal phase and the coexisting equilibrium fluid. Values obtained for σ(0,bcc) are on the order of a few k(B)T. Their values are not correlated to any of the electrostatic interaction parameters but rather show a systematic decrease with increasing size polydispersity and a lower value for the mixture as compared to the pure components. At the same time, σ(0) also shows an approximately linear correlation to the entropy of freezing. The equilibrium interfacial free energy of strictly monodisperse charged spheres may therefore be still greater.

The raspberry model for hydrodynamic interactions revisited. I. Periodic arrays of spheres and dumbbells

The so-called "raspberry" model refers to the hybrid lattice-Boltzmann and Langevin molecular dynamics scheme for simulating the dynamics of suspensions of colloidal particles, originally developed by Lobaskin and Dünweg [New J. Phys. 6, 54 (2004)], wherein discrete surface points are used to achieve fluid-particle coupling. This technique has been used in many simulation studies on the behavior of colloids. However, there are fundamental questions with regards to the use of this model. In this paper, we examine the accuracy with which the raspberry method is able to reproduce Stokes-level hydrodynamic interactions when compared to analytic expressions for solid spheres in simple-cubic crystals. To this end, we consider the quality of numerical experiments that are traditionally used to establish these properties and we discuss their shortcomings. We show that there is a discrepancy between the translational and rotational mobility reproduced by the simple raspberry model and present a way to numerically remedy this problem by adding internal coupling points. Finally, we examine a non-convex shape, namely, a colloidal dumbbell, and show that the filled raspberry model replicates the desired hydrodynamic behavior in bulk for this more complicated shape. Our investigation is continued in de Graaf et al. [J. Chem. Phys. 143, 084108 (2015)], wherein we consider the raspberry model in the confining geometry of two parallel plates.

Rheology and dynamics of colloidal superballs

Recent advances in colloidal synthesis make it possible to generate a wide array of precisely controlled, non-spherical particles. This provides a unique opportunity to probe the role that particle shape plays in the dynamics of colloidal suspensions, particularly at higher volume fractions, where particle interactions are important. We examine the role of particle shape by characterizing both the bulk rheology and micro-scale diffusion in a suspension of pseudo-cubic silica superballs. Working with these well-characterized shaped colloids, we can disentangle shape effects in the hydrodynamics of isolated particles from shape-mediated particle interactions. We find that the hydrodynamic properties of isolated superballs are marginally different from comparably sized hard spheres. However, shape-mediated interactions modify the suspension microstructure, leading to significant differences in the self-diffusion of the superballs. While this excluded volume interaction can be captured with a rescaling of the superball volume fraction, we observe qualitative differences in the shear thickening behavior of moderately concentrated superball suspensions that defy simple rescaling onto hard sphere results. This study helps to define the unknowns associated with the effects of shape on the rheology and dynamics of colloidal solutions.

Active dipole clusters: From helical motion to fission

The structure of a finite particle cluster is typically determined by total energy minimization. Here we consider the case where a cluster of soft-sphere dipoles becomes active, i.e., when the individual particles exhibit an additional self-propulsion along their dipole moments. We numerically solve the overdamped equations of motion for soft-sphere dipoles in a solvent. Starting from an initial metastable dipolar cluster, the self-propulsion generates a complex cluster dynamics. The final cluster state has in general a structure widely different to the initial one, the details depend on the model parameters and on the protocol of how the self-propulsion is turned on. The center of mass of the cluster moves on a helical path, the details of which are governed by the initial cluster magnetization. An instantaneous switch to a high self-propulsion leads to fission of the cluster. However, fission does not occur if the self-propulsion is increased slowly to high strengths. Our predictions can be verified through experiments with self-phoretic colloidal Janus particles and for macroscopic self-propelled dipoles in a highly viscous solvent.

Can the self-propulsion of anisotropic microswimmers be described by using forces and torques?

The self-propulsion of artificial and biological microswimmers (or active colloidal particles) has often been modelled by using a force and a torque entering into the overdamped equations for the Brownian motion of passive particles. This seemingly contradicts the fact that a swimmer is force-free and torque-free, i.e. that the net force and torque on the particle vanish. Using different models for mechanical and diffusiophoretic self-propulsion, we demonstrate here that the equations of motion of microswimmers can be mapped onto those of passive particles with the shape-dependent grand resistance matrix and formally external effective forces and torques. This is consistent with experimental findings on the circular motion of artificial asymmetric microswimmers driven by self-diffusiophoresis. The concept of effective self-propulsion forces and torques significantly facilitates the understanding of the swimming paths, e.g. for a microswimmer under gravity. However, this concept has its limitations when the self-propulsion mechanism of a swimmer is disturbed either by another particle in its close vicinity or by interactions with obstacles, such as a wall.

Jammed elastic shells - a 3D experimental soft frictionless granular system

We present a new experimental system of monodisperse, soft, frictionless, fluorescent labeled elastic shells for the characterization of structure, universal scaling laws and force networks in 3D jammed matter. The elastic shells in a jammed packing are deformed in such a way that at each contact one of the shells buckles with a dimple and the other remain spherical, closely resembling overlapping spheres. Using confocal microscopy, we obtained 3D stacks of images of shells at different volume fractions which were subsequently processed in ImageJ software to find their coordinates. The determination of 3D coordinates involved three steps: locating the edges of shells in all 2D slices, analyzing their shape and subsequently finding their 2D coordinates, and finally determining their 3D centers by grouping the corresponding 2D coordinates. From this analysis routine we obtained particle coordinates with sub-pixel accuracy. In a contact pair we also identified the shell that underwent buckling forming a dimple by analyzing the intensity profile of a line that connects the centers of particle pairs. The amorphous structure of the packing was analyzed as a function of distance to the jamming threshold by investigating the radial distribution function, bond order parameters, contact numbers and the number of dimples per particle (buckling number), which is a unique property of this system. We find that the power law scaling of the contact number with excess volume fraction deviated from theoretical and computer simulation predictions. In addition, the buckling number also showed a similar scaling as that of the contact number with distance to the jamming transition.

The passive diffusion of Leptospira interrogans

Motivated by recent experimental measurements, the passive diffusion of the bacterium Leptospira interrogans is investigated theoretically. By approximating the cell shape as a straight helix and using the slender-body-theory approximation of Stokesian hydrodynamics, the resistance matrix of Leptospira is first determined numerically. The passive diffusion of the helical cell is then obtained computationally using a Langevin formulation which is sampled in time in a manner consistent with the experimental procedure. Our results are in excellent quantitative agreement with the experimental results with no adjustable parameters.

Brownian motion of arbitrarily shaped particles in two dimensions

We implement microfabricated boomerang particles with unequal arm lengths as a model for nonsymmetric particles and study their Brownian motion in a quasi-two-dimensional geometry by using high-precision single-particle motion tracking. We show that because of the coupling between translation and rotation, the mean squared displacements of a single asymmetric boomerang particle exhibit a nonlinear crossover from short-time faster to long-time slower diffusion, and the mean displacements for fixed initial orientation are nonzero and saturate out at long times. The measured anisotropic diffusion coefficients versus the tracking point position indicate that there exists one unique point, i.e., the center of hydrodynamic stress (CoH), at which all coupled diffusion coefficients vanish. This implies that in contrast to motion in three dimensions where the CoH exists only for high-symmetry particles, the CoH always exists for Brownian motion in two dimensions. We develop an analytical model based on Langevin theory to explain the experimental results and show that among the six anisotropic diffusion coefficients only five are independent because the translation-translation coupling originates from the translation-rotation coupling. Finally, we classify the behavior of two-dimensional Brownian motion of arbitrarily shaped particles into four groups based on the particle shape symmetry group and discussed potential applications of the CoH in simplifying understanding of the circular motions of microswimmers.

Casimir effect in swimmer suspensions

We show that the Casimir effect can emerge in microswimmer suspensions. In principle, two effects conspire against the development of Casimir effects in swimmer suspensions. First, at low Reynolds number, the force on any closed volume vanishes, but here the relevant effect is the drag by the flow produced by the swimmers, which can be finite. Second, the fluid velocity and the pressure are linear on the swimmer force dipoles, and averaging over the swimmer orientations would lead to a vanishing effect. However, being that the suspension is a discrete system, the noise terms of the coarse-grained equations depend on the density, which itself fluctuates, resulting in effective nonlinear dynamics. Applying the tools developed for other nonequilibrium systems to general coarse-grained equations for swimmer suspensions, the Casimir drag is computed on immersed objects, and it is found to depend on the correlation function between the rescaled density and dipolar density fields. By introducing a model correlation function with medium-range order, explicit expressions are obtained for the Casimir drag on a body. When the correlation length is much larger than the microscopic cutoff, the average drag is independent of the correlation length, with a range that depends only on the size of the immersed bodies.