Active Brownian particles display self-propelled movement, which can be modelled as arising from a one-body force. Although their interparticle interactions are purely repulsive, for strong self propulsion the swimmers phase separate into dilute and
dense phases. We describe in detail a recent theory (Phys. Rev. E 100, 052604 (2019); Phys. Rev. Lett. 128, 26802 (2019)) for such motility induced phase-separation. Starting from the continuity equation and the force density balance, the description is based on four superadiabatic contributions to the internal force density. Here the superadiabatic forces are due to the flow in the system and they act on top of the adiabatic forces that arise from the equilibrium free energy. Phase coexistence is described by bulk state functions and agrees quantitatively with Brownian dynamics simulation results from the literature. We describe in detail all analytical steps to fully resolve the spatial and orientational dependence of the one-body density and current. The decomposition into angular Fourier series leads to coupling of total density, polarization and all higher modes. We describe the power functional approach, including the kinematic dependence of the superadiabatic force fields and the quiet life effect that pushes particles from fast to slow regions, and hence induces the phase separation.
Detailed control over the motion of colloidal particles is relevant in many applications in colloidal science such as lab-on-a-chip devices. Here, we use an external magnetic field to assemble paramagnetic colloidal spheres into colloidal rods of sev
eral lengths. The rods reside above a square magnetic pattern and are transported via modulation of the direction of the external magnetic field. The rods behave like bipeds walking above the pattern. Depending on their length, the bipeds perform topologically distinct classes of protected walks above the pattern. We demonstrate that it is possible to design parallel polydirectional modulation loops of the external field that command up to six classes of bipeds to walk on distinct predesigned paths. We use such parallel polydirectional loops to induce the collision of reactant bipeds, their polymerization addition reaction to larger bipeds, the separation of product bipeds from the educts, the sorting of different product bipeds, and also the parallel writing of a word consisting of several different letters.
We present a fundamental classification of forces relevant in nonequilibrium structure formation under collective flow in Brownian many-body systems. The internal one-body force field is systematically split into contributions relevant for the spatia
l structure and for the coupled motion. We demonstrate that both contributions can be obtained straightforwardly in computer simulations, and present a power functional theory that describes all types of forces quantitatively. Our conclusions and methods are relevant for flow in inertial systems, such as molecular liquids and granular media.
Single and double paramagnetic colloidal particles are placed above a magnetic square pattern and are driven with an external magnetic field precessing around a high symmetry direction of the pattern. The external magnetic field and that of the patte
rn confine the colloids into lanes parallel to a lattice vector of the pattern. The precession of the external field causes traveling minima of the magnetic potential along the direction of the lanes. At sufficiently high frequencies of modulation only the doublets respond to the external field and move in direction of the traveling minima along the lanes, while the single colloids cannot follow and remain static. We show how the doublets can induce a coordinated motion of the single colloids building colloidal trains made of a chain of several single colloids transported by doublets.
The excluded area between a pair of two-dimensional hard particles with given relative orientation is the region in which one particle cannot be located due to the presence of the other particle. The magnitude of the excluded area as a function of th
e relative particle orientation plays a major role in the determination of the bulk phase behaviour of hard particles. We use principal component analysis to identify the different types of excluded area corresponding to randomly generated two-dimensional hard particles modeled as non-self-intersecting polygons and star lines (line segments radiating from a common origin). Only three principal components are required to have an excellent representation of the value of the excluded area as a function of the relative particle orientation. Independently of the particle shape, the minimum value of the excluded area is always achieved when the particles are antiparallel to each other. The property that affects the value of the excluded area most strongly is the elongation of the particle shape. Principal component analysis identifies four limiting cases of excluded areas with one to four global minima at equispaced relative orientations. We study selected particle shapes using Monte Carlo simulations.
When an external field drives a colloidal system out of equilibrium, the ensuing colloidal response can be very complex and obtaining a detailed physical understanding often requires case-by-case considerations. In order to facilitate systematic anal
ysis, here we present a general iterative scheme for the determination of the unique external force field that yields a prescribed inhomogeneous stationary or time-dependent flow in an overdamped Brownian many-body system. The computer simulation method is based on the exact one-body force balance equation and allows to specifically tailor both gradient and rotational velocity contributions, as well as to freely control the one-body density distribution. Hence compressibility of the flow field can be fully adjusted. The practical convergence to a unique external force field demonstrates the existence of a functional map from both velocity and density to external force field, as predicted by the power functional variational framework. In equilibrium, the method allows to find the conservative force field that generates a prescribed target density profile, and hence implements the Mermin-Evans classical density functional map from density distribution to external potential. The conceptual tools developed here enable one to gain detailed physical insight into complex flow behaviour, as we demonstrate in prototypical situations.
We identify a structural one-body force field that sustains spatial inhomogeneities in nonequilibrium overdamped Brownian many-body systems. The structural force is perpendicular to the local flow direction, it is free of viscous dissipation, it is m
icroscopically resolved in both space and and time, and it can stabilize density gradients. From the time evolution in the exact (Smoluchowski) low-density limit, Brownian dynamics simulations and a novel power functional approximation, we obtain a quantitative understanding of viscous and structural forces, including memory and shear migration.
We study the percolation properties for a system of functionalized colloids on patterned substrates via Monte Carlo simulations. The colloidal particles are modeled as hard disks with three equally-distributed attractive patches on their perimeter. W
e describe the patterns on the substrate as circular potential wells of radius $R_p$ arranged in a regular square or hexagonal lattice. We find a nonmonotonic behavior of the percolation threshold (packing fraction) as a function of $R_p$. For attractive wells, the percolation threshold is higher than the one for clean (non-patterned) substrates if the circular wells are non-overlapping and can only be lower if the wells overlap. For repulsive wells we find the opposite behavior. In addition, at high packing fractions the formation of both structural and bond defects suppress percolation. As a result, the percolation diagram is reentrant with the non-percolated state occurring at very low and intermediate densities.
Topological insulators insulate in the bulk but exhibit robust conducting edge states protected by the topology of the bulk material. Here, we design a colloidal topological insulator and demonstrate experimentally the occurrence of edge states in a
classical particle system. Magnetic colloidal particles travel along the edge of two distinct magnetic lattices. We drive the colloids with a uniform external magnetic field that performs a topologically non-trivial modulation loop. The loop induces closed orbits in the bulk of the magnetic lattices. At the edge, where both lattices merge, the colloids perform skipping orbits trajectories and hence edge-transport. We also observe paramagnetic and diamagnetic colloids moving in opposite directions along the edge between two inverted patterns; the analogue of a quantum spin Hall effect in topological insulators. We present a new, robust, and versatile way of transporting colloidal particles, enabling new pathways towards lab on a chip applications.
Calculating one-body density profiles in equilibrium via particle-based simulation methods involves counting of events of particle occurrences at (histogram-resolved) space points. Here we investigate an alternative method based on a histogram of the
local force density. Via an exact sum rule the density profile is obtained with a simple spatial integration. The method circumvents the inherent ideal gas fluctuations. We have tested the method in Monte Carlo, Brownian Dynamics and Molecular Dynamics simulations. The results carry a statistical uncertainty smaller than that of the standard, counting, method, reducing therefore the computation time.
