Researcher profile

Frank Smallenburg

Frank Smallenburg contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

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

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.

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.

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.

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.

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.

preprint2020arXiv

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

Ohm&#39;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&#39;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.

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.

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.

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 &#34;optimal&#34; patch size.

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&#39;s Conjecture, it might not be the icosahedra that matter for glasses, but rather the tetrahedra inside them.