Designing self-assembling kinetics with differentiable statistical physics models.

The inverse problem of designing component interactions to target emergent structure is fundamental to numerous applications in biotechnology, materials science, and statistical physics. Equally important is the inverse problem of designing emergent kinetics, but this has received considerably less attention. Using recent advances in automatic differentiation, we show how kinetic pathways can be precisely designed by directly differentiating through statistical physics models, namely free energy calculations and molecular dynamics simulations. We consider two systems that are crucial to our understanding of structural self-assembly: bulk crystallization and small nanoclusters. In each case, we are able to assemble precise dynamical features. Using gradient information, we manipulate interactions among constituent particles to tune the rate at which these systems yield specific structures of interest. Moreover, we use this approach to learn nontrivial features about the high-dimensional design space, allowing us to accurately predict when multiple kinetic features can be simultaneously and independently controlled. These results provide a concrete and generalizable foundation for studying nonstructural self-assembly, including kinetic properties as well as other complex emergent properties, in a vast array of systems.

Using active colloids as machines to weave and braid on the micrometer scale.

Controlling motion at the microscopic scale is a fundamental goal in the development of biologically inspired systems. We show that the motion of active, self-propelled colloids can be sufficiently controlled for use as a tool to assemble complex structures such as braids and weaves out of microscopic filaments. Unlike typical self-assembly paradigms, these structures are held together by geometric constraints rather than adhesive bonds. The out-of-equilibrium assembly that we propose involves precisely controlling the 2D motion of active colloids so that their path has a nontrivial topology. We demonstrate with proof-of-principle Brownian dynamics simulations that, when the colloids are attached to long semiflexible filaments, this motion causes the filaments to braid. The ability of the active particles to provide sufficient force necessary to bend the filaments into a braid depends on a number of factors, including the self-propulsion mechanism, the properties of the filament, and the maximum curvature in the braid. Our work demonstrates that nonequilibrium assembly pathways can be designed using active particles.

Designing allostery-inspired response in mechanical networks.

Recent advances in designing metamaterials have demonstrated that global mechanical properties of disordered spring networks can be tuned by selectively modifying only a small subset of bonds. Here, using a computationally efficient approach, we extend this idea to tune more general properties of networks. With nearly complete success, we are able to produce a strain between any two target nodes in a network in response to an applied source strain on any other pair of nodes by removing only ∼1% of the bonds. We are also able to control multiple pairs of target nodes, each with a different individual response, from a single source, and to tune multiple independent source/target responses simultaneously into a network. We have fabricated physical networks in macroscopic 2D and 3D systems that exhibit these responses. This work is inspired by the long-range coupled conformational changes that constitute allosteric function in proteins. The fact that allostery is a common means for regulation in biological molecules suggests that it is a relatively easy property to develop through evolution. In analogy, our results show that long-range coupled mechanical responses are similarly easy to achieve in disordered networks.

Scaling ansatz for the jamming transition.

We propose a Widom-like scaling ansatz for the critical jamming transition. Our ansatz for the elastic energy shows that the scaling of the energy, compressive strain, shear strain, system size, pressure, shear stress, bulk modulus, and shear modulus are all related to each other via scaling relations, with only three independent scaling exponents. We extract the values of these exponents from already known numerical or theoretical results, and we numerically verify the resulting predictions of the scaling theory for the energy and residual shear stress. We also derive a scaling relation between pressure and residual shear stress that yields insight into why the shear and bulk moduli scale differently. Our theory shows that the jamming transition exhibits an emergent scale invariance, setting the stage for the potential development of a renormalization group theory for jamming.

Reply to the 'Comment on "Spatial structure of states of self stress in jammed systems"' by E. Lerner, Soft Matter, 2017, 13, DOI.

Lerner's theoretical analysis neglects a normalizing factor which distinguishes states of self stress from the stress response to a unit dipolar force along a bond. This factor leads to different spatial profiles upon ensemble averaging.

Divergence of Voronoi Cell Anisotropy Vector: A Threshold-Free Characterization of Local Structure in Amorphous Materials.

Characterizing structural inhomogeneity is an essential step in understanding the mechanical response of amorphous materials. We introduce a threshold-free measure based on the field of vectors pointing from the center of each particle to the centroid of the Voronoi cell in which the particle resides. These vectors tend to point in toward regions of high free volume and away from regions of low free volume, reminiscent of sinks and sources in a vector field. We compute the local divergence of these vectors, where positive values correspond to overpacked regions and negative values identify underpacked regions within the material. Distributions of this divergence are nearly Gaussian with zero mean, allowing for structural characterization using only the moments of the distribution. We explore how the standard deviation and skewness vary with the packing fraction for simulations of bidisperse systems and find a kink in these moments that coincides with the jamming transition.

Pinning Susceptibility: The Effect of Dilute, Quenched Disorder on Jamming.

We study the effect of dilute pinning on the jamming transition. Pinning reduces the average contact number needed to jam unpinned particles and shifts the jamming threshold to lower densities, leading to a pinning susceptibility, χ_{p}. Our main results are that this susceptibility obeys scaling form and diverges in the thermodynamic limit as χ_{p}∝|ϕ-ϕ_{c}^{∞}|^{-γ_{p}} where ϕ_{c}^{∞} is the jamming threshold in the absence of pins. Finite-size scaling arguments yield these values with associated statistical (systematic) errors γ_{p}=1.018±0.026(0.291) in d=2 and γ_{p}=1.534±0.120(0.822) in d=3. Logarithmic corrections raise the exponent in d=2 to close to the d=3 value, although the systematic errors are very large.

Spatial structure of states of self stress in jammed systems.

States of self stress, organizations of internal forces in many-body systems that are in equilibrium with an absence of external forces, can be thought of as the constitutive building blocks of the elastic response of a material. In overconstrained disordered packings they have a natural mathematical correspondence with the zero-energy vibrational modes in underconstrained systems. While substantial attention in the literature has been paid to diverging length scales associated with zero- and finite-energy vibrational modes in jammed systems, less is known about the spatial structure of the states of self stress. In this work we define a natural way in which a unique state of self stress can be associated with each bond in a disordered spring network derived from a jammed packing, and then investigate the spatial structure of these bond-localized states of self stress. This allows for an understanding of how the elastic properties of a system would change upon changing the strength or even existence of any bond in the system.

The Principle of Independent Bond-Level Response: Tuning by Pruning to Exploit Disorder for Global Behavior.

We introduce a principle unique to disordered solids wherein the contribution of any bond to one global perturbation is uncorrelated with its contribution to another. Coupled with sufficient variability in the contributions of different bonds, this "independent bond-level response" paves the way for the design of real materials with unusual and exquisitely tuned properties. To illustrate this, we choose two global perturbations: compression and shear. By applying a bond removal procedure that is both simple and experimentally relevant to remove a very small fraction of bonds, we can drive disordered spring networks to both the incompressible and completely auxetic limits of mechanical behavior.

Disordered surface vibrations in jammed sphere packings.

We study the vibrational properties near a free surface of disordered spring networks derived from jammed sphere packings. In bulk systems, without surfaces, it is well understood that such systems have a plateau in the density of vibrational modes extending down to a frequency scale ω*. This frequency is controlled by ΔZ = 〈Z〉 - 2d, the difference between the average coordination of the spheres and twice the spatial dimension, d, of the system, which vanishes at the jamming transition. In the presence of a free surface we find that there is a density of disordered vibrational modes associated with the surface that extends far below ω*. The total number of these low-frequency surface modes is controlled by ΔZ, and the profile of their decay into the bulk has two characteristic length scales, which diverge as ΔZ(-1/2) and ΔZ(-1) as the jamming transition is approached.

A computational toolbox for the assembly yield of complex and heterogeneous structures.

The self-assembly of complex structures from a set of non-identical building blocks is a hallmark of soft matter and biological systems, including protein complexes, colloidal clusters, and DNA-based assemblies. Predicting the dependence of the equilibrium assembly yield on the concentrations and interaction energies of building blocks is highly challenging, owing to the difficulty of computing the entropic contributions to the free energy of the many structures that compete with the ground state configuration. While these calculations yield well known results for spherically symmetric building blocks, they do not hold when the building blocks have internal rotational degrees of freedom. Here we present an approach for solving this problem that works with arbitrary building blocks, including proteins with known structure and complex colloidal building blocks. Our algorithm combines classical statistical mechanics with recently developed computational tools for automatic differentiation. Automatic differentiation allows efficient evaluation of equilibrium averages over configurations that would otherwise be intractable. We demonstrate the validity of our framework by comparison to molecular dynamics simulations of simple examples, and apply it to calculate the yield curves for known protein complexes and for the assembly of colloidal shells.

Finite-size scaling at the jamming transition.

We present an analysis of finite-size effects in jammed packings of N soft, frictionless spheres at zero temperature. There is a 1/N correction to the discrete jump in the contact number at the transition so that jammed packings exist only above isostaticity. As a result, the canonical power-law scalings of the contact number and elastic moduli break down at low pressure. These quantities exhibit scaling collapse with a nontrivial scaling function, demonstrating that the jamming transition can be considered a phase transition. Scaling is achieved as a function of N in both two and three dimensions, indicating an upper critical dimension of 2.

Self-assembly-based posttranslational protein oscillators.

Recent advances in synthetic posttranslational protein circuits are substantially impacting the landscape of cellular engineering and offer several advantages compared to traditional gene circuits. However, engineering dynamic phenomena such as oscillations in protein-level circuits remains an outstanding challenge. Few examples of biological posttranslational oscillators are known, necessitating theoretical progress to determine realizable oscillators. We construct mathematical models for two posttranslational oscillators, using few components that interact only through reversible binding and phosphorylation/dephosphorylation reactions. Our designed oscillators rely on the self-assembly of two protein species into multimeric functional enzymes that respectively inhibit and enhance this self-assembly. We limit our analysis to within experimental constraints, finding (i) significant portions of the restricted parameter space yielding oscillations and (ii) that oscillation periods can be tuned by several orders of magnitude using recent advances in computational protein design. Our work paves the way for the rational design and realization of protein-based dynamic systems.

Enhanced diffusion by binding to the crosslinks of a polymer gel.

Creating a selective gel that filters particles based on their interactions is a major goal of nanotechnology, with far-reaching implications from drug delivery to controlling assembly pathways. However, this is particularly difficult when the particles are larger than the gel's characteristic mesh size because such particles cannot passively pass through the gel. Thus, filtering requires the interacting particles to transiently reorganize the gel's internal structure. While significant advances, e.g., in DNA engineering, have enabled the design of nano-materials with programmable interactions, it is not clear what physical principles such a designer gel could exploit to achieve selective permeability. We present an equilibrium mechanism where crosslink binding dynamics are affected by interacting particles such that particle diffusion is enhanced. In addition to revealing specific design rules for manufacturing selective gels, our results have the potential to explain the origin of selective permeability in certain biological materials, including the nuclear pore complex.

Single-molecule electrophoresis of beta-hairpin peptides by electrical recordings and Langevin dynamics simulations.

We used single-channel electrical recordings and Langevin molecular dynamics simulations to explore the electrophoretic translocation of various beta-hairpin peptides across the staphylococcal alpha-hemolysin (alphaHL) protein pore at single-molecule resolution. The beta-hairpin peptides, which varied in their folding properties, corresponded to the C terminal residues of the B1 domain of protein G. The translocation time was strongly dependent on the electric force and was correlated with the folding features of the beta-hairpin peptides. Highly unfolded peptides entered the pore in an extended conformation, resulting in fast single-file translocation events. In contrast, the translocation of the folded beta-hairpin peptides occurred more slowly. In this case, the beta-hairpin peptides traversed the alphaHL pore in a misfolded or fully folded conformation. This study demonstrates that the interaction between a polypeptide and a beta-barrel protein pore is dependent on the folding features of the polypeptide.

Erratum: Collective dynamics of soft active particles [Phys. Rev. E 91, 032706 (2015)].

This corrects the article DOI: 10.1103/PhysRevE.91.032706.

Collective dynamics of soft active particles.

We present a model of soft active particles that leads to a rich array of collective behavior found also in dense biological swarms of bacteria and other unicellular organisms. Our model uses only local interactions, such as Vicsek-type nearest-neighbor alignment, short-range repulsion, and a local boundary term. Changing the relative strength of these interactions leads to migrating swarms, rotating swarms, and jammed swarms, as well as swarms that exhibit run-and-tumble motion, alternating between migration and either rotating or jammed states. Interestingly, although a migrating swarm moves slower than an individual particle, the diffusion constant can be up to three orders of magnitude larger, suggesting that collective motion can be highly advantageous, for example, when searching for food.

Contact nonlinearities and linear response in jammed particulate packings.

Packings of frictionless athermal particles that interact only when they overlap experience a jamming transition as a function of packing density. Such packings provide the foundation for the theory of jamming. This theory rests on the observation that, despite the multitude of disordered configurations, the mechanical response to linear order depends only on the distance to the transition. We investigate the validity and utility of such measurements that invoke the harmonic approximation and show that, despite particles coming in and out of contact, there is a well-defined linear regime in the thermodynamic limit.

Phonon dispersion and elastic moduli of two-dimensional disordered colloidal packings of soft particles with frictional interactions.

Particle tracking and displacement covariance matrix techniques are employed to investigate the phonon dispersion relations of two-dimensional colloidal glasses composed of soft, thermoresponsive microgel particles whose temperature-sensitive size permits in situ variation of particle packing fraction. Bulk, B, and shear, G, moduli of the colloidal glasses are extracted from the dispersion relations as a function of packing fraction, and variation of the ratio G/B with packing fraction is found to agree quantitatively with predictions for jammed packings of frictional soft particles. In addition, G and B individually agree with numerical predictions for frictional particles. This remarkable level of agreement enabled us to extract an energy scale for the interparticle interaction from the individual elastic constants and to derive an approximate estimate for the interparticle friction coefficient.