Turning statistical physics models into materials design engines

Inverse design of a pyrochlore lattice of DNA origami through model-driven experiments

Sophisticated statistical mechanics approaches and human intuition have demonstrated the possibility of self-assembling complex lattices or finite-size constructs. However, attempts so far have mostly only been successful in silico and often fail in experiment because of unpredicted traps associated with kinetic slowing down (gelation, glass transition) and competing ordered structures. Theoretical predictions also face the difficulty of encoding the desired interparticle interaction potential with the experimentally available nano- and micrometer-sized particles. To overcome these issues, we combine SAT assembly (a patchy-particle interaction design algorithm based on constrained optimization) with coarse-grained simulations of DNA nanotechnology to experimentally realize trap-free self-assembly pathways. We use this approach to assemble a pyrochlore three-dimensional lattice, coveted for its promise in the construction of optical metamaterials, and characterize it with small-angle x-ray scattering and scanning electron microscopy visualization.

Surrogate Modeling of the Relative Entropy for Inverse Design Using Smolyak Sparse Grids

Relative entropy minimization, a statistical-mechanics approach for finding potential energy functions that produce target structural ensembles, has proven to be a powerful strategy for the inverse design of nanoparticle self-assembly. For a given target structure, the gradient of the relative entropy with respect to the adjustable parameters of the potential energy function is computed by performing a simulation, and then these parameters are updated using iterative gradient-based optimization. Small parameter updates per iteration and many iterations can be required for numerical stability, but this incurs considerable computational expense because a new simulation must be performed to reevaluate the gradient at each iteration. Here, we investigate the use of surrogate modeling to decouple the process of minimizing the relative entropy from the computationally demanding process of determining its gradient. We approximate the relative-entropy gradient using Chebyshev polynomial interpolation on Smolyak sparse grids. Our approach potentially increases the robustness and computational efficiency of using the relative entropy for inverse design, primarily for physically informed potential energy functions that have a small number of adjustable parameters.

Engineering Azeotropy to Optimize the Self-Assembly of Colloidal Mixtures

The goal of inverse self-assembly is to design interparticle interactions capable of assembling the units into a desired target structure. The effective assembly of complex structures often requires the use of multiple components, each new component increasing the thermodynamic degrees of freedom and, hence, the complexity of the self-assembly pathway. In this work we explore the possibility to use azeotropy, i.e., a special thermodynamic condition where the system behaves effectively as a one-component system, as a way to control the self-assembly of an arbitrary number of components. Exploiting the mass-balance equations, we show how to select patchy particle systems that exhibit azeotropic points along the desired self-assembly pathway. As an example we map the phase diagram of a binary mixture that, by design, fully assembles into cubic (and only cubic) diamond crystal via an azeotropic point. The ability to explicitly include azeotropic points in artificial designs reveals effective pathways for the self-assembly of complex structures.

Rationalizing Euclidean Assemblies of Hard Polyhedra from Tessellations in Curved Space

Entropic self-assembly is governed by the shape of the constituent particles, yet a priori prediction of crystal structures from particle shape alone is nontrivial for anything but the simplest of space-filling shapes. At the same time, most polyhedra are not space filling due to geometric constraints, but these constraints can be relaxed or even eliminated by sufficiently curving space. We show using Monte Carlo simulations that the majority of hard Platonic solids self-assemble entropically into space-filling crystals when constrained to the surface volume of a 3-sphere. As we gradually decrease curvature to "flatten" space and compare the local morphologies of crystals assembling in curved and flat space, we show that the Euclidean assemblies can be categorized as either remnants of tessellations in curved space (tetrahedra and dodecahedra) or nontessellation-based assemblies caused by large-scale geometric frustration (octahedra and icosahedra).

Nonequilibrium design strategies for functional colloidal assemblies

We use a nonequilibrium variational principle to optimize the steady-state, shear-induced interconversion of self-assembled nanoclusters of DNA-coated colloids. Employing this principle within a stochastic optimization algorithm allows us to identify design strategies for functional materials. We find that far-from-equilibrium shear flow can significantly enhance the flux between specific colloidal states by decoupling trade-offs between stability and reactivity required by systems in equilibrium. For isolated nanoclusters, we find nonequilibrium strategies for amplifying transition rates by coupling a given reaction coordinate to the background shear flow. We also find that shear flow can be made to selectively break detailed balance and maximize probability currents by coupling orientational degrees of freedom to conformational transitions. For a microphase consisting of many nanoclusters, we study the flux of colloids hopping between clusters. We find that a shear flow can amplify the flux without a proportional compromise on the microphase structure. This approach provides a general means of uncovering design principles for nanoscale, autonomous, functional materials driven far from equilibrium.

Application of colloidal photonic crystals in study of organoids

As alternative disease models, other than 2D cell lines and patient-derived xenografts, organoids have preferable in vivo physiological relevance. However, both endogenous and exogenous limitations impede the development and clinical translation of these organoids. Fortunately, colloidal photonic crystals (PCs), which benefit from favorable biocompatibility, brilliant optical manipulation, and facile chemical decoration, have been applied to the engineering of organoids and have achieved the desirable recapitulation of the ECM niche, well-defined geometrical onsets for initial culture, in situ multiphysiological parameter monitoring, single-cell biomechanical sensing, and high-throughput drug screening with versatile functional readouts. Herein, we review the latest progress in engineering organoids fabricated from colloidal PCs and provide inputs for future research.

Template Design for Complex Block Copolymer Patterns Using a Machine Learning Method

This study represents the first attempt to address the inverse design problem of the guiding template for directed self-assembly (DSA) patterns using solely machine learning methods. By formulating the problem as a multi-label classification task, the study shows that it is possible to predict templates without requiring any forward simulations. A series of neural network (NN) models, ranging from the basic two-layer convolutional neural network (CNN) to the large NN models (32-layer CNN with 8 residual blocks), have been trained using simulated pattern samples generated by thousands of self-consistent field theory (SCFT) calculations; a number of augmentation techniques, especially suitable for predicting morphologies, have been also proposed to enhance the performance of the NN model. The exact match accuracy of the model in predicting the template of simulated patterns was significantly improved from 59.8% for the baseline model to 97.1% for the best model of this study. The best model also demonstrates an excellent generalization ability in predicting the template for human-designed DSA patterns, while the simplest baseline model is ineffective in this task.

Optimization of non-equilibrium self-assembly protocols using Markov state models

The promise of self-assembly to enable the bottom-up formation of materials with prescribed architectures and functions has driven intensive efforts to uncover rational design principles for maximizing the yield of a target structure. Yet, despite many successful examples of self-assembly, ensuring kinetic accessibility of the target structure remains an unsolved problem in many systems. In particular, long-lived kinetic traps can result in assembly times that vastly exceed experimentally accessible timescales. One proposed solution is to design non-equilibrium assembly protocols in which system parameters change over time to avoid such kinetic traps. Here, we develop a framework to combine Markov state model (MSM) analysis with optimal control theory to compute a time-dependent protocol that maximizes the yield of the target structure at a finite time. We present an adjoint-based gradient descent method that, in conjunction with MSMs for a system as a function of its control parameters, enables efficiently optimizing the assembly protocol. We also describe an interpolation approach to significantly reduce the number of simulations required to construct the MSMs. We demonstrate our approach with two examples; a simple semi-analytic model for the folding of a polymer of colloidal particles, and a more complex model for capsid assembly. Our results show that optimizing time-dependent protocols can achieve significant improvements in the yields of selected structures, including equilibrium free energy minima, long-lived metastable structures, and transient states.

Hierarchical Materials from High Information Content Macromolecular Building Blocks: Construction, Dynamic Interventions, and Prediction

Hierarchical materials that exhibit order over multiple length scales are ubiquitous in nature. Because hierarchy gives rise to unique properties and functions, many have sought inspiration from nature when designing and fabricating hierarchical matter. More and more, however, nature's own high-information content building blocks, proteins, peptides, and peptidomimetics, are being coopted to build hierarchy because the information that determines structure, function, and interfacial interactions can be readily encoded in these versatile macromolecules. Here, we take stock of recent progress in the rational design and characterization of hierarchical materials produced from high-information content blocks with a focus on stimuli-responsive and "smart" architectures. We also review advances in the use of computational simulations and data-driven predictions to shed light on how the side chain chemistry and conformational flexibility of macromolecular blocks drive the emergence of order and the acquisition of hierarchy and also on how ionic, solvent, and surface effects influence the outcomes of assembly. Continued progress in the above areas will ultimately usher in an era where an understanding of designed interactions, surface effects, and solution conditions can be harnessed to achieve predictive materials synthesis across scale and drive emergent phenomena in the self-assembly and reconfiguration of high-information content building blocks.

Bayesian modeling of pattern formation from one snapshot of pattern

Partial differential equations (PDEs) have been widely used to reproduce patterns in nature and to give insight into the mechanism underlying pattern formation. Although many PDE models have been proposed, they rely on the pre-request knowledge of physical laws and symmetries, and developing a model to reproduce a given desired pattern remains difficult. We propose a method, referred to as Bayesian modeling of PDEs (BM-PDEs), to estimate the best dynamical PDE for one snapshot of a objective pattern under the stationary state without ground truth. We apply BM-PDEs to nontrivial patterns, such as quasicrystals (QCs), a double gyroid, and Frank-Kasper structures. We also generate three-dimensional dodecagonal QCs from a PDE model. This is done by using the estimated parameters for the Frank-Kasper A15 structure, which closely approximates the local structures of QCs. Our method works for noisy patterns and the pattern synthesized without the ground-truth parameters, which are required for the application toward experimental data.

SAT-assembly: a new approach for designing self-assembling systems

We propose a general framework for solving inverse self-assembly problems, i.e. designing interactions between elementary units such that they assemble spontaneously into a predetermined structure. Our approach uses patchy particles as building blocks, where the different units bind at specific interaction sites (the patches), and we exploit the possibility of having mixtures with several components. The interaction rules between the patches is determined by transforming the combinatorial problem into a Boolean satisfiability problem (SAT) which searches for solutions where all bonds are formed in the target structure. Additional conditions, such as the non-satisfiability of competing structures (e.g. metastable states) can be imposed, allowing to effectively design the assembly path in order to avoid kinetic traps. We demonstrate this approach by designing and numerically simulating a cubic diamond structure from four particle species that assembles without competition from other polymorphs, including the hexagonal structure.

Inverse design of soft materials via a deep learning-based evolutionary strategy

Colloidal self-assembly—the spontaneous organization of colloids into ordered structures—has been considered key to produce next-generation materials. However, the present-day staggering variety of colloidal building blocks and the limitless number of thermodynamic conditions make a systematic exploration intractable. The true challenge in this field is to turn this logic around and to develop a robust, versatile algorithm to inverse design colloids that self-assemble into a target structure. Here, we introduce a generic inverse design method to efficiently reverse-engineer crystals, quasicrystals, and liquid crystals by targeting their diffraction patterns. Our algorithm relies on the synergetic use of an evolutionary strategy for parameter optimization, and a convolutional neural network as an order parameter, and provides a way forward for the inverse design of experimentally feasible colloidal interactions, specifically optimized to stabilize the desired structure.

Deposition control of model glasses with surface-mediated orientational order

We introduce a minimal model of solid-forming anisotropic molecules that displays, in thermal equilibrium, surface orientational order without bulk orientational order. The model reproduces the nonequilibrium behavior of recent experiments in which a bulk nonequilibrium structure grown by deposition contains regions of orientational order characteristic of the surface equilibrium. This order is deposited, in general, in a nonuniform way because of the emergence of a growth-poisoning mechanism that causes equilibrated surfaces to grow slower than non-equilibrated surfaces. We use evolutionary methods to design oscillatory protocols able to grow nonequilibrium structures with uniform order, demonstrating the potential of protocol design for the fabrication of this class of materials.

Open-source platform for block polymer formulation design using particle swarm optimization

Facile exploration of large design spaces is critical to the development of new functional soft materials, including self-assembling block polymers, and computational inverse design methodologies are a promising route to initialize this task. We present here an open-source software package coupling particle swarm optimization (PSO) with an existing open-source self-consistent field theory (SCFT) software for the inverse design of self-assembling block polymers to target bulk morphologies. To lower the barrier to use of the software and facilitate exploration of novel design spaces, the underlying SCFT calculations are seeded with algorithmically generated initial fields for four typical morphologies: lamellae, network phases, cylindrical phases, and spherical phases. In addition to its utility within PSO, the initial guess tool also finds generic applicability for stand-alone SCFT calculations. The robustness of the software is demonstrated with two searches for classical phases in the conformationally symmetric diblock system, as well as one search for the Frank-Kasper [Formula: see text] phase in conformationally asymmetric diblocks. The source code for both the initial guess generation and the PSO wrapper is publicly available.

Thermodynamic stability versus kinetic accessibility: Pareto fronts for programmable self-assembly

A challenge in designing self-assembling building blocks is to ensure the target state is both thermodynamically stable and kinetically accessible. These two objectives are known to be typically in competition, but it is not known how to simultaneously optimize them. We consider this problem through the lens of multi-objective optimization theory: we develop a genetic algorithm to compute the Pareto fronts characterizing the tradeoff between equilibrium probability and folding rate, for a model system of small polymers of colloids with tunable short-ranged interaction energies. We use a coarse-grained model for the particles' dynamics that allows us to efficiently search over parameters, for systems small enough to be enumerated. For most target states there is a tradeoff when the number of types of particles is small, with medium-weak bonds favouring fast folding, and strong bonds favouring high equilibrium probability. The tradeoff disappears when the number of particle types reaches a value m*, that is usually much less than the total number of particles. This general approach of computing Pareto fronts allows one to identify the minimum number of design parameters to avoid a thermodynamic-kinetic tradeoff. However, we argue, by contrasting our coarse-grained model's predictions with those of Brownian dynamics simulations, that particles with short-ranged isotropic interactions should generically have a tradeoff, and avoiding it in larger systems will require orientation-dependent interactions.

Neuroevolutionary Learning of Particles and Protocols for Self-Assembly

Within simulations of molecules deposited on a surface we show that neuroevolutionary learning can design particles and time-dependent protocols to promote self-assembly, without input from physical concepts such as thermal equilibrium or mechanical stability and without prior knowledge of candidate or competing structures. The learning algorithm is capable of both directed and exploratory design: it can assemble a material with a user-defined property, or search for novelty in the space of specified order parameters. In the latter mode it explores the space of what can be made, rather than the space of structures that are low in energy but not necessarily kinetically accessible.

From predictive modelling to machine learning and reverse engineering of colloidal self-assembly

An overwhelming diversity of colloidal building blocks with distinct sizes, materials and tunable interaction potentials are now available for colloidal self-assembly. The application space for materials composed of these building blocks is vast. To make progress in the rational design of new self-assembled materials, it is desirable to guide the experimental synthesis efforts by computational modelling. Here, we discuss computer simulation methods and strategies used for the design of soft materials created through bottom-up self-assembly of colloids and nanoparticles. We describe simulation techniques for investigating the self-assembly behaviour of colloidal suspensions, including crystal structure prediction methods, phase diagram calculations and enhanced sampling techniques, as well as their limitations. We also discuss the recent surge of interest in machine learning and reverse-engineering methods. Although their implementation in the colloidal realm is still in its infancy, we anticipate that these data-science tools offer new paradigms in understanding, predicting and (inverse) design of novel colloidal materials.

Inverse design of equilibrium cluster fluids applied to a physically informed model

Inverse design strategies have proven highly useful for the discovery of interaction potentials that prompt self-assembly of a variety of interesting structures. However, often the optimized particle interactions do not have a direct relationship to experimental systems. In this work, we show that Relative Entropy minimization is able to discover physically meaningful parameter sets for a model interaction built from depletion attraction and electrostatic repulsion that yield self-assembly of size-specific clusters. We then explore the sensitivity of the optimized interaction potentials with respect to deviations in the underlying physical quantities, showing that clustering behavior is largely preserved even as the optimized parameters are perturbed.

Inverse design of nonequilibrium steady states: A large-deviation approach

The design of small scale nonequilibrium steady states (NESS) is a challenging, open ended question. While similar equilibrium problems are tractable using standard thermodynamics, a generalized description for nonequilibrium systems is lacking, making the design problem particularly difficult. Here we show we can exploit the large-deviation behavior of a Brownian particle and design a variety of geometrically complex steady-state density distributions and flux field flows. We achieve this design target from direct knowledge of the joint large-deviation functional for the empirical density and flow, and a "relaxation" algorithm on the desired target states via adjustable force field parameters. We validate the method by replicating analytical results, and demonstrate its capacity to yield complex prescribed targets, such as rose-curve or polygonal shapes on the plane. We consider this dynamical fluctuation approach a first step towards the design of more complex NESS where general frameworks are otherwise still lacking.

Machine Learning Predictions of Block Copolymer Self-Assembly

Directed self-assembly of block copolymers is a key enabler for nanofabrication of devices with sub-10 nm feature sizes, allowing patterning far below the resolution limit of conventional photolithography. Among all the process steps involved in block copolymer self-assembly, solvent annealing plays a dominant role in determining the film morphology and pattern quality, yet the interplay of the multiple parameters during solvent annealing, including the initial thickness, swelling, time, and solvent ratio, makes it difficult to predict and control the resultant self-assembled pattern. Here, machine learning tools are applied to analyze the solvent annealing process and predict the effect of process parameters on morphology and defectivity. Two neural networks are constructed and trained, yielding accurate prediction of the final morphology in agreement with experimental data. A ridge regression model is constructed to identify the critical parameters that determine the quality of line/space patterns. These results illustrate the potential of machine learning to inform nanomanufacturing processes.

Simulating sticky particles: A Monte Carlo method to sample a stratification

Many problems in materials science and biology involve particles interacting with strong, short-ranged bonds that can break and form on experimental timescales. Treating such bonds as constraints can significantly speed up sampling their equilibrium distribution, and there are several methods to sample probability distributions subject to fixed constraints. We introduce a Monte Carlo method to handle the case when constraints can break and form. More generally, the method samples a probability distribution on a stratification: a collection of manifolds of different dimensions, where the lower-dimensional manifolds lie on the boundaries of the higher-dimensional manifolds. We show several applications of the method in polymer physics, self-assembly of colloids, and volume calculation in high dimensions.

Designing Patchy Interactions to Self-Assemble Arbitrary Structures

One of the fundamental goals of nanotechnology is to exploit selective and directional interactions between molecules to design particles that self-assemble into desired structures, from capsids, to nanoclusters, to fully formed crystals with target properties (e.g., optical, mechanical, etc.). Here, we provide a general framework which transforms the inverse problem of self-assembly of colloidal crystals into a Boolean satisfiability problem for which solutions can be found numerically. Given a reference structure and the desired number of components, our approach produces designs for which the target structure is an energy minimum, and also allows us to exclude solutions that correspond to competing structures. We demonstrate the effectiveness of our approach by designing model particles that spontaneously nucleate milestone structures such as the cubic diamond, the pyrochlore, and the clathrate lattices.

Self-assembly and entropic effects in pear-shaped colloid systems. II. Depletion attraction of pear-shaped particles in a hard-sphere solvent

We consider depletion effects of a pear-shaped colloidal particle in a hard-sphere solvent for two different model realizations of the pear-shaped colloidal particle. The two models are the pear hard Gaussian overlap (PHGO) particles and the hard pears of revolution (HPR). The motivation for this study is to provide a microscopic understanding for the substantially different mesoscopic self-assembly properties of these pear-shaped colloids, in dense suspensions, that have been reported in the previous studies. This is done by determining their differing depletion attractions via Monte Carlo simulations of PHGO and HPR particles in a pool of hard spheres and comparing them with excluded volume calculations of numerically obtained ideal configurations on the microscopic level. While the HPR model behaves as predicted by the analysis of excluded volumes, the PHGO model showcases a preference for splay between neighboring particles, which can be attributed to the special non-additive characteristics of the PHGO contact function. Lastly, we propose a potentially experimentally realizable pear-shaped particle model, the non-additive hard pear of revolution model, which is based on the HPR model but also features non-additive traits similar to those of PHGO particles to mimic their depletion behavior.

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.

Open-source code for self-consistent field theory calculations of block polymer phase behavior on graphics processing units

Self-consistent field theory (SCFT) is a powerful approach for computing the phase behavior of block polymers. We describe a fast version of the open-source Polymer Self-Consistent Field (PSCF) code that takes advantage of the massive parallelization provided by a graphical processing unit (GPU). Benchmarking double-precision calculations indicate up to 30× reduction in time to converge SCFT calculations of various diblock copolymer phases when compared to the Fortran CPU version of PSCF using the same algorithms, with the speed-up increasing with increasing unit cell size for the diblock polymer problems examined here. Where double-precision accuracy is not needed, single-precision calculations can provide speed-up of up to 60× in convergence time. These improvements in speed within an open-source format open up new vistas for SCFT-driven block polymer materials discovery by the community at large.

Customizing topographical templates for aperiodic nanostructures of block copolymers via inverse design

The limited complexity of self-assembled nanostructures of block copolymers seriously impedes their potential utility in the semiconductor industry. Therefore, the customizability of complex nanostructures has been a long-standing goal for the utilization of directed self-assembly in nanolithography. Herein, we integrated an advanced inverse design algorithm with a well-developed theoretical model to deduce inverse solutions of topographical templates to direct the self-assembly of block copolymers into reproducible target structures. The deduced templates were optimized by finely tuning the input parameters of the inverse design algorithm and through symmetric operation as well as nanopost elimination. More importantly, our developed algorithm has the capability to search inverse solutions of topographical templates for aperiodic nanostructures over exceptionally large areas. These results reveal design rules for guiding templates for the device-oriented nanostructures of block copolymers with prospective applications in nanolithography.

Nonequilibrium phase diagrams for actomyosin networks

Living cells dynamically modulate the local morphologies of their actin networks to perform biological functions, including force transduction, intracellular transport, and cell division. A major challenge is to understand how diverse structures of the actin cytoskeleton are assembled from a limited set of molecular building blocks. Here we study the spontaneous self-assembly of a minimal model of cytoskeletal materials, consisting of semiflexible actin filaments, crosslinkers, and molecular motors. Using coarse-grained simulations, we demonstrate that by changing concentrations and kinetics of crosslinkers and motors, as well as filament lengths, we can generate three distinct structural phases of actomyosin assemblies: bundled, polarity-sorted, and contracted. We introduce new metrics to distinguish these structural phases and demonstrate their functional roles. We find that the binding kinetics of motors and crosslinkers can be tuned to optimize contractile force generation, motor transport, and mechanical response. By quantitatively characterizing the relationships between the modes of cytoskeletal self-assembly, the resulting structures, and their functional consequences, our work suggests new principles for the design of active materials.

Optimizing self-consistent field theory block copolymer models with X-ray metrology

A block copolymer self-consistent field theory (SCFT) model is used for direct analysis of experimental X-ray scattering data obtained from thin films of polystyrene-b-poly(methyl methacrylate) (PS-b-PMMA) made from directed self-assembly. In a departure from traditional approaches, which reconstruct the real space structure using simple geometric shapes, we build on recent work that has relied on physics-based models to determine shape profiles and extract thermodynamic processing information from the scattering data. More specifically, an SCFT model, coupled to a covariance matrix adaptation evolutionary strategy (CMAES), is used to find the set of simulation parameters for the model that best reproduces the scattering data. The SCFT model is detailed enough to capture the essential physics of the copolymer self-assembly, but sufficiently simple to rapidly produce structure profiles needed for interpreting the scattering data. The ability of the model to produce a matching scattering profile is assessed, and several improvements are proposed in order to more accurately recreate the experimental observations. The predicted parameters are compared to those extracted from model fits via additional experimental methods and with predicted parameters from direct particle-based simulations of the same model, which incorporate the effects of fluctuations. The Flory-Huggins interaction parameter for PS-b-PMMA is found to be in agreement with reported ranges for this material. These results serve to strengthen the case for relying on physics-based models for direct analysis of scattering and light signal based experiments.

Machine learning and data science in soft materials engineering

In many branches of materials science it is now routine to generate data sets of such large size and dimensionality that conventional methods of analysis fail. Paradigms and tools from data science and machine learning can provide scalable approaches to identify and extract trends and patterns within voluminous data sets, perform guided traversals of high-dimensional phase spaces, and furnish data-driven strategies for inverse materials design. This topical review provides an accessible introduction to machine learning tools in the context of soft and biological materials by 'de-jargonizing' data science terminology, presenting a taxonomy of machine learning techniques, and surveying the mathematical underpinnings and software implementations of popular tools, including principal component analysis, independent component analysis, diffusion maps, support vector machines, and relative entropy. We present illustrative examples of machine learning applications in soft matter, including inverse design of self-assembling materials, nonlinear learning of protein folding landscapes, high-throughput antimicrobial peptide design, and data-driven materials design engines. We close with an outlook on the challenges and opportunities for the field.

Size scaling of failure strength with fat-tailed disorder in a fiber bundle model

We investigate the size scaling of the macroscopic fracture strength of heterogeneous materials when microscopic disorder is controlled by fat-tailed distributions. We consider a fiber bundle model where the strength of single fibers is described by a power law distribution over a finite range. Tuning the amount of disorder by varying the power law exponent and the upper cutoff of fibers' strength, in the limit of equal load sharing an astonishing size effect is revealed: For small system sizes the bundle strength increases with the number of fibers, and the usual decreasing size effect of heterogeneous materials is restored only beyond a characteristic size. We show analytically that the extreme order statistics of fibers' strength is responsible for this peculiar behavior. Analyzing the results of computer simulations we deduce a scaling form which describes the dependence of the macroscopic strength of fiber bundles on the parameters of microscopic disorder over the entire range of system sizes.

Optimising minimal building blocks for addressable self-assembly

Addressable structures are characterised by the set of unique components from which they are built and by the specific location that each component occupies. For an addressable structure to self-assemble, its constituent building blocks must be encoded with sufficient information to define their positions with respect to each other and to enable them to navigate to those positions. DNA, with its vast scope for encoding specific interactions, has been successfully used to synthesise addressable systems of several hundred components. In this work we examine the complementary question of the minimal requirements for building blocks to undergo addressable self-assembly driven by a controlled temperature quench. Our testbed is an idealised model of cubic particles patterned with attractive interactions. We introduce a scheme for optimising the interactions using a variant of basin-hopping and a negative design principle. The designed building blocks are tested dynamically in simple target structures to establish how their complexity affects the limits of reliable self-assembly.

Swarm Intelligence Platform for Multiblock Polymer Inverse Formulation Design

Multiblock polymers (MBPs) show great promise as a platform for synthesizing materials with highly customized nanostructures. In these systems, the subtle balance of chain configurational entropy and interactions between dissimilar block chemistries emerge a rich palette of self-assembled structural motifs that may be leveraged to impart novel functional, mechanical, or optical properties to the final material. Unfortunately, extensive study of MBP self-assembly has yielded few reliable heuristics for intuitively navigating the enormous molecular design space made available by modern polymer synthesis techniques. In order to make progress, new methods must be developed that allow researchers to efficiently screen molecular designs for achieving a desired state of self-assembly. Here, we introduce a platform for the automated discovery of tailored MBP formulations based on the Particle Swarm Optimization (PSO) method and a linear multiblock chain parametrization that enables continuous optimization of chain architecture. We apply the method to thin-film blends of linear ABC triblock polymers subject to lateral confinement, allowing all polymer and blend parameters to freely optimize in search of a prespecified, nontrivial target morphology. While we focus on pattern selection as a proof of principle, any computable equilibrium property can be optimized using the methods described here.

Digital Alchemy for Materials Design: Colloids and Beyond

Starting with the early alchemists, a holy grail of science has been to make desired materials by modifying the attributes of basic building blocks. Building blocks that show promise for assembling new complex materials can be synthesized at the nanoscale with attributes that would astonish the ancient alchemists in their versatility. However, this versatility means that making a direct connection between building-block attributes and bulk structure is both necessary for rationally engineering materials and difficult because building block attributes can be altered in many ways. Here we show how to exploit the malleability of the valence of colloidal nanoparticle "elements" to directly and quantitatively link building-block attributes to bulk structure through a statistical thermodynamic framework we term "digital alchemy". We use this framework to optimize building blocks for a given target structure and to determine which building-block attributes are most important to control for self-assembly, through a set of novel thermodynamic response functions, moduli, and susceptibilities. We thereby establish direct links between the attributes of colloidal building blocks and the bulk structures they form. Moreover, our results give concrete solutions to the more general conceptual challenge of optimizing emergent behaviors in nature and can be applied to other types of matter. As examples, we apply digital alchemy to systems of truncated tetrahedra, rhombic dodecahedra, and isotropically interacting spheres that self-assemble diamond, fcc, and icosahedral quasicrystal structures, respectively. Although our focus is on colloidal systems, our methods generalize to any building blocks with adjustable interactions.