Sire: An interoperability engine for prototyping algorithms and exchanging information between molecular simulation programs.

Sire is a Python/C++ library that is used both to prototype new algorithms and as an interoperability engine for exchanging information between molecular simulation programs. It provides a collection of file parsers and information converters that together make it easier to combine and leverage the functionality of many other programs and libraries. This empowers researchers to use sire to write a single script that can, for example, load a molecule from a PDBx/mmCIF file via Gemmi, perform SMARTS searches via RDKit, parameterize molecules using BioSimSpace, run GPU-accelerated molecular dynamics via OpenMM, and then display the resulting dynamics trajectory in a NGLView Jupyter notebook 3D molecular viewer. This functionality is built on by BioSimSpace, which uses sire's molecular information engine to interconvert with programs such as GROMACS, NAMD, Amber, and AmberTools for automated molecular parameterization and the running of molecular dynamics, metadynamics, and alchemical free energy workflows. Sire comes complete with a powerful molecular information search engine, plus trajectory loading and editing, analysis, and energy evaluation engines. This, when combined with an in-built computer algebra system, gives substantial flexibility to researchers to load, search for, edit, and combine molecular information from multiple sources and use that to drive novel algorithms by combining functionality from other programs. Sire is open source (GPL3) and is available via conda and at a free Jupyter notebook server at https://try.openbiosim.org. Sire is supported by the not-for-profit OpenBioSim community interest company.

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.

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, as well as '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 DNA. We also comment on the prospect of growing macroscopic structures in this manner.

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 two and three dimensions. 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.

Microscopic evidence for liquid-liquid separation in supersaturated CaCO3 solutions.

Recent experimental observations of the onset of calcium carbonate (CaCO3) mineralization suggest the emergence of a population of clusters that are stable rather than unstable as predicted by classical nucleation theory. This study uses molecular dynamics simulations to probe the structure, dynamics, and energetics of hydrated CaCO3 clusters and lattice gas simulations to explore the behavior of cluster populations before nucleation. Our results predict formation of a dense liquid phase through liquid-liquid separation within the concentration range in which clusters are observed. Coalescence and solidification of nanoscale droplets results in formation of a solid phase, the structure of which is consistent with amorphous CaCO3. The presence of a liquid-liquid binodal enables a diverse set of experimental observations to be reconciled within the context of established phase-separation mechanisms.

Uncovering the intrinsic size dependence of hydriding phase transformations in nanocrystals.

A quantitative understanding of nanocrystal phase transformations would enable more efficient energy conversion and catalysis, but has been hindered by difficulties in directly monitoring well-characterized nanoscale systems in reactive environments. We present a new in situ luminescence-based probe enabling direct quantification of nanocrystal phase transformations, applied here to the hydriding transformation of palladium nanocrystals. Our approach reveals the intrinsic kinetics and thermodynamics of nanocrystal phase transformations, eliminating complications of substrate strain, ligand effects and external signal transducers. Clear size-dependent trends emerge in nanocrystals long accepted to be bulk-like in behaviour. Statistical mechanical simulations show these trends to be a consequence of nanoconfinement of a thermally driven, first-order phase transition: near the phase boundary, critical nuclei of the new phase are comparable in size to the nanocrystal itself. Transformation rates are then unavoidably governed by nanocrystal dimensions. Our results provide a general framework for understanding how nanoconfinement fundamentally impacts broad classes of thermally driven solid-state phase transformations relevant to hydrogen storage, catalysis, batteries and fuel cells.

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.

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.

Dynamic order-disorder in atomistic models of structural glass formers.

The glass transition is the freezing of a liquid into a solid state without evident structural order. Although glassy materials are well characterized experimentally, the existence of a phase transition into the glass state remains controversial. Here, we present numerical evidence for the existence of a novel first-order dynamical phase transition in atomistic models of structural glass formers. In contrast to equilibrium phase transitions, which occur in configuration space, this transition occurs in trajectory space, and it is controlled by variables that drive the system out of equilibrium. Coexistence is established between an ergodic phase with finite relaxation time and a nonergodic phase of immobile molecular configurations. Thus, we connect the glass transition to a true phase transition, offering the possibility of a unified picture of glassy phenomena.

Decoupling of exchange and persistence times in atomistic models of glass formers.

With molecular dynamics simulations of a fluid mixture of classical particles interacting with pairwise additive Weeks-Chandler-Andersen potentials, we consider the time series of particle displacements and thereby determine the distributions for local persistence times and local exchange times. These basic characterizations of glassy dynamics are studied over a range of supercooled conditions and were shown to have behaviors, most notably decoupling, similar to those found in kinetically constrained lattice models of structural glasses. Implications are noted.