Stretchy and disordered: Toward understanding fracture in soft network materials via mesoscopic computer simulations.

Soft network materials exist in numerous forms ranging from polymer networks, such as elastomers, to fiber networks, such as collagen. In addition, in colloidal gels, an underlying network structure can be identified, and several metamaterials and textiles can be considered network materials as well. Many of these materials share a highly disordered microstructure and can undergo large deformations before damage becomes visible at the macroscopic level. Despite their widespread presence, we still lack a clear picture of how the network structure controls the fracture processes of these soft materials. In this Perspective, we will focus on progress and open questions concerning fracture at the mesoscopic scale, in which the network architecture is clearly resolved, but neither the material-specific atomistic features nor the macroscopic sample geometries are considered. We will describe concepts regarding the network elastic response that have been established in recent years and turn out to be pre-requisites to understand the fracture response. We will mostly consider simulation studies, where the influence of specific network features on the material mechanics can be cleanly assessed. Rather than focusing on specific systems, we will discuss future challenges that should be addressed to gain new fundamental insights that would be relevant across several examples of soft network materials.

Modeling the cholesteric pitch of apolar cellulose nanocrystal suspensions using a chiral hard-bundle model.

Cellulose nanocrystals (CNCs) are naturally sourced elongated nanocolloids that form cholesteric phases in water and apolar solvents. It is well accepted that CNCs are made of bundles of crystalline microfibrils clustered side-by-side, and there is growing evidence that each individual microfibril is twisted. Yet, the origin of the chiral interactions between CNCs remains unclear. In this work, CNCs are described with a simple model of chiral hard splinters, enabling the prediction of the pitch using density functional theory and Monte Carlo simulations. The predicted pitch P compares well with experimental observations in cotton-based CNC dispersions in apolar solvents using surfactants but also with qualitative trends caused by fractionation or tip sonication in aqueous suspensions. These results suggest that the bundle shape induces an entropy-driven chiral interaction between CNCs, which is the missing link in explaining how chirality is transferred from the molecular scale of cellulose chains to the cholesteric order.

Sharing the Load: Stress Redistribution Governs Fracture of Polymer Double Networks.

The stress response of polymer double networks depends not only on the properties of the constituent networks but also on the interactions arising between them. Here, we demonstrate, via coarse-grained simulations, that both their global stress response and their microscopic fracture mechanics are governed by load sharing through these internetwork interactions. By comparing our results with affine predictions, where stress redistribution is by definition homogeneous, we show that stress redistribution is highly inhomogeneous. In particular, the affine prediction overestimates the fraction of broken chains by almost an order of magnitude. Furthermore, homogeneous stress distribution predicts a single fracture process, while in our simulations, fracture of sacrificial chains takes place in two steps governed by load sharing within a network and between networks, respectively. Our results thus provide a detailed microscopic picture of how inhomogeneous stress redistribution after rupture of chains governs the fracture of polymer double networks.

The role of temperature in the rigidity-controlled fracture of elastic networks.

We study the influence of thermal fluctuations on the fracture of elastic networks, via simulations of the uniaxial extension of central-force spring networks with varying rigidity. Studying their failure response, both at the macroscopic and microscopic level, we find that an increase in temperature corresponds to a more homogeneous stress (re)distribution and induces thermally activated failure of springs. As a consequence, the material strength decreases upon increasing temperature, the microscopic damage spreads over a larger area and a more ductile fracture process is observed. These effects are modulated by network rigidity and can therefore be tuned via the network connectivity and the rupture threshold of the springs. Knowledge of the interplay between temperature and rigidity improves our understanding of the fracture of elastic network materials, such as (biological) polymer networks, and can help to refine design principles for tough soft materials.

Connectivity and plasticity determine collagen network fracture.

Collagen forms the structural scaffold of connective tissues in all mammals. Tissues are remarkably resistant against mechanical deformations because collagen molecules hierarchically self-assemble in fibrous networks that stiffen with increasing strain. Nevertheless, collagen networks do fracture when tissues are overloaded or subject to pathological conditions such as aneurysms. Prior studies of the role of collagen in tissue fracture have mainly focused on tendons, which contain highly aligned bundles of collagen. By contrast, little is known about fracture of the orientationally more disordered collagen networks present in many other tissues such as skin and cartilage. Here, we combine shear rheology of reconstituted collagen networks with computer simulations to investigate the primary determinants of fracture in disordered collagen networks. We show that the fracture strain is controlled by the coordination number of the network junctions, with less connected networks fracturing at larger strains. The hierarchical structure of collagen fine-tunes the fracture strain by providing structural plasticity at the network and fiber level. Our findings imply that low connectivity and plasticity provide protective mechanisms against network fracture that can optimize the strength of biological tissues.

Athermal Fracture of Elastic Networks: How Rigidity Challenges the Unavoidable Size-Induced Brittleness.

By performing extensive simulations with unprecedentedly large system sizes, we unveil how rigidity influences the fracture of disordered materials. We observe the largest damage in networks with connectivity close to the isostatic point and when the rupture thresholds are small. However, irrespective of network and spring properties, a more brittle fracture is observed upon increasing system size. Differently from most of the fracture descriptors, the maximum stress drop, a proxy for brittleness, displays a universal nonmonotonic dependence on system size. Based on this uncommon trend it is possible to identify the characteristic system size L^{*} at which brittleness kicks in. The more the disorder in network connectivity or in spring thresholds, the larger L^{*}. Finally, we speculate how this size-induced brittleness is influenced by thermal fluctuations.

Hard Competition: Stabilizing the Elusive Biaxial Nematic Phase in Suspensions of Colloidal Particles with Extreme Lengths.

We use computer simulations to study the existence and stability of a biaxial nematic N_{b} phase in systems of hard polyhedral cuboids, triangular prisms, and rhombic platelets, characterized by a long (L), medium (M), and short (S) particle axis. For all three shape families, we find stable N_{b} states provided the shape is not only close to the so-called dual shape with M=sqrt[LS] but also sufficiently anisotropic with L/S>9,11,14,23 for rhombi, (two types of) triangular prisms, and cuboids, respectively, corresponding to anisotropies not considered before. Surprisingly, a direct isotropic-N_{b} transition does not occur in these systems due to a destabilization of N_{b} by a smectic (for cuboids and prisms) or a columnar (for platelets) phase at small L/S or by an intervening uniaxial nematic phase at large L/S. Our results are confirmed by a density functional theory provided the third virial coefficient is included and a continuous rather than a discrete (Zwanzig) set of particle orientations is taken into account.

Laser Speckle Strain Imaging reveals the origin of delayed fracture in a soft solid.

Stresses well below the critical fracture stress can lead to highly unpredictable delayed fracture after a long period of seemingly quiescent stability. Delayed fracture is a major threat to the lifetime of materials, and its unpredictability makes it difficult to prevent. This is exacerbated by the lack of consensus on the microscopic mechanisms at its origin because unambiguous experimental proof of these mechanisms remains absent. We present an experimental approach to measure, with high spatial and temporal resolution, the local deformations that precipitate crack nucleation. We apply this method to study delayed fracture in an elastomer and find that a delocalized zone of very small strains emerges as a consequence of strongly localized damage processes. This prefracture deformation zone grows exponentially in space and time, ultimately culminating in the nucleation of a crack and failure of the material as a whole. Our results paint a microscopic picture of the elusive origins of delayed fracture and open the way to detect damage well before it manifests macroscopically.

Connectedness percolation of hard deformed rods.

Nanofiller particles, such as carbon nanotubes or metal wires, are used in functional polymer composites to make them conduct electricity. They are often not perfectly straight cylinders but may be tortuous or exhibit kinks. Therefore we investigate the effect of shape deformations of the rod-like nanofillers on the geometric percolation threshold of the dispersion. We do this by using connectedness percolation theory within a Parsons-Lee type of approximation, in combination with Monte Carlo integration for the average overlap volume in the isotropic fluid phase. We find that a deviation from a perfect rod-like shape has very little effect on the percolation threshold, unless the particles are strongly deformed. This demonstrates that idealized rod models are useful even for nanofillers that superficially seem imperfect. In addition, we show that for small or moderate rod deformations, the universal scaling of the percolation threshold is only weakly affected by the precise particle shape.

Density functional theory and simulations of colloidal triangular prisms.

Nanopolyhedra form a versatile toolbox to investigate the effect of particle shape on self-assembly. Here we consider rod-like triangular prisms to gauge the effect of the cross section of the rods on liquid crystal phase behavior. We also take this opportunity to implement and test a previously proposed version of fundamental measure density functional theory (0D-FMT). Additionally, we perform Monte Carlocomputer simulations and we employ a simpler Onsager theory with a Parsons-Lee correction. Surprisingly and disappointingly, 0D-FMT does not perform better than the Tarazona and Rosenfeld's version of fundamental measure theory (TR-FMT). Both versions of FMT perform somewhat better than the Parsons-Lee theory. In addition, we find that the stability regime of the smectic phase is larger for triangular prisms than for spherocylinders and square prisms.

Phase diagram of binary colloidal rod-sphere mixtures from a 3D real-space analysis of sedimentation-diffusion equilibria.

Self-assembly of binary particle systems offers many new opportunities for materials science. Here, we studied sedimentation equilibria of silica rods and spheres, using quantitative 3D confocal microscopy. We determined not only pressure, density and order parameter profiles, but also the experimental phase diagram exhibiting a stable binary smectic liquid-crystalline phase (Sm2). Using computer simulations we confirmed that the Sm2-phase can be stabilized by entropy alone, which opens up the possibility of combining new materials properties at a wide array of length scales.

Entropy-driven formation of chiral nematic phases by computer simulations.

Predicting the macroscopic chiral behaviour of liquid crystals from the microscopic chirality of the particles is highly non-trivial, even when the chiral interactions are purely entropic in nature. Here we introduce a novel chiral hard-particle model, namely particles with a twisted polyhedral shape and obtain a stable fully entropy-driven cholesteric phase by computer simulations. By slightly modifying the triangular base of the particle, we are able to switch from a left-handed prolate (calamitic) to a right-handed oblate (discotic) cholesteric phase using the same right-handed twisted particle model. Furthermore, we show that not only prolate and oblate chiral nematic phases, but also other novel entropy-driven phases, namely chiral blue phases, chiral nematic phases featuring both twist and splay deformations, chiral biaxial nematic phases with one of the axes twisted, can be obtained by varying particle biaxiality and chirality. Our results allow to identify general guidelines for the stabilization of these phases.

Phase diagrams of charged colloidal rods: Can a uniaxial charge distribution break chiral symmetry?

We construct phase diagrams for charged rodlike colloids within the second-virial approximation as a function of rod concentration, salt concentration, and colloidal charge. Besides the expected isotropic-nematic transition, we also find parameter regimes with a coexistence between a nematic and a second, more highly aligned nematic phase including an isotropic-nematic-nematic triple point and a nematic-nematic critical point, which can all be explained in terms of the twisting effect. We compute the Frank elastic constants to see if the twist elastic constant can become negative, which would indicate the possibility of a cholesteric phase spontaneously forming. Although the twisting effect reduces the twist elastic constant, we find that it always remains positive. In addition, we find that for finite aspect-ratio rods the twist elastic constant is also always positive, such that there is no evidence of chiral symmetry breaking due to a uniaxial charge distribution.

Cholesterics of colloidal helices: predicting the macroscopic pitch from the particle shape and thermodynamic state.

Building a general theoretical framework to describe the microscopic origin of macroscopic chirality in (colloidal) liquid crystals is a long-standing challenge. Here, we combine classical density functional theory with Monte Carlo calculations of virial-type coefficients to obtain the equilibrium cholesteric pitch as a function of thermodynamic state and microscopic details. Applying the theory to hard helices, we observe both right- and left-handed cholesteric phases that depend on a subtle combination of particle geometry and system density. In particular, we find that entropy alone can even lead to a (double) inversion in the cholesteric sense of twist upon changing the packing fraction. We show how the competition between single-particle properties (shape) and thermodynamics (local alignment) dictates the macroscopic chiral behavior. Moreover, by expanding our free-energy functional, we are able to assess, quantitatively, Straley's theory of weak chirality, which is used in several earlier studies. Furthermore, by extending our theory to different lyotropic and thermotropic liquid-crystal models, we analyze the effect of an additional soft interaction on the chiral behavior of the helices. Finally, we provide some guidelines for the description of more complex chiral phases, like twist-bend nematics. Our results provide new insights into the role of entropy in the microscopic origin of this state of matter.

Entropy-driven formation of large icosahedral colloidal clusters by spherical confinement.

Icosahedral symmetry, which is not compatible with truly long-range order, can be found in many systems, such as liquids, glasses, atomic clusters, quasicrystals and virus-capsids. To obtain arrangements with a high degree of icosahedral order from tens of particles or more, interparticle attractive interactions are considered to be essential. Here, we report that entropy and spherical confinement suffice for the formation of icosahedral clusters consisting of up to 100,000 particles. Specifically, by using real-space measurements on nanometre- and micrometre-sized colloids, as well as computer simulations, we show that tens of thousands of hard spheres compressed under spherical confinement spontaneously crystallize into icosahedral clusters that are entropically favoured over the bulk face-centred cubic crystal structure. Our findings provide insights into the interplay between confinement and crystallization and into how these are connected to the formation of icosahedral structures.