Source author record

Joost de Graaf

Joost de Graaf appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2026arXiv

The Geometry behind the Glass Transition and Frictional Jamming in Systems of Two-Dimensional Hard Disks

The relation between dynamics and structure in systems of Brownian bidisperse 2D hard disks with arrested dynamics is examined using numerical simulations. Surprisingly, the suspensions show dynamic arrest at an area fraction of ϕ {\approx} 0.777 over a wide range of disk-size ratios. This is in close agreement with the experimental findings of [Nat. Mater. 18, 1118 (2019)] (ϕ {\approx} 0.776) for a quasi-2D colloidal suspension of spheres with large-to-small size ratio of approximately 1.4. Intriguingly, this also matches a jamming transition (ϕ {\approx} 0.773 to 0.777) found experimentally in a 2D bidisperse granular packing of disks for a similar aspect ratio [Powders and Grains (2025)]. Adopting a geometric viewpoint allows for the identification of the floret pentagonal tiling (ϕ {\approx} 0.777343), which is comprised of congruent (elongated) pentagonal tiles. In analogy to foam and tissue models, it is the presence of this congruent reference state that induces a dynamical transition in the disordered fluid of the disks. That is, the change in dynamics observed both in a uniformly compressed granular medium and a colloidal one at finite temperature are implied to be caused by the same topological mechanism. Extending this reasoning to a honeycomb lattice could also explain the experimentally observed onset of caging dynamics (ϕ {\approx} 0.6) in a similar bidisperse, quasi-2D system of colloidal spheres [Nature 587, 225 (2020)]. An outlook is provided on how this concept may be applied to 3D: removing one in four particles from a face-centered-cubic arrangement leads to a square bipyramidal honeycomb} with associated volume fraction η {\approx} 0.55536. The ideas put forward here thus also provide a possible origin for random loose packing of (frictional) hard spheres in 3D and can potentially shed light onto the nature of the glass transition.

preprint2021arXiv

Activity-induced microswimmer interactions and cooperation in one-dimensional environments

Cooperative motion in biological microswimmers is crucial for their survival as it facilitates adhesion to surfaces, formation of hierarchical colonies, efficient motion, and enhanced access to nutrients. Synthetic microswimmers currently lack truly cooperative behavior that originates from activity-induced interactions. Here, we demonstrate that catalytic microswimmers show a variety of cooperative behaviors along one-dimensional paths. We show that their speed increases with the number of swimmers, while the activity induces a preferred distance between swimmers. Using a minimal model, we ascribe this behavior to an effective activity-induced potential that stems from a competition between chemical and hydrodynamic coupling. These interactions further induce active self-assembly into trains as well as compact chains that can elongate, break-up, become immobilized and remobilized. We identify the crucial role that environment morphology and swimmer directionality play on these highly dynamic chain behaviors. These activity-induced interactions open the door towards exploiting cooperation for increasing the efficiency of, as well as provide temporal and spatial control over, microswimmer motion, thereby enabling them to perform intricate tasks inside complex environments.

preprint2019arXiv

Slip length dependent propulsion speed of catalytic colloidal swimmers near walls

Catalytic colloidal swimmers that propel due to self-generated fluid flows exhibit strong affinity for surfaces. We here report experimental measurements of significantly different velocities of such microswimmers in the vicinity of substrates made from different materials. We find that velocities scale with the solution contact angle $θ$ on the substrate, which in turn relates to the associated hydrodynamic substrate slip length, as $V\propto(\cosθ+1)^{-3/2}$. We show that such dependence can be attributed to osmotic coupling between swimmers and substrate. Our work points out that hydrodynamic slip at the wall, though often unconsidered, can significantly impact the self-propulsion of catalytic swimmers.

preprint2016arXiv

Clustering of Self-Propelled Triangles with Surface Roughness

Self-propelled particles can spontaneously form dense phases from a dilute suspension in a process referred to as motility-induced phase separation. The properties of the out-of-equilibrium structures that are formed are governed by the specifics of the particle interactions and the strength of the activity. Thus far, most studies into the formation of these structures have focused on spherical colloids, dumbbells, and rod-like particles endowed with various interaction potentials. Only a few studies have examined the collective behavior of more complex particle shapes. Here, we increase the geometric complexity and use Molecular Dynamics simulations to consider the structures formed by triangular self-propelled particles with surface roughness. These triangles either move towards their apex or towards their base, i.e., they possess a polarity. We find that apex-directed triangles cluster more readily, more stably, and have a smoother cluster interface than their base-directed counterparts. A difference between the two polarities is in line with the results of [H.H. Wensink, et al., Phys. Rev. E 89, 010302(R) (2014)], however, we obtain the reversed result when it comes to clustering, namely that apex-directed particles cluster more readily. We further show that reducing the surface roughness negatively impacts the stability of the base-directed structures, suggesting that their formation is in large part due to surface roughness. Our results lay a solid foundation for future experimental and computational studies into the effect of roughness on the collective dynamics of swimmers.

preprint2016arXiv

Ionic Screening and Dissociation are Crucial for Understanding Chemical Self-Propulsion in Water

Water is a polar solvent and hence supports the bulk dissociation of itself and its solutes into ions, and the re-association of these ions into neutral molecules in a dynamic equilibrium, e.g., ${\rm H_2O_2}\leftrightharpoons {\rm H^+~+~HO_2^-}$. Using continuum theory, we study the influence of these reactions on the self-propulsion of colloids driven by surface chemical reactions (chemical swimmers) in aqueous solution. The association-dissociation reactions are here shown to have a strong influence on the swimmers' behaviour, and must therefore be included in future modelling. In particular, such bulk reactions permit charged swimmers to propel electrophoretically even if all species involved in the surface reactions are neutral. The bulk reactions also significantly modify the predicted speed of chemical swimmers propelled by ionic currents, by up to an order of magnitude. For swimmers whose surface reactions produce both anions and cations (ionic self-diffusiophoresis), the bulk reactions produce an additional reactive screening length, analogous to the Debye length in electrostatics. This in turn leads to an inverse relationship between swimmer radius and swimming speed, which could provide an alternative explanation for recent experimental observations on Pt-polystyrene Janus swimmers [S. Ebbens et al., Phys. Rev. E 85, 020401 (2012)]. We also use our continuum theory to investigate the effect of the Debye screening length itself, going beyond the infinitely thin limit approximation used by previous analytical theories. We identify significant departures from this limiting behavior for micron-sized swimmers under typical experimental conditions, and find that the limiting behavior fails entirely for nanoscale swimmers.

preprint2016arXiv

Lattice-Boltzmann Hydrodynamics of Anisotropic Active Matter

A plethora of active matter models exist that describe the behavior of self-propelled particles (or swimmers), both with and without hydrodynamics. However, there are few studies that consider shape-anisotropic swimmers and include hydrodynamic interactions. Here, we introduce a simple method to simulate self-propelled colloids interacting hydrodynamically in a viscous medium using the lattice-Boltzmann technique. Our model is based on raspberry-type viscous coupling and a force/counter-force formalism which ensures that the system is force free. We consider several anisotropic shapes and characterize their hydrodynamic multipolar flow field. We demonstrate that shape-anisotropy can lead to the presence of a strong quadrupole and octupole moments, in addition to the principle dipole moment. The ability to simulate and characterize these higher-order moments will prove crucial for understanding the behavior of model swimmers in confining geometries.

preprint2016arXiv

Limiting Spurious Flow in Simulations of Electrokinetic Phenomena

Electrokinetic transport phenomena can strongly influence the behaviour of macromolecules and colloidal particles in solution, with applications in, e.g., DNA translocation through nanopores, electro-osmotic flow in nanocapillaries, and electrophoresis of charged macromolecules. Numerical simulations are an important tool to investigate these electrokinetic phenomena, but are often plagued by spurious fluxes and spurious flows that can easily exceed physical fluxes and flows. Here, we present a method that reduces one of these spurious currents, spurious flow, by several orders of magnitude. We demonstrate the effectiveness and generality of our method for both electrokinetic lattice-Boltzmann and finite-element-method based algorithms by simulating a charged sphere in an electrolyte solution, and flow through a nanopore. We also show that previous attempts to suppress these spurious currents introduce other sources of error.

preprint2016arXiv

Nanoparticle Translocation through Conical Nanopores: A Finite Element Study of Electrokinetic Transport

Recent years have seen a surge of interest in nanopores because such structures show a strong potential for characterizing nanoparticles, proteins, DNA, and even single molecules. These systems have been extensively studied in experiment as well as by all-atom and coarse-grained simulations, with a strong focus on DNA translocation. However, the equally interesting problem of particle characterization using nanopores has received far less attention. Here, we theoretically investigate the translocation of a nanoparticle through a conical nanocapillary. We use a model based on numerically solving the coupled system of electrokinetic continuum equations, which we introduce in detail. Based on our findings, we formulate basic guidelines for obtaining the maximum current signal during the translocation event, which should be transferable to other nanopore geometries. In addition, the dependence of the signal strength on particle properties, such as surface charge and size, is evaluated. Finally, we identify conditions under which the translocation is prevented by the formation of a strong electro-osmotic barrier and show that the particle may even become trapped at the pore orifice, without imposing an external hydrostatic pressure difference.

preprint2016arXiv

Selective Trapping of DNA using Glass Microcapillaries

We show experimentally that a cheap glass microcapillary can accumulate λ-phage DNA at its tip and deliver the DNA into the capillary using a combination of electro-osmotic flow, pressure-driven flow, and electrophoresis. We develop an efficient simulation model for this phenomenon based on the electrokinetic equations and the finite-element method. Using our model, we explore the large parameter space of the trapping mechanism by varying the salt concentration, the capillary surface charge, the applied voltage, the pressure difference, and the mobility of the analyte molecules. Our simulation results show that this system can be tuned to capture a wide range of analyte molecules, such as DNA or proteins, based on their electrophoretic mobility. Our method for separation and pre-concentration of analytes has implications for the development of low-cost lab-on-a-chip devices.

preprint2016arXiv

The Efficiency of Self-Phoretic Propulsion Mechanisms with Surface Reaction Heterogeneity

We consider the efficiency of self-phoretic colloidal particles (swimmers) as a function of the heterogeneity in the surface reaction rate. The set of fluid, species, and electrostatic continuity equations is solved analytically using a linearization and numerically using a finite-element method. To compare spherical swimmers of different size and with heterogeneous catalytic conversion rates, a 'swimmer efficiency' functional $η$ is introduced. It is proven, that in order to obtain maximum swimmer efficiency the reactivity has to be localized at the pole(s). Our results also shed light on the sensitivity of the propulsion speed to details of the surface reactivity, a property that is notoriously hard to measure. This insight can be utilized in the design of new self-phoretic swimmers.

preprint2016arXiv

The Raspberry Model for Hydrodynamic Interactions Revisited. II. The Effect of Confinement

The so-called 'raspberry' model refers to the hybrid lattice-Boltzmann (LB) and Langevin molecular dynamics scheme for simulating the dynamics of suspensions of colloidal particles, originally developed by [V. Lobaskin and B. Dünweg, New J. Phys. 6, 54 (2004)], wherein discrete surface points are used to achieve fluid-particle coupling. In this paper, we present a follow up to our study of the effectiveness of the raspberry model in reproducing hydrodynamic interactions in the Stokes regime for spheres arranged in a simple-cubic crystal [L. Fischer, et al., J. Chem. Phys. 143, 084108 (2015)]. Here, we consider the accuracy with which the raspberry model is able to reproduce such interactions for particles confined between two parallel plates. To this end, we compare our LB simulation results to established theoretical expressions and finite-element calculations. We show that there is a discrepancy between the translational and rotational mobility when only surface coupling points are used, as also found in Part I of our joint publication. We demonstrate that adding internal coupling points to the raspberry, can be used to correct said discrepancy in confining geometries as well. Finally, we show that the raspberry model accurately reproduces hydrodynamic interactions between a spherical colloid and planar walls up to roughly one LB lattice spacing.

preprint2016arXiv

Understanding the Onset of Oscillatory Swimming in Microchannels

Self-propelled colloids (swimmers) in confining geometries follow trajectories determined by hydrodynamic interactions with the bounding surfaces. However, typically these interactions are ignored or truncated to lowest order. We demonstrate that higher-order hydrodynamic moments cause rod-like swimmers to follow oscillatory trajectories in quiescent fluid between two parallel plates, using a combination of lattice-Boltzmann simulations and far-field calculations. This behavior occurs even far from the confining walls and does not require lubrication results. We show that a swimmer's hydrodynamic quadrupole moment is crucial to the onset of the oscillatory trajectories. This insight allows us to develop a simple model for the dynamics near the channel center based on these higher hydrodynamic moments, and suggests opportunities for trajectory-based experimental characterization of swimmers' hydrodynamic properties.

preprint2015arXiv

Diffusiophoretic Self-Propulsion for Partially Catalytic Spherical Colloids

Colloidal spheres with a partial platinum surface coating perform auto-phoretic motion when suspended in hydrogen peroxide solution. We present a theoretical analysis of the self-propulsion velocity of these particles using a continuum multi-component, self-diffusiophoretic model. With this model as a basis, we show how the slip-layer approximation can be derived and in which limits it holds. First, we consider the differences between the full multi-component model and the slip-layer approximation. Then the slip model is used to demonstrate and explore the sensitive nature of the particle's velocity on the details of the molecule-surface interaction. We find a strong asymmetry in the dependence of the colloid's velocity as a function of the level of catalytic coating, when there is a different interaction between the solute and solvent molecules and the inert and catalytic part of the colloid, respectively. The direction of motion can even be reversed by varying the level of the catalytic coating. Finally, we investigate the robustness of these results with respect to variations in the reaction rate near the edge between the catalytic and inert parts of the particle. Our results are of significant interest to the interpretation of experimental results on the motion of self-propelled particles.

preprint2015arXiv

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 [V. Lobaskin and B. 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 the 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 [J. de Graaf, et al., J. Chem. Phys. 143, 084107 (2015)], wherein we consider the raspberry model in the confining geometry of two parallel plates.

preprint2014arXiv

Phase Behavior of a Family of Truncated Hard Cubes

In continuation of our work in [A.P. Gantapara et al., Phys. Rev. Lett. 111, 015501 (2013)], we investigate here the thermodynamic phase behavior of a family of truncated hard cubes, for which the shape evolves smoothly from a cube via a cuboctahedron to an octahedron. We used Monte Carlo simulations and free-energy calculations to establish the full phase diagram. This phase diagram exhibits a remarkable richness in crystal and mesophase structures, depending sensitively on the precise particle shape. In addition, we examined in detail the nature of the plastic crystal (rotator) phases that appear for intermediate densities and levels of truncation. Our results allow us to probe the relation between phase behavior and building-block shape and to further the understanding of rotator phases. Furthermore, the phase diagram presented here should prove instrumental for guiding future experimental studies on similarly-haped nanoparticles and the creation of new materials.

preprint2012arXiv

Electrostatic Interactions between Janus Particles

In this paper we study the electrostatic properties of `Janus' spheres with unequal charge densities on both hemispheres. We introduce a method to compare primitive-model Monte Carlo simulations of the ionic double layer with predictions of (mean-field) nonlinear Poisson-Boltzmann theory. We also derive practical DLVO-like expressions that describe the Janus-particle pair interactions by mean-field theory. Using a large set of parameters, we are able to probe the range of validity of the Poisson-Boltzmann approximation, and thus of DLVO-like theories, for such particles. For homogeneously charged spheres this range corresponds well to the range that was predicted by field-theoretical studies of homogeneously charged flat surfaces. Moreover, we find similar ranges for colloids with a Janus-type charge distribution. The techniques and parameters we introduce show promise for future studies of an even wider class of charged-patterned particles.

preprint2012arXiv

Phase diagram of colloidal hard superballs: from cubes via spheres to octahedra

The phase diagram of colloidal hard superballs, of which the shape interpolates between cubes and octahedra via spheres, is determined by free-energy calculations in Monte Carlo simulations. We discover not only a stable face-centered cubic (fcc) plastic crystal phase for near-spherical particles, but also a stable body-centered cubic (bcc) plastic crystal close to the octahedron shape, and in fact even coexistence of these two plastic crystals with a substantial density gap. The plastic fcc and bcc crystals are, however, unstable in the cube and octahedron limit, suggesting that the rounded corners of superballs play an important role in stablizing the rotator phases. In addition, we observe a two-step melting phenomenon for hard octahedra, in which the Minkowski crystal melts into a metastable bcc plastic crystal before melting into the fluid phase.

preprint2011arXiv

Dense Regular Packings of Irregular Non-Convex Particles

We present a new numerical scheme to study systems of non-convex, irregular, and punctured particles in an efficient manner. We employ this method to analyze regular packings of odd-shaped bodies, not only from a nanoparticle but also both from a computational geometry perspective. Besides determining close-packed structures for many shapes, we also discover a new denser configuration for Truncated Tetrahedra. Moreover, we consider recently synthesized nanoparticles and colloids, where we focus on the excluded volume interactions, to show the applicability of our method in the investigation of their crystal structures and phase behavior. Extensions to the presented scheme include the incorporation of soft particle-particle interactions, the study of quasicrystalline systems, and random packings.

preprint2009arXiv

A Triangular Tessellation Scheme for the Adsorption Free Energy at the Liquid-Liquid Interface: Towards Non-Convex Patterned Colloids

We introduce a new numerical technique, namely triangular tessellation, to calculate the free energy associated with the adsorption of a colloidal particle at a flat interface. The theory and numerical scheme presented here are sufficiently general to handle non-convex patchy colloids with arbitrary surface patterns characterized by a wetting angle, e.g., amphiphilicity. We ignore interfacial deformation due to capillary, electrostatic, or gravitational forces, but the method can be extended to take such effects into account. It is verified that the numerical method presented is accurate and sufficiently stable to be applied to more general situations than presented in this paper. The merits of the tessellation method prove to outweigh those of traditionally used semi-analytic approaches, especially when it comes to generality and applicability.

preprint2008arXiv

Curvature dependence of the electrolytic liquid-liquid interfacial tension

The interfacial tension of a liquid droplet surrounded by another liquid in the presence of microscopic ions is studied as a function of the droplet radius. An analytical expression for the interfacial tension is obtained within a linear Poisson-Boltzmann theory and compared with numerical results from non-linear Poisson-Boltzmann theory. The excess liquid-liquid interfacial tension with respect to the pure, salt-free liquid-liquid interfacial tension is found to decompose into a curvature-independent part due to short-ranged interfacial effects and a curvature-dependent electrostatic contribution. Several curvature-dependent regimes of different scaling of the electrostatic excess interfacial tension are identified. Symmetry relations of the interfacial tension upon swapping droplet and bulk liquid are found to hold in the low-curvature limit, which, e.g., lead to a sign change of the excess Tolman length. For some systems a low-curvature expansion up to second order turns out to be applicable if and only if the droplet size exceeds the Debye screening length in the droplet, independent of the Debye length in the bulk.

preprint2008arXiv

Spontaneous Charging and Crystallization of Water Droplets in Oil

We study the spontaneous charging and the crystallization of spherical micron-sized water-droplets dispersed in oil by numerically solving, within a Poisson-Boltzmann theory in the geometry of a spherical cell, for the density profiles of the cations and anions in the system. We take into account screening, ionic Born self-energy differences between oil and water, and partitioning of ions over the two media. We find that the surface charge density of the droplet as induced by the ion partitioning is significantly affected by the droplet curvature and by the finite density of the droplets. We also find that the salt concentration and the dielectric constant regime in which crystallization of the water droplets is predicted is enhanced substantially compared to results based on the planar oil-water interface, thereby improving quantitative agreement with recent experiments.

preprint2008arXiv

Stability of additive-free water-in-oil emulsions

We calculate ion distributions near a planar oil-water interface within non-linear Poisson-Boltzmann theory, taking into account the Born self-energy of the ions in the two media. For unequal self-energies of cations and anions, a spontaneous charge separation is found such that the water and oil phase become oppositely charged, in slabs with a typical thickness of the Debye screening length in the two media. From the analytical solutions, the corresponding interfacial charge density and the contribution to the interfacial tension is derived, together with an estimate for the Yukawa-potential between two spherical water droplets in oil. The parameter regime is explored where the plasma coupling parameter exceeds the crystallization threshold, i.e. where the droplets are expected to form crystalline structures due to a strong Yukawa repulsion, as recently observed experimentally. Extensions of the theory that we discuss briefly include numerical calculations on spherical water droplets in oil, and analytical calculations of the linear PB-equation for a finite oil-water interfacial width.