BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Frank Smallenburg

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

ResearcherUnclaimed source record
cond-mat.softcond-mat.stat-mechphysics.comp-phphysics.flu-dyn

Catalog footprint

What is connected

24works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

24 published item(s)

preprint2026arXiv

Characterizing the cage state of glassy systems and its sensitivity to frozen boundaries

Understanding the role that structure plays in the dynamical arrest observed in glassy systems remains an open challenge. Over the last decade, machine learning (ML) strategies have emerged as an important tool for probing this structure-dynamics relationship, particularly for predicting heterogeneous glassy dynamics from local structure. A recent advancement is the introduction of the cage state, a structural quantity that captures the average positions of particles while rearrangements are forbidden. During the caging regime, linear models trained on the cage state have been shown to outperform more complex ML methods trained on initial configurations only. In this paper, we explore the properties associated with the cage state in more detail to better understand why it serves as such an effective predictor for the dynamics. Specifically, we examine how the cage state in a binary hard-sphere mixture is influenced by both packing fraction and boundary conditions. Our results reveal that, as the system approaches the glassy regime, the cage state becomes increasingly influenced by long-range structural effects. This influence is evident both in its predictive power for particle dynamics and in the internal structure of the cage state, suggesting that the CS might be associated with some form of an amorphous growing structural length scale.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Machine learning short-ranged many-body interactions in colloidal systems using descriptors based on Voronoi cells

Machine learning (ML) strategies are opening the door to faster computer simulations, allowing us to simulate more realistic colloidal systems. Since the interactions in colloidal systems are often highly many-body, stemming from e.g. depletion and steric interactions, one of the challenges for these algorithms is capturing the many-body nature of these interactions. In this paper, we introduce a new ML-based strategy for fitting many-body interactions in colloidal systems where the many-body interaction is highly local. To this end, we develop Voronoi-based descriptors for capturing the local environment and fit the effective potential using a simple neural network. To test this algorithm, we consider a simple two-dimensional model for a colloid-polymer mixture, where the colloid-colloid interactions and colloid-polymer interactions are hard-disk like, while the polymers themselves interact as ideal gas particles. We find that a Voronoi-based description is sufficient to accurately capture the many-body nature of this system. Moreover, we find that the Pearson correlation function alone is insufficient to determine the predictive power of the network emphasizing the importance of additional metrics when assessing the quality of ML-based potentials.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Comparing machine learning techniques for predicting glassy dynamics

In the quest to understand how structure and dynamics are connected in glasses, a number of machine learning based methods have been developed that predict dynamics in supercooled liquids. These methods include both increasingly complex machine learning techniques, and increasingly sophisticated descriptors used to describe the environment around particles. In many cases, both the chosen machine learning technique and choice of structural descriptors are varied simultaneously, making it hard to quantitatively compare the performance of different machine learning approaches. Here, we use three different machine learning algorithms -- linear regression, neural networks, and GNNs -- to predict the dynamic propensity of a glassy binary hard-sphere mixture using as structural input a recursive set of order parameters recently introduced by Boattini et al. [Phys. Rev. Lett. 127, 088007 (2021)]. As we show, when these advanced descriptors are used, all three methods predict the dynamics with nearly equal accuracy. However, the linear regression is orders of magnitude faster to train making it by far the method of choice.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Efficient event-driven simulations of hard spheres

Hard spheres are arguably one of the most fundamental model systems in soft matter physics, and hence a common topic of simulation studies. Event-driven simulation methods provide an efficient method for studying the phase behavior and dynamics of hard spheres under a wide range of different conditions. Here, we examine the impact of several optimization strategies for speeding up event-driven molecular dynamics of hard spheres, and present a light-weight simulation code that outperforms existing simulation codes over a large range of system sizes and packing fractions. The presented differences in simulation speed, typically a factor of five to ten, save significantly on both CPU time and energy consumption, and may be a crucial factor for studying slow processes such as crystal nucleation and glassy dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Self-assembly of dodecagonal and octagonal quasicrystals in hard spheres on a plane

Quasicrystals are fascinating structures, characterized by strong positional order but lacking the periodicity of a crystal. In colloidal systems, quasicrystals are typically predicted for particles with complex or highly specific interactions, which makes experimental realization difficult. Here, we propose an ideal colloidal model system for quasicrystal formation: binary mixtures of hard spheres sedimented onto a flat substrate. Using computer simulations, we explore both the close-packing and spontaneous self-assembly of these systems over a wide range of size ratios and compositions. Surprisingly, we find that this simple, effectively two-dimensional model systems forms not only a variety of crystal phases, but also two quasicrystal phases: one dodecagonal and one octagonal. Intriguingly, the octagonal quasicrystal consists of three different tiles, whose relative concentrations can be continuously tuned via the composition of the binary mixture. Both phases form reliably and rapidly over a significant part of parameter space, making hard spheres on a plane an ideal model system for exploring quasicrystal self-assembly on the colloidal scale.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Correlation between plastic rearrangements and local structure in a cyclically driven glass

The correlation between local structure and the propensity for structural rearrangements has been widely investigated in glass forming liquids and glasses. In this paper we use the excess two-body entropy $S_2$ and tetrahedrality $\n_{tet}$ as the per-particle local structural order parameters to explore such correlations in a three-dimensional model glass subjected to cyclic shear deformation. We first show that for both liquid configurations and the corresponding inherent structures, local ordering increases upon lowering temperature, signaled by a decrease in the two-body entropy and an increase in tetrahedrality. When the inherent structures, or glasses, are periodically sheared athermally, they eventually reach absorbing states for small shear amplitudes, which do not change from one cycle to the next. Large strain amplitudes result in the formation of shear bands, within which particle motion is diffusive. We show that in the steady state, there is a clear difference in the local structural environment of particles that will be part of plastic rearrangements during the next shear cycle and that of particles that are immobile. In particular, particles with higher $S_2$ and lower $n_{tet}$ are more likely to go through rearrangements irrespective of the average energies of the configurations and strain amplitude. For high shear, we find very distinctive local order outside the mobile shear band region, where almost $30\%$ of the particles are involved in icosahedral clusters, contrasting strongly with the fraction of $<5\%$ found inside the shear band.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Climbing two hills is faster than one: collective barrier-crossing by colloids driven through a microchannel

Ohm's law is one of the most central transport rules stating that the total resistance of sequential single resistances is additive. While this rule is most commonly applied to electronic circuits, it also applies to other transport phenomena such as the flow of colloids or nanoparticles through channels containing multiple obstacles, as long as these obstacles are sufficiently far apart. Here we explore the breakdown of Ohm's law for fluids of repulsive colloids driven over two energetic barriers in a microchannel, using real-space microscopy experiments, particle-resolved simulations, and dynamical density functional theory. If the barrier separation is comparable to the particle correlation length, the resistance is highly non-additive, such that the resistance added by the second barrier can be significantly higher or lower than that of the first. Surprisingly, in some cases the second barrier can even add a {\it negative} resistance, such that two identical barriers are easier to cross than a single one. We explain this counterintuitive observation in terms of the structuring of particles trapped between the barriers.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Infinite-pressure phase diagram of binary mixtures of (non)additive hard disks

One versatile route to the creation of two-dimensional crystal structures on the nanometer to micrometer scale is the self-assembly of colloidal particles at an interface. Here, we explore the crystal phases that can be expected from the self-assembly of mixtures of spherical particles of two different sizes, which we map to (additive or non-additive) hard-disk mixtures. We map out the infinite-pressure phase diagram for these mixtures, using Floppy Box Monte Carlo simulations to systematically sample candidate crystal structures with up to 12 disks in the unit cell. As a function of the size ratio and number ratio of the two species of particles, we find a rich variety of periodic crystal structures. Additionally, we identify random tiling regions to predict random tiling quasicrystal stability ranges. Increasing non-additivity both gives rise to additional crystal phases and broadens the stability regime for crystal structures involving a large number of large-small contacts, including random tilings. Our results provide useful guidelines for controlling the self-assembly of colloidal particles at interfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Modeling of many-body interactions between elastic spheres through symmetry functions

Simple models for spherical particles with a soft shell have been shown to self-assemble into numerous crystal phases and even quasicrystals. However, most of these models rely on a simple pairwise interaction, which is usually a valid approximation only in the limit of small deformations, i.e. low densities. In this work, we consider a many-body yet simple model for the evaluation of the elastic energy associated with the deformation of a spherical shell. The resulting energy evaluation, however, is relatively expensive for direct use in simulations. We significantly reduce the associated numerical cost by fitting the potential using a set of symmetry functions. We propose a method for selecting a suitable set of symmetry functions that capture the most relevant features of the particle environment in a systematic manner. The fitted interaction potential is then used in Monte Carlo simulations to draw the phase diagram of the system in two dimensions. The system is found to form both a fluid and a hexagonal crystal phase.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Rotational and translational dynamics in dense fluids of patchy particles

We explore the effect of directionality on rotational and translational relaxation in glassy systems of patchy particles. Using molecular dynamics simulations we analyze the impact of two distinct patch geometries, one that enhance local icosahedral structure and one which does not strongly affect local order. We find that in nearly all investigated cases, rotational relaxation takes place on a much faster time scale than translational relaxation. By comparing to a simplified dynamical Monte Carlo model, we illustrate that rotational diffusion can be qualitatively explained as purely local motion within a fixed environment, which is not coupled strongly to the cage-breaking dynamics required for translational relaxation. Nonetheless, icosahedral patch placement has a profound effect on the local structure of the system, resulting in a dramatic slowdown at low temperatures which is strongest at an intermediate "optimal" patch size.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Tetrahedrality dictates dynamics in hard spheres

Glasses are ubiquitous amorphous solids that remain one of the big mysteries in condensed matter. Despite the vast body of literature on glasses, a unifying approach to link the structure and dynamics of glasses is still missing. A growing set of evidence, however, indicates the microscopic local geometry as a key ingredient. This originated from the seminal work of Frank, who conjectured that glasses may be the result of the local tendency of liquids to form icosahedral structures, which are not capable of globally filling space regularly. Here, we show that, for a fundamental glass model, dynamics can be fully understood by simply counting the number of tetrahedra. Both globally and locally, these local structures directly predict dynamical slowdown. After more than 60 years of Frank's Conjecture, it might not be the icosahedra that matter for glasses, but rather the tetrahedra inside them.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Close packing of rods on spherical surfaces

We study the optimal packing of short, hard spherocylinders confined to lie tangential to a spherical surface, using simulated annealing and molecular dynamics simulations. For clusters of up to twelve particles, we map out the changes in the geometry of the closest-packed configuration as a function of the aspect ratio $L/D$, where $L$ is the cylinder length and $D$ the diameter of the rods. We find a rich variety of cluster structures. For larger clusters, we find that the best-packed configurations up to around 100 particles are highly dependent on the exact number of particles and aspect ratio. For even larger clusters, we find largely disordered clusters for very short rods ($L/D = 0.25$), while slightly longer rods ($L/D = 0.5$ or $1$) prefer a global baseball-like geometry of smectic-like domains, similar to the behavior of large-scale nematic shells. Our results provide predictions for experimentally realizable systems of colloidal rods trapped at the interface of emulsion droplets.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Hard rectangles near curved hard walls: tuning the sign of the Tolman length

Combining analytic calculations, computer simulations, and classical density functional theory we determine the interfacial tension of orientable two-dimensional hard rectangles near a curved hard wall. Both a circular cavity holding the particles and a hard circular obstacle surrounded by particles are considered. We focus on moderate bulk densities (corresponding to area fractions up to 50 percent) where the bulk phase is isotropic and vary the aspect ratio of the rectangles and the curvature of the wall. The Tolman length, which gives the leading curvature correction of the interfacial tension, is found to change sign at a finite density, which can be tuned via the aspect ratio of the rectangles.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Flow of colloidal solids and fluids through constrictions: dynamical density functional theory versus simulation

Using both dynamical density functional theory and particle-resolved Brownian dynamics simulations, we explore the flow of two-dimensional colloidal solids and fluids driven through a linear channel with a geometric constriction. The flow is generated by a constant external force acting on all colloids. The initial configuration is equilibrated in the absence of flow and then the external force is switched on instantaneously. Upon starting the flow, we observe four different scenarios: a complete blockade, a monotonic decay to a constant particle flux (typical for a fluid), a damped oscillatory behaviour in the particle flux, and a long-lived stop-and-go behaviour in the flow (typical for a solid). The dynamical density functional theory describes all four situations but predicts infinitely long undamped oscillations in the flow which are always damped in the simulations. We attribute the mechanisms of the underlying stop-and-go flow to symmetry conditions on the flowing solid. Our predictions are verifiable in real-space experiments on magnetic colloidal monolayers which are driven through structured microchannels and can be exploited to steer the flow throughput in microfluidics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Phase diagram of the ST2 model of water

We evaluate the free energy of the fluid and crystal phases for the ST2 potential [F.H. Stillinger and A. Rahman, J. Chem. Phys. 60, 1545 (1974)] with reaction field corrections for the long-range interactions. We estimate the phase coexistence boundaries in the temperature-pressure plane, as well as the gas-liquid critical point and gas-liquid coexistence conditions. Our study frames the location of the previously identified liquid-liquid critical point relative to the crystalline phase boundaries, and opens the way for exploring crystal nucleation in a model where the metastable liquid-liquid critical point is computationally accessible.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Swim pressure on walls with curves and corners

The concept of swim pressure quantifies the average force exerted by microswimmers on confining walls in non-equilibrium. Here we explore how the swim pressure depends on the wall curvature and on the presence of sharp corners in the wall. For active Brownian particles at high dilution, we present a coherent framework which describes the force and torque on passive particles of arbitrary shape, in the limit of large particles compared to the persistence length of the swimmer trajectories. The resulting forces can be used to derive, for example, the activity-induced depletion interaction between two disks, as well as to optimize the shape of a tracer particle for high swimming velocity. Our predictions are verifiable in experiments on passive obstacles exposed to a bath of bacteria or artificial microswimmers.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Tuning the liquid-liquid transition by modulating the hydrogen bond angular flexibility in a model for water

We propose a simple extension of the well known ST2 model for water [F.H. Stillinger and A. Rahman, J. Chem. Phys. {\bf 60}, 1545 (1974)] that allows for a continuous modification of the hydrogen bond angular flexibility. We show that the bond flexibility affects the relative thermodynamic stability of the liquid and of the hexagonal (or cubic) ice. On increasing flexibility, the liquid-liquid critical point, which in the original ST2 model is located in the no-man's land (i. e. the region where ice is the thermodynamically stable phase) progressively moves to a temperature where the liquid is more stable than ice. Our study definitively proves that the liquid-liquid transition in ST2 is a genuine phenomenon, of high relevance in all tetrahedral network-forming liquids, including water.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Gelling by Heating

We introduce a simple model, a binary mixture of patchy particles, which has been designed to form a gel upon heating. Due to the specific nature of the particle interactions, notably the number and geometry of the patches as well as their interaction energies, the system is a fluid both at high and at low temperatures, whereas at intermediate temperatures the system forms a solid-like disordered open network structure, i.e. a gel. Using molecular dynamics we investigate the static and dynamic properties of this system.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Liquids more stable than crystals

All liquids (except helium due to quantum effects) crystallize at low temperatures, forming ordered structures. The competition between disorder, which stabilizes the liquid phase, and energy, which favors the ordered crystalline structure, inevitably turns in favor of the latter when temperature is lowered and the entropic contribution to the free energy becomes progressively less and less relevant. The "liquid" state survives at low temperatures only as a glass, an out-of-equilibrium arrested state of matter. This textbook description holds inevitably for atomic and molecular systems, where the interaction between particles is set by quantum mechanical laws. The question remains whether the same physics hold for colloidal particles, where inter-particle interactions are usually short-ranged and tunable. Here we show that for patchy colloids with limited valence, conditions can be found for which the disordered liquid phase is stable all the way down to the zero-temperature limit. Our results offer interesting cues for understanding the stability of gels and the glass forming ability of atomic and molecular network glasses.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Phase Diagram of a Reentrant Gel of Patchy Particles

We study the phase diagram of a binary mixture of patchy particles which has been designed to form a reversible gel. For this we perform Monte Carlo and molecular dynamics simulations to investigate the thermodynamics of such a system and compare our numerical results with predictions based on the analytical parameter-free Wertheim theory. We explore a wide range of the temperature-density-composition space that defines the three-dimensional phase diagram of the system. As a result, we delimit the region of thermodynamic stability of the fluid. We find that for a large region of the phase diagram the Wertheim theory is able to give a quantitative description of the system. For higher densities our simulations show that the system is crystallizing into a BCC structure. Finally we study the relaxation dynamics of the system by means of the density and temperature dependences of the diffusion coefficient. We show that there exists a density range where the system passes reversibly from a gel to a fluid upon both heating and cooling, encountering neither demixing nor phase separation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Understanding tetrahedral liquids through patchy colloids

We investigate the structural properties of a simple model for tetrahedral patchy colloids in which the patch width and the patch range can be tuned independently. For wide bond angles, a fully bonded network can be generated by standard Monte Carlo or molecular dynamics simulations of the model, providing a neat method for generating defect-free random tetrahedral networks. This offers the possibility of focusing on the role of the patch angular width on the structure of the fully bonded network. The analysis of the fully bonded configurations as a function of the bonding angle shows how the bonding angle controls the system compressibility, the strength of the pre-peak in the structure factor and ring size distribution. Comparison with models of liquid water and silica allows us to find the best mapping between these continuous potentials and the colloidal one. Building on previous studies focused on the connection between angular range and crystallization, the mapping makes it possible to shed new light on the glass-forming ability of network-forming tetrahedral liquids.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Vacancy-stabilized crystalline order in hard cubes

We examine the effect of vacancies on the phase behavior and structure of systems consisting of hard cubes using event-driven molecular dynamics and Monte Carlo simulations. We find a first-order phase transition between a fluid and a simple cubic crystal phase that is stabilized by a surprisingly large number of vacancies, reaching a net vacancy concentration of ~6.4% near bulk coexistence. Remarkably, we find that vacancies increase the positional order in the system. Finally, we show that the vacancies are delocalized and therefore hard to detect.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Measuring Colloidal Forces from Particle Position Deviations inside an Optical Trap

We measure interaction forces between pairs of charged PMMA colloidal particles suspended in a relatively low-polar medium (5 $\lesssim ε\lesssim$ 8) directly from the deviations of particle positions inside two time-shared optical traps. The particles are confined to optical point traps; one is held in a stationary trap and the other particle is brought closer in small steps while tracking the particle positions using confocal microscopy. From the observed particle positions inside the traps we calculate the interparticle forces using an ensemble-averaged particle displacement-force relationship. The force measurements are confirmed by independent measurements of the different parameters using electrophoresis and a scaling law for the liquid-solid phase transition. When increasing the salt concentration by exposing the sample to UV light, the force measurements agree well with the classical DLVO theory assuming a constant surface potential. On the other hand, when adding tetrabutylammonium chloride (TBAC) to vary the salt concentration, surface charge regulation seems to play an important role.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Phase diagrams of colloidal spheres with a constant zeta-potential

We study suspensions of colloidal spheres with a constant zeta-potential within Poisson-Boltzmann theory, quantifying the discharging of the spheres with increasing colloid density and decreasing salt concentration. We use the calculated renormalized charge of the colloids to determine their pairwise effective screened-Coulomb repulsions. Bulk phase diagrams in the colloid concentration-salt concentration representation follow, for various zeta-potentials, by a mapping onto published fits of phase boundaries of point-Yukawa systems. Although the resulting phase diagrams do feature face-centered cubic (fcc) and body-centered cubic (bcc) phases, they are dominated by the (re-entrant) fluid phase due to the colloidal discharging with increasing colloid concentration and decreasing salt concentration.

0 citations0 reviewsNo code linkedOpen access