The role of pressure in inverse design for assembly

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.

Molecularly Informed Field Theories from Bottom-up Coarse-Graining

Polymer formulations possessing mesostructures or phase coexistence are challenging to simulate using atomistic particle-explicit approaches due to the disparate time and length scales, while the predictive capability of field-based simulations is hampered by the need to specify interactions at a coarser scale (e.g., χ-parameters). To overcome the weaknesses of both, we introduce a bottom-up coarse-graining methodology that leverages all-atom molecular dynamics to molecularly inform coarser field-theoretic models. Specifically, we use relative-entropy coarse-graining to parametrize particle models that are directly and analytically transformable into statistical field theories. We demonstrate the predictive capability of this approach by reproducing experimental aqueous poly(ethylene oxide) (PEO) cloud-point curves with no parameters fit to experimental data. This synergistic approach to multiscale polymer simulations opens the door to de novo exploration of phase behavior across a wide variety of polymer solutions and melt formulations.

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.

BoltzmaNN: Predicting effective pair potentials and equations of state using neural networks

Neural networks (NNs) are employed to predict equations of state from a given isotropic pair potential using the virial expansion of the pressure. The NNs are trained with data from molecular dynamics simulations of monoatomic gases and liquids, sampled in the NVT ensemble at various densities. We find that the NNs provide much more accurate results compared to the analytic low-density limit estimate of the second virial coefficient and the Carnahan-Starling equation of state for hard sphere liquids. Furthermore, we design and train NNs for computing (effective) pair potentials from radial pair distribution functions, g(r), a task that is often performed for inverse design and coarse-graining. Providing the NNs with additional information on the forces greatly improves the accuracy of the predictions since more correlations are taken into account; the predicted potentials become smoother, are significantly closer to the target potentials, and are more transferable as a result.

Surveying the free energy landscape of clusters of attractive colloidal spheres

Controlling the assembly of colloidal particles into specific structures has been a long-term goal of the soft materials community. Much can be learned about the process of self-assembly by examining the early stage assembly into clusters. For the simple case of hard spheres with short-range attractions, the rigid clusters of N particles (where N is small) have been enumerated theoretically and tested experimentally. Less is known, however, about how the free energy landscapes are altered when the inter-particle potential is long-ranged. In this work, we demonstrate how adaptive biasing in molecular simulations may be used to pinpoint shifts in the stability of colloidal clusters as the inter-particle potential is varied. We also discuss the generality of our techniques and strategies for application to related molecular systems.

Grand canonical inverse design of multicomponent colloidal crystals

Inverse design methods are powerful computational approaches for creating colloidal systems which self-assemble into a target morphology by reverse engineering the Hamiltonian of the system. Despite this, these optimization procedures tend to yield Hamiltonians which are too complex to be experimentally realized. An alternative route to complex structures involves the use of several different components, however, conventional inverse design methods do not explicitly account for the possibility of phase separation into compositionally distinct structures. Here, we present an inverse design scheme for multicomponent colloidal systems by combining active learning with a method to directly compute their ground state phase diagrams. This explicitly accounts for phase separation and can locate stable regions of Hamiltonian parameter space which grid-based surveys are prone to miss. Using this we design low-density, binary structures with Lennard-Jones-like pairwise interactions that are simpler than in the single component case and potentially realizable in an experimental setting. This reinforces the concept that ground states of simple, multicomponent systems might be rich with previously unappreciated diversity, enabling the assembly of non-trivial structures with only few simple components instead of a single complex one.