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

Stephen Whitelam

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

ResearcherUnclaimed source record
cond-mat.stat-mechcond-mat.softBiomoleculesMachine LearningNeural and Evolutionary ComputingBiological Physicscond-mat.mes-hallcond-mat.mtrl-sciphysics.chem-phphysics.comp-phphysics.app-phSubcellular Processes

Catalog footprint

What is connected

39works
12topics
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

39 published item(s)

preprint2026arXiv

Nonlinear thermodynamic computing out of equilibrium

We present the design for a thermodynamic computer that can perform arbitrary nonlinear calculations in or out of equilibrium. Simple thermodynamic circuits, fluctuating degrees of freedom in contact with a thermal bath and confined by a quartic potential, display an activity that is a nonlinear function of their input. Such circuits can therefore be regarded as thermodynamic neurons, and can serve as the building blocks of networked structures that act as thermodynamic neural networks, universal function approximators whose operation is powered by thermal fluctuations. We simulate a digital model of a thermodynamic neural network, and show that its parameters can be adjusted by genetic algorithm to perform nonlinear calculations at specified observation times, regardless of whether the system has attained thermal equilibrium. This work expands the field of thermodynamic computing beyond the regime of thermal equilibrium, enabling fully nonlinear computations, analogous to those performed by classical neural networks, at specified observation times.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Learning stochastic dynamics and predicting emergent behavior using transformers

We show that a neural network originally designed for language processing can learn the dynamical rules of a stochastic system by observation of a single dynamical trajectory of the system, and can accurately predict its emergent behavior under conditions not observed during training. We consider a lattice model of active matter undergoing continuous-time Monte Carlo dynamics, simulated at a density at which its steady state comprises small, dispersed clusters. We train a neural network called a transformer on a single trajectory of the model. The transformer, which we show has the capacity to represent dynamical rules that are numerous and nonlocal, learns that the dynamics of this model consists of a small number of processes. Forward-propagated trajectories of the trained transformer, at densities not encountered during training, exhibit motility-induced phase separation and so predict the existence of a nonequilibrium phase transition. Transformers have the flexibility to learn dynamical rules from observation without explicit enumeration of rates or coarse-graining of configuration space, and so the procedure used here can be applied to a wide range of physical systems, including those with large and complex dynamical generators.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Training neural networks using Metropolis Monte Carlo and an adaptive variant

We examine the zero-temperature Metropolis Monte Carlo algorithm as a tool for training a neural network by minimizing a loss function. We find that, as expected on theoretical grounds and shown empirically by other authors, Metropolis Monte Carlo can train a neural net with an accuracy comparable to that of gradient descent, if not necessarily as quickly. The Metropolis algorithm does not fail automatically when the number of parameters of a neural network is large. It can fail when a neural network's structure or neuron activations are strongly heterogenous, and we introduce an adaptive Monte Carlo algorithm, aMC, to overcome these limitations. The intrinsic stochasticity and numerical stability of the Monte Carlo method allow aMC to train deep neural networks and recurrent neural networks in which the gradient is too small or too large to allow training by gradient descent. Monte Carlo methods offer a complement to gradient-based methods for training neural networks, allowing access to a distinct set of network architectures and principles.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Correspondence between neuroevolution and gradient descent

We show analytically that training a neural network by conditioned stochastic mutation or neuroevolution of its weights is equivalent, in the limit of small mutations, to gradient descent on the loss function in the presence of Gaussian white noise. Averaged over independent realizations of the learning process, neuroevolution is equivalent to gradient descent on the loss function. We use numerical simulation to show that this correspondence can be observed for finite mutations,for shallow and deep neural networks. Our results provide a connection between two families of neural-network training methods that are usually considered to be fundamentally different.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Improving the accuracy of nearest-neighbor classification using principled construction and stochastic sampling of training-set centroids

A conceptually simple way to classify images is to directly compare test-set data and training-set data. The accuracy of this approach is limited by the method of comparison used, and by the extent to which the training-set data cover configuration space. Here we show that this coverage can be substantially increased using coarse graining (replacing groups of images by their centroids) and stochastic sampling (using distinct sets of centroids in combination). We use the MNIST and Fashion-MNIST data sets to show that a principled coarse-graining algorithm can convert training images into fewer image centroids without loss of accuracy of classification of test-set images by nearest-neighbor classification. Distinct batches of centroids can be used in combination as a means of stochastically sampling configuration space, and can classify test-set data more accurately than can the unaltered training set. On the MNIST and Fashion-MNIST data sets this approach converts nearest-neighbor classification from a mid-ranking- to an upper-ranking member of the set of classical machine-learning techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Evolutionary reinforcement learning of dynamical large deviations

We show how to calculate the likelihood of dynamical large deviations using evolutionary reinforcement learning. An agent, a stochastic model, propagates a continuous-time Monte Carlo trajectory and receives a reward conditioned upon the values of certain path-extensive quantities. Evolution produces progressively fitter agents, eventually allowing the calculation of a piece of a large-deviation rate function for a particular model and path-extensive quantity. For models with small state spaces the evolutionary process acts directly on rates, and for models with large state spaces the process acts on the weights of a neural network that parameterizes the model's rates. This approach shows how path-extensive physics problems can be considered within a framework widely used in machine learning.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Learning to grow: control of material self-assembly using evolutionary reinforcement learning

We show that neural networks trained by evolutionary reinforcement learning can enact efficient molecular self-assembly protocols. Presented with molecular simulation trajectories, networks learn to change temperature and chemical potential in order to promote the assembly of desired structures or choose between competing polymorphs. In the first case, networks reproduce in a qualitative sense the results of previously-known protocols, but faster and with higher fidelity; in the second case they identify strategies previously unknown, from which we can extract physical insight. Networks that take as input the elapsed time of the simulation or microscopic information from the system are both effective, the latter more so. The evolutionary scheme we have used is simple to implement and can be applied to a broad range of examples of experimental self-assembly, whether or not one can monitor the experiment as it proceeds. Our results have been achieved with no human input beyond the specification of which order parameter to promote, pointing the way to the design of synthesis protocols by artificial intelligence.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Long-Range Exciton Diffusion in Two-Dimensional Assemblies of Cesium Lead Bromide Perovskite Nanocrystals

Förster Resonant Energy Transfer (FRET)-mediated exciton diffusion through artificial nanoscale building block assemblies could be used as a new optoelectronic design element to transport energy. However, so far nanocrystal (NC) systems supported only diffusion length of 30 nm, which are too small to be useful in devices. Here, we demonstrate a FRET-mediated exciton diffusion length of 200 nm with 0.5 cm2/s diffusivity through an ordered, two-dimensional assembly of cesium lead bromide perovskite nanocrystals (PNC). Exciton diffusion was directly measured via steady-state and time-resolved photoluminescence (PL) microscopy, with physical modeling providing deeper insight into the transport process. This exceptionally efficient exciton transport is facilitated by PNCs high PL quantum yield, large absorption cross-section, and high polarizability, together with minimal energetic and geometric disorder of the assembly. This FRET-mediated exciton diffusion length matches perovskites optical absorption depth, opening the possibility to design new optoelectronic device architectures with improved performances, and providing insight into the high conversion efficiencies of PNC-based optoelectronic devices.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Minimal physical requirements for crystal growth self-poisoning

Self-poisoning is a kinetic trap that can impair or prevent crystal growth in a wide variety of physical settings. Here we use dynamic mean-field theory and computer simulation to argue that poisoning is ubiquitous because its emergence requires only the notion that a molecule can bind in two (or more) ways to a crystal; that those ways are not energetically equivalent; and that the associated binding events occur with sufficiently unequal probability. If these conditions are met then the steady-state growth rate is in general a non-monotonic function of the thermodynamic driving force for crystal growth, which is the characteristic of poisoning. Our results also indicate that relatively small changes of system parameters could be used to induce recovery from poisoning.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the microscopic origin and macroscopic implications of lane formation in mixtures of oppositely-driven particles

Colloidal particles of two types, driven in opposite directions, can segregate into lanes [Vissers et al. Soft Matter 7, 2352 (2011)]. This phenomenon can be reproduced by two-dimensional Brownian dynamics simulations of model particles [Dzubiella et al. Phys. Rev. E 65, 021402 (2002)]. Here we use computer simulation to assess the generality of lane formation with respect to variation of particle type and dynamical protocol. We find that laning results from rectification of diffusion on the scale of a particle diameter: oppositely-driven particles must, in the time taken to encounter each other in the direction of the drive, diffuse in the perpendicular direction by about one particle diameter. This geometric constraint implies that the diffusion constant of a particle, in the presence of those of the opposite type, grows approximately linearly with Peclet number, a prediction confirmed by our numerics over a range of model parameters. Such environment-dependent diffusion is statistically similar to an effective interparticle attraction; consistent with this observation, we find that oppositely-driven non-attractive colloids display features characteristic of the simplest model system possessing both interparticle attractions and persistent motion, the driven Ising lattice gas [Katz, Leibowitz, Spohn, J. Stat. Phys. 34, 497 (1984)]. These features include long-ranged correlations in the disordered regime, and a critical regime characterized by a change in slope of the particle current with Peclet number and by fluctuations that grow with system size. By analogy, we suggest that lane formation in the driven colloid system is in the macroscopic limit a phase transition, but that macroscopic phase separation would not occur in finite time upon starting from disordered initial conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Selective gas capture via kinetic trapping

Conventional approaches to the capture of CO_2 by metal-organic frameworks focus on equilibrium conditions, and frameworks that contain little CO_2 in equilibrium are often rejected as carbon-capture materials. Here we use a statistical mechanical model, parameterized by quantum mechanical data, to suggest that metal-organic frameworks can be used to separate CO_2 from a typical flue gas mixture when used under {\em nonequilibrium} conditions. The origin of this selectivity is an emergent gas-separation mechanism that results from the acquisition by different gas types of different mobilities within a crowded framework. The resulting distribution of gas types within the framework is in general spatially and dynamically heterogeneous. Our results suggest that relaxing the requirement of equilibrium can substantially increase the parameter space of conditions and materials for which selective gas capture can be effected.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Crystallization and arrest mechanisms of model colloids

We performed dynamic simulations of spheres with short-range attractive interactions for many values of interaction strength and range. Fast crystallization occurs in a localized region of this parameter space, but the character of crystallization pathways is not uniform within this region. Pathways range from one-step, in which a crystal nucleates directly from a gas, to two-step, in which substantial liquid-like clusters form and only subsequently become crystalline. Crystallization can fail because of slow nucleation from either gas or liquid, or because of dynamic arrest caused by strong interactions. Arrested states are characterized by the formation of networks of face-sharing tetrahedra that can be detected by a local common neighbor analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hierarchical assembly may be a way to make large information-rich structures

Self-assembly in the laboratory can now yield `information-rich' nanostructures in which each component is of a distinct type and has a defined spatial position. Ensuring the thermodynamic stability of such structures requires inter-component interaction energies to increase logarithmically with structure size, in order to counter the entropy gained upon mixing component types in solution. However, self-assembly in the presence of strong interactions results in general in kinetic trapping, so suggesting a limit to the size of an (equilibrium) structure that can be self-assembled from distinguishable components. Here we study numerically a two-dimensional hierarchical assembly scheme already considered in experiment. We show that this scheme is immune to the kinetic traps associated with strong `native' interactions (interactions designed to stabilize the intended structure), and so, in principle, offers a way to make large information-rich structures. In this scheme the size of an assembled structure scales exponentially with the stage of assembly, and assembly can continue as long as random motion is able to bring structures into contact. The resulting superstructure could provide a template for building in the third dimension. The chief drawback of this scheme is that it is particularly susceptible to kinetic traps that result from `non-native' interactions (interactions not required to stabilize the intended structure); the scale on which such a scheme can be realized therefore depends upon how effectively this latter kind of interaction can be suppressed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Predicting the outcome of the growth of binary solids far from equilibrium

The growth of multicomponent structures in simulations and experiments often results in kinetically trapped, nonequilibrium objects. In such cases we have no general theoretical framework for predicting the outcome of the growth process. Here we use computer simulations to study the growth of two-component structures within a simple lattice model. We show that kinetic trapping happens for many choices of growth rate and inter-component interaction energies, and that qualitatively distinct kinds of kinetic trapping are found in different regions of parameter space. In a region in which the low-energy structure is an `antiferromagnet' or `checkerboard', we show that the grown nonequilibrium structure displays a component-type stoichiometry that is different to the equilibrium one but is insensitive to growth rate and solution conditions. This robust nonequilibrium stoichiometry can be predicted via a mapping to the jammed random tiling of dimers studied by Flory, a finding that suggests a way of making defined nonequilibrium structures in experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Ramachandran number: an order parameter for protein geometry

Three-dimensional protein structures usually contain regions of local order, called secondary structure, such as $α$-helices and $β$-sheets. Secondary structure is characterized by the local rotational state of the protein backbone, quantified by two dihedral angles called $ϕ$ and $ψ$. Particular types of secondary structure can generally be described by a single (diffuse) location on a two-dimensional plot drawn in the space of the angles $ϕ$ and $ψ$, called a Ramachandran plot. By contrast, a recently-discovered nanomaterial made from peptoids, structural isomers of peptides, displays a secondary-structure motif corresponding to two regions on the Ramachandran plot [Mannige et al., Nature 526, 415 (2015)]. In order to describe such `higher-order' secondary structure in a compact way we introduce here a means of describing regions on the Ramachandran plot in terms of a single Ramachandran number, ${\mathcal{R}}$, which is a structurally meaningful combination of $ϕ$ and $ψ$. We show that the potential applications of ${\mathcal{R}}$ are numerous: it can be used to describe the geometric content of protein structures, and can be used to draw diagrams that reveal, at a glance, the frequency of occurrence of regular secondary structures and disordered regions in large protein datasets. We propose that ${\mathcal{R}}$ might be used as an order parameter for protein geometry for a wide range of applications.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Common physical framework explains phase behavior and dynamics of atomic, molecular and polymeric network-formers

We show that the self-assembly of a diverse collection of building blocks can be understood within a common physical framework. These building blocks, which form periodic honeycomb networks and nonperiodic variants thereof, range in size from atoms to micron-scale polymers, and interact through mechanisms as different as hydrogen bonds and covalent forces. A combination of statistical mechanics and quantum mechanics shows that one can capture the physics that governs assembly of these networks by resolving only the geometry and strength of building block interactions. The resulting framework reproduces a broad range of phenomena seen experimentally, including periodic and nonperiodic networks in thermal equilibrium, and nonperiodic supercooled and glassy networks away from equilibrium. Our results show how simple `design criteria' control assembly of a wide variety of networks, and suggest that kinetic trapping can be a useful way of making functional assemblies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Emergent rhombus tilings from molecular interactions with $M$-fold rotational symmetry

We show that model molecules with particular rotational symmetries can self-assemble into network structures equivalent to rhombus tilings. This assembly happens in an emergent way, in the sense that molecules spontaneously select irregular 4-fold local coordination from a larger set of possible local binding geometries. The existence of such networks can be rationalized by simple geometrical arguments, but the same arguments do not guarantee networks' spontaneous self-assembly. This class of structures must in certain regimes of parameter space be able to reconfigure into networks equivalent to triangular tilings.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Examples of molecular self-assembly at surfaces

The self-assembly of molecules at surfaces can be caused by a range of physical mechanisms. Assembly can be driven by intermolecular forces, or molecule-surface forces, or both; it can result in structures that are in equilibrium or that are kinetically trapped. Here we review examples of self-assembly at surfaces that have been studied within the User program of the Molecular Foundry at Lawrence Berkeley National Laboratory, focusing on a physical understanding of what causes patterns seen in experiment. Some apparently disparate systems can be described in similar physical terms, indicating that simple factors -- such as the geometry and energy scale of intermolecular binding -- are key to understanding the self-assembly of those systems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Growth of equilibrium structures built from a large number of distinct component types

We use simple analytic arguments and lattice-based computer simulations to study the growth of structures made from a large number of distinct component types. Components possess 'designed' interactions, chosen to stabilize an equilibrium target structure in which each component type has a defined spatial position, and 'undesigned' interactions that allow components to bind in a compositionally-disordered way. We find that high-fidelity growth of the equilibrium target structure can happen in the presence of substantial attractive undesigned interactions, as long as the energy scale of the set of designed interactions is chosen appropriately. This observation may help explain why equilibrium DNA 'brick' structures self-assemble even if undesigned interactions are not suppressed [Ke et al. Science 338, 1177 (2012)]. We also find that high-fidelity growth of the target structure is most probable when designed interactions are drawn from a distribution that is as narrow as possible. We use this result to suggest how to choose complementary DNA sequences in order to maximize the fidelity of multicomponent self-assembly mediated by complementary DNA interactions. We also comment on the prospect of growing macroscopic structures in this manner

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Modeling sequence-specific polymers using anisotropic coarse-grained sites allows quantitative comparison with experiment

Certain sequences of peptoid polymers (synthetic analogs of peptides) assemble into bilayer nanosheets via a nonequilibrium assembly pathway of adsorption, compression, and collapse at an air-water interface. As with other large-scale dynamic processes in biology and materials science, understanding the details of this supramolecular assembly process requires a modeling approach that captures behavior on a wide range of length and time scales, from those on which individual sidechains fluctuate to those on which assemblies of polymers evolve. Here we demonstrate that a new coarse-grained modeling approach is accurate and computationally efficient enough to do so. Our approach uses only a minimal number of coarse-grained sites, but retains independently fluctuating orientational degrees of freedom for each site. These orientational degrees of freedom allow us to accurately parameterize both bonded and nonbonded interactions, and to generate all-atom configurations with sufficient accuracy to perform atomic scattering calculations and to interface with all-atom simulations. We have used this approach to reproduce all available experimental X-ray scattering spectra (for stacked nanosheets, and for peptoids adsorbed at air-water interfaces and in solution), in order to resolve the microscopic, real-space structures responsible for these Fourier-space features. By interfacing with all-atom simulations, we have also laid the foundations for future multiscale simulations of sequence-specific polymers that communicate in both directions across scales.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Self-assembly at a nonequilibrium critical point

We use analytic theory and computer simulation to study patterns formed during the growth of two-component assemblies in 2D and 3D. We show that these patterns undergo a nonequilibrium phase transition, at a particular growth rate, between mixed and demixed arrangements of component types. This finding suggests that principles of nonequilibrium statistical mechanics can be used to predict the outcome of multicomponent self-assembly, and suggests an experimental route to the self-assembly of multicomponent structures of a qualitatively defined nature.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Competing thermodynamic and dynamic factors select molecular assemblies on a gold surface

Controlling the self-assembly of surface-adsorbed molecules into nanostructures requires understanding physical mechanisms that act across multiple length and time scales. By combining scanning tunneling microscopy with hierarchical ab initio and statistical mechanical modeling of 1,4-substituted benzenediamine (BDA) molecules adsorbed on a gold (111) surface, we demonstrate that apparently simple nanostructures are selected by a subtle competition of thermodynamics and dynamics. Of the collection of possible BDA nanostructures mechanically stabilized by hydrogen bonding, the interplay of intermolecular forces, surface modulation, and assembly dynamics select at low temperature a particular subset: low free energy oriented linear chains of monomers, and high free energy branched chains.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Selective nucleation in porous media

Geometrical arguments suggest that pore-mediated nucleation happens in general in a two-step fashion, the first step being nucleation within the pore, the second being nucleation from the filled pore into solution [Page & Sear, Phys. Rev. Lett. 97, 65701 (2006)]. The free energy barriers controlling the two steps of this process show opposite dependencies on pore size, implying that for given thermodynamic conditions there exists a pore size for which nucleation happens fastest. Here we show, within the two- and three-dimensional Ising lattice gas, that this preferred pore size tracks the size of the bulk critical nucleus, up to a numerical prefactor. This observation suggests a simple prescription for directing nucleation to certain locations within heterogeneous porous media.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Control of pathways and yields of protein crystallization through the interplay of nonspecific and specific attractions

We use computer simulation to study crystal-forming model proteins equipped with interactions that are both orientationally specific and nonspecific. Distinct dynamical pathways of crystal formation can be selected by tuning the strengths of these interactions. When the nonspecific interaction is strong, liquidlike clustering can precede crystallization; when it is weak, growth can proceed via ordered nuclei. Crystal yields are in certain parameter regimes enhanced by the nonspecific interaction, even though it promotes association without local crystalline order. Our results suggest that equipping nanoscale components with weak nonspecific interactions (such as depletion attractions) can alter both their dynamical pathway of assembly and optimize the yield of the resulting material.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Do hierarchical structures assemble best via hierarchical pathways?

Hierarchically structured natural materials possess functionalities unattainable to the same components organized or mixed in simpler ways. For instance, the bones and teeth of mammals are far stronger and more durable than the mineral phases from which they are derived because their constituents are organized hierarchically from the molecular scale to the macroscale. Making similarly functional synthetic hierarchical materials will require an understanding of how to promote the self-assembly of structure on multiple lengthscales, without falling foul of numerous possible kinetic traps. Here we use computer simulation to study the self-assembly of a simple hierarchical structure, a square crystal lattice whose repeat unit is a tetramer. Although the target material is organized hierarchically, it self-assembles most reliably when its dynamic assembly pathway consists of the sequential addition of monomers to a single structure. Hierarchical dynamic pathways via dimer and tetramer intermediates are also viable modes of assembly, but result in general in lower yield: these intermediates have a stronger tendency than monomers to associate in ways not compatible with the target structure. In addition, assembly via tetramers results in a kinetic trap whereby material is sequestered in trimers that cannot combine to form the target crystal. We use analytic theory to relate dynamical pathways to the presence of equilibrium phases close in free energy to the target structure, and to identify the thermodynamic principles underpinning optimum self-assembly in this model: 1) make the free energy gap between the target phase and the most stable fluid phase of order kT, and 2) ensure that no other dense phases (liquids or close-packed solids of monomers or oligomers) or fluids of incomplete building blocks fall within this gap.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Patterning a surface so as to speed nucleation from solution

Motivated by the question of how to pattern a surface in order to best speed nucleation from solution, we build on the work of Page and Sear [Phys. Rev. Lett. 97, 65701 (2006)] and calculate rates and free energy profiles for nucleation in the 3d Ising model in the presence of cuboidal pores. Pores of well-chosen aspect ratio can dramatically speed nucleation relative to a planar surface made of the same material, while badly-chosen pores provide no such enhancement. For a given pore, the maximum nucleation rate is achieved when one of its two horizontal dimensions attains a critical length, largely irrespective of the other dimension (provided that the latter is large enough). This observation implies that patterning a surface in a raster-like fashion is a better strategy for speeding nucleation than e.g. scoring long grooves in it.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Random and ordered phases of off-lattice rhombus tiles

We study the covering of the plane by non-overlapping rhombus tiles, a problem well-studied only in the limiting case of dimer coverings of regular lattices. We go beyond this limit by allowing tiles to take any position and orientation on the plane, to be of irregular shape, and to possess different types of attractive interactions. Using extensive numerical simulations we show that at large tile densities there is a phase transition from a fluid of rhombus tiles to a solid packing with broken rotational symmetry. We observe self-assembly of broken-symmetry phases, even at low densities, in the presence of attractive tile-tile interactions. Depending on tile shape and interactions the solid phase can be random, possessing critical orientational fluctuations, or crystalline. Our results suggest strategies for controlling tiling order in experiments involving `molecular rhombi'.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Self-assembly of multicomponent structures in and out of equilibrium

Theories of phase change and self-assembly often invoke the idea of a `quasiequilibrium', a regime in which the nonequilibrium association of building blocks results nonetheless in a structure whose properties are determined solely by an underlying free energy landscape. Here we study a prototypical example of multicomponent self-assembly, a one-dimensional fiber grown from red and blue blocks. If the equilibrium structure possesses compositional correlations different from those characteristic of random mixing, then it cannot be generated without error at any finite growth rate: there is no quasiequilibrium regime. However, by exploiting dynamic scaling, structures characteristic of equilibrium at one point in phase space can be generated, without error, arbitrarily far from equilibrium. Our results thus suggest a `nonperturbative' strategy for multicomponent self-assembly in which the target structure is, by design, not the equilibrium one.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Approximating the dynamical evolution of systems of strongly-interacting overdamped particles

We describe collective-move Monte Carlo algorithms designed to approximate the overdamped dynamics of self-assembling nanoscale components equipped with strong, short-ranged and anisotropic interactions. Conventional Monte Carlo simulations comprise sequential moves of single particles, proposed and accepted so as to satisfy detailed balance. Under certain circumstances such simulations provide an approximation of overdamped dynamics, but the accuracy of this approximation can be poor if e.g. particle-particle interactions vary strongly with distance or angle. The twin requirements of simulation efficiency (trial moves of appreciable scale are needed to ensure reasonable sampling) and dynamical fidelity (true in the limit of vanishingly small trial moves) then become irreconcilable. As a result, single-particle moves can underrepresent important collective modes of relaxation, such as self-diffusion of particle clusters. However, one way of using Monte Carlo simulation to mimic real collective modes of motion, retaining the ability to make trial moves of reasonable scale, is to make explicit moves of collections of particles. We will outline ways of doing so by iteratively linking particles to their environment. Linking criteria can be static, conditioned upon properties of the current state of a system, or dynamic, conditioned upon energy changes resulting from trial virtual moves of particles. We argue that the latter protocol is better-suited to approximating real dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Design rules for the self-assembly of a protein crystal

Theories of protein crystallization based on spheres that form close-packed crystals predict optimal assembly within a `slot' of second virial coefficients and enhanced assembly near the metastable liquid-vapor critical point. However, most protein crystals are open structures stabilized by anisotropic interactions. Here, we use theory and simulation to show that assembly of one such structure is not predicted by the second virial coefficient or enhanced by the critical point. Instead, good assembly requires that the thermodynamic driving force be on the order of the thermal energy and that interactions be made as nonspecific as possible without promoting liquid-vapor phase separation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Electrostatics and aggregation: how charge can turn a crystal into a gel

The crystallization of proteins or colloids is often hindered by the appearance of aggregates of low fractal dimension called gels. Here we study the effect of electrostatics upon crystal and gel formation using an analytic model of hard spheres bearing point charges and short range attractive interactions. We find that the chief electrostatic free energy cost of forming assemblies comes from the entropic loss of counterions that render assemblies charge-neutral. Because there exists more accessible volume for these counterions around an open gel than a dense crystal, there exists an electrostatic entropic driving force favoring the gel over the crystal. This driving force increases with increasing sphere charge, but can be counteracted by increasing counterion concentration. We show that these effects cannot be fully captured by pairwise-additive macroion interactions of the kind often used in simulations, and we show where on the phase diagram to go in order to suppress gel formation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Limit of validity of Ostwald's rule of stages in a statistical mechanical model of crystallization

We have only rules of thumb with which to predict how a material will crystallize, chief among which is Ostwald's rule of stages. It states that the first phase to appear upon transformation of a parent phase is the one closest to it in free energy. Although sometimes upheld, the rule is without theoretical foundation and is not universally obeyed, highlighting the need for microscopic understanding of crystallization controls. Here we study in detail the crystallization pathways of a prototypical model of patchy particles. The range of crystallization pathways it exhibits is richer than can be predicted by Ostwald's rule, but a combination of simulation and analytic theory reveals clearly how these pathways are selected by microscopic parameters. Our results suggest strategies for controlling self-assembly pathways in simulation and experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Microscopic implications of S-DNA

Recent experiments [J. van Mameren et al. PNAS 106, 18231 (2009)] provide a detailed spatial picture of overstretched DNA, showing that under certain conditions the two strands of the double helix separate at about 65 pN. It was proposed that this observation rules out the existence of an elongated, hybridized form of DNA ('S-DNA'). Here we argue that the S-DNA picture is consistent with the observation of unpeeling during overstretching. We demonstrate that assuming the existence of S-DNA does not imply DNA overstretching to consist of the complete or near-complete conversion of the molecule from B- to S-form. Instead, this assumption implies in general a more complex dynamic coexistence of hybridized and unhybridized forms of DNA. We argue that such coexistence can rationalize several recent experimental observations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Nonclassical assembly pathways of anisotropic particles

Advances in synthetic methods have spawned an array of nanoparticles and bio-inspired molecules of diverse shapes and interaction geometries. Recent experiments indicate that such anisotropic particles exhibit a variety of 'nonclassical' self-assembly pathways, forming ordered assemblies via intermediates that do not share the architecture of the bulk material. Here we apply mean field theory to a prototypical model of interacting anisotropic particles, and find a clear thermodynamic impetus for nonclassical ordering in certain regimes of parameter space. In other parameter regimes, by contrast, assembly pathways are selected by dynamics. This approach suggests a means of predicting when anisotropic particles might assemble in a manner more complicated than that assumed by classical nucleation theory.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Self-assembly of amphiphilic peanut-shaped nanoparticles

We use computer simulation to investigate the self-assembly of Janus-like amphiphilic peanut-shaped nanoparticles, finding phases of clusters, bilayers and micelles in accord with ideas of packing familiar from the study of molecular surfactants. However, packing arguments do not explain the hierarchical self-assembly dynamics that we observe, nor the coexistence of bilayers and faceted polyhedra. This coexistence suggests that experimental realizations of our model can achieve multipotent assembly of either of two competing ordered structures.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Transformation from spots to waves in a model of actin pattern formation

Actin networks in certain single-celled organisms exhibit a complex pattern-forming dynamics that starts with the appearance of static spots of actin on the cell cortex. Spots soon become mobile, executing persistent random walks, and eventually give rise to traveling waves of actin. Here we describe a possible physical mechanism for this distinctive set of dynamic transformations, by equipping an excitable reaction-diffusion model with a field describing the spatial orientation of its chief constituent (which we consider to be actin). The interplay of anisotropic actin growth and spatial inhibition drives a transformation at fixed parameter values from static spots to moving spots to waves.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

There and (slowly) back again: Entropy-driven hysteresis in a model of DNA overstretching

When pulled along its axis, double-stranded DNA elongates abruptly at a force of about 65 pN. Two physical pictures have been developed to describe this overstretched state. The first proposes that strong forces induce a phase transition to a molten state consisting of unhybridized single strands. The second picture instead introduces an elongated hybridized phase, called S-DNA, structurally and thermodynamically distinct from standard B-DNA. Little thermodynamic evidence exists to discriminate directly between these competing pictures. Here we show that within a microscopic model of DNA we can distinguish between the dynamics associated with each. In experiment, considerable hysteresis in a cycle of stretching and shortening develops as temperature is increased. Since there are few possible causes of hysteresis in a system whose extent is appreciable in only one dimension, such behavior offers a discriminating test of the two pictures of overstretching. Most experiments are performed upon nicked DNA, permitting the detachment (`unpeeling') of strands. We show that the long-wavelength progression of the unpeeled front generates hysteresis, the character of which agrees with experiment only if we assume the existence of S-DNA. We also show that internal melting (distinct from unpeeling) can generate hysteresis, the degree of which is strongly dependent upon the nonextensive loop entropy of single-stranded DNA.

0 citations0 reviewsNo code linkedOpen access