Achiral symmetry breaking and positive Gaussian modulus lead to scalloped colloidal membranes

Dynamics of a self-interacting sheet in shear flow

Solution processing of 2D materials such as graphene is important for applications thereof, yet a complete fundamental understanding of how 2D materials behave dynamically in solution is lacking. Here, we extend previous work by Silmore et al., Soft Matter, 2021, 17(18), 4707-4718 by adding short-ranged Lennard-Jones interactions to 2D sheets in shear flow. We find that the addition of these interactions allows for a rich landscape of conformations which depend on the balance between shear strength, bending rigidity, and interaction strength as well as the initial configuration of the sheet. We explore this conformational space and classify sheets as flat, tumbling, 1D folded, or 2D folded based on their conformational properties. We use kinetic and energetic arguments to explain why sheets adopt certain conformations within the folded regime. Finally, we calculate the stresslet and find that, even in the absence of thermal fluctuations and multiple sheet interactions, shear-thinning followed by shear-thickening behavior can appear.

Chiral fluid membranes with orientational order and multiple edges

We carry out Monte Carlo simulations on fluid membranes with orientational order and multiple edges in the presence and absence of external forces. The membrane resists bending and has an edge tension, the orientational order couples with the membrane surface normal through a cost for tilting, and there is a chiral liquid crystalline interaction. In the absence of external forces, a membrane initialized as a vesicle will form a disk at low chirality, with the directors forming a smectic-A phase with alignment perpendicular to the membrane surface except near the edge. At large chirality a catenoid-like shape or a trinoid-like shape is formed, depending on the number of edges in the initial vesicle. This shape change is accompanied by cholesteric ordering of the directors and multiple π walls connecting the membrane edges and wrapping around the membrane neck. If the membrane is initialized instead in a cylindrical shape and stretched by an external force, it maintains a nearly cylindrical shape but additional liquid crystalline phases appear. For large tilt coupling and low chirality, a smectic-A phase forms where the directors are normal to the surface of the membrane. For lower values of the tilt coupling, a nematic phase appears at zero chirality with the average director oriented perpendicular to the long axis of the membrane, while for nonzero chirality a cholesteric phase appears. The π walls are tilt walls at low chirality and transition to twist walls as chirality is increased. We construct a continuum model of the director field to explain this behavior.

The thermodynamic origins of chiral twist in monolayer assemblies of rod-like colloids

The propagation of chirality across scales is a common but poorly understood phenomenon in soft matter. Here, using computer simulations, we study twisted monolayer assemblies formed by both chiral and achiral rod-like particles in the presence of non-adsorbing polymer and characterise the thermodynamic driving forces responsible for the twisting. We observe assemblies with both like and inverted chirality relative to the rods and show that the preferred twist is already determined during the initial stage of the self-assembly. Depending on the geometry of the constituent rods, the chiral twist is regulated by either the entropy gain of the polymer, or of the rods, or both. This can include important contributions from changes in both the surface area and volume of the monolayer and from rod fluctuations perpendicular to the monolayer. These findings can deepen our understanding of why chirality propagates and of how to control it.

Spontaneous organization and phase separation of skyrmions in chiral active matter

Skyrmions are topologically protected vortex-like excitations that hold promise for applications such as information processing and electron manipulation. Here we combine theoretical analysis and numerical simulations to show that skyrmions can spontaneously emerge in chiral active matter without external confinements or regulation. Strikingly, these activity-driven skyrmions can either self-organize into a periodic, stable square lattice consisting of half Néel skyrmions and antiskyrmions, where the in-plane flows display an antiferromagnetic vortex array, or undergo phase separation between skyrmions with different topological numbers. We identify that the emerging skyrmion dynamics stems from the competition between the chiral and polar coherence length scales dictated by the interplay of intrinsic chirality, polarity, and elasticity in the system. Our results reveal unanticipated topological excitations, self-organization, and phase separation in non-equilibrium systems and also suggest a potential way towards engineering complicated bespoke skyrmionic structures through manipulating active matter.

Controlling the shape and topology of two-component colloidal membranes

Changes in the geometry and topology of self-assembled membranes underlie diverse processes across cellular biology and engineering. Similar to lipid bilayers, monolayer colloidal membranes have in-plane fluid-like dynamics and out-of-plane bending elasticity. Their open edges and micrometer-length scale provide a tractable system to study the equilibrium energetics and dynamic pathways of membrane assembly and reconfiguration. Here, we find that doping colloidal membranes with short miscible rods transforms disk-shaped membranes into saddle-shaped surfaces with complex edge structures. The saddle-shaped membranes are well approximated by Enneper's minimal surfaces. Theoretical modeling demonstrates that their formation is driven by increasing the positive Gaussian modulus, which in turn, is controlled by the fraction of short rods. Further coalescence of saddle-shaped surfaces leads to diverse topologically distinct structures, including shapes similar to catenoids, trinoids, four-noids, and higher-order structures. At long timescales, we observe the formation of a system-spanning, sponge-like phase. The unique features of colloidal membranes reveal the topological transformations that accompany coalescence pathways in real time. We enhance the functionality of these membranes by making their shape responsive to external stimuli. Our results demonstrate a pathway toward control of thin elastic sheets' shape and topology-a pathway driven by the emergent elasticity induced by compositional heterogeneity.

Thermally fluctuating, semiflexible sheets in simple shear flow

We perform Brownian dynamics simulations of semiflexible colloidal sheets with hydrodynamic interactions and thermal fluctuations in shear flow. As a function of the ratio of bending rigidity to shear energy (a dimensionless quantity we denote S) and the ratio of bending rigidity to thermal energy, we observe a dynamical transition from stochastic flipping to crumpling and continuous tumbling. This dynamical transition is broadened by thermal fluctuations, and the value of S at which it occurs is consistent with the onset of chaotic dynamics found for athermal sheets. The effects of different dynamical conformations on rheological properties such as viscosity and normal stress differences are also quantified. Namely, the viscosity in a dilute dispersion of sheets is found to decrease with increasing shear rate (shear-thinning) up until the dynamical crumpling transition, at which point it increases again (shear-thickening), and non-zero first normal stress differences are found that exhibit a local maximum with respect to temperature at large S (small shear rate). These results shed light on the dynamical behavior of fluctuating 2D materials dispersed in fluids and should greatly inform the design of associated solution processing methods.

Multiphase Coexistences in Rod-Polymer Mixtures

Using recently derived analytical equations of state for hard rod dispersions, we predict the phase behavior of athermal rod-polymer mixtures with free volume theory. The rods are modeled as hard spherocylinders, while the nonadsorbing polymer chains are described as penetrable hard spheres. It is demonstrated that all of the different types of phase states that are stable for pure colloidal rod dispersions can coexist with any combination of these phases if polymers are added, depending on the concentrations, rod aspect ratio, and polymer-rod size ratio. This includes novel two-, three-, and four-phase coexistences and isostructural coexistences between dilute and concentrated phases of the same kind, even for the more ordered (liquid) crystal phases. This work provides insight into the conditions at which particular multiphase coexistences are expected for well-defined model colloidal rod-polymer mixtures. We provide a quantitative map detailing the various types of isostructural coexistences, which confirms an early qualitative hypothesis by Bolhuis et al. ( J. Chem. Phys. 107, 1997 1551).

Axisymmetric membranes with edges under external force: buckling, minimal surfaces, and tethers

We use theory and numerical computation to determine the shape of an axisymmetric fluid membrane with a resistance to bending and constant area. The membrane connects two rings in the classic geometry that produces a catenoidal shape in a soap film. In our problem, we find infinitely many branches of solutions for the shape and external force as functions of the separation of the rings, analogous to the infinite family of eigenmodes for the Euler buckling of a slender rod. Special attention is paid to the catenoid, which emerges as the shape of maximal allowable separation when the area is less than a critical area equal to the planar area enclosed by the two rings. A perturbation theory argument directly relates the tension of catenoidal membranes to the stability of catenoidal soap films in this regime. When the membrane area is larger than the critical area, we find additional cylindrical tether solutions to the shape equations at large ring separation, and that arbitrarily large ring separations are possible. These results apply for the case of vanishing Gaussian curvature modulus; when the Gaussian curvature modulus is nonzero and the area is below the critical area, the force and the membrane tension diverge as the ring separation approaches its maximum value. We also examine the stability of our shapes and analytically show that catenoidal membranes have markedly different stability properties than their soap film counterparts.

Deformation and orientational order of chiral membranes with free edges

Motivated by experiments on colloidal membranes composed of chiral rod-like viruses, we use Monte Carlo methods to simulate these systems and determine the phase diagram for the liquid crystalline order of the rods and the membrane shape. We generalize the Lebwohl-Lasher model for a nematic with a chiral coupling to a curved surface with edge tension and a resistance to bending, and include an energy cost for tilting of the rods relative to the local membrane normal. The membrane is represented by a triangular mesh of hard beads joined by bonds, where each bead is decorated by a director. The beads can move, the bonds can reconnect and the directors can rotate at each Monte Carlo step. When the cost of tilt is small, the membrane tends to be flat, with the rods only twisting near the edge for low chiral coupling, and remaining parallel to the normal in the interior of the membrane. At high chiral coupling, the rods twist everywhere, forming a cholesteric state. When the cost of tilt is large, the emergence of the cholesteric state at high values of the chiral coupling is accompanied by the bending of the membrane into a saddle shape. Increasing the edge tension tends to flatten the membrane. These results illustrate the geometric frustration arising from the inability of a surface normal to have twist.

Buckling, crumpling, and tumbling of semiflexible sheets in simple shear flow

As 2D materials such as graphene, transition metal dichalcogenides, and 2D polymers become more prevalent, solution processing and colloidal-state properties are being exploited to create advanced and functional materials. However, our understanding of the fundamental behavior of 2D sheets and membranes in fluid flow is still lacking. In this work, we perform numerical simulations of athermal semiflexible sheets with hydrodynamic interactions in shear flow. For sheets initially oriented near the flow-vorticity plane, we find buckling instabilities of different mode numbers that vary with bending stiffness and can be understood with a quasi-static model of elasticity. For different initial orientations, chaotic tumbling trajectories are observed. Notably, we find that sheets fold or crumple before tumbling but do not stretch again upon applying greater shear.

Chiral molecules on curved colloidal membranes

Colloidal membranes, self assembled monolayers of aligned rod like molecules, offer a template for designing membranes with definite shapes and curvature, and possibly new functionalities in the future. Often the constituent rods, due to their molecular chirality, are tilted with respect to the membrane normal. Spatial patterns of this tilt on curved membranes result from a competition among depletion forces, nematic interactions, molecular chirality and boundary effects. We present a covariant theory for the tilt pattern on minimal surfaces, like helicoids and catenoids, which have been generated in the laboratory only recently. We predict several non-uniform tilt patterns, some of which are consistent with experimental observations and some, which are yet to be discovered.

Shapes of fluid membranes with chiral edges

We carry out Monte Carlo simulations of a colloidal fluid membrane with a free edge and composed of chiral rodlike viruses. The membrane is modeled by a triangular mesh of beads connected by bonds in which the bonds and beads are free to move at each Monte Carlo step. Since the constituent viruses are experimentally observed to twist only near the membrane edge, we use an effective energy that favors a particular sign of the geodesic torsion of the edge. The effective energy also includes the membrane bending stiffness, edge bending stiffness, and edge tension. We find three classes of membrane shapes resulting from the competition of the various terms in the free energy: branched shapes, chiral disks, and vesicles. Increasing the edge bending stiffness smooths the membrane edge, leading to correlations among the membrane normals at different points along the edge. The normalized power spectrum for edge displacements shows a peak with increasing preferred geodesic torsion. We also consider membrane shapes under an external force by fixing the distance between two ends of the membrane and finding the shape for increasing values of the distance between the two ends. As the distance increases, the membrane twists into a ribbon, with the force eventually reaching a plateau.

Force-Induced Formation of Twisted Chiral Ribbons

We demonstrate that an achiral stretching force transforms disk-shaped colloidal membranes composed of chiral rods into twisted ribbons with handedness opposite the preferred twist of the rods. Using an experimental technique that enforces torque-free boundary conditions we simultaneously measure the force-extension curve and the ribbon shape. An effective theory that accounts for the membrane bending energy and uses geometric properties of the edge to model the internal liquid crystalline degrees of freedom explains both the measured force-extension curve and the force-induced twisted shape.

All twist and no bend makes raft edges splay: Spontaneous curvature of domain edges in colloidal membranes

Using theory and experiments, we study the interface between two immiscible domains in a colloidal membrane composed of rigid rods of different lengths. Geometric considerations of rigid rod packing imply that a domain of sufficiently short rods in a background membrane of long rods is more susceptible to twist than the inverse structure, a long-rod domain in a short-rod membrane. The midplane tilt at the interdomain edge forces splay, which, in turn, manifests as spontaneous edge curvature with energetics controlled by the length asymmetry of constituent rods. A thermodynamic model of such tilt-curvature coupling at interdomain edges explains a number of experimental observations, including annularly shaped long-rod domains, and a nonmonotonic dependence of edge twist on domain radius. Our work shows how coupling between orientational and compositional degrees of freedom in two-dimensional fluids gives rise to complex shapes of fluid domains, analogous to shape transitions in 3D fluid vesicles.

Direct Observation of Topological Defects in Striped Block Copolymer Discs and Polymersomes

Topology and defects are of fundamental importance for ordered structures on all length scales. Despite extensive research on block copolymer self-assembly in solution, knowledge about topological defects and their effect on nanostructure formation has remained limited. Here, we report on the self-assembly of block copolymer discs and polymersomes with a cylinder line pattern on the surface that develops specific combinations of topological defects to satisfy the Euler characteristics for closed spheres as described by Gauss-Bonnet theorem. The dimension of the line pattern allows the direct visualization of defect emergence, evolution, and annihilation. On discs, cylinders either form end-caps that coincide with λ^+1/2 disclinations or they bend around τ^+1/2 disclinations in 180° turns (hairpin loops). On polymersomes, two λ^+1/2 defects connect into three-dimensional (3D) Archimedean spirals, while two τ^+1/2 defects form 3D Fermat spirals. Electron tomography reveals two complementary line patterns on the inside and outside of the polymersome membrane, where λ^+1/2 and τ^+1/2 disclinations always eclipse on opposing sides ("defect communication"). Attractive defects are able to annihilate with each other into +1 disclinations and stabilize anisotropic polymersomes with sharp tips through screening of high-energy curvature. This study fosters our understanding of the behavior of topological defects in self-assembled polymer materials and aids in the design of polymersomes with preprogrammed shapes governed by synthetic block length and topological rules.

Modeling the dynamics of phospholipids in the fluid phase of liposomes

We present the derivation of a new model to describe neutron spin echo spectroscopy and quasi-elastic neutron scattering data on liposomes. We compare the new model with existing approaches and benchmark it with experimental data. The analysis indicates the importance of including all major contributions in the modeling of the intermediate scattering function. Simultaneous analysis of the experimental data on lipids with full contrast and tail contrast matched samples reveals highly confined lipid tail motion. A comparison of their dynamics demonstrates the statistical independence of tail-motion and height-height correlation of the membrane. A more detailed analysis indicates that the lipid tail relaxation is confined to a potential with cylindrical symmetry, in addition to the undulation and diffusive motion of the liposome. Despite substantial differences in the chemistry of the fatty acid tails, the observation indicates a universal behavior. The analysis of partially deuterated systems confirms the strong contribution of the lipid tail to the intermediate scattering function. Within the time range from 5 to 100 ns, the intermediate scattering function can be described by the height-height correlation function. The existence of the fast-localized tail motion and the contribution of slow translational diffusion of liposomes determine the intermediate scattering function for t < 5 ns and t > 100 ns, respectively. Taking into account the limited time window lowers the bending moduli by a factor of 1.3 (DOPC) to 2 (DMPC) compared to the full range.

Conformational switching of chiral colloidal rafts regulates raft-raft attractions and repulsions

Membrane-mediated particle interactions depend both on the properties of the particles themselves and the membrane environment in which they are suspended. Experiments have shown that chiral rod-like inclusions dissolved in a colloidal membrane of opposite handedness assemble into colloidal rafts, which are finite-sized reconfigurable droplets consisting of a large but precisely defined number of rods. We systematically tune the chirality of the background membrane and find that, in the achiral limit, colloidal rafts acquire complex structural properties and interactions. In particular, rafts can switch between 2 chiral states of opposite handedness, which alters the nature of the membrane-mediated raft-raft interactions. Rafts with the same chirality have long-ranged repulsions, while those with opposite chirality acquire attractions with a well-defined minimum. Both attractive and repulsive interactions are qualitatively explained by a continuum model that accounts for the coupling between the membrane thickness and the local tilt of the constituent rods. These switchable interactions enable assembly of colloidal rafts into intricate higher-order architectures, including stable tetrameric clusters and "ionic crystallites" of counter-twisting domains organized on a binary square lattice. Furthermore, the properties of individual rafts, such as their sizes, are controlled by their complexation with other rafts. The emergence of these complex behaviors can be rationalized purely in terms of generic couplings between compositional and orientational order of fluids of rod-like elements. Thus, the uncovered principles might have relevance for conventional lipid bilayers, in which the assembly of higher-order structures is also mediated by complex membrane-mediated interactions.

Self assembly of cyclic polygon shaped fluid colloidal membranes through pinning

2D fluid monolayer membranes of rod-like viruses spontaneously form in a mixture of rods and polymers through depletion attraction. The rods are uniformly oriented within the bulk and twist in a zone around the membrane edge. Surprisingly, we find that cyclic polygonal shaped colloidal membranes form when polymers are added to a mixture of long and short-thick rods with the long and short-thick rods forming the faceted core and lobes of the polygon, respectively. We demonstrate that the origin of this anisotropic shape lies in the phenomenon of spreading of one liquid over another in the presence of disorder. As a membrane of short-thick rods spreads over another of longer rods, the edge bound rods untwist to become part of the newly formed two-rod interface. However, a small fraction of rods fail to untwist as the two rod interface forms and act as mobile pinning centers. Capillary flow of short-thick rods drives all the pinning centers to a single location in the composite membrane which now acts like a junction. This pinning junction inhibits complete engulfing of one membrane by the other. Repeated sequential events like this then lead to formation of multiple junctions and the overall cyclic polygon topology. We find that pinning junctions are weakly cross-linked in nature instead of being topological defects. We outline the necessary and sufficient constraints on the nature of rods to obtain stable out of equilibrium cyclic polygon membranes. Our results show a unique counter-intuitive scenario where defects lead to self-assembly of ordered structures.

Topological bounds of bending energy for lipid vesicles

The Helfrich bending energy plays an important role in providing a mechanism for the conformation of a lipid vesicle in theoretical biophysics, which is governed by the principle of energy minimization over configurations of appropriate topological characteristics. We will show that the presence of a quantity called the spontaneous curvature obstructs the existence of a minimizer of the Helfrich energy over the set of embedded ring tori. In addition, despite the well-realized knowledge that lipid vesicles may present themselves in a variety of shapes of complicated topology, there is a lack of topological bounds for the Helfrich energy. To overcome these difficulties, we consider a general scale-invariant anisotropic curvature energy that extends the Canham elastic bending energy developed in modeling a biconcave-shaped red blood cell. We will show that, up to a rescaling of the generating radii, there is a unique minimizer of the energy over the set of embedded ring tori, in the entire parameter regime, which recovers the Willmore minimizer in its Canham isotropic limit. We also show how elevated anisotropy favors energetically a clear transition from spherical-, to ellipsoidal-, and then to biconcave-shaped surfaces, for a lipid vesicle. We then establish some genus-dependent topological lower and upper bounds for the anisotropic energy. Finally, we derive the shape equation of the generalized bending energy, which extends the well-known Helfrich shape equation.

Filamentous phages as building blocks for reconfigurable and hierarchical self-assembly

Filamentous bacteriophages such as fd-like viruses are monodisperse rod-like colloids that have well defined properties of diameter, length, rigidity, charge and chirality. Engineering these viruses leads to a library of colloidal rods, which can be used as building blocks for reconfigurable and hierarchical self-assembly. Their condensation in an aqueous solution with additive polymers, which act as depletants to induce attraction between the rods, leads to a myriad of fluid-like micronic structures ranging from isotropic/nematic droplets, colloid membranes, achiral membrane seeds, twisted ribbons, π-wall, pores, colloidal skyrmions, Möbius anchors, scallop membranes to membrane rafts. These structures, and the way that they shape-shift, not only shed light on the role of entropy, chiral frustration and topology in soft matter, but also mimic many structures encountered in different fields of science. On the one hand, filamentous phages being an experimental realization of colloidal hard rods, their condensation mediated by depletion interactions constitutes a blueprint for the self-assembly of rod-like particles and provides a fundamental foundation for bio- or material-oriented applications. On the other hand, the chiral properties of the viruses restrict the generalities of some results but vastly broaden the self-assembly possibilities.

Curvature instability of chiral colloidal membranes on crystallization

Buckling and wrinkling instabilities are failure modes of elastic sheets that are avoided in the traditional material design. Recently, a new paradigm has appeared where these instabilities are instead being utilized for high-performance applications. Multiple approaches such as heterogeneous gelation, capillary stresses, and confinement have been used to shape thin macroscopic elastic sheets. However, it remains a challenge to shape two-dimensional self-assembled monolayers at colloidal or molecular length scales. Here, we show the existence of a curvature instability that arises during the crystallization of finite-sized monolayer membranes of chiral colloidal rods. While the bulk of the membrane crystallizes, its edge remains fluid like and exhibits chiral ordering. The resulting internal stresses cause the flat membrane to buckle macroscopically and wrinkle locally. Our results demonstrate an alternate pathway based on intrinsic stresses instead of the usual external ones to assemble non-Euclidean sheets at the colloidal length scale.

Buckling and wrinkling are instabilities which involve thin elastic sheets and are well-investigated phenomena at the macroscale. Here Saikia et al. investigate curvature instabilities at the colloidal lengthscale in quasi-2D monolayers of rod-like viruses across the fluid-crystal phase transition.

Determining the Gaussian Modulus and Edge Properties of 2D Materials: From Graphene to Lipid Bilayers

The dominant deformation behavior of two-dimensional materials (bending) is primarily governed by just two parameters: bending rigidity and the Gaussian modulus. These properties also set the energy scale for various important physical and biological processes such as pore formation, cell fission and generally, any event accompanied by a topological change. Unlike the bending rigidity, the Gaussian modulus is, however, notoriously difficult to evaluate via either experiments or atomistic simulations. In this Letter, recognizing that the Gaussian modulus and edge tension play a nontrivial role in the fluctuations of a 2D material edge, we derive closed-form expressions for edge fluctuations. Combined with atomistic simulations, we use the developed approach to extract the Gaussian modulus and edge tension at finite temperatures for both graphene and various types of lipid bilayers. Our results possibly provide the first reliable estimate of this elusive property at finite temperatures and appear to suggest that earlier estimates must be revised. In particular, we show that, if previously estimated properties are employed, the graphene-free edge will exhibit unstable behavior at room temperature. Remarkably, in the case of graphene, we show that the Gaussian modulus and edge tension even change sign at finite temperatures.