Source author record

Laura Filion

Laura Filion 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

14works
3topics
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

14 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

Modeling of effective interactions between ligand coated nanoparticles through symmetry functions

Ligand coated nanoparticles are complex objects consisting of a metallic or semiconductor core with organic ligands grafted on their surface. These organic ligands provide stability to a nanoparticle suspension. In solutions, the effective interactions between such nanoparticles are mediated through a complex interplay of interactions between the nanoparticle cores, the surrounding ligands and the solvent molecules. While it is possible to compute these interactions using fully atomistic molecular simulations, such computations are too expensive for studying self-assembly of a large number of nanoparticles. The problem can be made tractable by removing the degrees of freedom associated with the ligand chains and solvent molecules and using the potentials of mean force (PMF) between nanoparticles. In general, the functional dependence of the PMFs on the inter-particle distance is unknown and can be quite complex. In this article, we present a method to model the two-body and three-body potentials of mean force between ligand coated nanoparticles through a linear combination of symmetry functions. The method is quite general and can be extended to model interactions between different types of macromolecules.

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.

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

Predicting the phase behavior of mixtures of active spherical particles

An important question in the field of active matter is whether or not it is possible to predict the phase behavior of these systems. Here, we study the phase coexistence of binary mixtures of torque-free active Brownian particles, for both systems with purely repulsive interactions and systems with attractions. Using Brownian dynamics simulations, we show that phase coexistences can be predicted quantitatively for these systems by measuring the pressure and "reservoir densities". Specifically, in agreement with previous literature, we find that the coexisting phases are in mechanical equilibrium, i.e. the two phases have the same pressure. Importantly, we also demonstrate that the coexisting phases are in chemical equilibrium by bringing each phase into contact with particle reservoirs, and showing that for each species these reservoirs are characterized by the same density for both phases. Using this requirement of mechanical and chemical equilibrium we accurately construct the phase boundaries from properties which can be measured purely from the individual coexisting phases. This result highlights that torque-free active Brownian systems follow simple coexistence rules, thus shedding new light on their thermodynamics.

preprint2019arXiv

Unsupervised learning for local structure detection in colloidal systems

We introduce a simple, fast, and easy to implement unsupervised learning algorithm for detecting different local environments on a single-particle level in colloidal systems. In this algorithm, we use a vector of standard bond-orientational order parameters to describe the local environment of each particle. We then use a neural-network-based autoencoder combined with Gaussian mixture models in order to autonomously group together similar environments. We test the performance of the method on snapshots of a wide variety of colloidal systems obtained via computer simulations, ranging from simple isotropically interacting systems, to binary mixtures, and even anisotropic hard cubes. Additionally, we look at a variety of common self-assembled situations such as fluid-crystal and crystal-crystal coexistences, grain boundaries, and nucleation. In all cases, we are able to identify the relevant local environments to a similar precision as "standard", manually-tuned and system-specific, order parameters. In addition to classifying such environments, we also use the trained autoencoder in order to determine the most relevant bond orientational order parameters in the systems analyzed.

preprint2016arXiv

Removing grain boundaries from three-dimensional colloidal crystals using active dopants

Using computer simulations we explore how grain boundaries can be removed from three-dimensional colloidal crystals by doping with a small fraction of active colloids. We show that for sufficient self-propulsion, the system is driven into a crystal-fluid coexistence. In this phase separated regime, the active dopants become mobile and spontaneously gather at the grain boundaries. The resulting surface melting and recrystallization of domains result in the motion of the grain boundaries over time and lead to the formation of a large single crystal. However, when the self-propulsion is too low to cause a phase separation, we observe no significant enhancement of grain growth.

preprint2016arXiv

Vapour-Liquid Coexistence of an Active Lennard-Jones fluid

We study a three-dimensional system of self-propelled Lennard-Jones particles using Brownian Dynamics simulations. Using recent theoretical results for active matter, we calculate the pressure and report equations of state for the system. Additionally, we chart the vapour-liquid coexistence and show that the coexistence densities can be well described using simple power laws. Lastly, we demonstrate that our out-of-equilibrium system shows deviations from both the law of rectilinear diameters and the law of corresponding states.

preprint2012arXiv

A parameter-free, solid-angle based, nearest-neighbor algorithm

We propose a parameter-free algorithm for the identification of nearest neighbors. The algorithm is very easy to use and has a number of advantages over existing algorithms to identify nearest- neighbors. This solid-angle based nearest-neighbor algorithm (SANN) attributes to each possible neighbor a solid angle and determines the cutoff radius by the requirement that the sum of the solid angles is 4π. The algorithm can be used to analyze 3D images, both from experiments as well as theory, and as the algorithm has a low computational cost, it can also be used "on the fly" in simulations. In this paper, we describe the SANN algorithm, discuss its properties, and compare it to both a fixed-distance cutoff algorithm and to a Voronoi construction by analyzing its behavior in bulk phases of systems of carbon atoms, Lennard-Jones particles and hard spheres as well as in Lennard-Jones systems with liquid-crystal and liquid-vapor interfaces.

preprint2012arXiv

Predicting Patchy Particle Crystals: Variable Box Shape Simulations and Evolutionary Algorithms

We consider several patchy particle models that have been proposed in literature and we investigate their candidate crystal structures in a systematic way. We compare two different algorithms for predicting crystal structures: (i) an approach based on Monte Carlo simulations in the isobaric-isothermal ensemble and (ii) an optimization technique based on ideas of evolutionary algorithms. We show that the two methods are equally successful and provide consistent results on crystalline phases of patchy particle systems.

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.

preprint2010arXiv

Crystal nucleation of hard spheres using molecular dynamics, umbrella sampling and forward flux sampling: A comparison of simulation techniques

Over the last number of years several simulation methods have been introduced to study rare events such as nucleation. In this paper we examine the crystal nucleation rate of hard spheres using three such numerical techniques: molecular dynamics, forward flux sampling and a Bennett-Chandler type theory where the nucleation barrier is determined using umbrella sampling simulations. The resulting nucleation rates are compared with the experimental rates of Harland and Van Megen [J. L. Harland and W. van Megen, Phys. Rev. E 55, 3054 (1997)], Sinn et al. [C. Sinn et al., Prog. Colloid Polym. Sci. 118, 266 (2001)] and Schatzel and Ackerson [K. Schatzel and B.J. Ackerson, Phys. Rev. E, 48, 3766 (1993)] and the predicted rates for monodisperse and 5% polydisperse hard spheres of Auer and Frenkel [S. Auer and D. Frenkel, Nature 409, 1020 (2001)]. When the rates are examined in long-time diffusion units, we find agreement between all the theoretically predicted nucleation rates, however, the experimental results display a markedly different behaviour for low supersaturation. Additionally, we examined the pre-critical nuclei arising in the molecular dynamics, forward flux sampling, and umbrella sampling simulations. The structure of the nuclei appear independent of the simulation method, and in all cases, the nuclei contain on average significantly more face-centered-cubic ordered particles than hexagonal-close-packed ordered particles.