Source author record

Robert L. Jack

Robert L. Jack 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

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

53 published item(s)

preprint2025arXiv

Non-equilibrium pathways between cluster morphologies in active phase separation: necking, rupture and cavitation

We investigate the dynamical pathways of a geometric phase transition in a two-dimensional active lattice gas undergoing motility-induced phase separation. The transition is between metastable morphologies of the liquid cluster: a system-spanning "slab" and a compact "droplet". We generate trajectories of this transition in both directions using forward flux sampling. We find that the droplet-to-slab transition always follows a similar mechanism to its equilibrium counterpart, but the reverse (slab-to-droplet) transition depends on rare non-equilibrium fluctuations. At low Peclet numbers the equilibrium and non-equilibrium pathways compete, while at high Peclet numbers the equilibrium pathway is entirely suppressed, and the only allowed mechanism involves a large vapour bubble. We discuss the implications of these findings for active matter systems more generally.

preprint2022arXiv

Convergence of random-weight sequential Monte Carlo methods

We investigate the properties of a sequential Monte Carlo method where the particle weight that appears in the algorithm is estimated by a positive, unbiased estimator. We present broadly-applicable convergence results, including a central limit theorem, and we discuss their relevance for applications in statistical physics. Using these results, we show that the resampling step reduces the impact of the randomness of the weights on the asymptotic variance of the estimator. In addition, we explore the limits of convergence of the sequential Monte Carlo method, with a focus on almost sure convergence. We construct an example algorithm where we can prove convergence in probability, but which does not converge almost surely, even in the non-random-weight case.

preprint2022arXiv

Cycle counts and affinities in stochastic models of non-equilibrium systems

For non-equilibrium systems described by finite Markov processes, we consider the number of times that a system traverses a cyclic sequence of states (a cycle). The joint distribution of the number of forward and backward instances of any given cycle is described by universal formulae which depend on the cycle affinity, but are otherwise independent of system details. We discuss the similarities and differences of this result to fluctuation theorems, and generalize the result to families of cycles, relevant under coarse-graining. Finally, we describe the application of large deviation theory to this cycle counting problem.

preprint2022arXiv

Direct Imaging of Contacts and Forces in Colloidal Gels

Colloidal dispersions are prized as model systems to understand basic properties of materials, and are central to a wide range of industries from cosmetics to foods to agrichemicals. Among the key developments in using colloids to address challenges in condensed matter is to resolve the particle coordinates in 3D, allowing a level of analysis usually only possible in computer simulation. However in amorphous materials, relating mechanical properties to microscopic structure remains problematic. This makes it rather hard to understand, for example, mechanical failure. Here we address this challenge by studying the contacts and the forces between particles, as well as their positions. To do so, we use a colloidal model system (an emulsion) in which the interparticle forces and local stress can be linked to the microscopic structure. We demonstrate the potential of our method to reveal insights into the failure mechanisms of soft amorphous solids by determining local stress in a colloidal gel. In particular, we identify `force chains' of load-bearing droplets, and local stress anisotropy, and investigate their connection with locally rigid packings of the droplets.

preprint2022arXiv

Disordered collective motion in dense assemblies of persistent particles

We explore the emergence of nonequilibrium collective motion in disordered non-thermal active matter when persistent motion and crowding effects compete, using simulations of a two-dimensional model of size polydisperse self-propelled particles. In stark contrast with monodisperse systems, we find that polydispersity stabilizes a homogeneous active liquid at arbitrary large persistence times, characterized by remarkable velocity correlations and irregular turbulent flows. For all persistence values, the active fluid undergoes a nonequilibrium glass transition at large density. This is accompanied by collective motion, whose nature evolves from near-equilibrium spatially heterogeneous dynamics at small persistence, to a qualitatively different intermittent dynamics when persistence is large. This latter regime involves a complex time evolution of the correlated displacement field

preprint2022arXiv

Irreversibility and biased ensembles in active matter: Insights from stochastic thermodynamics

Active systems evade the rules of equilibrium thermodynamics by constantly dissipating energy at the level of their microscopic components. This energy flux stems from the conversion of a fuel, present in the environment, into sustained individual motion. It can lead to collective effects without any equilibrium equivalent, such as a phase separation for purely repulsive particles, or a collective motion (flocking) for aligning particles. Some of these effects can be rationalized by using equilibrium tools to recapitulate nonequilibrium transitions. An important challenge is then to delineate systematically to which extent the character of these active transitions is genuinely distinct from equilibrium analogs. We review recent works that use stochastic thermodynamics tools to identify, for active systems, a measure of irreversibility comprising a coarse-grained or informatic entropy production. We describe how this relates to the underlying energy dissipation or thermodynamic entropy production, and how it is influenced by collective behavior. Then, we review the possibility to construct thermodynamic ensembles out-of-equilibrium, where trajectories are biased towards atypical values of nonequilibrium observables. We show that this is a generic route to discovering unexpected phase transitions in active matter systems, which can also inform their design.

preprint2022arXiv

Mechanical analysis of a dynamical phase transition for particles in a channel

We analyse biased ensembles of trajectories for a two-dimensional system of particles, evolving by Langevin dynamics in a channel geometry. This bias controls the degree of particle clustering. On biasing to large clustering, we observe a dynamical phase transition where the particles break symmetry and accumulate at one of the walls. We analyse the mechanical properties of this symmetry-broken state using the Irving-Kirkwood stress tensor. The biased ensemble is characterised by body forces which originate in random thermal noises, but have finite averages in the presence of the bias. We discuss the connection of these forces to Doob's transform and optimal control theory.

preprint2022arXiv

Perpendicular and Parallel Phase Separation in Two Species Driven Diffusive Lattice Gases

We study three different lattice models in which two species of diffusing particles are driven in opposite directions by an electric field. We focus on dynamical phase transitions that involve phase separation into domains that may be parallel or perpendicular to a driving field. In all cases, the perpendicular state appears for weak driving, consistent with previous work. For strong driving, we introduce two models that support the parallel state. In one model, this state occurs because of the inclusion of dynamical rules that enhance lateral diffusion during collisions; in the other, it is a result of a nearest-neighbour attractive/repulsive interaction between particles of the same/opposite species. We discuss the connections between these results and the behaviour found in off-lattice systems, including laning and freezing by heating.

preprint2022arXiv

Rectification in a mixture of active and passive particles subject to a ratchet potential

We study by simulation a mixture of active (run-and-tumble) and passive (Brownian) particles with repulsive exclusion interactions in one dimension, subject to a ratchet (smoothed sawtooth) potential. Such a potential is known to rectify active particles at one-body level, creating a net current in the `easy direction'. This is the direction in which one encounters the lower maximum force en route to the top of a potential barrier. The exclusion constraint results in single-file motion, so the mean velocities of active and passive particles are identical; we study the effects of activity level, Brownian diffusivity, particle size, initial sequence of active and passive particles, and active/passive concentration ratio on this mean velocity (i.e., the current per particle). We show that in some parameter regimes the sign of the current is reversed. This happens when the passive particles are at high temperature and so would cross barriers relatively easily, and without rectification, except that they collide with `cold' active ones, which would otherwise be localized near the potential minima. In this case, the reversed current arises because hot passive particles push cold active ones preferentially in the direction with the lower spatial separation between the bottom and top of the barrier. A qualitatively similar mechanism operates in a mixture containing passive particles of two very different temperatures, although there is no quantitative mapping between that case and the systems studied here.

preprint2022arXiv

Tracer dynamics in one dimensional gases of active or passive particles

We consider one-dimensional systems comprising either active run-and-tumble particles (RTPs) or passive Brownian random walkers. These particles are either noninteracting or have hardcore exclusions. We study the dynamics of a single tracer particle embedded in such a system - this tracer may be either active or passive, with hardcore exclusion from environmental particles. In an active hardcore environment, both active and passive tracers show long-time subdiffusion: displacements scale as $t^{1/4}$ with a density-dependent prefactor that is independent of tracer type, and differs from the corresponding result for passive-in-passive subdiffusion. In an environment of noninteracting active particles, the passive-in-passive results are recovered at low densities for both active and passive tracers, but transient caging effects slow the tracer motion at higher densities, delaying the onset of any $t^{1/4}$ regime. For an active tracer in a passive environment, we find more complex outcomes, which depend on details of the dynamical discretization scheme. We interpret these results by studying the density distribution of environmental particles around the tracer. In particular, sticking of environment particles to the tracer cause it to move more slowly in noninteracting than in interacting active environments, while the anomalous behaviour of the active-in-passive cases stems from a `snowplough' effect whereby a large pile of diffusive environmental particles accumulates in front of a RTP tracer during a ballistic run.

preprint2020arXiv

Assessing the structural heterogeneity of supercooled liquids through community inference

We present an information-theoretic approach inspired by distributional clustering to assess the structural heterogeneity of particulate systems. Our method identifies communities of particles that share a similar local structure by harvesting the information hidden in the spatial variation of two- or three-body static correlations. This corresponds to an unsupervised machine learning approach that infers communities solely from the particle positions and their species. We apply this method to three models of supercooled liquids and find that it detects subtle forms of local order, as demonstrated by a comparison with the statistics of Voronoi cells. Finally, we analyze the time-dependent correlation between structural communities and particle mobility and show that our method captures relevant information about glassy dynamics.

preprint2020arXiv

Collective motion in large deviations of active particles

We analyse collective motion that occurs during rare (large deviation) events in systems of active particles, both numerically and analytically. We discuss the associated dynamical phase transition to collective motion, which occurs when the active work is biased towards larger values, and is associated with alignment of particles' orientations. A finite biasing field is needed to induce spontaneous symmetry breaking, even in large systems. Particle alignment is computed exactly for a system of two particles. For many-particle systems, we analyse the symmetry breaking by an optimal-control representation of the biased dynamics, and we propose a fluctuating hydrodynamic theory that captures the emergence of polar order in the biased state.

preprint2020arXiv

Ergodicity and large deviations in physical systems with stochastic dynamics

In ergodic physical systems, time-averaged quantities converge (for large times) to their ensemble-averaged values. Large deviation theory describes rare events where these time averages differ significantly from the corresponding ensemble averages. It allows estimation of the probabilities of these events, and their mechanisms. This theory has been applied to a range of physical systems, where it has yielded new insights into entropy production, current fluctuations, metastability, transport processes, and glassy behaviour. We review some of these developments, identifying general principles. We discuss a selection of dynamical phase transitions, and we highlight some connections between large-deviation theory and optimal control theory.

preprint2020arXiv

Giant leaps and long excursions: fluctuation mechanisms in systems with long-range memory

We analyse large deviations of time-averaged quantities in stochastic processes with long-range memory, where the dynamics at time t depends itself on the value q_t of the time-averaged quantity. First we consider the elephant random walk and a Gaussian variant of this model, identifying two mechanisms for unusual fluctuation behaviour, which differ from the Markovian case. In particular, the memory can lead to large deviation principles with reduced speeds, and to non-analytic rate functions. We then explain how the mechanisms operating in these two models are generic for memory-dependent dynamics and show other examples including a non-Markovian symmetric exclusion process.

preprint2020arXiv

Inference, prediction and optimization of non-pharmaceutical interventions using compartment models: the PyRoss library

PyRoss is an open-source Python library that offers an integrated platform for inference, prediction and optimisation of NPIs in age- and contact-structured epidemiological compartment models. This report outlines the rationale and functionality of the PyRoss library, with various illustrations and examples focusing on well-mixed, age-structured populations. The PyRoss library supports arbitrary structured models formulated stochastically (as master equations) or deterministically (as ODEs) and allows mid-run transitioning from one to the other. By supporting additional compartmental subdivision ad libitum, PyRoss can emulate time-since-infection models and allows medical stages such as hospitalization or quarantine to be modelled and forecast. The PyRoss library enables fitting to epidemiological data, as available, using Bayesian parameter inference, so that competing models can be weighed by their evidence. PyRoss allows fully Bayesian forecasts of the impact of idealized NPIs by convolving uncertainties arising from epidemiological data, model choice, parameters, and intrinsic stochasticity. Algorithms to optimize time-dependent NPI scenarios against user-defined cost functions are included. PyRoss's current age-structured compartment framework for well-mixed populations will in future reports be extended to include compartments structured by location, occupation, use of travel networks and other attributes relevant to assessing disease spread and the impact of NPIs. We argue that such compartment models, by allowing social data of arbitrary granularity to be combined with Bayesian parameter estimation for poorly-known disease variables, could enable more powerful and robust prediction than other approaches to detailed epidemic modelling. We invite others to use the PyRoss library for research to address today's COVID-19 crisis, and to plan for future pandemics.

preprint2019arXiv

Correction of coarse-graining errors by a two-level method: application to the Asakura-Oosawa model

We present a method that exploits self-consistent simulation of coarse-grained and fine-grained models, in order to analyse properties of physical systems. The method uses the coarse-grained model to obtain a first estimate of the quantity of interest, before computing a correction by analysing properties of the fine system. We illustrate the method by applying it to the Asakura-Oosawa (AO) model of colloid-polymer mixtures. We show that the liquid-vapour critical point in that system is affected by three-body interactions which are neglected in the corresponding coarse-grained model. We analyse the size of this effect and the nature of the three-body interactions. We also analyse the accuracy of the method, as a function of the associated computational effort.

preprint2019arXiv

Dynamical phase coexistence in the Fredrickson-Andersen model

We analyse a first-order dynamical phase transition that takes place in the Fredrickson--Andersen (FA) model. We construct a two-dimensional spin system whose thermodynamic properties reproduce the dynamical large deviations of the FA model and we analyse this system numerically, comparing our results with finite-size scaling theory. This allows us to rationalise recent results for the FA model, including the exponential divergence of its susceptibility at phase coexistence. We also discuss a simple interfacial model that reproduces quantitatively the behaviour of the FA model at coexistence.

preprint2019arXiv

Large deviations and optimal control forces for hard particles in one dimension

We analyse large deviations of the dynamical activity in one-dimensional systems of diffusing hard particles. Using an optimal-control representation of the large-deviation problem, we analyse effective interaction forces which can be added to the system, to aid sampling of biased ensembles of trajectories. We find several distinct regimes, as a function of the activity and the system size: we present approximate analytical calculations that characterise the effective interactions in several of these regimes. For high activity the system is hyperuniform and the interactions are long-ranged and repulsive. For low activity, there is a near-equilibrium regime described by macroscopic fluctuation theory, characterised by long-ranged attractive forces. There is also a far-from-equilibrium regime in which one of the interparticle gaps becomes macroscopic and the interactions depend strongly on the size of this gap. We discuss the extent to which transition path sampling of these ensembles is improved by adding suitable control forces.

preprint2016arXiv

Absence of dissipation in trajectory ensembles biased by currents

We consider biased ensembles of trajectories associated with large deviations of currents in equilibrium systems. The biased ensembles are characterised by non-zero currents and lack the time-reversal symmetry of the equilibrium state, but we show that they retain a generalised time-reversal symmetry, involving a spatial transformation that inverts the current. This means that these ensembles lack dissipation. Hence, they differ significantly from non-equilibrium steady states where currents are induced by external forces. One consequence of this result is that maximum entropy assumptions (MaxEnt/MaxCal), widely used for modelling thermal systems away from equilibrium, have quite unexpected implications, including apparent superfluid behaviour in a classical model of shear flow.

preprint2016arXiv

Coarse-grained depletion potentials for anisotropic colloids: application to lock-and-key systems

When a colloid is mixed with a depletant such as a non-adsorbing polymer, one observes attractive effective interactions between the colloidal particles. If these particles are anisotropic, analysis of these effective interactions is challenging in general. We present a method for inference of approximate (coarse-grained) effective interaction potentials between such anisotropic particles. Using the example of indented (lock-and-key) colloids, we show how numerical solutions can be used to integrate out the (hard sphere) depletant, leading to a depletion potential that accurately characterises the effective interactions. The accuracy of the method is based on matching of contributions to the second virial coefficient of the colloids. The simplest version of our method yields a piecewise-constant effective potential; we also show how this scheme can be generalised to other functional forms, where appropriate.

preprint2016arXiv

Ideal bulk pressure of active Brownian particles

The extent to which active matter might be described by effective equilibrium concepts like temperature and pressure is currently being discussed intensely. Here we study the simplest model, an ideal gas of non-interacting active Brownian particles. While the mechanical pressure exerted onto confining walls has been linked to correlations between particles' positions and their orientations, we show that these correlations are entirely controlled by boundary effects. We also consider a definition of local pressure, which describes interparticle forces in terms of momentum exchange between different regions of the system. We present three pieces of analytical evidence which indicate that such a local pressure exists, and we show that its bulk value differs from the mechanical pressure exerted on the walls of the system. We attribute this difference to the fact that the local pressure in the bulk does not depend on boundary effects, contrary to the mechanical pressure. We carefully examine these boundary effects using a channel geometry, and we show a virial formula for the pressure correctly predicts the mechanical pressure even in finite channels. However, this result no longer holds in more complex geometries, as exemplified for a channel that includes circular obstacles.

preprint2016arXiv

Population dynamics method with a multi-canonical feedback control

We discuss the GiardinĂ -Kurchan-Peliti population dynamics method for evaluating large deviations of time averaged quantities in Markov processes [Phys. Rev. Lett. \textbf{96}, 120603 (2006)]. This method exhibits systematic errors which can be large in some circumstances, particularly for systems with weak noise, with many degrees of freedom, or close to dynamical phase transitions. We show how these errors can be mitigated by introducing control forces within the algorithm. These forces are determined by an iteration-and-feedback scheme, inspired by multicanonical methods in equilibrium sampling. We demonstrate substantially improved results in a simple model and we discuss potential applications to more complex systems.

preprint2016arXiv

Structure of inactive states of a binary Lennard-Jones mixture

We study the structure of inactive states in a prototypical model glass, the Kob-Andersen binary Lennard-Jones mixture. These inactive states are obtained by transition path sampling and are at dynamical phase coexistence with an active equilibrium state. Configurations in the inactive states are kinetically stable and are located in deeper basins of the energy landscape than their active counterparts. By analyzing trajectory-to-trajectory fluctuations within the inactive state, we assess correlations between kinetic stability, energy and other structural properties. We show that measures of local order associated to stable local packings and bond-orientational order are weakly correlated with energy and kinetic stability. We discuss what kinds of structural measurement might capture the relevant dynamical features of the inactive state.

preprint2016arXiv

The melting of stable glasses is governed by nucleation-and-growth dynamics

We discuss the microscopic mechanisms by which low-temperature amorphous states, such as ultrastable glasses, transform into equilibrium fluids, after a sudden temperature increase. Experiments suggest that this process is similar to the melting of crystals, thus differing from the behaviour found in ordinary glasses. We rationalize these observations using the physical idea that the transformation process takes place very close to a `hidden' equilibrium first-order phase transition, which is observed in systems of coupled replicas. We illustrate our views using simulation results for a simple two-dimensional plaquette spin model, which is known to exhibit a range of glassy behaviour. Our results suggest that nucleation-and-growth dynamics, as found near ordinary first-order transitions, is also the correct theoretical framework to analyse the melting of ultrastable glasses. Our approach provides a unified understanding of multiple experimental observations, such as propagating melting fronts, large kinetic stability ratios, and `giant' dynamic lengthscales.

preprint2015arXiv

Dynamical phase transitions in one-dimensional hard-particle systems

We analyse a one-dimensional model of hard particles, within ensembles of trajectories that are conditioned (or biased) to atypical values of the time-averaged dynamical activity. We analyse two phenomena that are associated with these large deviations of the activity: phase separation (at low activity) and the formation of hyperuniform states (at high activity). We consider a version of the model which operates at constant volume, and a version at constant pressure. In these non-equilibrium systems, differences arise between the two ensembles, because of the extra freedom available to the constant-pressure system, which can change its total density. We discuss the relationships between different ensembles, mechanical equilibrium, and the probability cost of rare density fluctuations.

preprint2015arXiv

Effective interactions and large deviations in stochastic processes

We discuss the relationships between large deviations in stochastic systems, and "effective interactions" that induce particular rare events. We focus on the nature of these effective interactions in physical systems with many interacting degrees of freedom, which we illustrate by reviewing several recent studies. We describe the connections between effective interactions, large deviations at "level 2.5", and the theory of optimal control. Finally, we discuss possible physical applications of variational results associated with those theories.

preprint2015arXiv

Evidence for a disordered critical point in a glass-forming liquid

Using computer simulations of an atomistic glass-forming liquid, we investigate the fluctuations of the overlap between a fluid configuration and a quenched reference system. We find that large fluctuations of the overlap develop as temperature decreases, consistent with the existence of the random critical point that is predicted by effective field theories. We discuss the scaling of fluctuations near the presumed critical point, comparing the observed behaviour with that of the random-field Ising model. We argue that this critical point directly reveals the existence of an interfacial tension between amorphous metastable states, a quantity relevant both for equilibrium relaxation and for nonequilibrium melting of stable glass configurations.

preprint2015arXiv

Overlap and activity glass transitions in plaquette spin models with hierarchical dynamics

We consider thermodynamic and dynamic phase transitions in plaquette spin models of glasses. The thermodynamic transitions involve coupled (annealed) replicas of the model. We map these coupled-replica systems to a single replica in a magnetic field, which allows us to analyse the resulting phase transitions in detail. For the triangular plaquette model (TPM), we find for the coupled-replica system a phase transition between high- and low-overlap phases, occuring at a coupling eps*(T), which vanishes in the low-temperature limit. Using computational path sampling techniques, we show that a single TPM also displays space-time transitions between active and inactive dynamical phases. These first-order dynamical transitions occur at a critical counting field s_c(T)>=0 that appears to vanish at zero temperature, in a manner reminiscent of the thermodynamic overlap transition. In order to extend the ideas to three dimensions we introduce the square pyramid model which also displays both overlap and activity transitions. We discuss a possible common origin of these various phase transitions, based on long-lived (metastable) glassy states.

preprint2015arXiv

Phase transition for quenched coupled replicas in a plaquette spin model of glasses

We study a three-dimensional plaquette spin model whose low temperature dynamics is glassy, due to localised defects and effective kinetic constraints. While the thermodynamics of this system is smooth at all temperatures, we show that coupling it to a second system with a fixed (quenched) configuration can lead to a phase transition, at finite coupling. The order parameter is the overlap between the copies, and the transition is between phases of low and high overlap. We find critical points whose properties are consistent with random-field Ising universality. We analyse the interfacial free energy cost between the high- and low-overlap states that coexist at (and below) the critical point, and we use this cost as the basis for a finite-size scaling analysis. We discuss these results in the context of mean-field and dynamical facilitation theories of the glass transition.

preprint2015arXiv

Unusual liquid phases for indented colloids with depletion interactions

We study indented spherical colloids, interacting via depletion forces. These systems exhibit liquid-vapor phase transitions whose properties are determined by a combination of strong "lock-and-key" bonds and weaker non-specific interactions. As the propensity for lock-and-key binding increases, the critical point moves to significantly lower density, and the coexisting phases change their structure. In particular, the liquid phase is "porous", exhibiting large percolating voids. The properties of this system depend strongly on the topological structure of an underlying bond network: we comment on the implications of this fact for the assembly of equilibrium states with controlled porous structures.

preprint2014arXiv

Geometrical interpretation of fluctuating hydrodynamics in diffusive systems

We discuss geometric formulations of hydrodynamic limits in diffusive systems. Specifically, we describe a geometrical construction in the space of density profiles --- the Wasserstein geometry --- which allows the deterministic hydrodynamic evolution of the systems to be related to steepest descent of the free energy, and show how this formulation can be related to most probable paths of mesoscopic dissipative systems. The geometric viewpoint is also linked to fluctuating hydrodynamics of these systems via a saddle point argument.

preprint2014arXiv

Hyperuniformity and phase separation in biased ensembles of trajectories for diffusive systems

We analyse biased ensembles of trajectories for diffusive systems. In trajectories biased either by the total activity or the total current, we use fluctuating hydrodynamics to show that these systems exhibit phase transtions into `hyperuniform' states, where large-wavelength density fluctuations are strongly suppressed. We illustrate this behaviour numerically for a system of hard particles in one dimension and we discuss how it appears in simple exclusion processes. We argue that these diffusive systems generically respond very strongly to any non-zero bias, so that homogeneous states with "normal" fluctuations (finite compressibility) exist only when the bias is very weak.

preprint2014arXiv

Information-theoretic measurements of coupling between structure and dynamics in glass-formers

We analyse the connections between structure and dynamics in two model glass-formers, using the mutual information between an initial configuration and the ensuing dynamics to compare the predictive value of different structural observables. We consider the predictive power of normal modes, locally favoured structures, and coarse-grained measurements of local energy and density. The mutual information allows the influence of the liquid structure on the dynamics to be analysed quantitatively as a function of time, showing that normal modes give the most useful predictions on short time scales while local energy and density are most strongly predictive at long times.

preprint2014arXiv

Self-assembly and crystallisation of indented colloids at a planar wall

We report experimental and simulation studies of the structure of a monolayer of indented ("lock and key") colloids, on a planar surface. On adding a non-absorbing polymer with prescribed radius and volume fraction, depletion interactions are induced between the colloids, with controlled range and strength. For spherical particles, this leads to crystallisation, but the indented colloids crystallise less easily than spheres, in both simulation and experiment. Nevertheless, simulations show that indented colloids do form plastic (rotator) crystals. We discuss the conditions under which this occurs, and the possibilities of lower-symmetry crystal states. We also comment on the kinetic accessibility of these states.

preprint2014arXiv

Self-assembly of colloidal polymers via depletion-mediated lock and key binding

We study the depletion-induced self-assembly of indented colloids. Using state-of-the-art Monte Carlo simulation techniques that treat the depletant particles explicitly, we demonstrate that colloids assemble by a lock-and-key mechanism, leading to colloidal polymerization. The morphology of the chains that are formed depends sensitively on the size of the colloidal indentation, with smaller values additionally permitting chain branching. In contrast to the case of spheres with attractive patches, Wertheim's thermodynamic perturbation theory fails to provide a fully quantitative description of the polymerization transition. We trace this failure to a neglect of packing effects and we introduce a modified theory that accounts better for the shape of the colloids, yielding improved agreement with simulation.

preprint2014arXiv

The Statistical Mechanics of Dynamic Pathways to Self-assembly

We describe some of the important physical characteristics of the `pathways', i.e. dynamical processes, by which molecular, nanoscale and micron-scale self-assembly occurs. We highlight the fact that there exist features of self-assembly pathways that are common to a wide range of physical systems, even though those systems may be different in respect of their microscopic details. We summarize some existing theoretical descriptions of self-assembly pathways, and highlight areas -- notably, the description of self-assembly pathways that occur `far' from equilibrium -- that are likely to become increasingly important.

preprint2013arXiv

Controlling crystal self-assembly using a real-time feedback scheme

We simulate crystallisation of hard spheres with short-ranged attractive potentials, as a model self-assembling system. We show how measurements of correlation and response functions during assembly can be used to tune the interaction parameters as assembly proceeds, in order to obtain high-quality crystals. The method we use is independent of details of the interaction potential and of the structure of the final crystal - we propose that it could be applied to a wide range of self-assembling systems.

preprint2013arXiv

Counting metastable states in a kinetically constrained model using a patch repetition analysis

We analyse metastable states in the East model, using a recently-proposed patch-repetition analysis based on time-averaged density profiles. The results reveal a hierarchy of states of varying lifetimes, consistent with previous studies in which the metastable states were identified and used to explain the glassy dynamics of the model. We establish a mapping between these states and configurations of systems of hard rods, which allows us to analyse both typical and atypical metastable states. We discuss connections between the complexity of metastable states and large-deviation functions of dynamical quantities, both in the context of the East model and more generally in glassy systems.

preprint2013arXiv

Dynamical correlations in a glass-former with randomly pinned particles

The effects of randomly pinning particles in a model glass-forming fluid are studied, with a focus on the dynamically heterogeneous relaxation in the presence of pinning. We show how four-point dynamical correlations can be analysed in real space, allowing direct extraction of a length scale that characterises dynamical heterogeneity. In the presence of pinning, the relaxation time of the glassy system increases by up to two decades, but there is almost no increase in either the four-point correlation length or the strength of the four-point correlations. We discuss the implications of these results for theories of the glass transition.

preprint2013arXiv

Dynamical phase transitions in supercooled liquids: interpreting measurements of dynamical activity

We study dynamical phase transitions in a model supercooled liquid. These transitions occur in ensembles of trajectories that are biased towards low (or high) dynamical activity. We compare two different measures of activity that were introduced in recent papers and we find that they are anti-correlated with each other. To interpret this result, we show that the two measures couple to motion on different length and time scales. We find that inactive states with very slow structural relaxation nevertheless have increased molecular motion on short scales. We discuss these results in terms of the potential energy landscape of the system and in terms of the liquid structure in active/inactive states.

preprint2013arXiv

Investigating amorphous order in stable glasses by random pinning

We use a random pinning procedure to investigate stable glassy states associated with large deviations of the activity in a model glass-former. We pin particles both from active (equilibrium) configurations and from stable (inactive) glassy states. By comparing the distribution of the overlap between states that share pinned particles, we infer that the inactive states are characterised by a structural length scale that is comparable with the system size. This is a manifestation of amorphous order in these glassy states, which helps to explain their stability.

preprint2013arXiv

Large deviations of the dynamical activity in the East model: analysing structure in biased trajectories

We consider large deviations of the dynamical activity in the East model. We bias this system to larger than average activity and investigate the structure that emerges. To best characterise this structure, we exploit the fact that there are effective interactions that would reproduce the same behaviour in an equilibrium system. We combine numerical results with linear response theory and variational estimates of these effective interactions, giving the first insights into such interactions in a many-body system, across a wide range of biases. The system exhibits a hierarchy of responses to the bias, remaining quasi-equilibrated on short length scales, but deviating far from equilibrium on large length scales. We discuss the connection between this hierarchy and the hierarchical aging behaviour of the system.

preprint2013arXiv

The effect of boundary adaptivity on hexagonal ordering and bistability in circularly confined quasi hard discs

The behaviour of materials under spatial confinement is sensitively dependent on the nature of the confining boundaries. In two dimensions, confinement within a hard circular boundary inhibits the hexagonal ordering observed in bulk systems at high density. Using colloidal experiments and Monte Carlo simulations, we investigate two model systems of quasi hard discs under circularly symmetric confinement. The first system employs an adaptive circular boundary, defined experimentally using holographic optical tweezers. We show that deformation of this boundary allows, and indeed is required for, hexagonal ordering in the confined system. The second system employs a circularly symmetric optical potential to confine particles without a physical boundary. We show that, in the absence of a curved wall, near perfect hexagonal ordering is possible. We propose that the degree to which hexagonal ordering is suppressed by a curved boundary is determined by the `strictness' of that wall.

preprint2012arXiv

Quantifying reversibility in a phase-separating lattice gas: an analogy with self-assembly

We present dynamic measurements of a lattice gas during phase separation, which we use as an analogy for self-assembly of equilibrium ordered structures. We use two approaches to quantify the degree of 'reversibility' of this process: firstly, we count events in which bonds are made and broken; secondly, we use correlation-response measurements and fluctuation-dissipation ratios to probe reversibility during different time intervals. We show how correlation and response functions can be related directly to microscopic (ir)reversibility and we discuss time-dependence and observable- dependence of these measurements, including the role of fast and slow degrees of freedom during assembly.

preprint2012arXiv

Space-time phase transitions in the East model with a softened kinetic constraint

The East model has a dynamical phase transition between an active (fluid) and inactive (glass) state. We show that this phase transition generalizes to "softened" systems where constraint violations are allowed with small but finite probabilities. Moreover, we show that the first order coexistence line separating the active and inactive phases terminates in a finite-temperature space-time critical point. Implications of these results for equilibrium dynamics are discussed.

preprint2011arXiv

Analyzing mechanisms and microscopic reversibility of self-assembly

We use computer simulations to investigate self-assembly in a system of model chaperonin proteins, and in an Ising lattice gas. We discuss the mechanisms responsible for rapid and efficient assembly in these systems, and we use measurements of dynamical activity and assembly progress to compare their propensities for kinetic trapping. We use the analytic solution of a simple minimal model to illustrate the key features associated with such trapping, paying particular attention to the number of ways that particles can misbind. We discuss the relevance of our results for the design and control of self-assembly in general.

preprint2011arXiv

Mechanisms of kinetic trapping in self-assembly and phase transformation

In self-assembly processes, kinetic trapping effects often hinder the formation of thermodynamically stable ordered states. In a model of viral capsid assembly and in the phase transformation of a lattice gas, we show how simulations in a self-assembling steady state can be used to identify two distinct mechanisms of kinetic trapping. We argue that one of these mechanisms can be adequately captured by kinetic rate equations, while the other involves a breakdown of theories that rely on cluster size as a reaction coordinate. We discuss how these observations might be useful in designing and optimising self-assembly reactions.

preprint2011arXiv

Predicting the self-assembly of a model colloidal crystal

We investigate the self-assembly (crystallisation) of particles with hard cores and isotropic, square-well interactions, using a Monte Carlo scheme to simulate overdamped Langevin dynamics. We measure correlation and response functions during the early stages of assembly, and we analyse the results using fluctuation-dissipation theorems, aiming to predict which systems will self-assemble successfully and which will get stuck in disordered states. The early-time correlation and response measurements are made before significant crystallisation has taken place, indicating that dynamical measurements are valuable in measuring a system's propensity for kinetic trapping.

preprint2011arXiv

Preparation and relaxation of very stable glassy states of a simulated liquid

We prepare metastable glassy states in a model glass-former made of Lennard-Jones particles by sampling biased ensembles of trajectories with low dynamical activity. These trajectories form an inactive dynamical phase whose `fast' vibrational degrees of freedom are maintained at thermal equilibrium by contact with a heat bath, while the `slow' structural degrees of freedom are located in deep valleys of the energy landscape. We examine the relaxation to equilibrium and the vibrational properties of these metastable states. The glassy states we prepare by our trajectory sampling method are very stable to thermal fluctuations and also more mechanically rigid than low-temperature equilibrated configurations.

preprint2011arXiv

Random pinning in glassy spin models with plaquette interactions

We use a random pinning procedure to study amorphous order in two glassy spin models. On increasing the concentration of pinned spins at constant temperature, we find a sharp crossover (but no thermodynamic phase transition) from bulk relaxation to localisation in a single state. At low temperatures, both models exhibit scaling behaviour. We discuss the growing length and time scales associated with amorphous order, and the fraction of pinned spins required to localize the system in a single state. These results, obtained for finite dimensional interacting models, provide a theoretical scenario for the effect of random pinning that differs qualitatively from previous approaches based either on mean-field, mode-coupling, or renormalization group reatments.

preprint2010arXiv

Finite-temperature critical point of a glass transition

We generalize the simplest kinetically constrained model of a glass-forming liquid by softening kinetic constraints, allowing them to be violated with a small finite rate. We demonstrate that this model supports a first-order dynamical (space-time) phase transition, similar to those observed with hard constraints. In addition, we find that the first-order phase boundary in this softened model ends in a finite-temperature dynamical critical point, which we expect to be present in natural systems. We discuss links between this critical point and quantum phase transitions, showing that dynamical phase transitions in $d$ dimensions map to quantum transitions in the same dimension, and hence to classical thermodynamic phase transitions in $d+1$ dimensions. We make these links explicit through exact mappings between master operators, transfer matrices, and Hamiltonians for quantum spin chains.

preprint2009arXiv

Metastable states and space-time phase transitions in a spin-glass model

We study large deviations of the dynamical activity in the random orthogonal model (ROM). This is a fully connected spin-glass model with one-step replica symmetry breaking behaviour, consistent with the random first-order transition scenario for structural glasses. We show that this model displays dynamical (space-time) phase-transitions between active and inactive phases, as demonstrated by singularities in large deviation functions. We argue that such transitions are generic in systems with long-lived metastable states.

preprint2009arXiv

Second-order dynamic transition in a p=2 spin-glass model

We consider the dynamics of a disordered p-spin model with p=2, analyzing the dynamics within Ruelle's thermodynamic formalism, We use an indicator of the dynamical activity to construct the relevant dynamical Gibbs ensemble. We prove that the dynamics in the low-temperature (spin glass) phase of the model take place at a second-order phase transition between dynamically active and inactive trajectories. We also show that the same behaviour is found in a related model of a three-dimensional ferromagnet.