Equilibrium mechanisms of self-limiting assembly

2D capsid formation within an oscillatory energy landscape: orderly self-assembly depends on the interplay between a dynamic potential and intrinsic relaxation times

Multiple dissipative self-assembly protocols designed to create novel structures or to reduce kinetic traps have recently emerged. Specifically, temporal oscillations of particle interactions have been shown effective at both aims, but investigations thus far have focused on systems of simple colloids or their binary mixtures. In this work, we expand our understanding of the effect of temporally oscillating interactions to a two-dimensional coarse-grained viral capsid-like model that undergoes a self-limited assembly. This model includes multiple intrinsic relaxation times due to the internal structure of the capsid subunits and, under certain interaction regimes, proceeds via a two-step nucleation mechanism. We find that oscillations much faster than the local intrinsic relaxation times can be described via a time averaged inter-particle potential across a wide range of interaction strengths, while oscillations much slower than these relaxation times result in structures that adapt to the attraction strength of the current half-cycle. Interestingly, oscillation periods similar to these relaxation times shift the interaction window over which orderly assembly occurs by enabling error correction during the half-cycles with weaker attractions. Our results provide fundamental insights to non-equilibrium self-assembly on temporally variant energy landscapes.

Programmable Crowding and Tunable Phases in a Binary Mixture of Colloidal Particles under Light-Driven Thermal Convection

We employ photothermally driven self-assembly of colloidal particles to design microscopic structures with programmable size and tunable order. The experimental system is based on a binary mixture of "plasmonic heater" gold nanoparticles and "assembly building block" microparticles. Photothermal heating of the gold nanoparticles under visible light causes a natural convection flow that efficiently assembles the microscale building block particles (diameter 1-10 μm) into a monolayer. We identify the onset of active Brownian motion of colloidal particles under this convective flow by varying the conditions of light intensity, gold nanoparticle concentration, and sample height. We realize a crowded assembly of microparticles around the center of illumination and show that the size of the particle crowd can be programmed using patterned light illumination. In a binary mixture of gold nanoparticles and polystyrene microparticles, we demonstrate the formation of rapid and large-scale crystalline monolayers, covering an area of 0.88 mm2 within 10 min. We find that the structural order of the assembly can be tuned by varying the surface charge of the nanoparticles and the size of the microparticles, giving rise to the formation of different phases-colloidal crystals, crowds, and gels. Using Monte Carlo simulations, we explain how the phases emerge from the interplay between hydrodynamic and electrostatic interactions, as well as the assembly kinetics. Our study demonstrates the promise of self-assembly with programmable shapes and structural order under nonequilibrium conditions using an accessible setup comprising only binary mixtures and LED light.

Economical routes to size-specific assembly of self-closing structures

Programmable self-assembly has seen an explosion in the diversity of synthetic crystalline materials, but developing strategies that target "self-limiting" assemblies has remained a challenge. Among these, self-closing structures, in which the local curvature defines the finite global size, are prone to polymorphism due to thermal bending fluctuations, a problem that worsens with increasing target size. Here, we show that assembly complexity can be used to eliminate this source of polymorphism in the assembly of tubules. Using many distinct components, we prune the local density of off-target geometries, increasing the selectivity of the tubule width and helicity to nearly 100%. We further show that by reducing the design constraints to target either the pitch or the width alone, fewer components are needed to reach complete selectivity. Combining experiments with theory, we reveal an economical limit, which determines the minimum number of components required to create arbitrary assembly sizes with full selectivity.

Programming patchy particles for materials assembly design

Direct design of complex functional materials would revolutionize technologies ranging from printable organs to novel clean energy devices. However, even incremental steps toward designing functional materials have proven challenging. If the material is constructed from highly complex components, the design space of materials properties rapidly becomes too computationally expensive to search. On the other hand, very simple components such as uniform spherical particles are not powerful enough to capture rich functional behavior. Here, we introduce a differentiable materials design model with components that are simple enough to design yet powerful enough to capture complex materials properties: rigid bodies composed of spherical particles with directional interactions (patchy particles). We showcase the method with self-assembly designs ranging from open lattices to self-limiting clusters, all of which are notoriously challenging design goals to achieve using purely isotropic particles. By directly optimizing over the location and interaction of the patches on patchy particles using gradient descent, we dramatically reduce the computation time for finding the optimal building blocks.

Nanocrystal Assemblies: Current Advances and Open Problems

We explore the potential of nanocrystals (a term used equivalently to nanoparticles) as building blocks for nanomaterials, and the current advances and open challenges for fundamental science developments and applications. Nanocrystal assemblies are inherently multiscale, and the generation of revolutionary material properties requires a precise understanding of the relationship between structure and function, the former being determined by classical effects and the latter often by quantum effects. With an emphasis on theory and computation, we discuss challenges that hamper current assembly strategies and to what extent nanocrystal assemblies represent thermodynamic equilibrium or kinetically trapped metastable states. We also examine dynamic effects and optimization of assembly protocols. Finally, we discuss promising material functions and examples of their realization with nanocrystal assemblies.

Limits of economy and fidelity for programmable assembly of size-controlled triply periodic polyhedra

We propose and investigate an extension of the Caspar-Klug symmetry principles for viral capsid assembly to the programmable assembly of size-controlled triply periodic polyhedra, discrete variants of the Primitive, Diamond, and Gyroid cubic minimal surfaces. Inspired by a recent class of programmable DNA origami colloids, we demonstrate that the economy of design in these crystalline assemblies-in terms of the growth of the number of distinct particle species required with the increased size-scale (e.g., periodicity)-is comparable to viral shells. We further test the role of geometric specificity in these assemblies via dynamical assembly simulations, which show that conditions for simultaneously efficient and high-fidelity assembly require an intermediate degree of flexibility of local angles and lengths in programmed assembly. Off-target misassembly occurs via incorporation of a variant of disclination defects, generalized to the case of hyperbolic crystals. The possibility of these topological defects is a direct consequence of the very same symmetry principles that underlie the economical design, exposing a basic tradeoff between design economy and fidelity of programmable, size controlled assembly.

Narcissistic self-sorting in Zn(II) porphyrin derived semiconducting nanostructures

The narcissistic self-sorted phenomenon is explicitly attributed to the structural similarities in organic molecules. Although such relevant materials are rarely explored, self-sorted structures from macrocyclic π-conjugated-based p- and n-type organic semiconductors facilitate the increase of exciton dissociation and charge separation in bulk heterojunction solar cells. Herein, we report two extended π-conjugated derivatives consisting of zinc-porphyrin-linked benzothiadiazole acting as an acceptor (PB) and anthracene as a donor (PA). Despite having the same porphyrin π-conjugated core in PA and PB, variations in donor and acceptor moieties make the molecular packing form one-dimensional (1D) self-assembled nanofibers via H- and J-type aggregates. Interestingly, a dissimilar aggregate of PA and PB exists as a mixture (PA + PB), promoting narcissistic self-sorted structures. Electrochemical impedance investigation reveals that the electronic characteristics of self-sorting assemblies are influenced by the difference in electrostatic potentials for PA and PB, resulting in a transitional electrical conductivity of 0.14 S cm-1. Therefore, the design of such materials for the fabrication of effective photovoltaics is promoted by these extraordinary self-sorted behaviors in comparable organic π-conjugated molecules.

Electro-Sorting Create Heterogeneity: Constructing A Multifunctional Janus Film with Integrated Compositional and Microstructural Gradients for Guided Bone Regeneration

Biology remains the envy of flexible soft matter fabrication because it can satisfy multiple functional needs by organizing a small set of proteins and polysaccharides into hierarchical systems with controlled heterogeneity in composition and microstructure. Here, it is reported that controlled, mild electronic inputs (<10 V; <20 min) induce a homogeneous gelatin-chitosan mixture to undergo sorting and bottom-up self-assembly into a Janus film with compositional gradient (i.e., from chitosan-enriched layer to chitosan/gelatin-contained layer) and tunable dense-porous gradient microstructures (e.g., porosity, pore size, and ratio of dense to porous layers). This Janus film performs is shown multiple functions for guided bone regeneration: the integration of compositional and microstructural features confers flexible mechanics, asymmetric properties for interfacial wettability, molecular transport (directional growth factor release), and cellular responses (prevents fibroblast infiltration but promotes osteoblast growth and differentiation). Overall, this work demonstrates the versatility of electrofabrication for the customized manufacturing of functional gradient soft matter.

Hierarchical assembly is more robust than egalitarian assembly in synthetic capsids

Self-assembly of complex and functional materials remains a grand challenge in soft material science. Efficient assembly depends on a delicate balance between thermodynamic and kinetic effects, requiring fine-tuning affinities and concentrations of subunits. By contrast, we introduce an assembly paradigm that allows large error-tolerance in the subunit affinity and helps avoid kinetic traps. Our combined experimental and computational approach uses a model system of triangular subunits programmed to assemble into T = 3 icosahedral capsids comprising 60 units. The experimental platform uses DNA origami to create monodisperse colloids whose three-dimensional geometry is controlled to nanometer precision, with two distinct bonds whose affinities are controlled to kBT precision, quantified in situ by static light scattering. The computational model uses a coarse-grained representation of subunits, short-ranged potentials, and Langevin dynamics. Experimental observations and modeling reveal that when the bond affinities are unequal, two distinct hierarchical assembly pathways occur, in which the subunits first form dimers in one case and pentamers in another. These hierarchical pathways produce complete capsids faster and are more robust against affinity variation than egalitarian pathways, in which all binding sites have equal strengths. This finding suggests that hierarchical assembly may be a general engineering principle for optimizing self-assembly of complex target structures.

Supersolid Phases of Bosonic Particles in a Bubble Trap

Confinement can have a considerable effect on the behavior of particle systems and is therefore an effective way to discover new phenomena. A notable example is a system of identical bosons at low temperature under an external field mimicking an isotropic bubble trap, which constrains the particles to a portion of space close to a spherical surface. Using path integral Monte Carlo simulations, we examine the spatial structure and superfluid fraction in two emblematic cases. First, we look at soft-core bosons, finding the existence of supersolid cluster arrangements with polyhedral symmetry; we show how different numbers of clusters are stabilized depending on the trap radius and the particle mass, and we characterize the temperature behavior of the cluster phases. A detailed comparison with the behavior of classical soft-core particles is provided too. Then, we examine the case, of more immediate experimental interest, of a dipolar condensate on the sphere, demonstrating how a quasi-one-dimensional supersolid of clusters is formed on a great circle for realistic values of density and interaction parameters. Crucially, this supersolid phase is only slightly disturbed by gravity. We argue that the predicted phases can be revealed in magnetic traps with spherical-shell geometry, possibly even in a lab on Earth. Our results pave the way for future simulation studies of correlated quantum systems in curved geometries.

Thermal stability and secondary aggregation of self-limiting, geometrically frustrated assemblies: Chain assembly of incommensurate polybricks

In geometrically frustrated assemblies, equilibrium self-limitation manifests in the form of a minimum in the free energy per subunit at a finite, multisubunit size which results from the competition between the elastic costs of frustration within an assembly and the surface energy at its boundaries. Physical realizations-from ill-fitting particle assemblies to self-twisting protein superstructures-are capable of multiple mechanisms of escaping the cumulative costs of frustration, resulting in unlimited equilibrium assembly, including elastic modes of "shape flattening" and the formation of weak, defective bonds that screen intra-assembly stresses. Here we study a model of one-dimensional chain assembly of incommensurate "polybricks" and determine its equilibrium assembly as a function of temperature, concentration, degree of shape frustration, elasticity, and interparticle binding, notably focusing on how weakly cohesive, defective bonds give rise to strongly temperature-dependent assembly. Complex assembly behavior derives from the competition between multiple distinct local minima in the free-energy landscape, including self-limiting chains, weakly bound aggregates of self-limiting chains, and strongly bound, elastically defrustrated assemblies. We show that this scenario, in general, gives rise to anomalous multiple aggregation behavior, in which disperse subunits (stable at low concentration and high temperature) first exhibit a primary aggregation transition to self-limiting chains (at intermediate concentration and temperature) which are ultimately unstable to condensation into unlimited assembly of finite-chains through weak binding beyond a secondary aggregation transition (at low temperature and high concentration). We show that window of stable self-limitation is determined both by the elastic costs of frustration in the assembly as well as energetic and entropic features of intersubunit binding.

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).

Geometric frustration of hard-disk packings on cones

Conical surfaces pose an interesting challenge to crystal growth: A crystal growing on a cone can wrap around and meet itself at different radii. We use a disk-packing algorithm to investigate how this closure constraint can geometrically frustrate the growth of single crystals on cones with small opening angles. By varying the crystal seed orientation and cone angle, we find that-except at special commensurate cone angles-crystals typically form a seam that runs along the axial direction of the cone, while near the tip, a disordered particle packing forms. We show that the onset of disorder results from a finite-size effect that depends strongly on the circumference and not on the seed orientation or cone angle. This finite-size effect occurs also on cylinders, and we present evidence that on both cylinders and cones, the defect density increases exponentially as circumference decreases. We introduce a simple model for particle attachment at the seam that explains the dependence on the circumference. Our findings suggest that the growth of single crystals can become frustrated even very far from the tip when the cone has a small opening angle. These results may provide insights into the observed geometry of conical crystals in biological and materials applications.

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.

Self-assembly coupled to liquid-liquid phase separation

Liquid condensate droplets with distinct compositions of proteins and nucleic acids are widespread in biological cells. While it is known that such droplets, or compartments, can regulate irreversible protein aggregation, their effect on reversible self-assembly remains largely unexplored. In this article, we use kinetic theory and solution thermodynamics to investigate the effect of liquid-liquid phase separation on the reversible self-assembly of structures with well-defined sizes and architectures. We find that, when assembling subunits preferentially partition into liquid compartments, robustness against kinetic traps and maximum achievable assembly rates can be significantly increased. In particular, both the range of solution conditions leading to productive assembly and the corresponding assembly rates can increase by orders of magnitude. We analyze the rate equation predictions using simple scaling estimates to identify effects of liquid-liquid phase separation as a function of relevant control parameters. These results may elucidate self-assembly processes that underlie normal cellular functions or pathogenesis, and suggest strategies for designing efficient bottom-up assembly for nanomaterials applications.

Isotropic and Anisotropic Monolayer Structures in RF Discharge Plasma

We present the results of an experimental and analytical study of the structural and dynamic properties of a monolayer consisting of dust grains in an electrostatic trap in an RF discharge plasma. The possibility of forming a monolayer with an isotropic distribution for interparticle distance and kinetic energy of particles in the structure has been experimentally shown. Isotropy has crucial importance for the study of various processes in such systems, including the kinetics of phase transitions, the formation of directed flows, wave propagation, and others.

Grand canonical Brownian dynamics simulations of adsorption and self-assembly of SAS-6 rings on a surface

The Spindle Assembly Abnormal Protein 6 (SAS-6) forms dimers, which then self-assemble into rings that are critical for the nine-fold symmetry of the centriole organelle. It has recently been shown experimentally that the self-assembly of SAS-6 rings is strongly facilitated on a surface, shifting the reaction equilibrium by four orders of magnitude compared to the bulk. Moreover, a fraction of non-canonical symmetries (i.e., different from nine) was observed. In order to understand which aspects of the system are relevant to ensure efficient self-assembly and selection of the nine-fold symmetry, we have performed Brownian dynamics computer simulation with patchy particles and then compared our results with the experimental ones. Adsorption onto the surface was simulated by a grand canonical Monte Carlo procedure and random sequential adsorption kinetics. Furthermore, self-assembly was described by Langevin equations with hydrodynamic mobility matrices. We find that as long as the interaction energies are weak, the assembly kinetics can be described well by coagulation-fragmentation equations in the reaction-limited approximation. By contrast, larger interaction energies lead to kinetic trapping and diffusion-limited assembly. We find that the selection of nine-fold symmetry requires a small value for the angular interaction range. These predictions are confirmed by the experimentally observed reaction constant and angle fluctuations. Overall, our simulations suggest that the SAS-6 system works at the crossover between a relatively weak binding energy that avoids kinetic trapping and a small angular range that favors the nine-fold symmetry.

Anticooperative Supramolecular Oligomerization Mediated by V-Shaped Monomer Design and Unconventional Hydrogen Bonds

After more than three decades of extensive investigations on supramolecular polymers, strategies for self-limiting growth still remain challenging. Herein, we exploit a new V-shaped monomer design to achieve anticooperatively formed oligomers with superior robustness and high luminescence. In toluene, the monomer-oligomer equilibrium is shifted to the monomer side, enabling the elucidation of the molecular packing modes and the resulting (weak) anticooperativity. Steric effects associated with an antiparallel staircase organization of the dyes are proposed to outcompete aromatic and unconventional B-F⋅⋅⋅H-N/C interactions, restricting the growth at the stage of oligomers. In methylcyclohexane (MCH), the packing modes and the anticooperativity are preserved; however, pronounced solvophobic and chain-enwrapping effects lead to thermally ultrastable oligomers. Our results shed light on understanding anticooperative effects and restricted growth in self-assembly.

Building blocks of non-Euclidean ribbons: size-controlled self-assembly via discrete frustrated particles

Geometric frustration offers a pathway to soft matter self-assembly with controllable finite sizes. While the understanding of frustration in soft matter assembly derives almost exclusively from continuum elastic descriptions, a current challenge is to understand the connection between microscopic physical properties of misfitting "building blocks" and emergent assembly behavior at the mesoscale. We present and analyze a particle-based description of what is arguably the best studied example for frustrated soft matter assembly, negative-curvature ribbon assembly, observed in both assemblies of chiral surfactants and shape-frustrated nanoparticles. Based on our particle model, known as saddle wedge monomers, we numerically test the connection between microscopic shape and interactions of the misfitting subunits and the emergent behavior at the supra-particle scale, specifically focussing on the propagation and relaxation of inter-particle strains, the emergent role of extrinsic shape on frustrated ribbons and the equilibrium regime of finite width selection. Beyond the intuitive role of shape misfit, we show that self-limitation is critically dependent on the finite range of cohesive interactions, with larger size finite assemblies requiring increasing short-range interparticle forces. Additionally, we demonstrate that non-linearities arising from discrete particle interactions alter self-limiting behavior due to both strain-softening in shape-flattened assembly and partial yielding of highly strained bonds, which in turn may give rise to states of hierarchical, multidomain assembly. Tracing the regimes of frustration-limited assembly to the specific microscopic features of misfitting particle shapes and interactions provides necessary guidance for translating the theory of size-programmable assembly into design of intentionally-frustrated colloidal particles.

Explicit demonstration of geometric frustration in chiral liquid crystals

Many solid materials and liquid crystals exhibit geometric frustration, meaning that they have an ideal local structure that cannot fill up space. For that reason, the global phase must be a compromise between the ideal local structure and geometric constraints. As an explicit example of geometric frustration, we consider a chiral liquid crystal confined in a long cylinder with free boundaries. When the radius of the tube is sufficiently small, the director field forms a double-twist configuration, which is the ideal local structure. However, when the radius becomes larger (compared with the natural twist of the liquid crystal), the double-twist structure cannot fill space, and hence the director field must transform into some other chiral structure that can fill space. This space-filling structure may be either (1) a cholesteric phase with single twist, or (2) a set of double-twist regions separated by a disclination, which can be regarded as the beginning of a blue phase. We investigate these structures using theory and simulations, and show how the relative free energies depend on the system size, the natural twist, and the disclination energy. As another example, we also study a cholesteric liquid crystal confined between two infinite parallel plates with free boundaries.

Geometrically programmed self-limited assembly of tubules using DNA origami colloids

Nature is replete with self-assembled materials that have one or more self-limited dimensions, including shells, tubules, and fibers. Despite significant advances in making nanometer- and micrometer-scale subunits, the programmable assembly of similar self-limiting architectures from synthetic components has remained largely out of reach. In this article, we create geometrically programmed subunits using DNA origami and study their assembly into tubules with a self-limited width. We show that the average self-limited dimension can be tuned by changing the local curvature encoded in a single subunit. Exploiting the programmability of our system, we further test the tradeoffs between fidelity and complexity embodied by two paradigms for self-limited assembly: self-closure through programmed curvature and addressable assembly through programmed specific interactions.

Self-assembly is one of the most promising strategies for making functional materials at the nanoscale, yet new design principles for making self-limiting architectures, rather than spatially unlimited periodic lattice structures, are needed. To address this challenge, we explore the tradeoffs between addressable assembly and self-closing assembly of a specific class of self-limiting structures: cylindrical tubules. We make triangular subunits using DNA origami that have specific, valence-limited interactions and designed binding angles, and we study their assembly into tubules that have a self-limited width that is much larger than the size of an individual subunit. In the simplest case, the tubules are assembled from a single component by geometrically programming the dihedral angles between neighboring subunits. We show that the tubules can reach many micrometers in length and that their average width can be prescribed through the dihedral angles. We find that there is a distribution in the width and the chirality of the tubules, which we rationalize by developing a model that considers the finite bending rigidity of the assembled structure as well as the mechanism of self-closure. Finally, we demonstrate that the distributions of tubules can be further sculpted by increasing the number of subunit species, thereby increasing the assembly complexity, and demonstrate that using two subunit species successfully reduces the number of available end states by half. These results help to shed light on the roles of assembly complexity and geometry in self-limited assembly and could be extended to other self-limiting architectures, such as shells, toroids, or triply periodic frameworks.

Polymorphic self-assembly of helical tubules is kinetically controlled

In contrast to most self-assembling synthetic materials, which undergo unbounded growth, many biological self-assembly processes are self-limited. That is, the assembled structures have one or more finite dimensions that are much larger than the size scale of the individual monomers. In many such cases, the finite dimension is selected by a preferred curvature of the monomers, which leads to self-closure of the assembly. In this article, we study an example class of self-closing assemblies: cylindrical tubules that assemble from triangular monomers. By combining kinetic Monte Carlo simulations, free energy calculations, and simple theoretical models, we show that a range of programmable size scales can be targeted by controlling the intricate balance between the preferred curvature of the monomers and their interaction strengths. However, their assembly is kinetically controlled-the tubule morphology is essentially fixed shortly after closure, resulting in a distribution of tubule widths that is significantly broader than the equilibrium distribution. We develop a simple kinetic model based on this observation and the underlying free-energy landscape of assembling tubules that quantitatively describes the distributions. Our results are consistent with recent experimental observations of tubule assembly from triangular DNA origami monomers. The modeling framework elucidates design principles for assembling self-limited structures from synthetic components, such as artificial microtubules that have a desired width and chirality.

Thermodynamic Size Control in Curvature-Frustrated Tubules: Self-Limitation with Open Boundaries

We use computational modeling to investigate the assembly thermodynamics of a particle-based model for geometrically frustrated assembly, in which the local packing geometry of subunits is incompatible with uniform, strain-free large-scale assembly. The model considers discrete triangular subunits that drive assembly toward a closed, hexagonal-ordered tubule, but have geometries that locally favor negative Gaussian curvature. We use dynamical Monte Carlo simulations and enhanced sampling methods to compute the free energy landscape and corresponding self-assembly behavior as a function of experimentally accessible parameters that control assembly driving forces and the magnitude of frustration. The results determine the parameter range where finite-temperature self-limiting assembly occurs, in which the equilibrium assembly size distribution is sharply peaked around a well-defined finite size. The simulations also identify two mechanisms by which the system can escape frustration and assemble to unlimited size, and determine the particle-scale properties of subunits that suppress unbounded growth.

Nanoparticle cluster formation mechanisms elucidated via Markov state modeling: Attraction range effects, aggregation pathways, and counterintuitive transition rates

Nanoparticle clusters are promising candidates for developing functional materials. However, it is still a challenging task to fabricate them in a predictable and controllable way, which requires investigation of the possible mechanisms underlying cluster formation at the nanoscale. By constructing Markov state models (MSMs) at the microstate level, we find that for highly dispersed particles to form a highly aggregated cluster, there are multiple coexisting pathways, which correspond to direct aggregation, or pathways that need to pass through partially aggregated, intermediate states. Varying the range of attraction between nanoparticles is found to significantly affect pathways. As the attraction range becomes narrower, compared to direct aggregation, some pathways that need to pass through partially aggregated intermediate states become more competitive. In addition, from MSMs constructed at the macrostate level, the aggregation rate is found to be counterintuitively lower with a lower free-energy barrier, which is also discussed.

Controlling cluster size in 2D phase-separating binary mixtures with specific interactions

By varying the concentration of molecules in the cytoplasm or on the membrane, cells can induce the formation of condensates and liquid droplets, similar to phase separation. Their thermodynamics, much studied, depends on the mutual interactions between microscopic constituents. Here, we focus on the kinetics and size control of 2D clusters, forming on membranes. Using molecular dynamics of patchy colloids, we model a system of two species of proteins, giving origin to specific heterotypic bonds. We find that concentrations, together with valence and bond strength, control both the size and the growth time rate of the clusters. In particular, if one species is in large excess, it gradually saturates the binding sites of the other species; the system then becomes kinetically arrested and cluster coarsening slows down or stops, thus yielding effective size selection. This phenomenology is observed both in solid and fluid clusters, which feature additional generic homotypic interactions and are reminiscent of the ones observed on biological membranes.

Anti-cooperative Self-Assembly with Maintained Emission Regulated by Conformational and Steric Effects

Herein, we present a strategy to enable a maintained emissive behavior in the self‐assembled state by enforcing an anti‐cooperative self‐assembly involving weak intermolecular dye interactions. To achieve this goal, we designed a conformationally flexible monomer unit 1 with a central 1,3‐substituted (diphenyl)urea hydrogen bonding synthon that is tethered to two BODIPY dyes featuring sterically bulky trialkoxybenzene substituents at the meso‐position. The competition between attractive forces (H‐bonding and aromatic interactions) and destabilizing effects (steric and competing conformational effects) limits the assembly, halting the supramolecular growth at the stage of small oligomers. Given the presence of weak dye–dye interactions, the emission properties of molecularly dissolved 1 are negligibly affected upon aggregation. Our findings contribute to broadening the scope of emissive supramolecular assemblies and controlled supramolecular polymerization.

Characterizing the free-energy landscapes of DNA origamis

We show how coarse-grained modelling combined with umbrella sampling using distance-based order parameters can be applied to compute the free-energy landscapes associated with mechanical deformations of large DNA nanostructures. We illustrate this approach for the strong bending of DNA nanotubes and the potentially bistable landscape of twisted DNA origami sheets. The homogeneous bending of the DNA nanotubes is well described by the worm-like chain model; for more extreme bending the nanotubes reversibly buckle with the bending deformations localized at one or two "kinks". For a twisted one-layer DNA origami, the twist is coupled to the bending of the sheet giving rise to a free-energy landscape that has two nearly-degenerate minima that have opposite curvatures. By contrast, for a two-layer origami, the increased stiffness with respect to bending leads to a landscape with a single free-energy minimum that has a saddle-like geometry. The ability to compute such landscapes is likely to be particularly useful for DNA mechanotechnology and for understanding stress accumulation during the self-assembly of origamis into higher-order structures.

Cumulative geometric frustration and superextensive energy scaling in a nonlinear classical XY-spin model

Geometric frustration results from a discrepancy between the locally favored arrangement of the constituents of a system and the geometry of the embedding space. Geometric frustration can be either noncumulative, which implies an extensive energy growth, or cumulative, which implies superextensive energy scaling and highly cooperative ground-state configurations which may depend on the dimensions of the system. Cumulative geometric frustration was identified in a variety of continuous systems including liquid crystals, filament bundles, and molecular crystals. However, a spin-lattice model which clearly demonstrates cumulative geometric frustration was lacking. In this paper we describe a nonlinear variation of the XY-spin model on a triangular lattice that displays cumulative geometric frustration. The model is studied numerically and analyzed in three distinct parameter regimes, which are associated with different energy minimizing configurations. We show that, despite the difference in the ground-state structure in the different regimes, in all cases the superextensive power-law growth of the frustration energy for small domains grows with the same universal exponent that is predicted from the structure of the underlying compatibility condition.

Virus Assembly Pathways Inside a Host Cell

Simple RNA viruses self-assemble spontaneously and encapsulate their genome into a shell called the capsid. This process is mainly driven by the attractive electrostatics interaction between the positive charges on capsid proteins and the negative charges on the genome. Despite its importance and many decades of intense research, how the virus selects and packages its native RNA inside the crowded environment of a host cell cytoplasm in the presence of an abundance of nonviral RNA and other anionic polymers has remained a mystery. In this paper, we perform a series of simulations to monitor the growth of viral shells and find the mechanism by which cargo-coat protein interactions can impact the structure and stability of the viral shells. We show that coat protein subunits can assemble around a globular nucleic acid core by forming nonicosahedral cages, which have been recently observed in assembly experiments involving small pieces of RNA. We find that the resulting cages are strained and can easily be split into fragments along stress lines. This suggests that such metastable nonicosahedral intermediates could be easily reassembled into the stable native icosahedral shells if the larger wild-type genome becomes available, despite the presence of a myriad of nonviral RNAs.

Tiling a tubule: how increasing complexity improves the yield of self-limited assembly

The ability to design and synthesize ever more complicated colloidal particles opens the possibility of self-assembling a zoo of complex structures, including those with one or more self-limited length scales. An undesirable feature of systems with self-limited length scales is that thermal fluctuations can lead to the assembly of nearby, off-target states. We investigate strategies for limiting off-target assembly by using multiple types of subunits. Using simulations and energetics calculations, we explore this concept by considering the assembly of tubules built from triangular subunits that bind edge to edge. While in principle, a single type of triangle can assemble into tubules with a monodisperse width distribution, in practice, the finite bending rigidity of the binding sites leads to the formation of off-target structures. To increase the assembly specificity, we introduce tiling rules for assembling tubules from multiple species of triangles. We show that the selectivity of the target structure can be dramatically improved by using multiple species of subunits, and provide a prescription for choosing the minimum number of subunit species required for near-perfect yield. Our approach of increasing the system's complexity to reduce the accessibility of neighboring structures should be generalizable to other systems beyond the self-assembly of tubules.

Cumulative geometric frustration in physical assemblies

Geometric frustration arises whenever the constituents of a physical assembly locally favor an arrangement that cannot be realized globally. Recently, such frustrated assemblies were shown to exhibit filamentation, size limitation, large morphological variations and other exotic response properties. While these unique characteristics can be shown to be a direct outcome of the geometric frustration, some geometrically frustrated systems do not exhibit any of the above phenomena. In this work we exploit the intrinsic approach to provide a framework for directly addressing the frustration in physical assemblies. The framework highlights the role of the compatibility conditions associated with the intrinsic fields describing the physical assembly. We show that the structure of the compatibility conditions determines the behavior of small assemblies and in particular predicts their superextensive energy growth exponent. We illustrate the use of this framework to several well-known frustrated assemblies.

Disorder in cellular packing can alter proliferation dynamics to regulate growth

The mechanisms by which an organ regulates its growth are not yet fully understood, especially when the cells are closely packed as in epithelial tissues. We explain growth arrest as a collective dynamical transition in coupled oscillators on disordered lattices. As the cellular morphologies become homogeneous over the course of development, the signals induced by cell-cell contact increase beyond a critical value that triggers coordinated cessation of the cell-cycle oscillators driving cell division. Thus, control of cell proliferation is causally related to the geometry of cellular packing.

Size-encoded hierarchical self-assembly of nanoparticles into chains and tubules

Hierarchical structures with sophisticated patterns allow the emergence of challenging properties. However, the highly cooperative and specific interactions needed for assembly spanning different length scales are typically lacking in inorganic nanoparticles (NPs). Here we show that size can be a common structural driving force for controlling hierarchical assembly of inorganic NPs into anisotropic superstructures. It involves first the self-limiting assembly of small CdS NPs into large supraparticles and their subsequent spontaneous organization into chains and tubules hundreds of nanometers long. Our quantitative calculations based on DLVO theory reveals an intrinsic size effect relating to the dimension change of assembly units in accordance with a negative cooperativity. It is shown that the size increase in building blocks creates an effective kinetic barrier contrast at different attachment sites due to the increase of interparticle electrostatic repulsion, switching the assembly from thermodynamically preferred 3D to kinetically favored 1D pathway. The size-encoded hierarchical assembly is accompanied by the ligand-controlled Oswald ripening process, which is responsible for the variation of hierarchical patterns from chains to tubules. The general principle in governing multistage inorganic NP ordering represents an important guideline toward the complex mesoscale structures that may surpass traditional materials in both design and functionality.

Vesicle shape transformations driven by confined active filaments

In active matter systems, deformable boundaries provide a mechanism to organize internal active stresses. To study a minimal model of such a system, we perform particle-based simulations of an elastic vesicle containing a collection of polar active filaments. The interplay between the active stress organization due to interparticle interactions and that due to the deformability of the confinement leads to a variety of filament spatiotemporal organizations that have not been observed in bulk systems or under rigid confinement, including highly-aligned rings and caps. In turn, these filament assemblies drive dramatic and tunable transformations of the vesicle shape and its dynamics. We present simple scaling models that reveal the mechanisms underlying these emergent behaviors and yield design principles for engineering active materials with targeted shape dynamics.

Probing the size-dependent polarizability of mesoscopic ionic clusters and their induced-dipole interactions

Mesoscopic clusters composed of oppositely charged particles are ubiquitous in synthetic and biological soft materials. The effective interaction between these clusters is influenced by their polarizability, that is, the ability of their constituent charges to re-arrange in response to an external electrical field. Here, using coarse-grained simulations, we show that the polarizability of electrically neutral ionic clusters decreases as the number of constituent charges increases and/or their Coulombic interaction strength increases for various ion valencies, ion densities, and degrees of cluster boundary hardness. For clusters of random ionomers and their counterions, their polarizability is shown to depend on the number of polymer chains. The variation of the cluster polarizability with the cluster size indicates that throughout the assembly, the induced-dipole interactions between the clusters may be reduced substantially as they acquire more charges while maintaining zero net charge. Under certain conditions, the induced-dipole interactions may become repulsive, as inferred from our simulations with a polarizable solvent. As a result, the dipole-induced related interactions can serve as a counterbalancing force that contributes to the self-limiting aggregation of charge-containing assemblies.

Solution of disordered microphases in the Bethe approximation

The periodic microphases that self-assemble in systems with competing short-range attractive and long-range repulsive (SALR) interactions are structurally both rich and elegant. Significant theoretical and computational efforts have thus been dedicated to untangling their properties. By contrast, disordered microphases, which are structurally just as rich but nowhere near as elegant, have not been as carefully considered. Part of the difficulty is that simple mean-field descriptions make a homogeneity assumption that washes away all of their structural features. Here, we study disordered microphases by exactly solving a SALR model on the Bethe lattice. By sidestepping the homogenization assumption, this treatment recapitulates many of the key structural regimes of disordered microphases, including particle and void cluster fluids as well as gelation. This analysis also provides physical insight into the relationship between various structural and thermal observables, between criticality and physical percolation, and between glassiness and microphase ordering.