Experimental soft-matter science

Sidney R. Nagel The James Franck and Enrico Fermi Institutes and the Department of Physics The University of Chicago, Chicago, IL 60637

(Dated: January 31, 2017)

Soft materials consist of basic units that are significantly larger than an atom but much smaller than the overall dimensions of the sample. The label “soft condensed matter” emphasizes that the large basic building blocks of these materials produce low elastic moduli that govern a material’s ability to withstand deformations. Aside from softness, there are many other properties that are also caused by the large size of the constituent building blocks. Soft matter is dissipative, disordered, far-from-equilibrium, non-linear, thermal and entropic, slow, observable, gravity-affected, patterned, non-local, interfa- cially elastic, memory-forming and active. This is only a partial list of how matter created from large component particles is distinct from “hard matter” composed of constituents at an atomic scale. Issues inherent in soft matter raise problems that are broadly important in diverse areas of science and require multiple modes of attack. For example, far-from-equilibrium be- havior is confronted in biology, chemistry, geophysics, astrophysics and nuclear physics. Similarly, issues dealing with disorder appear broadly throughout many branches of inquiry wherever rugged landscapes are invoked. This article reviews the discussions that occurred during a workshop held on January 30-31, 2016 in which opportunities in soft-matter experiment were surveyed. Soft matter has had an exciting history of discovery and continues to be a fertile ground for future research. PACS numbers: 47.57.-s, 47.54.-r, 82.35.-x, 05.70.Ln, 46.65.+g, 83., 47.

CONTENTS Acknowledgments 24

I. Introduction: 1 A. Participants in the workshop held in Arlington VA on A. Characteristics of soft condensed matter 1 January 30-31, 2016 on which this report is based 24 B. Soft-matter systems and the disciplines that study 1. Scientific Advisory Committee 24 them 5 2. Attendees 24

II. Colloids 6 References 24

III. Granular Materials, (non-Colloidal) Suspensions, Emulsions and Foams 7 I. INTRODUCTION: IV. Glasses and Jamming 9 A. Characteristics of soft condensed matter V. Packing 11 The label “soft condensed matter” does not do justice VI. Liquid Crystals 12 to the broad scope of this branch of research. “Softness” VII. Polymers and Rheology 13 emphasizes the low elastic moduli that govern a mate- rial’s ability to withstand deformations. However, that VIII. Fluids 14 is only one attribute and many other properties go hand- IX. Mechanical Meta-materials 15 in-hand with the malleability of a material. arXiv:1610.10013v2 [cond-mat.soft] 30 Jan 2017 To understand some of the issues, it is relevant to be X. Active and Adaptive Matter 17 reminded of the common cause of the emergence of low XI. Self Assembly 18 elastic moduli in these soft materials. In colloids, poly- mers, foams, granular matter, emulsions and other soft XII. Patterns and Structure Formation 18 matter, the atoms are organized on a mesoscopic scale XIII. Simulations and Big Data 20 into entities that are much larger than an atom but still much smaller than the overall size of the material. These XIV. Instrumentation 21 large entities interact with one another to produce the elastic response of the material. However, because the XV. Outreach 22 individual constituents are large, the moduli must neces- XVI. Concluding Remarks 23 sarily be small. This is essentially a dimensional-analysis letters to nature

Experiments that cleanly reveal the nature of interactions that this street is shed starting from the free end rather than from the between a deformable body and a ¯ow are dif®cult to perform, top, where the ¯ow ®rst encounters the ®lament. Despite the with complications including the control of three-dimensional ¯ow oscillations of the street, the ®lament remains a static, seemingly effects and the challenge of ¯ow visualization. Soap ®lm is a rigid body. Indeed, transient small disturbances (such as a small convenient experimental system that resembles two-dimensional displacement at the free end, or a thin vortex wake produced by a hydrodynamics in many respects3±7. In this experiment, we take cylinder placed upstream) do not destroy the stability of this state. advantage of the relative simplicity of this system, and use the However, if the external perturbation is suf®ciently large, the ¯owing soap ®lm as a laminar two-dimensional ¯ow tunnel. The ®lament jumps to a stable ¯apping state (Fig. 2b), whose persistence experimental set-up is shown in Fig. 1. The horizontal span of the is quite robust to perturbation. Unlike a simple pendulum, the tunnel is 8.5 cm, the typical ®lm thickness 3±4 mm, and ¯ow speeds undulation is well ®tted by a travelling harmonic wave with a 200±300 cm s-1 (see Methods). The experiment is performed 80 cm spatially varying envelope: below the ¯ow entrance, at the midline. The ®laments are 0.15 mm y x; t†ˆf x†sin 2pvt ‡ 2px=l† in diameter and 2±6 cm long. Their linear density is about 2:0 3 10 2 4 g cm 2 1. The rigidity of the ®lament is estimated to be Here, y x; t† is the horizontal displacement of the ®lament from the the order of 0.1 erg cm, tested in a steady soap ®lm where the centre-line, measured at a vertical distance x from the ®xed point for de¯ection is measured under a known force. The ¯uid wets the time t. f(x) is a spatial envelope function (increasing monotonically ®lament, with surface tension forces constraining the ®lament to lie from the ®xed point), v is the ¯apping frequency and l the always in the plane of the ®lm. The length of the ®lament can be wavelength. When the ®lament length is about 3 cm long and changed incrementally through the holder at the upstream end, ¯ow speed 280 cm s-1 (the corresponding Reynolds number is of which abuts the ®lm perpendicularly. the order of 104), we observed a typical ¯apping frequency of about At a ®xed ¯ow rate, there exists a critical ®lament length Lc, below 50 Hz and an amplitude (the total excursion of the free end) of which the ®lament is stretched straight and aligned with the ¯ow. about 1.5 cm. A thin vortex street is still evident, but is strongly Regardless of the magnitude of an externally imposed disturbance, modi®ed by the large-scale ¯apping motion. The successive small after a few oscillations, the ®lament quickly returns to this rest state. eddies produced by a single stroke are now of a single sign, unlike If, however, the length of the ®lament exceeds Lc, the system the alternating signs in the stretched-straight case. This local vortex becomes bi-stable, and two different stable dynamical states are observed. One remains the stretched-straight state, as when L , Lc. 2 Figure 2a shows such a state (exposure time 0.5 ms). Using a monochromatic light source, the photograph shows instantaneous interference patternsargument that reveal (Frenkel, small 2002). thickness The variations. bulk and A shear moduli narrow von KaÂrmaÂofn vortex a material `street' have8,9 (that the is, units a train [energy]/[volume]. of vortices of While alternating signs) canthe be interaction seen trailing energy behind between the ®lament. two constituent We note particles is perhaps not exceedingly different from that between individual atoms, the volume of each particle is large. Thus the moduli are orders of magnitude smaller than they would be for an atomicH or molecular crystal. Thus materials made up of such large basic units must have specialS attributes. However, even a cursory exam- ination suggests that there are many other properties, aside from low elastic moduli, that are also caused by the large size of the component building blocks and that would make such materials distinctive. DissipativeTT matter: Because the individual building blocks are large, they each have many low-energy degrees of freedom (e.g., vibrations) that can be easily activated when two particles collide or scrape past one another. This generatesF strong dissipation during flow or defor- mation, either due to frictional interactions or to the in- elasticityNN of collisions. Kinetic energy, for example due to sound vibrations, is rapidly absorbed into the internal degrees of freedom of the constituents. Thus granular or colloidal materials are highly damped systems. Disordered matter: Because the individual particles TT are large, it is prohibitively unlikely that any two of them will have precisely the same number of atoms in the same arrangement. Thus there is disorder in the size and shape of the constituents. (One counter-example is Figure 2 The two stableFIG. states 1 The of the oscillatory ®lament. When flapping the ®lament state length of a is filament greater than in a flowing matter made from viruses,L where the DNA regulatesthe critical the value, it cantwo-dimensional stay in either a stretched-straight soap film. state Flow or an structures, oscillatory state. modulated The by structure precisely.) Even this small amount of¯ow disorder is visualized usingthe interference moving filament, patterns under are a advectedsodium lamp. downstream.a, The ®lament is The flow is Figure 1 Experimental set-up.can Ahave vertical dramatic soap ®lm is fed effects. from the top For (H) with example, a solution ofin astretched granular straight invisualized the ¯ow. Von KausingÂrmaÂn typeinterference vortices are shed imaging. from its free From end. b (Zhang, The et al., 1.5% Dawn, a commercialmaterial brand of dishwashing consisting liquid of soap, hard in distilled grains, water. forces The willoscillatory predomi- ¯apping2000). state. Flow structures, modulated by the ®lament, are advected pressure head is ®xed usingnantly an over¯ow propagate mechanism. through The attached the nozzle material (S) regulates in pathsdownstream. where Interference fringes run in the ¯ow direction, indicating that the soap ®lm is the ¯ux. Two nylon threadsslightly under tension greater (T), spreadoverlap from between the nozzle, particles run parallel, is and present.under little Par- stretching and compression. c, The free end of the ®lament shows a `®gure of then converge into the receivingticles container that have (L). The asperities upper end of change the silk ®lament the form (F) is ®xed of theeight' interac- trajectory owingOne to the class existence of offar-from-equilibrium the travelling waves associated behavior with the ¯apping occurs when using a thin tube, whichtion is introduced between normal particles to the ®lm andplane. dominate frictionalmotion. contacts. continual energy input pushes the system into a dynamic But this is just the beginning. Because the energy to regime where new structures form; relaxation is insuffi- 836 overcome a barrier between two© nearby 2000 Macmillan configurations Magazines LtdcientlyNATURE rapid| toVOL allow 408 | 14 the DECEMBER unstressed 2000 | www.nature.com state to be an ap- is small, it is easy for even an initially ordered pack to propriate starting point for describing the characteristic lose all remnants of long-range order. Thus much of the structure and dynamics of the material. A classic exam- physics of dealing with these systems must confront is- ple of such a situation is the turbulent flow of a liquid. sues of glassiness: slow and non-exponential relaxations, When the boundary of a fluid is flexible, startlingly strik- rugged energy landscapes, heterogeneous dynamics (es- ing images emerge that illuminate the underlying dynam- pecially near the glass transition), and extra modes of ics. Figure 1 (Zhang et al., 2000) shows the interaction excitation not present in crystalline matter. of a thin filament with a two-dimensional flow. Pattern Non-linear matter: Because the component particles formation is often another consequence of energy input are large, it is easy to drive a soft material out of its overtaking the ability of the medium to relax. linear-response regime. Thus at low forces, it is pos- Another class of far-from-equilibrium behavior occurs sible to deform a soft material and even to change its when thermal energy is insufficient for exploring the rel- nature from that of a rigid solid into a rearranging fluid. evant phase space. The glass transition is perhaps the Descriptions of the excitations within a rigid material prototypical example of how a system develops such be- must often confront anharmonic effects especially when havior (Berthier and Biroli, 2011; Debenedetti and Still- the system is on the verge of a rearrangement event where inger, 2001; Ediger et al., 1996). In a supercooled liquid, non-linear behavior can be dominant. relaxation times become exceeding long as the tempera- Far-from-equilibrium matter: For the reasons men- ture is lowered. When these relaxation times exceed the tioned above, soft matter is often far from equilibrium. time that any experiment can practically explore, the sys- R EPORTS

glass transition (8, 26). To better characterize moving particles: For ⌬t*when␣2 is a by observing the samples in 3D for several the ␣-relaxation, we determined the distribution maximum, the fastest particles are precisely the hours. The distribution of cluster sizes for a of particle displacements P[⌬x(⌬t)] (Fig. 3). particles contributing to the tails of P[⌬x(⌬t*)], given volume fraction is broad (Fig. 5A); –␮ Although this distribution is gaussian for purely thus making ␣2 large (Fig. 3). We chose a P(nc) ϳ nc with ␮ϭ2.2 Ϯ 0.2 for the diffusive particles, it is expected to be consid- cutoff ⌬r*foragivensamplesuchthatover supercooled fluids, similar to the value ␮ϭ erably broader near the ␣-relaxation (11, 18, time, 5% of the particles had displacements ͉⌬ជr͉ 1.9 Ϯ 0.1 seen in simulations (14). An expo- 19). Deviations from a gaussian are quantified Ն ⌬r*(11, 14, 29), although at any given time, nent ␮Ͻ3impliesthatquantitieswhichde- 2 by a nongaussian parameter the exact fraction may not be 5%. On average, pend on ͗nc͘,suchasaverageclustersize,will ͗⌬x4͘ these particles had moved five times farther be dominated by the largest clusters; thus, struc- ␣ (⌬t) ϭ Ϫ1(1)than the ensemble of particles (⌬r*/͌͗x2͘ Ϸ tural relaxation occurs because of a small num- 2 3͗⌬x2͘2 5). The 5% most mobile particles were also ber of large clusters of cooperative fast parti- the simplest combination of the second and examined in simulations (14), but the cutoffs cles, rather than many individual fast particles fourth moments of a 1D P[⌬x(⌬t)], which is ⌬r*weretypicallylargerthantheparticlera- moving independently. It is likely that the dis- zero for a gaussian distribution (27). Broader dius a,whereasforourdata,⌬r*istypically tributions of cluster sizes are even broader than

distributions result in large values of ␣2. As 0.4a –0.8a.Thedifferencemaybeduetothe indicated because the largest clusters extend out in Fig. 1B, for supercooled fluids (open cir- binary size distribution or to the different par- of the viewing volume.

cles), ␣2 is largest for lag times correspond- ticle interaction potential used in the simula- The cluster size increased dramatically as ␾ ing to the end of the cage-trapping plateau in tions. To look for spatial correlations of these increased (Fig. 5B) (31), consistent with the the MSD. We see a clear indication of the fast particles, we constructed the 3D Delaunay increased size of the cooperatively rearranging approach to the glass transition in the rise of triangulation of the particle positions (30), regions of the Adam and Gibbs hypothesis (6).

the peak value of ␣2 as ␾ increases toward which provides the nearest neighbor connectiv- Moreover, there is a pronounced drop in the ␾G. Note that the magnitude of ␣2 may be ity, and we identified the clusters of connected average cluster size at ␾G Ϸ 0.58 (vertical increased due to the slight polydispersity (5% fast particles. dashed line in Fig. 5B), which signifies the by radius) of the particles (17). For the supercooled fluid, the fast particles onset of the colloidal glass transition. on July 5, 2016 A dramatic change in the behavior of ␣2 were strongly spatially correlated and exhibited The larger clusters are generally extended occurred at ␾Ϸ0.58; at lower ␾, ␣2 exhib- large extended clusters (Fig. 4A). This result is structures (Fig. 4); thus, we plot the number of ited a distinct peak near the ␣-relaxation, adramaticdemonstrationthatthe␣-relaxation particles in a cluster against the cluster’s radius

whereas at higher ␾, the peak in ␣2 was much in colloidal fluids occurs by means of cooper- of gyration Rg (Fig. 5C). The power law scaling broader but not as high. We identify this ative particle motion: when one particle moves, observed is indicative of fractal structure with a

sharp change as the glass transition and de- another particle moves by closely following the fractal dimension df ϭ 1.9 Ϯ 0.4 for all volume termine ␾G ϭ 0.58 Ϯ 0.01, in agreement with first (12, 13, 19, 20). We calculated the angles fractions, comparable to the preliminary value previous work (18, 21). For glasses (closed between displacement vectors of neighboring df Ϸ 1.75 seen in the simulations (14). We symbols), ␣2 drops at longer lag times, even fast particles; the distribution of these angles is further characterized the structure by measuring at lag times when the MSD begins to rise strongly peaked at 0° and shows that neighbor- the number of neighbors, Nf,ofeachfastparti- (19). The upturn in the MSD at longer lag ing particles move in parallel directions. More- cle. For supercooled fluids, the distribution

times for the glasses has been seen in other over, we found that the displacement vectors P(Nf)exhibitedabroadpeak,withϳ10% of the http://science.sciencemag.org/ experiments (18, 19, 28) and may be due to are more likely to point toward other fast par- particles having Nf Ն 7, indicating dense re- activated processes (1). ticles than elsewhere, confirming that the mo- gions of cooperating particles (32). A typical To study structural relaxations in super- tion is cooperative. P(Nf)forfluidsisshownbytheopensymbolsin cooled fluid samples, we examined the fastest We characterized the nature of these clusters Fig. 5D. 3

tem falls out of equilibrium and becomes a glass. While Downloaded from the glass transition occurs in nearly all materials, it plays a special role in polymer science where entanglements can help to suppress the ability to crystallize and freeze the system into far-from-equilibrium configurations. Many attempts have been made to describe glassy and other far-from-equilibrium materials as having a fictive tem- perature. However, it is not at all clear that a single Fig. 3. Distribution function P(⌬x) for ␾ϭ0.56, temperature can adequately describeat ⌬t* allϭ 1000 aspects s, corresponding of to the the peak in ␣2(⌬t) (Fig. 1B). The dashed line is the best fit dynamics of a glassy material. gaussian, and the solid line is a fit of a stretched Fig. 4. TheC locations of the fastest particles (large spheres) and the other particles (smaller spheres). exponential to the tails of the distribution [P ϳ The spheres are drawn smaller for clarity; the particles all have the same physical size, which is the Thermal and entropic matter: Because␤ the energy exp(Ϫ͉x/x0͉ ) with ␤ϭ0.8; we found that size of the large spheres shown in this figure. (A) “Supercooled” sample with ␾ϭ0.56, ⌬t* ϭ 0.8 Ͻ␤Ͻ1.5 for different choices of ␾ and ⌬t. 1000 s; the fastest particles had a displacement Ͼ0.67 ␮m. The red cluster contained 69 particles; necessary to deform a soft materialSmaller values is small, of ␤ coincide room- with larger values the light blue cluster contained 50 particles. (B) “Glassy” sample with ␾ϭ0.61, ⌬t* ϭ 720 s; the of ␣2]. The data within the dotted lines are the fastest particles had a displacement Ͼ0.33 ␮m. The largest cluster (red) contained 21 particles. The temperature thermal excitations areslowest often 95%; particles sufficient in the fastest to 5% have “speed” of a particle was determined over a time ⌬t* corresponding to the ␣-relaxation for (A) and rearrange the system. Indeed, for͉⌬ soft-matterx͉ Ͼ 0.2 ␮m. systems, the ␤-relaxation for (B); see text for details. the appropriate energy scale is often given not in units of www.sciencemag.org SCIENCE VOL 287 28 JANUARY 2000 629 eV , but in units of kBT . Just by waiting at room tem- perature under an applied stress, a soft-matter system creeps and deforms. For many of the same reasons that temperature plays a significant role, entropic considerations have an increased importance in determining the structure of soft matter. In many cases, ordered phases of matter are created due FIG. 2 Fastest particles in (A) a supercooled and (B) a to the role of entropy. Liquid crystals, for example, order glassy colloid. In (A) red (lower central) and blue (upper into novel phases precisely because creating order in one left) clusters contain 69 and 50 particles respectively. In (B) degree of freedom (such as molecular orientation) can al- the largest cluster (red at the upper left corner of the box) low many more configurations in others (such as transla- contains 21 particles. In both (A) and (B) the particles that are not moving as rapidly are shown with smaller size to aid tions). Understanding the variety of order created out of visualization. From (Weeks et al., 2000). (C) Motion in a entropic interactions has given rise to an entire industry. glass composed of ellipsoidal particles. Both rotational and Slow matter: The diffusion constant for a particle is translational motion can be determined. Colored (shaded) inversely proportional to its linear size so that the time particles in the image shows the spatial distributions of the to diffuse its own length, L, increases as L3. Therefore fastest 8% of translational motions of the ellipsoids. From large particles, such as colloids, diffuse much more slowly (Zheng et al., 2011). than ones of atomic-scale. Large particles also accelerate more slowly than less massive ones. These are just two reasons why, for example, colloids and granular materials are affected by Earth’s gravity, g. A container filled with have such long relaxation times. Even if a material could sand grains of size L ≈ 1mm and density ρ ≈ 1g/cm3 reach equilibrium, it would take much longer to do so if is one example where the room-temperature thermal the comprising particles are large. energy, kBT , is approximately 12 orders of magnitude 4 Observable matter: One advantage of being slow is that smaller than ρgL , the energy of an elementary rear- soft matter allows observation of processes that would rangement caused by raising one particle over its neigh- otherwise be too rapid to track in more conventional sit- bor. As with other far-from-equilibrium phenomena, the uations. Moreover, large size allows direct visualization question arises: if temperature is no longer sufficient to of behavior that otherwise could not be perceived with allow equilibration, how does one describe the state of a clarity. For example, the large individual particles and sandpile? What is the ensemble over which one should the slow dynamics of colloidal systems has made them average in order to obtain a meaningful description? ideal for observing slow relaxation in glass-forming fluids Patterned matter: Because soft matter is often far from as shown in Fig. 2 (Weeks et al., 2000; Zheng et al., 2011) equilibrium, it is prone to form a wide array of distinct and the vibrational excitations in disordered solids (Chen patterns. As mentioned above, turbulence occurs as a et al., 2010). On a different front, because the phenom- low-viscosity liquid is forced to flow at high speeds. As ena appearing in soft-matter systems is often ubiquitous, the energy cascades from the large length scales where counter-intuitive and easily observable, these systems are it is injected into the system down to a smaller (Kol- ideal for demonstrating scientific enquiry to a general mogorov) scale where it is dissipated, new structures are public and for engaging younger students just starting created, such as eddies and entangled vortices, that are to think about a career in research. not observed in the motionless ground state. Gravity-affected matter: Some, but not all, soft-matter Other patterns occur between the constituent phases systems are naturally far-from-equilibrium because they in a material and display intricate and varied forms in VOLUME 52, NUMBER 16 PHYSICAL REVIEW LETTERS 16 APRrL 1984

contrast, the aggregation occurs on much longer give an estimate of the Hausdorff dimension, time scales which can be controllably varied from D = 1.7 + 0.1. —30 min to several weeks. The internal structure of each cluster also exhi- We examine the structure of the clusters using a bits scale invariance. %e show this by recording transmission electron microscope (TEM), and en- the center coordinates of each individual gold parti- sure that no additional aggregation occurs as the cle and calculating the point-to-point correlation liquid is evaporated from the highly hydrophyllic functions of the images of the clusters, c (r) =— grid. However, awhile the aggregates are three- (p(r)p(0))/(p(r)). Here p(r) is 1 at a particle dimensional in solution, edge-on examination with center and 0 elsewhere, and the average is over po- the TEM shows that upon drying they collapse to sition and orientation of the cluster. In our form nearly flat, two-dimensional structures. By analysis, we account for the finite size of the clus- examining the full area of the TEM grid, clusters of ters and correct for edge effects by finding the widely varying sizes can be found and photo- smallest circle that just encloses each cluster. Then, graphed, representing a snapshot in time of the ag- centered on each individual gold particle we find gregating colloid. A typical picture of a relatively the largest circle which fits inside the outer circle, large aggregate is shown in Fig. 1. It exhibits the and include only other particles within the small cir- tenuous, ramified appearance characteristically ob- cle in the calculation of c(r). This ensures that for served for colloidal aggregates. all directions there is an equal probability of finding To illustrate the scale-invariant nature of the ag- a second particle a distance r from the first, but lim- gregates as they grow, and to obtain an estimate of its the maximum extent over which we can measure their Hausdorff dimension, we measure the number correlations to less than one quarter the cluster size. of particles in a cluster, N, as a function of the size For a two-dimensional projection of these frac- L. Since there are relatively few overlapping parti- tals, we expect6 c(r) —r for r « L, where cles, even in the largest cluster, N can be measured a = 2 —D. A representative series of correlation with high precision. As a measure of the size, we functions calculated for clusters of increasing size is take the geometric mean of the longest linear shown by the solid lines in Fig. 3. For the very dimension of the cluster and the length perpendicu- small clusters, the statistics are poor and the curves lar to that axis. A logarithmic plot of N versus L is are noisy. However, once the number of particles shown in Fig. 2 for nearly 100 different clusters on in the cluster becomes greater than about 1000, the the grid. The data are very well described by a correlation functions exhibit an extended linear re- power-law behavior, as expected for a fractal. 6 gime, with the slopes essentially independent of Similar results were obtained for pictures from oth- cluster size. The limiting slope, shown by the er grids, and linear least-squares fits to the data sets dashed line in Fig. 3, yields n= —0.23 +0.1, where 4 the error reflects the cluster to cluster variation. The resulting Hausdorff dimension is D = 1.77

0 2 3 4 Ln L

FIG. 2. N vs L, where the solid line is a least-squares fit to the data, with the slope giving D —1.75. Wean be converted to mass by multiplying by the mass of a single ' FIG.FIG. 41. TEMTEM imageimage of goldof colloidtypical aggregategold colloid showingaggregate. DLA gold particle, —10 g, while L can be converted to structure. From (Weitz and Oliveria, 1984). This cluster contains 4739 gold particles. nanorneters by multiplying by 14.5. 1434 et al., 1985). Elastic response and elementary excita- FIG. 3 Branching patterns seen in the limbs of trees on top tions are governed to a large extent by pressure variations and in dielectric breakdown of a Lichtenberg figure below. at the surface of a jammed solid (Silbert et al., 2009). Surface tension in a dense suspension governs many as- pects of rheology (Brown and Jaeger, 2012). Surfactants both space and time. Classic instabilities, such as the (large molecules with a polar head and a hydrophobic viscous-fingering instability (Saffman and Taylor, 1958) tail) change the surface tension at a fluid surface dra- or the Belousov-Zhabotinsky reaction (Belousov, 1959; matically. They can produce surface flows and create a Zhabotinsky, 1964), produce mesmerizing and detailed myriad array of structures. Fracture is the creation of features that progressively unfold in time. Another com- an internal surface in a material under stress. It too is mon pattern based on dilation symmetry for the pene- transformed in soft materials into a novel form of material tration of space is shown in Fig. 3. Nature appears to be failure. An effect commonly used in liquid crystals is to showing off the tricks that it can do! use a surface to anchor the orientation of the anisotropic Non-local matter: When patterns form, they cannot molecules. This provides one means to control the defect always be described by local rules. Thus, for example, a structures in the order parameter. branch in a diffusion-limited aggregation (DLA) pattern Memory-retaining matter: Matter that is malleable as shown in Fig. 4 can be shielded from further growth can be altered by repeated applied stresses. Because by other parts of the interface that are far away when the material would typically be disordered by such treat- measured along the interface itself (Witten and Sander, ment, this manipulation can lead to different internal 1981); when those segments circle back and produce an structures that are a signature of how the material has overhang they effectively prevent other particles from ap- been treated. That is, the material can retain a mem- proaching and attaching themselves to the initial pro- ory of how it was formed or previously manipulated. The tected branch. manner in which such memories can be encoded, read out Interfacial elastic matter: In hard-matter solids, in- and erased can vary from system to system. Sometimes ternal and external interfaces are exceedingly important the memory effects can be very subtle. for controlling the electronic properties of devices. Inter- Active matter: Because the individual ingredients of a faces, both internal and external, also play an important soft-matter material are large, it is comparatively easy role in soft-matter. Because soft matter is deformable, to have these large entities inject energy at a local scale. elastic deformation at a surface is transmitted efficiently This can be done, for example, by creating different sur- into the bulk. Flow in microfluidic channels and in Hele- face activity on different parts of a particle. When ex- Shaw cells are affected by surface topography (Ben-Jacob posed to an external optical beam or when placed in con- LETTER RESEARCH

It is well known that with increasing concentration, rod-like mole- of microtubule bundles in 3D active networks. Once created, a fracture cules undergo a transition to a nematic liquid-crystalline phase. To line terminated with a pair of oppositely charged disclination defects. explore this regime we created a flat 2D oil–water interface stabilized Eventually, even as the fracture self-healed, the defects remained with a surfactant, with BANs being dispersed in aqueous phase. Over unbound and streamed around until eventually annihilating with time extensile microtubule bundles adsorbed onto the PEG brush oppositely charged defects (Supplementary Video 6). The rates of formed by the surfactant molecules, eventually covering the entire defect creation and annihilation were balanced, creating steady-state surface with a dense liquid-crystalline monolayer of locally aligned streaming dynamics that persisted for many hours. These observa- bundles. The adsorbed 2D layer constituted an active microtubule tions exemplify how active nematics are fundamentally different from liquid-crystalline phase, characterized by fast streaming flows and defect equilibrium ones, in which fractures, internal streaming flows and unbinding (Supplementary Video 6). Further information about the spontaneous unbinding of defect pairs are never observed. nature of active microtubule liquid crystals can be garnered by exam- In a biological context, active fluids are frequently confined to the ining the structure and dynamics of their defects. In general, active liquid cytoplasm and it has been proposed that such confinement leads to crystals can have either nematic symmetry24,25 (as found in monolayers emergence of coherent flows that can enhance cellular transport, a of amoeboid cells or vertically shaken monolayers of granular rods26,27) phenomenon known as cytoplasmic streaming30,31. For this reason or polar symmetry (which is found in motility assays at high concentra- we encapsulated BANs in aqueous droplets emulsified in fluorinated tions5,28). Polar liquid crystals form vortex and aster defects5,29. oil. When squeezed between two surfaces, such water-in-oil active In contrast, we found that active microtubule liquid crystals form droplets exhibited an unforeseen emergent property: persistent auto- disclination defects of charge 1/2 or -1/2, implying the presence of nomous motility. This motility was highly robust, limited only by the nematic symmetry (Fig. 3a-c). This is expected because the basic building sample lifetime, which was typically a few days (Fig. 4a, Supplementary blocks of these materials are symmetrically extensile microtubule Video 7). Instead of moving along straight lines, the active emulsion bundles. In equilibrium liquid crystals, defects are largely static struc- droplets preferentially moved in periodic patterns, with their average tures whose presence is determined by either internal frustrations or velocities reaching up to ,1 mms21. In comparison, droplets without external boundary conditions. In contrast, defects in active micro- any chemical fuel did not exhibit any motion (Fig. 4b). tubule nematics exhibited unique spatiotemporal dynamics. They To elucidate the microscopic mechanism that drives droplet moti- REPORTS lity, we imaged the internal microtubule dynamics. For small droplets were created as uniformly alignedcoherent nematic motion, and giant domains density fluctuations are extend,g is the acceleration buckle due to gravity, varies between experimental system thus has all the physical and internally self-fracture (Fig. 3d),intimately similar related consequences to the of the dynamical orientation- G =4andG cascades=7.Wevariedthetotalnumberof(lessingredients than of an30 activem nematicm), (1– microtubule4). bundles extended and pushed against al order that develops in a sufficiently dense particles in the cell, Ntotal,between1500and The system is in a very dynamic steady state, grouping of self-driven objects with anisotropic 2820. Ntotal in each instance was counted by with particle motion (movie S1) organized in body shape. This has substantial implications for hand. The area fraction, f,occupiedbythepar- macroscopic swirls. Swirling motions do not biological pattern formation and movement ecol- ticles is the total projected area of all the rods necessarily imply the existence of giant number ogy (6): The coupling of density fluctuations to divided by the surface area of the cell. f varies fluctuations (20, 21); however, particle motions alignment of individuals will affect populations as from 35% to 66%. Our experimental system is in our system generate anomalously large fluc- addiverse as herds of cattle, swarms of locusts (7), similar to those used to study the phase behaviora tuations in density. Figure 1A shows a typical in- b schools of fish (8, 9), motile cells (10), and of inelastic spheres (16, 17). Galanis et al.(18) stantaneous configuration, and the Fig. 1B inset filaments driven by motor proteins (11–13). shook rods in a similar setup, albeit with much showsthe orientational correlation function G2(r)= We report here that persistent giant number less confinement in the vertical direction. The 〈cos2(qi − qj)〉,wherei,j run over pairs of particles fluctuations and the coupling of particle currents particles were imaged with a digital camera (19). separated by a distance r and oriented at angles qi to particle orientation arise in a far simpler driven The rods gain kinetic energy through frequent and qj with respect to a reference axis. The angle system, namely, an agitated monolayer of rodlike collisions with the floor and the ceiling of the brackets denote an average over all such pairs and particles shown in (14)toexhibitliquidcrystal- cell. Because the axes of the particles are almost about 150 images spaced 15 s apart in time. The line order. These fluctuations have also been always inclined to the horizontal, these collisions data in the inset show that the systems with Ntotal = observed in computer simulations of a simple impart or absorb momentum in the horizontal 2500 and Ntotal =2820displayquasi–long-ranged model of the flocking of apolar particles by Chaté plane. Collisions between particles conserve mo- nematic order, where G2(r)decaysasapowerof et al.(15). The rods we used were cut to a length mentum but also drive horizontal motion by con- the separation, r.Ontheotherhand,thesystem l = 4.6 ± 0.16 (SEM) mm from copper wire of verting vertical motion into motion in the plane. with Ntotal = 1500 shows only short-ranged diameter d = 0.8 mm. The ends of the rods were Interparticle collisions as well as particle-wall nematic order, with G2(r)decayingexponentially etched to give them the shape of a rolling collisions are inelastic, and all particle motion with r.Detailsofthecrossoverbetweenthesetwo pin. The rods were confined in a quasi-two- would cease within a few collision times if the behaviors can be found in [Supporting Online b dimensional cell 1 mm tall and with a circular vibrations were switched off. The momentum of Material (SOM) text]. Autocorrelations of the cross-section 13 cm in diameter. The cell was the system of rods is not conserved either, be- density field as well as of the orientation of a mounted in the horizontal plane on a perma- cause the walls of the cell can absorb or impart tagged particle decay to zero on much shorter time on January 1, 2017 nent magnet shaker and vibrated vertically at a momentum. The rods are apolar; that is, individ- scales (SOM text), so we expect these images to frequency f =200Hz,withanamplitude,A, ual particles do not have a distinct head and tail be statistically independent. To quantify the between 0.025 and 0.043 mm. The resultant that determine fore-aft orientation or direction of number fluctuations, we extracted from each dimensionless acceleration G =(4p2f 2A)/g,where motion and can form a true nematic phase. The image the number of particles in subsystems of 5

known example of the immense industrial impact of soft matter. Likewise, dissipation has been a bane to the

efficient processing of materials in industries that manip- ulate particulate matter. The cosmetic,http://science.sciencemag.org/ paint and food industriescdc rely on the rheology of structured fluids. Com- plex flows in polymers lead to dramatic sensitivityd of ma- terial structure and properties to processing conditions.

The different attributes affordedDownloaded from by soft matter raise important problems in basic science that need to be addressed from multiple perspectives. Some issues cut c across all areas of science; soft matter allows us to ad- dress common problems from a fresh perspective. For ex- ample, matter that is far-from-equilibrium is confronted in biology, chemistry, geophysics, astrophysics and nu-

B clear physics. Even a small increment in our understand- Fig. 1. Giant number fluctuations in active granular rods. (A)Asnapshotof of various sizes. Theing number could fluctuations have in each enormous subsystem are determined ramifications. Similarly, issues the nematic order assumed by the rods. There are 2820 particles (counted by from images taken every 15 s over a period of 40 min (19). The squares hand) in the cell (area fraction is 66%) being sinusoidally vibrated represent the systemdealing shown in (A). with It is adisorder dense system where appear the nematic broadly in many scientific perpendicular to the plane of the image, at a peak acceleration of G =5. order is well developed. The magnitude of the scaled number fluctuations The sparse region at the top between 10 and 11 o'clock is an instance of a decreases in moreendeavors. dilute systems, where Even the nematic if order not is weaker specifically (SOM addressing disorder, large density fluctuation. These takeseveralminutestorelaxandform text). Deviations from the central limit theorem result are still visible at an area elsewhere. (B)Themagnitudeofthenumberfluctuations(quantifiedbyDN fraction ≅ 58% (diamonds)wherever but not at rough an area fraction energy≅ 35% landscapes (circles). (Inset) are invoked (e.g., high- and normalized by pN)againstthemeannumberofparticles,forsubsystems The nematic-orderenergy correlation function physics as a functio andnofspatialseparation. string theory, biological cell differen- ffiffiffi tiation) much is to be learned by fully understanding the 106 6 JULY 2007 VOL 317 SCIENCE www.sciencemag.org energy landscape underlying glasses and the glass transi- tion. One can make similar claims about the deep issues Figureconfronted 4 | inMotile many of water-in-oil the other headings emulsion on that droplets.list. In a, Droplets containing extensilesome cases, microtubule these materials bundles have provided exhibit a novel spontaneous exper- autonomous motility when partiallyimental window compressed for studying between common chamber phenomena surfaces. that A droplet trajectory taken over could not be easily undertaken in atom-scale systems. Figure 3 Dynamics of 2D streamingFIG. 5 nematics A) Snapshot confined of nematic to order fluid assumed interfaces. by rods vi- a time interval of 33 min is overlaid onto a bright-field droplet image. Scale bar | brated perpendicular to the plane of the image showing giant a, Schematic illustrations of the nematicnumber fluctuations. director configurationFrom (Narayan et al.around, 2007). B) Active for a and b is 80 mm. b, In the absence of ATP, passive droplets exert no internal disclination defects of charge 1/2liquid (left) crystals and 2 exhibiting1/2 (right). buckling,b, c folding, Active and liquidinternal frac- forces,B. Soft-matter and the systems only and contribution the disciplines to that their study movement them is minor drift. ture. The two arrows indicate the positions of two topological crystals exhibit disclinations of bothcharges 1/2 associated (top) and with2 the1/2 disclination (bottom) defects charge, in the order c, FluorescenceIn addition to theimage types of of active behavior microtubule displayed, a bundles list that spontaneously adsorb indicating the presence of nematicparameter. order. Scale The scale bars, bar 15 ism 20m.µmd., From A sequence (Sanchez et of al., onto the oil–water interface. The resulting active liquid crystalline phase 2012). of the different classes of systems encompassed by the images demonstrates buckling, folding and internal fracture of a nematic exhibitslabel “soft streaming matter” is flows, also large. indicated Above, with we mentioned blue arrows. The red arrow indicates domain. The fracture line terminates with a pair of oppositely charged thea few direction sample systems: of instantaneous liquid crystals, droplet colloids, velocity. foams, The image is focused on the granular matter, suspensions and polymers. There are tact with a surrounding fluid, the different ends of the disclination defects (red arrow tracks 1/2 disclination; blue arrow tracks 21/2 dropletother classes surface of material that is thatin contact share a withstrong the intellectual coverslip. Scale bar, 100 mm. d, Image particle can undergo different interactions that will pro- base with soft matter. High on this list would be simple disclination). After the fracture linepel themself-heals, through the the disclinationsolvent in which pair they remains are immersed. of the droplet taken at a midplane indicates that the droplet interior is largely liquids. Liquids are disordered atomically; they are dis- unbound. Time lapse, 15 s; scale bar,When 20 theremm. are a large number of such interactive swim- devoid of microtubule bundles. Scale bar, 100 mm. sipative (viscosity is important); they are non-linear and mers, they can produce dramatic swarms of activity and can easily form observable patterns; they are soft (but structure (Marchetti et al., 2013). Two examples of ac- for a reason different15 from NOVEMBER that for soft matter 2012 based | on VOL 491 | NATURE | 433 tive systems are illustrated in Fig. 5 (Narayan et al., 2007; large constituent particles) with a vanishing static shear Sanchez et al., 2012).©2012 Macmillan Publishers Limited. All rights reserved modulus; their interfaces are easily deformed and deter- mine much of a fluid’s behavior; they can be pushed far This list shows an impressive number of ways in which from equilibrium either by constant energy injection, as matter created from large component particles is dis- in turbulent flow, or by super-cooling into a glass phase tinct from matter composed of simpler constituents at where thermal energy is insufficient to allow adequate an atomic scale. Softness is but one of many properties exploration of phase space. Thus, one home for fluid dy- distinguishing these materials from their “hard-matter” namics is definitely in the soft-matter community. counterparts. Indeed, softness is perhaps one of the least Another example that shares common features with interesting attributes listed above. soft matter is found in living systems. Biological matter Each of the characteristics mentioned can allow spe- is soft. It is often active with energy input on a mi- cial functionalities thus creating technological opportu- croscopic level. Much of biological activity, though pur- nities. The industry spawned by the new forms of liquid- posive, is also seemingly disordered and issues such as crystalline order is perhaps the most obvious and well- crowding and jamming arise. Living biological matter is 6 the ultimate far-from-equilibrium system. Soft matter, example of a colloid is milk. In a colloidal suspension the with its interest in issues tied up with disorder, patterns, particles are too small to settle appreciably under the structure formation, non-linear and far-from-equilibrium influence of gravity but are strongly affected by thermal behavior, is a close cousin to biological matter and is ar- energies. In other non-colloidal suspensions, the particles guably the closest intellectual discipline to biology that are larger; they can sink and are less prone to thermal exists in the physical sciences. rearrangements. An emerging activity with a home in the soft-matter Because they do not settle over time, colloidal suspen- community deals with meta-materials. The aim is to or- sions can mimic the properties of atomic-scale fluids but ganize structures at a macroscopic level in order to create at a much larger size and much slower rates as shown, for unique material properties and functionality. Included in example, in Fig. 2. They thus allow the observation of this endeavor are the intellectually rich fields of origami organizational principles that govern the many-body in- and kirigami based on the folding and tearing of two- teractions inherent in strongly interacting matter. Much dimensional sheets (see e.g., (Chen et al., 2016)). It also previous research has been devoted to understanding sim- includes other approaches to incorporate long-range in- ple colloids in equilibrium and at the glass transition teractions between local entities or to create global re- (Cates and Evans, 2000; Hunter and Weeks, 2012; Kle- sponse. Recent work has seen the incorporation of topo- man and Lavrentovich, 2003; Witten and Pincus, 2004). logical constraints to control excitations. Colloids have become such a good test system in part due Soft-matter research is highly multi-disciplinary. Be- to a variety of tools that have been invented to manip- cause of the enormous range of both materials and at- ulate and measure their properties with unprecedented tributes, it is not surprising that investigators come from precision. As described in the section on Instrumenta- many disciplines including physics, chemistry, materials tion below, these include optical tweezers, diffusion wave science, biology, geophysics, applied mathematics in ad- spectroscopy (DWS) and confocal microscopy. dition to chemical, biological, mechanical, and civil en- New physics can be revealed by adding complexity to gineering departments. The range of interests has many the colloids and to their inter-particle interactions and overlapping threads. It is essential to the health of the by pushing the systems to be far-from-equilibrium. The field to recognize that they all coexist under one roof. goal is to understand the really complex dynamic inter- Soft-matter science provides an excellent tool for out- actions of real systems. These include global collective reach and is accessible to a general audience because and emergent phenomena. Such behavior results from the phenomena that are investigated are so pertinent colloidal-scale interactions and can be sensitive to small to daily life: patterns, food, cosmetics, plastics, mod- changes in these interactions. For example, the field is ern display technology. The subject matter adapts itself poised to address the spontaneous assembly (or disas- readily to simple demonstrations, and as such can be ex- sembly) of structures, the emergence (or averaging out) hibited in museums and public lectures. In addition, the of large-scale fluctuations. Understanding some of these research often creates many stunning and memorable im- issues would be important for material processing. ages whose import is easily grasped. One forefront issue is to use colloidal systems to study the onset of non-equilibrium behavior. While, the small There is no ideal way to organize a report on experi- size of colloids often allows them to reach equilibrium, mental soft matter. While different materials are stud- the relaxation times are sufficiently long that the sys- ied, many issues transcend specific materials. These give tem cannot always keep up with the forcing. For exam- intellectual coherence to the field. The following sections ple, slow compression can push a colloidal system into will give an overview of some promising areas of research. a glassy or jammed phase where the system has fallen Some topics, such as biologically inspired matter, are not out of equilibrium with the thermal bath. These glassy given a section of their own but are referred to through- systems, with radial inter-particle interactions, provide a out the text. The following sections summarize the dis- minimal model for systems like metallic glasses and are cussions that occurred at a workshop held in January 30- an excellent model system with which to study the tran- 31, 2016 devoted to surveying current challenges facing sition from a fluid to a rigid solid and the transition from soft-matter experiment. The participants and discussion a thermal to an athermal system, topics discussed else- leaders are listed in the Appendix. where in this report. A fascinating aspect would be to use colloids as a tool to gain insight into non-equilibrium behavior. Can these systems reveal dynamics in unex- II. COLLOIDS pected experimental phenomena? Colloidal aggregates are often disordered. While some- A colloidal suspension is a mixture of microscopic insol- times a bane, disorder can also provide advantages: it uble particles, between approximately 1nm and 1µm in can make new material properties accessible that would diameter, that are dispersed in another medium such as otherwise be difficult to obtain in ordered matter. For a solvent fluid. Perhaps the most commonly encountered example, disorder can provide fault-tolerant tools. Also, 7 because of disorder, a given bond can have widely vary- between theory, computation and experiment. ing contributions to different global properties (Goodrich et al., 2015). Can novel materials be made from disor- dered colloidal aggregates that use this property? III. GRANULAR MATERIALS, (NON-COLLOIDAL) Another avenue to be pursued is to push the study of SUSPENSIONS, EMULSIONS AND FOAMS colloidal materials into new and extreme regimes. For ex- ample, these systems could be pushed to higher driving Granular physics is central to many aspects of soft- rates and into new physical environments. Their dynam- matter research. This subfield concerns materials com- ics would need to be studied at much faster time scales. posed of macroscopic “grains”, which are usually solid Colloidal suspensions can also be driven by external forc- particles but also may be liquid droplets as in emulsions ing or can be made “active” so that the individual par- or gas bubbles as in foams (Duran, 2012; Halsey and ticles have their own energy source. This suggests quali- Mehta, 2003; Jaeger et al., 1996). By contrast with the tatively new capabilities for how materials self organize. colloidal regime, the particles are large enough that ther- For example, it can provide an experimental platform to mal motion is usually negligible. Grain-grain interactions study flocking or swarming. Likewise nanoparticles can are often dominated by contact forces, but may also be display directional motion in response to external fields. mediated by an interstitial fluid or by electrostatic in- Biological matter has many mechanisms that teractions. There are strong scientific connections with exquisitely accomplish and control specific tasks. other areas highlighted in this report, including colloids, For example, a protein can interact with another mechanics, rheology, glasses and jamming, simulations molecule at one point on its surface which then allows and big data, pattern formation, and packing. Active the binding of a second molecule at a distant location. matter, human crowds, and traffic also have a granular Such “allosteric” interactions are an important way in nature but involve more complex “grains” with “social” which protein activity is regulated (Ribeiro and Ortiz, interactions. 2016). Colloidal systems may present an opportunity It is worth emphasizing at the outset the enormous to study allostery in a physical system and examine range of size and conditions where granular matter ap- what aspects of structure allow specific and localized pears. It can be found as sand on beaches, as screes long-range interactions. on mountain sides, as icebergs floating off the edge of Another bio-inspired example is to use E. coli as an glaciers, and, at the celestial scale, as planetary rings and inspiration for designing colloids (which, in analogy, we asteroids. It is also found ubiquitously in industrial pro- call “e-colloids”) that would report on chemical interac- cesses from the pharmaceutical to the civil-construction tions within a material and react to them in a specified industries. Understanding the science of granular mat- way. A major challenge is to design a more general chem- ter is important for many societal needs with implications istry that could engineer colloids with selected surface for soil remediation, mining and resource extraction, car- and bulk properties in order to react to, and report on, bon sequestration, erosion, sedimentation, siltation, and a variety of interactions. For example, a tissue spanning water purification. two distinct (immiscible) fluid phases, taking informa- Granular materials exhibit four prototypical proper- tion from one phase to cause a reaction in the other, ties: they are (i) disordered and heterogeneous, (ii) could form a natural platform for decision-making struc- highly dissipative, (iii) far-from-equilibrium with (iv) tures. This also offers the possibility of remote control or strongly nonlinear response to applied forces. Because robotic applications. The colloids can act as sensors and there is no understanding of what ensemble can be used amplify the signal at larger scales. Another goal would to obtain average properties, there is no adequate statis- be to design colloids that respond only when the collec- tical mechanics of a granular material. The usual tools tive arrangement of the particles is in a specific structure. and concepts we rely upon for crystalline materials, flu- Such reporter e-colloids would generalize one of the most ids, and even most other forms of soft matter, simply amazing mechanisms in the biological world and would do not apply. Friction, important for both statics and open up new ways of designing smart and programmable dynamics, is poorly understood. Except in the most ide- materials. alized situations, it is not understood how friction (and In closing this section, one should mention the fore- dissipation generally) changes the way a material reacts front issue of designing materials with specified rheo- to external perturbations. It has thus proven uncom- logical properties. For example, it would be useful to monly difficult to understand and manipulate the behav- have materials with a temperature-independent viscosity ior of a large collection of seemingly simple objects such or materials that behave in desired ways under extreme as sand grains. The emergent behavior can be surpris- shear. These are examples of an inverse problem; one ing, often with no parallel in other forms of matter as has starts with a desired property and then attempts to fig- been emphasized in various reviews. ure out how to make that property occur. This aspect of In many cases, we encounter granular material in a materials science necessarily requires a strong interaction static, but far-from-equilibrium, configuration. This is

Vol 435|23 June 2005|doi:10.1038/nature03805 LETTERS

Contact force measurements and stress-induced anisotropy in granular materials

T. S. Majmudar1 & R. P. Behringer1

Interparticle forces in granular media form an inhomogeneous strain) disks subjected to pure shear and isotropic compression. distribution of filamentary force chains. Understanding such Figure 1 shows a diagram of the experimental set-up and some forces and their spatial correlations, specifically in response to typical images; details of the set-up and the experimental procedure forces at the system boundaries1,2, represents a fundamental goal are described in Fig. 1 legend. Although previous approaches14–17 of granular mechanics. The problem is of relevance to civil have obtained contact forces using photoelastic techniques, they were engineering, geophysics and physics3–5, being important for the neither automated nor suitableforalargeenoughnumberof understanding of jamming, shear-induced yielding and mechan- ical response. Here we report measurements of the normal and tangential grain-scale forces inside a two-dimensional system of photoelastic disks that are subject to pure shear and isotropic compression. Various statistical measures show the underlying differences between these two stress states. These differences appear in the distributions of normal forces (which are more rounded for compression than shear), although not in the distri- butions of tangential forces (which are exponential in both cases). Sheared systems show anisotropy in the distributions of both the contact network and the contact forces. Anisotropy also occurs in the spatial correlations of forces, which provide a quantitative replacement for the idea of force chains. Sheared systems have 8 long-range correlations in the direction of force chains, whereas isotropically compressed systems have short-range correlations regardless of the direction. Under the action of external stresses, grains in dry granular materials form an inhomogeneous contact network, which carries most of the external load by way of force chains. The resultant network is different for shearing than for isotropic compression and 6–8 is history-dependent owing to friction. Previous experiments have FIG. 7 Liquid drop impact craters from a 3.1 mm water drop. Images show ring-shaped, solid pellet-shaped, asym- reported an exponential tail for the distribution of contact force metric residues as impact parameters are varied. From (Zhao magnitudes. This tail can be successfully predicted by many et al., 2015). 9–11 models with radically different mathematical structures and FIG. 6 Force chains in a two-dimensional granular material affect the material properties and can it be manipulated that is compressed uniaxially. The material is placed between to control a material’s behavior? microscopic assumptions. Testing the validity of these models Figure 1 | Experimental set-up and representativetwo crossed polarizers. data When. undera, Schematic stress, the granular par- Even a seemingly stationary granular material has ticles become birefringent and rotate the polarization of the requires that the predicted force distributions be verified by measure- diagram of the biaxial test cell. The biaxialtransmitted test light.apparatus The particles rests that are horizontally under stress are thus complicated dynamics if it is under an applied stress. For on a sheet of Plexiglas and is used to impartmade visible pure so that shear the stress and paths, isotropic from one side of the ma- example, it can show sub-threshold creep where local par- ments of full vectorial contact forces in the bulk of the sample. terial to the other, are made apparent. From (Majmudar and ticle rearrangements allow small reductions in the stress. compression. Motorized linear slides moveBehringer, two 2005). walls of the biaxial cell Thus, one needs to understand the nature of the linear It is also important to find other distinguishing signatures character- 21 regime for the different elastic moduli. The material can izing the nature of force chain networks under different boundary precisely and independently with a velocity of 0.024 cm s . The system is also give way in a large-scale avalanche – a prototypi- illuminated from below and a high-resolutionan example camera of a jammed captures state of matter. digital Because there cal catastrophic failure mode of a material along a slip conditions—an important goal of the present work. are external stresses (gravity or applied stress) and be- plane. The failure mechanisms depend on the formation images from above. Each image capturescause roughly thermal energy 250 is particles too negligible located to create a rear- history of the pile and how it was packed. The essence In the following experiments, we visualize internal stresses in each rangement event, the system is doomed to be stuck in of avalanche behavior can perhaps be captured in the around the centre of the cell, roughly 10%an amorphous of the total state that number is far from of its particles. lowest-energy collapse of simple experimental piles of sand or in simpli- 12,13 grain and by solving the full inverse photoelastic problem for ground state. The idea of jamming has been a produc- fied cellular-automata models. Clearly, however, richer The system is imaged through crossed circulartive way to polarizers. think about other For far-from-equilibrium each type sys- of sets of behavior can be found in real-world situations en- each disk, we obtain normal and tangential force components for load, incremental deformations are appliedtems (Liu in and a quasi-static Nagel, 2010). While manner, much effort has been countered, for example, in geophysical phenomena (and devoted to studying jamming in the simplest of cases, lab-scale counterparts). each contact between disks. We use this microscopic contact force beginning with a stress-free state. The shearedcurrent research states and arefuture created opportunities by involve focus- Just as Newtonian fluids routinely display a great ing on more complicated and realistic systems with more range of patterns, their complex cousins – granular ma- information to investigate differences in the distributions of contact compressing in one direction and expandinggeneral particle by an shape equal and interactions. amount Recent in theattention terials with their non-linear and dissipative interactions forces, and the force chain structure, arising from two different types other direction, with strains (1xx 1yy has focusedDL/L on) understanding ranging shear from jamming, 0 to the 0.042. flow of – likewise demonstrate a wealth of structures (such as ¼ ¼frictionalj particlesj under conditions of uniform shear (Bi oscillons (Umbanhowar et al., 1996), localized excita- of loads: pure shear and isotropic compression. We find that forces L (,40 cm) is the initial system length inet al. the, 2011).x or Newy issuesdirection. arise when Isotropically one considers spa- tions in thin granular layers) when subjected to vibra- compressed states are created by compressingtial heterogeneities in both and directions boundary effects. with The strains kinetics of tion. When a granular material flows, its flow profile is have distinctive angular distributions and spatial correlations jamming, involving the time evolution of a flowing sys- highly heterogeneous and the material has a non-local ranging from 0 to 0.016. The particles usedtem, is in a broad the new experiment challenge. are either rheology. Many remarkable phenomena take place such depending on the macroscopic preparation. In particular, forces In studying the mechanics of a dry (or indeed, even a as the segregation of grains according to their physical 0.8 cm or 0.9 cm in diameter and 0.6 cm inwet) height, granular medium, with a one Young’s can ask how modulus microstructure, of properties such as size, roughness or density. When a have long-range correlations in the direction of force chains for 4 MPa and a friction coefficient of 0.8. Thee.g., particle number shape or ratio roughness, of small affects the to dynamics large of projectile suddenly impacts a granular bed, as would oc- flow or the force transmission through a static pile. At cur in a meteor hitting the earth, the material gives way sheared systems, but are correlated over a much shorter range, disks is 4:1. b, Typical system size imagesa larger for scale, an isotropically a widely discussed, butcompressed unresolved, issue and flows in distinctive patterns as shown in Fig. 7. is the role of mesostructure (structure that can be iden- One phenomenon associated with flowing dry grains regardless of direction, for isotropically compressed systems. state (top) and a sheared state (bottom). ctified, An at example an intermediate of scale the between observed the particle stress size is the build-up of static electricity as individual parti- Our experimental system is a two-dimensional (2D) array of pattern for a single disk at the resolutionand ( the0.01 system cm dimension). per pixel) For example, used Fig. in 6these shows cles become charged and aggregate into larger clumps. the, existence of force chains in a two-dimensional granu- This can occur even if the particles interact with oth- approximately 2,500 bidisperse photoelastic (birefringent under studies. lar material under stress. How does such mesostructure ers made of the same material. Such charging is ubiqui-

1Department of Physics & Center for Nonlinear and Complex Systems, Duke University, Durham, North Carolina, 27708-0305, USA. 1079 © 2005 Nature Publishing Group MORPHOLOGY OF RAIN WATER CHANNELING IN … PHYS. REV. APPLIED 2, 044004 (2014) 2 IEEE TRANSACTIONS ON ROBOTICS 9

FIG. 8 Channeling of water in dry hydrophilic sandy soils at Fig. 2. Jamming gripper can achieve three separate gripping modes. (Left) different grain diameters. In the left panel, the grain diameter Static friction from surface contact. (Center) Geometric constraints from in- terlocking. (Right) Vacuum suction from an airtight seal. Normally, it would is 0.5 mm; in the right panel it is 1.0 mm. From (Wei et al., be unlikely that the interlocking or vacuum modes would be achieved without 2014). some additional contribution from friction.

tous, important for many industrial applications, puz- Optimal performance of a universal jamming gripper is main- 112 zling and poorly understood (P¨ahtz et al., 2010). In tained by resetting the gripper to a neutral state between gripping 113 addition, particle shape, heterogeneity, bonding (includ- tasks. Prior to the work presented here, this was accomplished 114 ing van der Waals forces, capillary bridges, or chemical by shaking the gripper, by kneading or massaging the gripper, 115 bonds) and many-body elasticity are important for un- or by pushing the gripper against some resetting apparatus that 116 derstanding the nature of granular flow. With all of these is mounted in the workspace, for example. We call this process 117 FIG. 2. Images (a)–(d) show deionized-water channeling in dry hydrophilic model sandy soils with bead diameters DFIG.varying 9 Universal from grippers based on jamming in granular ma- interactions, it is a challenge to design the properties at Fig. 1. Universal jamming gripper is able to grip a wide variety of objects manually resetting the gripper,andwithoutit,theabilityto 118 0.18 to 1 mm, at a rain rate of Q 14.5 cm=h and a raindrop impingingthe grain speed level of inU orderT to1 controlm=s. The a material’s color images behavior. in thewithoutterial. first grasp row A rubberplanning are the membrane,or sensory feedback. containing Multiple a objects granular can material,be gripped at ¼ ¼ once, as demonstrated here with salt and pepper shakers. grip subsequent objects degrades rapidly. We have found that 119 original images taken at steady states; the gray-scale images in the secondA vast row and are industrially obtained important by subtracting area has the to do background with can images be deformed taken so that its shape conforms to that of an ob- ject to be manipulated. The membrane can be evacuated so positive pressure can be used to replace this procedure with 120 granular suspensions in which the grains are dispersed prior to the rain from the ones taken at steady states, converting to gray scale, and then enhancing the contrast. The samplethat the packing granular in material becomes rigid due to the external ashortburstofairthatquicklyunjamsandresetsthegrip- 121 in a fluid. Such fluids have a complex rheology with each image is 26 cm wide, 25 cm high, and 0.8 cm thick. As labeled in the images, zwet is the infiltration depth of rain81 waterairThe pressure. in approach soils, Thisd thatis produces we propose a universal in this gripper paper that is to can use lift both a per. We also find that incorporating positive pressure improves 122 shear-thinning and shear-thickening regimes. Cornstarch large variety of objects with different shapes. Here the gripper 82 positive and negative pressure to modulate the jamming transi- the gripper’s speed, reliability, error tolerance, and placement 123 the water channel width, and d0 is the water channel separation. Thein saturated water (oobleck) hydraulic is a particularly conductivity dramaticκs [Eq. example(2)] for of theseis shown four picking samples up a salt and pepper shaker. From (Amend 83 tion in a universal jamming gripper. We design, manufacture, accuracy. In addition, the fast ejection that positive pressure 124 are determined as 73 cm=h, 204 cm=h, 567 cm=h, and 2300 cm=theh, respectively.latter regime; a person can run across a vat of the et al., 2012). 84 and test a prototype gripper that attaches to a commercial robot can provide enables the gripper to launch objects a significant 125 material but will sink if she pauses at any point. The 85 arm. Consisting of a single mass of granular material encased in distance—a capability that we call shooting,whichmayserveas 126 study of this extreme shear thickening has required the 86 an elastic membrane, the gripper can passively conform to the anewmethodforrobotstoextendtheirworkspaceandperform 127 development of new experimental tools. Other aspects Glass formation and jamming occurs when a material the separation between channels is very similar. The use of deionized-water and glycerol-water mixtures87 shape are of the used target as object, then vacuum-harden to grip it rigidly, tasks like sorting objects into bins in a factory or throwing away 128 of granular suspensions that are currently being investi- that can flow loses that ability as either the tempera- different-width sample cells excludes the effects of sample the rain water supply so that the effects of88 rainlater water using positive quality pressure to reverse this transition—releasing trash in a home. 129 gated have to do with sedimentation, clogging and ero- ture, pressure, density or some other control parameter is 89 the object and returning to a deformable state. An example of cell size. As an example, in a 26-cm-wide packing of 1-mm sion.can An also example be investigated. in Fig. 8 shows An the important morphology parameter of varied called (Berthier the and Biroli, 2011; Debenedetti and Still- In this paper, we develop a new universal jamming gripper 130 90 this gripper can be seen in Fig. 1. water channeling in sandy soils (Wei et al., 2014). inger, 2001; Ediger et al., 1996). It differs from crystal- that incorporates positive pressure. We quantify the gripper’s 131 hydrophilic glass beads, we see either one water channel saturated hydraulic conductivity is widely91 usedThis universal in these jamming gripper is an example of a passive uni- Active materials are emphasized throughout this re- lization, the more easily understood but not necessarily ability to grip objects of different shapes and sizes, as well as its 132 92 versal gripper that exploits the temperature-independent fluid- form near the center or twowater channels form close to each port.studies. Granular It matter quantitatively is an arena where shows the how challenge fast of a fluidmore cancommonly move found form of rigidity formation, which ability to tolerate errors in the location of the target object; we 133 93 like to solid-like phase transition of granular materials known as side; in a 55-cm-wide packing of 1-mm hydrophilic glass creatingthrough such materialsthe pore with spaces smart, in shape-changing a soil. For the and modeloccurs sandy when soils a liquid, cooled below its freezing tempera- test the gripper’s maximum speed and placement precision; we 134 94 jamming [26]–[31]. This gripper leverages three possible grip- reporting particles can be found. Soft robotics based on ture, nucleates crystallites that then grow to encompass test the gripper’s ability to grip multiple objects at once and to 135 beads, we see four to six water channels form with similar we use, the value of the saturated hydraulic95 ping conductivity modes for operation:κs 1) static friction from surface contact; granular materials has demonstrated exceptional success the entire system. The issues associated with rigidity in shoot objects of varying weight and shape. Our testing reveals 136 separation. can be determined as 96 2) geometric constraints from capture of the objectIEEE by interlock- as shown in Fig. 9. It provides a promising new way of amorphous solids are so common and relevant that the the capabilities and limitations of the gripper, and we compare 137 Figures 2(a)–2(d) show typical steady states of rain water using jamming to control motion and rigidity. 97 ing;study and of 3) all vacuum glass-forming suction when atomic an airtight and molecular seal is achieved systems on these with a manual reset gripper in order to isolate the perfor- 138 98 some portion of the object’s surface [25]. These three gripping There are an evolving set of tools thatρg are being2 used is aptly covered by the umbrella of soft-matter science. mance contribution from positive pressure. We demonstrate that 139 channeling in dry hydrophilic model sandy soils with κ K D ; 99 modes are illustrated2 in Fig. 2. The friction force results from to study granular material in alls its¼ forms.μ 0 These in- While windowð Þ glass is the most common example of dramatic improvements in performance are possible through the 140 varying bead diameters. The applied rainfall rate is clude magnetic resonance imaging, X-ray tomography,100 thea glass,slight ( nearly<0.5%) all volume liquids contraction can be made of the glassy membrane if only that addition of positive pressure, and we compare the performance 141 Q 14.5 cm=h, far smaller than the saturated hydraulic high-speed video and ultrasound that reveal structure,101 occursthey duringare cooled evacuation, rapidly which, enough. in turn, (Counter-examples causes a pinch force are to of a positive pressure jamming gripper with related grippers in 142 conductivity¼ of soils. The top row consists of raw images confocalwhere microscopyρ and thatμ are can measure the density contact forces, and thebig-102 viscositydevelop,ordinary normal water of theto and, the point of course, of contact. liquid Analytical helium calculations at atmo- the field. We conclude that this gripper has potential applications 143 data,applied machine-learning fluid, respectively, and evolutionaryD algorithmsis the soil that103 grainforspheric diameter, these values pressure.) and have been There previously are metallic, presented semi-conducting, [25]. in a variety of settings. 144 taken after a steady state is achieved. The bottom row shows 104 polymericBy achieving and one insulating or more glasses. of the three On larger gripping scales, modes, glass the Proof can design, analyze and identify novel phenomena in ex- periments.K0 is theDusty intrinsic plasmas now permeability offer a promising of the way105 soil.to jammingformation For a gripper random has been can grip studied many in different colloidal objects fluids, with and widely jam- gray-scale images obtained by subtracting the background II. DESIGN AND MANUFACTURE 145 images taken before the rain started from the raw ones in studyclose interactions packing in granularof spheres matter with in new a porosity ways. 106 of ϵvarying≈ming0.36 of shape,, courseK0 weight,has occurs and in fragility, suspensions including and objects macroscopic that are 4 107 traditionallysamples such challenging as foams for and other sandpiles. universal It grippers. has been For much exam- In its simplest form, a jamming gripper needs only to include 146 a value of 6.3 × 10− [3,20]. In theories of channeling, the the top row, converting the result to gray scale, and then 108 ple,debated we have how successfully the physics been of able the to glass grip transition a coin, a tetrahedron, overlaps some granular material that is contained in a flexible mem- 147 rain rate often enters as the dimensionless ratio Q=κ , and enhancing the contrast using a histogram equalization IV. GLASSES AND JAMMING 109 ahemisphere,arawegg,ajacktoy,andafoamearplug.Whenwith the physicss of jamming. Is the glass transition uni- brane in order to achieve its gripping behavior (the combination 148 experiments are performed at variable Q110. Here,mounted we to vary the robotQ, arm, the gripper functions entirely in open of ground coffee and a latex balloon has been found to work 149 algorithm provided by MATLAB. The three parameters that 111 loop—without grasp planning, vision, or sensory feedback. well [25]). No motors, cables, or linkages are required (just an 150 we are interested in are labeled on the processed gray-scale but we also shall vary the fluid viscosity as an alternative means of achieving a change in Q=κs. images of the figure. They are the infiltration depth zwet of the wetting front, the channel width d, and the channel When impinging on sandy soils, rain droplets create craters on the soil surface, whose size and shape depend separation d0. In the following subsections, we discuss the effects of rain conditions and sandy soil properties on these strongly on both the impinging speed and the soil grain parameters, respectively. diameter [21,22]. To determine if the craters affect the stability of the horizontal wetting front and, thus, control the locations of the water channels, we vary the raindrop A. Effects of raindrop impinging speed and rain rate impinging speed UT by adjusting the distance h between We first study the influence of the rain conditions, the capillary tips on the rain source and the soil surface at such as the raindrop impinging speed UT and rain rate two fixed rain rates, Q 14.5 cm=h and Q 96.0 cm=h. Q, using 1-mm dry hydrophilic glass bead packing. Both When a low rain rate¼ of Q 14.5 cm=¼h is applied, ¼

044004-3 10 versal (the same for a wide array of materials) or is it off by this inability to go to long times (although recent system-dependent so that the details of bonding (i.e., simulations have obtained systems with time scales that covalent versus ionic versus hydrogen bonding) matter? are comparable to those in experiment (Berthier et al., More provocatively, one can ask whether the glass transi- 2016)). tion in an associated liquid is the same as the freezing of When placed under stress, a glass will creep, yield and the dynamics in a colloidal sample or the yielding tran- flow. It is often attractive to think of these processes sition in a granular material. Such questions were the in terms of the motion of defects. While the concept of genesis of the jamming phase diagram (Liu and Nagel, a defect is clear in a crystalline material, an important 2010) where it was postulated that understanding e.g., current issue is to understand what constitutes a defect glasses would be relevant to other types of amorphous in an amorphous solid. It is far from obvious what those matter. Dynamic heterogeneities have been studied in defects might be in a glass, granular pile or other amor- many materials as the system dynamics becomes arrested phous material. This has led to a body of research at- (Berthier et al., 2011). tempting to identify from the local structure the location As an ordinary liquid is supercooled, the dynamics slow of soft spots where a particulate material is most likely to down until the time for relaxation becomes much longer undergo a rearrangement (Schoenholz et al., 2016). Much than any observer is prepared to wait. That is, the ma- remains to be done to understand soft-spot dynamics and terial falls out of equilibrium when the thermal energy how such “defects” interact with one another. in the sample is insufficient for the system to explore ad- Colloidal suspensions can be a way to study the physics equqtely the phase space in the accessible experimental of rigidity formation as a liquid state is cooled (or com- time. Glass formation is the prototypical example of how pressed) as would occur when a glass is formed out of a system can fall out of equilibrium, not because it is be- a supercooled liquid. As shown in Fig. 2, colloidal sys- ing constantly driven, but because its thermal energy is tems have allowed the visualization of local rearranging too small. Thus when studying a glass, one should be clusters that appear as a fluid loses its ability to flow aware that this state was reached, not by remaining in freely (Weeks et al., 2000). In recent years, colloidal sys- (quasi-) equilibrium while cooling through a transition, tems have been the most highly exploited experimental but by the system being forced to fall from equilibrium. system for testing the ideas that have emerged from com- In that sense, experimental low-temperature glasses are putational studies of jamming (Liu and Nagel, 2010). By similar to the jammed state of matter. studying the position and dynamics of individual parti- Unlike crystals, glasses exist in an incredibly large cles, it has been possible to observe the non-monotonic number of nearly equivalent states in an exceptionally correlation function as a function of compression, to mea- ramified, complex energy landscape. Dealing with so sure the spatial properties of localized normal modes, or many relevant energy minima is one of the central prob- to observe soft spots, the sites most likely to fail as the lems of statistical physics. As such, understanding glassy material is stressed (Caswell et al., 2013; Chen et al., behavior can open a window on many other physical 2011; Zhang et al., 2009). problems, ranging from cosmology and geology to chem- New challenges await in this area. These have to do istry and materials science, that deal with rough energy with studying how a particle’s (non-spherical) shape af- landscapes. It is useful to divide the studies of glasses fects jamming and glass formation and how more com- into two categories (i) those that study the glass transi- plicated inter-particle interactions change the transition tion and (ii) those that study the properties of the glass and the nature of the dynamics in the rigid phase. For itself and investigate how it is different from a crystal example, studies of colloids that have long-range or three- (Anderson et al., 2012). body interactions, have non-spherical shapes, or have One question that has plagued the glass-transition lit- frictional contacts can produce considerably different erature for decades is whether there is a true thermody- physics from the simple frictionless spheres with finite- namic phase transition that separates the liquid and the range repulsive interactions most often studied on com- glassy states as the temperature is lowered. From exper- puters. It would be of particular importance to under- iments alone, it has been impossible to answer this seem- stand the role of a small temperature on the jamming ingly simple question. One finds that as the temperature scenario. While some rearrangements would be thermally is lowered, relaxation times and viscosity rapidly increase allowed, they are far from sufficient to bring the entire by over 17 decades in a narrow temperature range (Ediger system to equilibrium. It is also important to consider et al., 1996); but it is impossible to tell if they diverge be- what other axes are relevant to the jamming phase dia- cause no matter how patient the experimenter, her exper- gram or whether the three canonical ones (temperature, iment invariably falls out of equilibrium at temperatures inverse density and shear stress) are sufficient. that are tens of degrees away from where the divergence As in so many of the other areas covered in this re- would appear to exist. All one can do is extrapolate port, the possibility of creating glasses and jammed mat- into this low-temperature, long-relaxation time region. ter out of active particles is the start of an entirely new Simulations on large particulate systems likewise are cut area of research. Active swimmers have correlated mo- REPORTS

thethe same same as as that that of of long long linear polymers of the way, vf(M)takesonaformidenticaltotheFlory perature, the similarities between between jammed jammed ball- ball- samesame type, type, but as M becomes very small, Tg in- form (Eq. 2) for Tg(M). Fitted to appropriate chains and glassy glassy polymers polymers are are pronounced pronounced enough enough creasescreases (19).). As As with with ball-chains, the behavior of regimes of the rf(M) data, we find for x = 7.5, toto suggest suggest that that the the jamming jamming idea idea captures captures much much of of linearlinear and and cyclic cyclic polymers polymers follows opposite trends vb = 2.3 and l(vb – ve) = 2.6; taking l = x = 7.5, thethe physics. physics. Because Because it it has has long long been been a a matter matter of of asas M becomes small. then ve =1.95.Forx =11.2,thebestfitgivesvb = debate whether glass transitions transitions in in polymers polymers and and in in The jamming phase diagram suggests a way to 2.5, and l(vb – ve) = 6; taking l = x = 11.2 gives simple glass-formers glass-formers are are fundamentally fundamentally similar, similar, it it is is directly compare glassy polymers with jammed ve =1.96. attractive to to think think that that jamming, jamming, whether whether of of grains grains chains.chains. Here, temperature T and specific packing Finally, chains with larger loop sizes pack less or chains, may provide a a unifying unifying connection. connection. volume v ≡ 1/r are orthogonal axes; jammed/glassy densely than chains with small x.Theanalogyof We have shown that long, floppy floppy chains chains pack pack statesstates occupy occupy the high-density, low-temperature vf(M)withTg(M)suggeststhatstiffpolymerswill intointo a a low-density low-density structure structure whose whose rigidity rigidity is is regionregion of of the the phase phase diagram, and fluid states occupy have a higher Tg than those that are more flexible. flexible. chiefly provided by by the the backbone. backbone. This This can can be be thethe rest. rest. The The glass glass transition transition is the transition from For vinyl polymers with rigid side groups that understood as the jamming jamming of of semi-rigid semi-rigid loops loops liquidliquid to to glass glass along along the the Taxis, whereas jamming is restrict chain flexibility, Tg is higher for for those those thatthat formed formed when when the the chains chains were were compacted. compacted. By By thethe transition transition from from unjammed unjammed to jammed states with large side groups than for those whose side invokinginvoking an an analogy analogy between between the the specific specific pack- pack- alongalong the v axis at T =0(17). In this sense, Tg(M) groups are small (20). inging volume volume and and the the glass glass transition transition temperature, temperature, andandvff((M))≡ 1/rff((M)shouldbeanalogousquantities. These observations do not prove that the poly- we have shown that the packing of of chains chains When we compare Tg(M)/Tg,∞ for linear and cyclic mer glass transition must be attributable to the parallels the polymer glass glass transition transition in in important important poly(dimethylsiloxane) (PDMS) and vf(M)/vf,∞ for jamming of chains, especially at finite temper- temper- respects.respects. If If these these two two phenomena phenomena are are indeed indeed thethe packing packing of of chains, chains, we we see that there is a marked atures and with realistic interactions. The role of closely connected, as as our our data data suggest, suggest, it it would would similaritysimilarity between the two (Fig. 4). temperature points to important distinctions be- be- be a beautiful illustration of of how how molecular molecular ge- ge- As shown earlier, in a packing of floppy linear tween polymers and ball-chains: (i) The stiffness stiffness ometry and symmetry, independent independent of of the the spe- spe- chains,chains, chain ends act as non-interacting, density- of a ball-chain is given solely by its construction, cific microscopic interactions, interactions, can can influence influence the the enhancingenhancing defects. We can then write the specific but a polymer becomes stiffer at lower T,and(ii),and(ii) structure of of condensed condensed matter. matter. packing volume vff inin the the floppy regime as the for the packing of macroscopic objects, thermal weighted sum of ve,, contributed contributed by ends, and motion is irrelevant, and the generalization of References and Notes Notes 1. T. T. Aste, Aste, D. D. Weaire, Weaire, TheThe Pursuit Pursuit of of Perfect Perfect Packing Packing (Taylor(Taylor on July 2, 2013 vvbb > vee,, due due to to the the bulk bulk temperature is an open question (21). Therefore, on July 2, 2013 & Francis, New York, York, ed. ed. 2, 2, 2008). 2008). the purely geometric jamming of chains cannot be be 2. R. R. Zallen, Zallen, TheThe Physics Physics of of Amorphous Amorphous Solids Solids 2l 2l exactly analogous to the glass transition in poly- (Wiley-Interscience,(Wiley-Interscience, New New York, York, 1998). 1998). vff M 1 − vb ve ðð Þ¼M þ M mers. Some aspects of temperature can perhaps be 3. J. J. D. D. Bernal, Bernal, NatureNature 183183,, 141 141 (1959). (1959). 4. G. G. D. D. Scott, Scott, NatureNature 194194,, 956 956 (1962). (1962). 2ll vb − ve simulated by varying the tapping strength used to vb − ð Þ 3 compact the packing (6, 22). Interaction can also 5. J. J. L. L. Finney, Finney, Proc.Proc. R. R. Soc. Soc. London London Ser. Ser. A A 319319,, 479 479 (1970). (1970). ¼ M ð Þ 6. J. J. B. B. Knight, Knight, C. C. G. G. Fandrich, Fandrich, C.-H. C.-H. Lau, Lau, H. H. M. M. Jaeger, Jaeger, be introduced, for instance, by making the chains S.S. R. R. Nagel, Nagel, Phys.Phys. Rev. Rev. E E Stat. Stat. Phys. Phys. Plasmas Plasmas Fluids Fluids Relat. Relat. Here ll ~ x isis the the extent extent that end influences slightly sticky with a thin coating of viscous oil.11 Interdiscip.Interdiscip. Top. Top. 5151,, 3957 3957 (1995). (1995). propagate into the bulk of the chain. Cast in this But even for hard particles and no explicit tem- 7. F. F. X. X. Villarruel, Villarruel, B. B. E. E. Lauderdale, Lauderdale, D. D. M. M. Mueth, Mueth, H. M. Jaeger, Phys.Phys. Rev. Rev. E E Stat. Stat. Phys. Phys. Plasmas Plasmas Fluids Fluids www.sciencemag.org Relat. Interdiscip. Top. 61, 6914 (2000). www.sciencemag.org tion that can lead to crowding and slow dynamics.Fig. 3. Tomographic Are Relat. Interdiscip. Top. 61, 6914 (2000). reconstructionreconstruction of of a a pack- 8. S. S. R. R. Williams, Williams, A. A. P. P. Philipse, Philipse, Phys.Phys. Rev. Rev. E E Stat. Stat. Nonlin. Nonlin. the properties of such systems the same as the conven- SoftSoft Matter Matter Phys. Phys. 6767,, 051301 051301 (2003). (2003). inging that that consists consists of of a a sin- sin- 9. A. A. Donev Donev etet al al.,., ScienceScience 303303,, 990 990 (2004). (2004). tional amorphous particle counterparts orglelongchain( are thereM new=4096). 10. W. W. Man Man etet al al.,., Phys.Phys. Rev. Rev. Lett. Lett. 9494,, 198001 198001 (2005). (2005). features that emerge? Are the flows in such((A active)Thefullreconstructed)Thefullreconstructed mate- 11. Z. Z. Zeravcic, Zeravcic, N. N. Xu, Xu, A. A. J. J. Liu, Liu, S. S. R. R. Nagel, Nagel, W. W. van van Saarloos, Saarloos, three-dimensionalthree-dimensional image. image. Europhys.Europhys. Lett. Lett. 8787,, 26001 26001 (2009). (2009). rials related to the excitations of their static((B counterparts)Asmallsliceofthe)Asmallsliceofthe 12. M. M. Mailman, Mailman, C. C. F. F. Schreck, Schreck, C. C. S. S. O O’’Hern,Hern, B. B. Chakraborty, Chakraborty, Phys. Rev. Lett. 102, 255501 (2009). and can one gain intuition into biological crowdingpacking, indicated by its by Phys. Rev. Lett. 102, 255501 (2009).

13. R. W. Cahn, Contemp. Phys. 21, 43 (1980). Downloaded from 13. R. W. Cahn, Contemp. Phys. 21, 43 (1980). Downloaded from analogy with jammed systems (Henkes et al.thethe, shaded shaded 2011)? box box in in (A), (A), 14. R. R. L. L. Jernigan, Jernigan, P. P. J. J. Flory, Flory, J.J. Chem. Chem. Phys. Phys. 5050,, 4178 4178 (1969). (1969). projected onto the hori- 15. I. I. Carmesin, Carmesin, K. K. Kremer, Kremer, MacromoleculesMacromolecules 2121,, 2819 2819 (1988). (1988). zontalzontal plane. plane. The The brighter 16. L. L. E. E. Silbert, Silbert, A. A. J. J. Liu, Liu, S. S. R. R. Nagel, Nagel, Phys.Phys. Rev. Rev. E E 7373,, aparticleis,thecloseritaparticleis,thecloserit 041304041304 (2006). (2006). is to the center plane of 17. A. A. J. J. Liu, Liu, S. S. R. R. Nagel, Nagel, NatureNature 396396,, 21 21 (1998). (1998). V. PACKING is to the center plane of 18. T. G. Fox, P. J. Flory, J. Appl. Phys. 21, 581 (1950). the slice. Highlighted 18. T. G. Fox, P. J. Flory, J. Appl. Phys. 21, 581 (1950). the slice. Highlighted 1 cm 19. S. S. J. J. Clarson, Clarson, K. K. Dodgson, Dodgson, J. J. A. A. Semlyen, Semlyen,PolymerPolymer2626,, 930 930 particles are arranged (1985).(1985). A material’s properties, such as whetherinin it three three isjammed near-minimal near-minimal loops. loops. 20. A. A. Eisenberg, Eisenberg, in in PhysicalPhysical Properties Properties of of Polymers Polymers,, or able to flow, are obviously dependent on how its con- FIG. 10 (A) X-ray tomography image of a “granular poly- J.J. E. E. Mark, Mark, Ed. Ed. (American (American Chemical Chemical Society, Society, Washington, Washington, Fig. 4. The glass transitionmer” temperature consisting of many1.2 equal-length segments of a connected DC, ed. 2, 2, 1993). 1993). stituent particles are assembled. When objects pack, one 21. S. S. F. F. Edwards, Edwards, D. D. V. V. Grinev, Grinev, Adv.Adv. Phys. Phys. 5151,, 1669 1669 (2002). (2002). TTgg((M)oflinearandcyclicPDMS(openand)oflinearandcyclicPDMS(openandchain of spheres. (B) Close-up image of red (shaded) region would like to be able to predict how theysolid will squares) assemble compared with the specific 22. E. E. R. R. Nowak, Nowak, J. J. B. B. Knight, Knight, E. E. Ben-Naim, Ben-Naim, H. H. M. M. Jaeger, Jaeger, solid squares) comparedin with (A). the specific The connecting links are visible from which it is possi- S.S. R. R. Nagel, Nagel, Phys.Phys. Rev. Rev. E E Stat. Stat. Phys. Phys. Plasmas Plasmas Fluids Fluids Relat. Relat. f, and what the packing density and overallpacking properties volume ofvf(M) ≡ 1/rf(M)oflinear f, Interdiscip. Top. 57, 1971 (1998). packing volume vf(M) ≡ 1/ f(M)oflinear v Interdiscip. Top. 57, 1971 (1998). v

ble to determine which spheres are attached together in each / / and cyclic x =11.2ball-chains(openand 1.0 23. We thank P. Eng for his assistance and T. Witten for

1.0

and cyclic x =11.2ball-chains(openand f 23. We thank P. Eng for his assistance and T. Witten for the material will be. One property of paramount impor- f

chain. From (Zouv et al., 2009). fruitful discussions. This work is supported by grants solid circles). Both are normalized by their v solid circles). Both are normalized by their glass fruitful discussions. This work is supported by grants tance for many applications is the thresholdasymptotic to materialM → values T and v . glass fromfrom NSF NSF Materials Materials Research Research Science Science and and Engineering Engineering asymptotic M → ∞ values Tg,∞ and vf,∞. transition CenterCenter (DMR-0820054) (DMR-0820054) and and the the U.S. U.S. Department Department of of failure from different applied stresses. To makePDMS dataheadway, are taken from (19). (Inset) The g, ,

PDMS data are taken from (19). (Inset) The g, , 0.8 EnergyEnergy (DOE) (DOE) (DE-FG02-03ER46088). (DE-FG02-03ER46088). GSECARS GSECARS is is one needs to know what inputs about the((TT particles,,vv)planeofthejammingphasediagram,)planeofthejammingphasediagram, and 0.8 supported by grants from NSF (EAR-0622171) and / T jamming supported by grants from NSF (EAR-0622171) and / T jamming with trajectories for glasscertainly transition and plays ag central role; frictionless particles simu- DOE (DE-FG02-94ER14466). Use of the Advanced with trajectories for glass transition and g DOE (DE-FG02-94ER14466). Use of the Advanced about the preparation conditions are necessary to pro- T jammed T jamming. temperature jammed Photon Source is supported by DOE, Office of jamming. lated on a computer do not displaytemperature nearly the range of Photon Source is supported by DOE, Office of duce accurate predictions. specific volume Science, Office of Basic Energy Sciences, under density variation that is found in laboratoryspecific volume packings. It Science, Office of Basic Energy Sciences, under 0.6 contractcontract no. no. DE-AC02-06CH11357. DE-AC02-06CH11357. The importance of this problem becomes clear when 0 1 2 3 4 5 is essential that all10 the interactions10 10 between10 objects10 (10e.g., 1 June 2009; accepted 1 September 2009 one considers the canonical example of two containers M (no. of monomers) 1 June 2009; accepted 1 September 2009 friction, cohesion, charging)M be(no. understood.of monomers) Monodis- 10.1126/science.1177114 filled with granular matter with different preparation or perse spheres have been fairly well studied. However, processing histories. As originally poured, a packing of 410 there is great16 OCTOBER latitude 2009 to create VOL 326moreSCIENCE exotic ingredients.www.sciencemag.org spheres will have one density but after just a few taps or It is now possible to create designer particles of differ- shakes (or by more gentle sedimentation) the grains set- ent shapes by using various chemical syntheses or, at a tle to a density that is over 10% greater than its original macroscopic scales, by using 3D printing. One can have value in the rapidly poured system (Onoda and Liniger, particles that are polydisperse, non-spherical, polymeric 1990). Particles with different shapes will have very dif- or with concave surfaces. Likewise, boundary conditions, ferent densities (Glotzer and Solomon, 2007). finite-size effects, external fields (e.g., gravity) are impor- There are also more subtle effects; packings of spherical tant. An example is the polymer-chain packing visual- grains that were prepared by different tapping protocols ized with X-ray tomography shown in Fig. 10. One goal to have the same density can be told apart by the way in is to be able to predict reliably the structure, density and which they respond to further perturbations (Josserand failure properties of such a packing. et al., 2000). This is a form of material memory that For many purposes, one wants packings with tailored is hidden in the packing structure. In another example, properties such as extreme loose or close packing. Thus the forces measured at different points along the bottom one wants to enhance or frustrate high packing densi- surface of a sandpile (due to the weight of the sand above ties. The role of packing becomes even more important it) varies with preparation conditions (Bouchaud et al., when one uses shapes that are more complicated than 2001; Geng et al., 2001; Wittmer et al., 1996). Can small simple spheres. One goal has been to find the parti- perturbations (such as from thermal expansion) disrupt cle shape that gives the highest-density packing (Donev the subtle structure of such packings? et al., 2004; Roth and Jaeger, 2016; Torquato and Still- It has proven notoriously difficult to characterize ex- inger, 2010). Another goal is to determine the limit of perimental packings. It requires tools that can acquire material parameters that can be obtained from different dynamic images of static and actively driven packings protocols. In the inverse problem, one wants to predict over a large range of scales deep into the bulk. Magnetic what shape particles should have to create a packing with resonance imaging (Nakagawa et al., 1993) or X-ray to- a particular property. One successful example used an mography has severe limitations in terms of resolution, evolutionary algorithm to predict a particle shape that speed, and size. Ultrasonic imaging (Han et al., 2016) matched a particular desired rheological feature (Miskin has issues with contrast but may be useful in some exper- and Jaeger, 2013). imental situations. When one goes to nanoscale particles, there are challenges for obtaining real-time imaging. To create reproducible packings, one must control the constituent particles’ shape and size. Moreover, friction VI. LIQUID CRYSTALS 12

Liquid-crystals exhibit some properties that resemble those of a liquid while others that resemble those of a solid. This behavior stems from the fact that the large constituent molecules in a liquid crystal are anisotropic and are arranged with an intermediate degree of order in some directions while they appear amorphous in oth- ers. Thus liquid crystals may flow while still retaining some molecular order. There are many different types of liquid crystals and, because the order in these materials is highly anisotropic, many liquid crystals have unusual and highly useful optical properties that can be manip- ulated to great effect in display technology. This has led to a > $100 billion/year industry (Yokoyama, 2012) and many new phases of liquid crystals are discovered each year. With this abundance of opportunity and new ma- 1 terial, it is not surprising that many fundamental issues FIG. 11 Two type- 2 strings in a liquid crystal crossing each need to be addressed. other and reconnecting. These defects have been taken as analogs of cosmic strings in the early universe. From (Chuang The symmetries and director order parameters that de- et al., 1991). scribe liquid crystals lead to a variety of topological de- fects that mediate the material’s structure. For example, when the director is pinned at two points with incompat- One can also extend the range of building blocks for ible alignments, a topological defect may be required to liquid-crystal materials. For example, one could enlarge mediate the transition between the two regions of space. the particle-size range by using anisotropic colloids. This Defect lines can even be tied into knots (Tkalec et al., leads to a complex space of possible outcomes because 2011). Experiment is essential for uncovering many of aspects, such as shape, flexibility and chirality, of the these subtle and intriguing structures and one can visu- individual particles matter. These shapes can lead to alize such topologically complex fields by making use of new possibilities for self-organization. the optical birefringence of these anisotropic materials. Liquid crystals also provide a unique platform for an- The development of large-scale atomistic molecular swering fundamental questions about how nature can cre- simulations could facilitate molecular design and the ma- ate novel forms of order, not only in condensed matter terial synthesis necessary to create materials with novel but also in more general contexts. Indeed, some of the characteristics. For example, shape-changing liquid crys- defects that have been discovered in liquid crystals, such tals would provide new mechanisms for collective reor- as those in Fig. 11, have been cited as similar to those ganization; by altering the structure one could induce a that appear in the fabric of the early universe (Chuang change of phase in the material. Another option for trig- et al., 1991). Thus liquid crystals not only have had im- gering large-scale reorganization from small changes in portant technological applications but also have provided input, is to use lyotropic molecules (molecules that cre- an exceptional resource for posing and investigating fun- ate liquid-crystal phases when in contact with a solvent). damental questions about collective behavior. In liquid crystals, surfaces often determine the orienta- A recurrent theme throughout this report, and of par- tion of the director field in its vicinity. Not surprisingly, ticular importance here, is the structure and dynamics therefore, interfaces play a central role in determining when a material is far from equilibrium. One can in- the properties of a liquid crystal. An important area for vestigate the classification of non-equilibrium states and future research is clearly to control surface interactions study the kinetics of defect motion. There are orders of so that the material will respond in a desired manner to magnitude in the amplification of response imparted by an external perturbation. New probes of surfaces and the long-range liquid-crystal organization. Moreover, the greater control over the chemistry at different scales is response can often be rapid. Although this amplification therefore desired for novel capabilities and applications. is at the core of many applications, its origin is not well In order to advance these goals, it is necessary to de- understood. velop new tools and instrumentation. One important ef- Another aspect occurs if one combines the polarity and fort would be to create advanced imaging capabilities so flow of liquid crystals with activity, so that individual that information can be obtained over a wide range of particles have their own energy source. In such cases, new length scales – from the microscopic level of the individ- regimes of behavior and organization can be expected as ual molecules to the macroscopic level of the order pa- seen in Fig. 5. The organization of such active liquid rameter fields. Another challenge would be the synthesis crystals cannot yet be predicted. Experiments would be of active liquid crystals. Here one would need to develop particularly revealing for this novel form of matter. a general, robust and efficient way to drive materials out 13 of equilibrium. cling to reinforced tires. Understanding supramolecular properties raises problems in statistical mechanics. These include how to design charge distributions on brushes and VII. POLYMERS AND RHEOLOGY create stimuli-responsive brushes for catalysis. Many material properties are determined by disorder Polymers are high-molecular-weight molecules that are including defects in networks, polydispersity in molec- concatenated sequences of monomers units that may ular weight and heterogeneity in charge distributions. or may not be identical. They are the basis of syn- Fluctuations influence material fabrication and proper- thetic plastics and natural biopolymers. These long- ties. They appear as phase transitions are approached, chain molecules, with a very high aspect ratio of length in the wettability of a surface, in the persistence length to diameter, are difficult to crystallize both due to steric of chains, in structural color, and potentially in the per- impurities in the repeat units and to their entanglements; meability of membranes. Understanding disorder is thus materials made from them are typically disordered and a continuing challenge. have unique properties due to their length and often- A problem common throughout much of soft matter flexible structure. Of course, there is a large variation in is to develop a better understanding of non-linear dy- the choice of the component monomers with a rich and namics. Understanding non-linear interactions from a flexible set of interactions and novel functional behav- macroscopic perspective is particularly important in en- ior. Block copolymers, which are polymers constructed tangled polymers where new models and measurements of long blocks of two distinct monomer units, can intro- are required. This has practical importance for process- duce structure on a tens of nanometer size scale. Poly- ing polymeric materials. When materials do not behave electrolytes in solvents such as water can lead to charged as expected for simple flow fields, such as shear or ex- polymers. For some reviews on polymers see chapters in tensional flow, it is particularly important to understand (Cates and Evans, 2000; Kleman and Lavrentovich, 2003; their behavior in mixed fields. This is especially relevant Witten and Pincus, 2004). to manufacturing. In the case of additive manufacturing, It has been often quoted that “There is a great future cooperative effects related to rheology of filled polymers, in plastics.” Given that the entire plastics industry re- where flow is required at high volume fractions, become lies on understanding the properties of polymers, much important. is known about certain classes of these molecules. How- Non-linear dynamics also plays a distinct role in the ever, given the wide variety of polymeric matter and in- complex rheology of active materials. Understanding teractions between polymers and with solvents, there are them fully could lead to opportunities in reconfigurable many outstanding issues. materials. In addition, a material can be trained so that There are non-trivial implications of topology and ar- it retains a memory of how it was processed. Glasses, chitecture for creating new materials. Star polymers have including glassy polymers, display memory in their relax- three or more arms radiating from a single node. There ation behavior. One intriguing question is whether such is an emerging focus on stars in which different arms are memory formation can be used for practical purposes. chemically different from one another (miktoarms) and As emphasized elsewhere in this report, interfaces are complex brush conformations (including bottlebrush). crucial for determining the elastic properties of soft mat- Such architectures change the thermodynamics of the ter. This is also true for polymeric materials. For ex- melts and have implications for the self-assembly of block ample, interfaces are important for charge transport in copolymers and polyelectrolytes where charge position- electronic polymers, where the polymer/electrode inter- ing, particularly important in batteries, is carefully con- face must be better understood in the presence of disor- trolled. Here, formulating frameworks to describe ther- der. An important goal is to create techniques that would modynamics of charge-containing systems is particularly allow interfacial processes to be visualized and character- important. Liquid-crystal elastomers are materials that ized. This would yield insight into problems such as stick- merge the orientational order of a liquid crystal into a slip and adhesion. One overarching question is to deter- cross-linked polymer network. These materials can dis- mine how much structural or chemical control is needed play extremely large, yet reversible, changes of dimension over the polymeric building blocks in order to create col- which is desirable in applications that require actuation lective behaviors of interest. It is unclear whether the (Warner and Terentjev, 2003). interfacial chemistry is sufficient to enable this control. Supramolecular polymers, in which the monomers are Synthetic polymers have a very broad range of mechan- held together by non-covalent bonds, present further ical properties and are easy to process. However they are opportunities to formulate new materials because they difficult to prepare so as to be absolutely monodisperse. can undergo reversible transitions between monomer and Building function into synthetic polymers is achieved polymer configurations. Likewise, in vitrimers polymer by adding chemical functionality to a monomer or ap- networks can be made, broken and reformed due to dy- pending a side chain with a specific functionality. Most namic covalent bonding. Applications range from recy- synthetic polymers are soluble in organic solvents, al- View Article Online

over the apparent one (r is a number larger than unity). In McCarthy experimentally observed that below a critical den- other words, the solid surface energy can be seen as multiplied sity of defects (i.e. below a critical roughness), there is indeed a by the factor r, which yields: serious deterioration of the water-repellent properties.11 cos h* 5 r cos h (1) 2.2. Examples of water-repellent materials where h is the Young contact angle, fixed by the chemical More than 200 plants, and many insects are (at least partially) natures of the solid, liquid and vapour. water-repellent (to protect themselves against water). The Eqn. (1) predicts that the contact angle on a hydrophobic morphology of the plant surfaces was studied comprehensively by Barthlott and Neinhuis, and many different designs were material (h . 90u) will increase with the roughness (h* . h). 2 This looks like a simple and attractive solution for inducing reported. However, it seems that the most efficient ones (in superhydrophobicity: the rougher the material, the higher the term of contact angle) consist of two hierarchical structures: contact angle. However, this is not that simple, for two typically bumps of about 10 mm, and submicronic microfibers. reasons: firstly, contact angles generally spread in quite a large Such structures decorate for example the surface of the lotus interval, contrasting with eqn. (1) which predicts a unique leaf, the archetype of a natural water-repellent surface. angle. This interval, often referred to as the contact angle These surfaces are thus very rough, which favours robust 12 hysteresis, is responsible for the sticking of drops, an effect in fakir states. Extending the idea of hierarchical structures contradiction with water repellency. In a Wenzel state, the naturally leads to fractal surfaces, which were achieved by the 3 contact angle hysteresis will be very large: trying to remove a Kao group, and were indeed found to be superhydrophobic. liquid makes it contact itself (owing to the fraction left in the Contact angles as high as 174u were measured on these surfaces, textures), which yields a low ‘‘receding’’ contact angle—values with a corresponding hysteresis smaller than 5u—yielding as low as 40u were reported, making this state hydrophilic-like amazing non-stick properties for water drops. Many techni- in the receding stage.5 The second reason which makes it ques for achieving disordered materials were proposed since, 13 impossible to reach high values of h*, as expected from eqn. (1) as reported by Nakajima et al. for r large and h . 90u, can be guessed quite easily: for very On the other hand, microfabrication techniques recently rough hydrophobic materials, the energy stored for following promoted much more regular structures such as pillars (as the solid surface is much larger than the energy associated with sketched in Fig. 1), and it was shown that such textures can 11,14 the air pockets sketched in Fig. 1b.5–8 also induce superhydrophobicity. Fig. 2 shows a milli- In this state (first suggested by Cassie and Baxter), the liquid metric water drop sitting on such a substrate, whose colours only contacts the solid through the top of the asperities, on a originate from the regularity of the structures. This was also 9 15 fraction that we denote as ws. If only air were present between observed for spherical regularly spaced microbeads, and it the solid and the liquid (as for a water drop on a very hot suggests more generally the possibility of taking advantage of plate), the ‘‘contact angle’’ would be 180u: the smaller ws, the the microtextures for inducing other properties than the water- Published on 22 April 2005. Downloaded 11/07/2016 23:09:11. closer to this extreme situation, and thus the higher the repellency alone. hydrophobicity. More precisely, the contact angle h* of such a Since pillar surfaces can be of low roughness, the Wenzel ‘‘fakir’’ drop (Fig. 1b) is an average between the angles on the state might be preferred. However, it turns out that very often, solid (of cosine cosh), and on the air (of cosine 21), respec- fakir drops are observed, in spite of a higher surface energy: 5,8,16,17 tively weighed by the fractions ws and 1 2 ws, which yields: the Cassie state can be metastable. This is observed in Fig. 3, where the two states are coexisting on the same cos h* 521 + ws (cos h + 1) (2) substrate, for which the density of pillars (of diameter 2 mm and height 12 mm) is about 1%. Two drops of same volume

For h 5 110u and ws 5 10%, we find that h* is about 160u. In this case, 90% of the drop base contacts air! This makes it understandable that the corresponding hysteresis is observed to be very low (typically around 5 to 10u), as first reported by 14 Johnson and Dettre:10 the liquid has very little interactions with its substrate. Hence, this statethough will be some, the like (only) polyethylene repellent oxide, are soluble in wa- ter. On the other hand, biological polymers, like pep- one, since it achieves both a largetides, contact are truly angle monodisperse, and a and small assemble into configu- rations where the three-dimensional arrangement of the hysteresis (this can be observed further,chains imparts in Fig. a specific 3). functionality to the molecule. h* monotonously increases as wTheys decreases, often bind and suggesting must be dispersed that in water to main- tain their well-defined shape. ws should be made as small as possible.There is But much reducing to be gainedw ifs thealso synthetic and biolog- makes the roughness decrease, soical that worlds we can reach be merged. the The critical synthetic polymers, with their ease of processing, could6,7 serve to protect the bio- FIG. 12 Millimeter-size water drop on a hydrophobic sur- roughness rc below which the Wenzellogical state polymer is favoured. in environmentsThe other than water. One face textured with regularly-shaped micropillars. The regular quantity rc is easily deduced fromcould the intersectionenvision synthetic-biological of eqn. (1) micellar structures in spacing of the pillars creates colored stripes. From (Callies organic or aqueous media where drugs could be hidden and Qu´er´e,2005) and (2), and is found to be (ws within2 1)/cos the organich + corews for, which delivery. isConversely, a biolog- generally close to 21/cosh (since weical will construct often could have: be storedws % within1). a polymeric matrix for subsequent use. One can imagine implants where the the soft-matter community. For example, one needs only For h 5 120u (a high value for thestructural Young integrity angle, of obtained a polymer could on be coupledFig. 2 withMillimetricto supercool water a liquid drop to on see a hydrophobic a strong increase surface in its vis- textured with fluorinated substrates), the fakir statethe biocompatibility will thus be of peptides favoured or proteins.regularly-spaced There are cosity micropillars. and a transition The into a texture glass phase. also One induces can also structural still few examples where the¨ merging of these two worlds drive a flowing liquid so hard that it becomes turbulent for roughness factors larger thanhas 2. been Conversely, effectively accomplished; Oner and the potentialcolours. is great. and new structures, such as eddies and vortex lines, form. To end this section with rheology in its title, men- In each case, the liquid is in a state that is far from tion should be made of shear thickening and shear thin- equilibrium and displays non-linear response. Even after 56 | Soft Matter, 2005, 1, 55–61 ning suspensions. It is still debated to what extent shearThisdecades journal of intense is ß The effort, Royal our understanding Society of of the Chemistry glass 2005 thickening is due to hydro-clusters (Cheng et al., 2011) transition and of turbulent flow are unsatisfactory. and to what extent it is a jamming phenomenon that re- A rich set of questions await an understanding in fluids quires an understanding of the role of surfaces (Brown that are forced into complex situations as can occur at and Jaeger, 2012). Magneto- and electro-rheological flu- a fluid interface. Examples occur when a fluid interacts ids also have rheological properties that can be varied. with a surface (such as a super-hydrophobic surface (Cal- These suspensions contain polarizable particles that re- lies and Qu´er´e,2005)) shown in Fig. 12 or when several spond to the application of an external field. For exam- fluids coexist at a line (such as at contact line on a surface ple, the application of a magnetic field can alter the rhe- (Guo et al., 2013; Kaz et al., 2012; Style et al., 2013)). ology of a magneto-rheological fluid by changing its vis- Free interfaces present the added complication that they cosity (Wereley, 2013). It has been pointed out that there can evolve with the internal motion of the liquid as shown are intricate interactions between the fluid interfaces, the in Fig. 1 and thus provide dynamically evolving, non- suspension particles, and the applied field (Orellana and linear, geometric constraints on the flow. This creates Jaeger, 2013). a highly coupled problem that can often lead to effects such as the formation of surface singularities where two liquid drops coalesces or one drop breaks apart (Eggers, VIII. FLUIDS 1997; Shi et al., 1994). Some of this singular behavior has been worked out in recent years but certainly much The subject of fluids has a long and rich history that remains to be done especially with non-Newtonian fluids. has impacted and drawn from many disciplines. Histor- Beyond Newton: Much previous research has dealt ically, fluid dynamics was one of the frontiers where the with simple “Newtonian” fluids (those that can be ad- mechanics of continuous media was unraveled and contin- equately described by a viscosity, a density and a surface ues to be a cornerstone for various branches of chemistry tension). The future of fluid dynamics in soft matter calls and physics. Fluid flows are essential for creating many for going beyond this restriction and understanding non- mundane marvels of our daily existence and for produc- Newtonian fluids and the effects that complex molecules ing phenomena observed in more esoteric environments can exert on the bulk and surface properties of a liquid. encountered in geophysics, astrophysics and engineering. It is well known that a variation of the interfacial ten- There is an important overlap between researchers in sion will set up complex flows along a surface. These fluids and those in soft matter. Fluids are involved in Marangoni flows can lead to many counter-intuitive phe- many aspects of soft matter: colloids, foams, liquid crys- nomena such as the tears of wine often observed at dinner tals, polymers, self-assembly, active matter, biomaterials, parties. There, above the wine’s surface, a thin band of and microfluidics. Perhaps more significantly, as empha- wine (composed of alcohol and water, two liquids with sized in the introduction, many scientific issues that are different surface tensions) drains down the side of the important for understanding fluids have resonance within glass in the form of tears. An example is shown in Fig. 13 ARTICLES NATURE PHYSICS DOI: 10.1038/NPHYS2560

ac a Micro-bubbles

b Smoke

25 mm de t = 300 ms t = 155 ms t = 105 ms t = 60 ms

222 A. E. Hosoi and J. W. M. Bush

b

f

Figure 3. Photo of the tears and ridges: an oblique perspective of the meniscus region between the reservoir (left) and the adjoining thin film (right) on a plate inclined at 4 relative to the horizontal. Note the three descending tears, the surface deformations corresponding to the ridges, and the forking on the large central ridge which is beginning to exhibit dendritic structure. The horizontal 15 g scale is 8 cm and the fluid is 65% methanol and 35% water. h

25.0 mm 18.8 mm 12.5 mm 9.4 mm

Figure 4. A plan view illustrating the reservoir (bottom), the thin film region (middle) and the FIG.incipient 13 tear line Tears (top). Kalliroscope of wine. reveals The convection image throughout. shows The the ridges tear are theline four is white at triangular structures at the lower right. Each ridge has a clearly visible dark spine which corresponds theto a line top, of surface a thin convergence film region and downwelling. in the The middle, fine-scale convection and four is marked ridges by light at thestreaks bottom further up the right. plate, and From the irregular (Hosoi convection and inBush, the deep 2001). reservoir can be seen in the lower right corner. The horizontal scale is 6 cm and the fluid is 65% methanol and 35% water. c The ridges were found to align for any non-zero slope angle, and to exist even when the bounding wall was vertical, which is clear evidence that the convection is (HosoiMarangoni and rather thanBush, Rayleigh–B 2001).enard:´ The if the convective problem instability becomes were buoyancy com- driven, it would disappear on a vertical plate. A series of experiments was also pounded when large surfactant molecules are introduced at a fluid interface. The surfactants alter the local sur- face tension creating fluid motion, which in turn causes FIG. 14 Vortex loop in water created by an accelerated air- the surfactants to flow. Research in this area continues foil with the topology of a knot. From (Kleckner and Irvine, Figureto be vibrant. 1 The The creation coupled problem of ofvortices surfactants inter-with designed2013). shape and topology. acting with the flows in the liquid leaves a great deal to 50 mm be understood| especially when the boundary geometry a, Theof the conventional medium is involved. method for generatingassembled a vortex structure, ring, have in unusual which elastic a properties. burst of In fluidAdding is forced large active through molecules an into orifice. the bulk ofb a liquid,, A vortexthese ringmaterials, in theair behavior visualized and response with depend on as shown in Fig. 5B, creates internally driven motion in how the many components are assembled rather than on smoke.which topologicalc, A vortex defects ring in the in order water parameter traced inter- bythe a properties line of of ultrafine the individual gas components bubbles, themselves. act directly with the flows (Sanchez et al., 2012). An- The advantage of such a design strategy is that one can d 15 whichother show example finer of topological core considerations details than interacting smokeattempt or dye. to created,e materials, A vortex with a ring set of candesired prop- with fluid dynamics is in the case of vortex formation. erties in mind that may be difficult to find in nature. alternativelyMuch as magnetic be field generated lines, vortex lines as the either starting extend An vortex example thatof a has suddenly become popular accelerated, is to create mate- to the edge of the system or wrap around until they join rials with negative Poisson ratios (i.e., auxetic response

speciallyback on themselves. designed This wing. is similar For to the a complex wing knot- within the which trailing compressing edge a material angled along inward, one direction de- 1 10

ting created out of defect lines in liquid crystals men- creases the other two dimensions as well) (Greaves et al., t 2 the starting vortex moves in the opposite of the direction of wing motion 0 ¬ tioned previously (Tkalec et al., 2011). Such vortex lines 2011). While negative Poisson ratio materials do exist r can become highly entangled and it can be asked whether naturally, they are less common than the normal approx- U f, Thea fluid starting can evolve so vortex that the vortexis a result lines reconnect of conservation with imately incompressible of circulation—the ones. It is an open question bound about = ∗ themselves allowing the tangle to relax to a simpler ge- how much freedom one has to design in multiple distinct t 5 circulationometry. Studies, around as shown a inwing Fig. 14, is show balanced how vortex byproperties the counter-rotating into a material. starting Linked rings vortex.rings cang, beh, tied A intorendering topological knotted of a wing structures tied that intoThe a knot disadvantage (g), of used meta-material to generate design is that a one can subsequently relax (Kleckner and Irvine, 2013). has to develop a means to guide the assembly of pieces Trefoil knot knottedThere vortex are many industries(h). that design fluids with or- into the determined structure. When the individual thogonal required functionalities (e.g., drilling mud, con- pieces themselves are microscopic, then this can be a 0 ditioner, or toothpaste) but these are fluids are typically challenge. However, if the pieces can be coaxed to self 0 20 40 60 80 designed empirically: “add things until it works”. There assemble in the proper manner, or if only their shape 20 3 wasare proved currently no universal. More design recently, principles to achieve thesebut solutions not their precise placementwere shown is important, to this diffi- be Re/10 (Re = U0Ch/)ν multiple functionalities in a single fluid. The creation of28 culty could be circumvented and a significant problem unstablesuch an overarching to linear set of rules perturbations is one outstanding chal- ;would however, be solved. With when the advent the of LIA 3D printing is tech- no lenge for the field. nology many different structures can be built at least longer valid, the interplay of global andat a moderately local induction small scale. One complicates challenge is to push Figure 2 Scaling of trefoil knot vortex loops. a, An overlay of vortex knots such technology to create much larger structures with | theIX. situation MECHANICAL considerably META-MATERIALS and the evolutionever finer features. of vortex Even if that knots were to remains become possi- taken at the same rescaled time. The vortex r.m.s. radii are ¯r 60, 45, 30 ble, it can still be asked whether such technology could = a matterA rapidly growing of theoretical research area deals with debate the design bothmass-produce in the material inviscid in sufficient quantity and to viscid be tech- and 22.5 mm for red, yellow, green and blue colouring, respectively. The and21 fabrication,22 of composite materials that, due to their nologically useful. 1 cases . Simulations of knotted vortices with large cores, based generating wing speed is 3.10 m s for all except the largest vortex knot, 29 1 on Navier–Stokes dynamics, suggest that they are short lived ; for which it is 2.15 m s . The biggest vortex (red) is slightly larger than the however, numerical integration of Biot–Savart vortex evolution imaging field of view, resulting in some clipping on the left edge. b, The suggests that finite core size may enhance stability21. Quantitatively, same knots as shown in a, scaled inversely proportional to the original capturing the details of even the evolution of a simple ring remains hydrofoil dimensions. c, A photograph of the four knot-generating wings. 30 a surprisingly resilient problem . Resolving such subtle questions d, The rescaled time, t⇤, of the first reconnection event. With the exception therefore requires both the generation of knotted flows in the of the smallest linked rings (which are affected by background flow from laboratory, and an effective means of three-dimensional (3D) the apparatus), we observe that the rescaled reconnection time is imaging of the evolution of the resulting flow with a high degree independent of the Reynolds number. of spatial and temporal resolution. The conventional method for making a vortex loop is to force a to create a knotted vortex23, but to our knowledge this has never burst of fluid out of an orifice (Fig. 1a–c and Supplementary Movie been demonstrated experimentally (our own attempts indicate that S1). It has been suggested that two perturbed rings could be collided the strong perturbations to the shape of a vortex resulting from

254 NATURE PHYSICS VOL 9 APRIL 2013 www.nature.com/naturephysics | | | © 2013 Macmillan Publishers Limited. All rights reserved PHYSICAL REVIEW LETTERS week ending PRL 99, 084301 (2007) 24 AUGUST 2007

tained through a mesh refinement study. The elastomeric ing mutually orthogonal ellipses above a nominal strain of stress-strain behavior was modeled as a nearly incompress- 0.04. This corresponds to the plateau region immediately ible neo-Hookean solid with a shear modulus of 3.25 MPa. after the departure from linearity. The pattern transforma- The sheet structures were modeled in two manners: (1) con- tion is a result of the compressive loading of the vertical sidering a representative volume element (RVE) of the interhole ligaments which undergo a buckling instability, domain with appropriate periodic boundary conditions thereby triggering the change to the new configuration. and containing small internal defects to aid the initiation This is a truly elastic event and is completely reversible of instabilities and (2) considering the full specimen, thus capturing any small geometry differences at the boundaries as well as any boundary effects on either constraining (top Experiments Simulations or bottom surfaces) or triggering (free side surfaces) the transformation. Experimental and numerical results for the nominal stress versus nominal strain behavior of the square array of circular holes, up to a compressive strain of 0.10, are shown in Fig. 1.Thebehaviorischaracterizedbyaninitial linear elastic behavior with a sudden departure from line- arity to a plateau stress. The departure from linearity is a result of a sudden transformation in the periodic pattern as =0.02∋ shown in the snapshots of deformed configurations at 0 MPa 0.001 different strains (Fig. 2). Figure 1 also provides a direct comparison between the numerical and experimental stress-strain results. There is excellent quantitative agree- ment between the full model calculation and the experi- ment. During the first loading in the experiment, the instability is abrupt and characterized by the sharp load drop seen in Fig. 1; subsequent experimental cycles do not exhibit the load drop and coincide more directly with the full model simulation result. =0.04∋ The experimental and numerical images of the deformed 0 MPa 0.06 configurations are shown in Fig. 2. A gradual and homoge- neous compression of the circular holes pattern during the linear elastic range of the deformation is replaced by a transformation to a strikingly different pattern of alternat-

=0.06∋ 0 MPa 0.11

16

A variety of routes to meta-material design are cur- rently being exploited to control the shape, mechanical properties or actuation of a material. One of these is based on the paradigm of origami that models three- dimensional structures by folding a two-dimensional sheet. (If the sheet has prescribed cuts this is called kirigami.) Because a fold along one crease constrains the ability of the sheet to distort elsewhere, an open theo- retical challenge is to develop the rules for folding that do, or do not, produce bending of the flat sections of the FIG. 15 (A) Thick sheet with circular holes. (B) This geom- =0.10∋ sheet. It is important to understand what are the degrees etry produces a negative Poisson ratio well into the non-linear of freedom accessible to a folded sheet (Silverberg et al., regime; as the material is compressed in one direction, it is 0 MPa 0.15 2015). From the experimental perspective with which compressed in the perpendicular direction as well. The ma- terial is birefringent under strain so that the stresses can be this report is predominantly concerned, one wishes to visualized. From (Mullin et al., 2007). create such objects efficiently from the designs that have FIG. 2 (color online). Experimental (left) and numerical been evolved theoretically or computationally.FIG. 1. Nominal More ex- stress vs nominal strain curves for the square (right) images of the square lattice at different levels of macro- citing would be to find ways in whicharray the showing material itself experimentalmaterial by and careful computational manipulation of the (full bonds. model and scopic strain: 2%, 4%, 6%, and 10%. The experimental samples can develop the desired folds to allow the emergence of The field of meta-materials is still in its infancy and unusual and desired behavior. RVE model) results.it The is impossible departure to tell from with linearitycertainty where is the research result will of were viewed through cross polars and the colors give a qualita- Using evolutionary computer algorithms,an elastic one buckling can sim- instabilitylead. It is possible in the microstructure to design materials that with triggers a desired a tive indicator of the regions of localized stress which provide a ulate what shapes of particles canpattern lead to a transformation. desired rhe- property either in its linear or in its non-linear response. useful comparison with the calculated principal stress levels. ological behavior. One wants to discover a microscopic Both behaviors could be useful in different contexts. Go- building block that will achieve a desired outcome. This ing beyond linearity poses an overarching challenge. The has been done in designing a shear-thickening granular material often undergoes catastrophic failure once the084301-2 material with an unusual stress-strain curve in which linear regime has been exceeded. Yet one would often the material gets progressively stiffer as the material is like to avoid failure but still make use of large-scale re- sheared (Miskin and Jaeger, 2013). How far can such sponse. A beautiful example of a material that has a large ideas be pushed to create novel and useful materials non-linear auxetic response is shown in Fig. 15. A square (O’Hern and Shattuck, 2013)? lattice of holes drilled into block will break the symme- Another tactic has been to design properties in net- try under uniaxial compression and produce a negative work structures. Pruning specific bonds in computer- Poisson ratio for large amplitude deformations (Florijn simulated networks derived from jammed packings allows et al., 2014; Mullin et al., 2007). the Poisson ratio to be varied at will (Goodrich et al., Another goal is to create materials that can transform 2015). Again, from an experimental point of view, the their shape and function in response to a stimulus. One issue is how the material can be coaxed into pruning the example is a self-folding surface such as occurs when an appropriate bonds itself without resorting to a simula- ultra-thin sheet encapsulates a liquid drop as shown in tion. Such results raise the issue about what are the Fig. 16. This large-scale motion takes advantage of the limits to which function can be programmed into a ma- ability of thin sheets to have large distortions without an terial. The Poisson ratio is a global property. Can similar enormous energy penalty (Paulsen et al., 2015). control be obtained at a local level, for example, as ex- By playing one set of forces off against another, one emplified by biological proteins with allosteric behavior? can add functionality to a material . When they are in There, binding a molecule to one specific site determines balance small shifts in steady-state conditions can lead whether a far-distant site will have the ability to bind to to one set of forces overcoming the others. If one member a third molecule. Can similar allosteric responses be de- buckles it can set off an avalanche of activity. One pos- signed into physical materials? The answer appears to be sibility is to pit surface tension against elastic forces at in the affirmative (Rocks et al., 2016). Likewise material an interface. Working near a phase transition is another failure how cracks propagate when a material is stressed way that a small change in a control parameter (e.g., to the point of failure also appears to be related to the pressure, temperature, shear stress) can cause a colossal average number of constraints in the material (Driscoll change in a material. One can design inhomogeneity into et al., 2016). Topological meta-materials also require un- a material to provide local control. derstanding how the degrees of freedom (allowed parti- Taking lessons from biology, such as in the case of al- cle motions) interplay with the inter-particle constraints lostery mentioned above, can lead to new forms of behav- (Kane and Lubensky, 2014; Nash et al., 2015; Paulose ior. Basic understanding of how a cell moves or muscles et al., 2015). These examples, based on understanding contract can provide ideas for how similar mechanisms the constraints in a network, may be just the opening can be employed in physical systems. On an even grander wedge of what kinds of functions can be imparted to a scale, one can conceive of active matter forming the en- 17

241 nm

FIG. 16 Side and top views of a circular polystyrene sheet wrapping a water drop immersed in silicone oil. Sheet thick- ness is 241 nm. Scale bar is 1 mm. From (Paulsen et al., 2015). tities that comprise a meta-material. This has been just a taste of what approaches can be used to add functionality to a material by assembling and manipulating the elements at a mesoscopic level. The possibilities seem limitless. FIG. 17 Bird flocks over the Rome train station. From (Cav- agna et al., 2010).

X. ACTIVE AND ADAPTIVE MATTER as shown in Fig. 17 or in schools of fish. The flocking rep- The study of active matter seeks to attain a fundamen- resents a non-equilibrium transition into a phase that has tal understanding of collective properties that emerge in coexistence of high- and low-density regions and corre- an ensemble of particles or agents each of which consumes lated motion of the agents (Cavagna and Giardina, 2014; energy to propel itself. Motion arises naturally in such Marchetti et al., 2013). Introducing active agents gener- systems and their understanding requires a broad per- ates a new degree of complexity and produces new fea- spective that unifies the global features of the flow with tures not observed in matter based on more conventional the microscopic propulsive mechanisms responsible for inter-particle interactions. the activity (Marchetti et al., 2013). We have referred Active liquid crystals represent a synergy between dif- to them throughout this report in the sections on col- ferent soft-matter efforts including fluid dynamics, rheol- loids, polymers, liquid crystals, fluids, mechanical meta- ogy, liquid-crystals, polymer physics, biophysics and ex- materials, self-assembly, pattern formation and granular treme mechanics including the study of fracture. One materials. The challenge is to determine the fundamental beautiful example of an active liquid crystal is shown in difference between driven non-equilibrium systems, such Fig. 5. Here the active liquid crystals show buckling, fold- as turbulent flows, where energy is typically injected at ing, internal fracture of a nematic domain and a pair of the macroscopic system size and then is dissipated at a topological defects. Another example of self-organization smaller level, and active matter where the energy is in- in a synthetic colloidal suspension is shown in Fig. 18B jected at the microscopic scale and proceeds to generate of a two-dimensional “living crystal” created from light- motion over large distances. (Figure 5A shows a case that activated artificial surfers (Palacci et al., 2013). perhaps lies between these two extremes. The entire ap- It is interesting to compare the generality of active sys- paratus is vibrated, but the individual particles couple tems with biological organisms, which are also inherently locally to this overall input of energy in different ways out of equilibrium. Biological systems usually employ a so that they seemingly have their own energy source.) limited number of structural and dynamical motifs in or- Active matter is inherently far from equilibrium in an der to build diverse functionalities. These include cell unexplored and interesting fashion so that much of our division, motility, and reproduction. By contrast, ac- intuition built upon experience with inactive systems is tive matter goes beyond biology by putting components bound to be inadequate. This is a new frontier between together in novel ways that had not previously been ob- different subfields of soft matter. served. It has less constraints. Studying active matter The most obvious examples of active matter come thus perhaps has the potential to teach us ways of form- from the realm of biology. On the molecular scale, self- ing structure beyond what biology has done. organization occurs in actin filaments or microtubules One important area is the study of active gels such and controls the structure of cells; on a much larger scale, as occurs when molecular motors using chemical energy the agents are the animals themselves and their activity drive mechanical motion in a network of cross-linked fil- leads to the stunning and mesmerizing flocks of starlings aments. These active gels are important in the dynamics 18

B! C! tals used in the semiconducting industry, by minimizing their free energy; kBT is important. Other structures A! form spontaneously. For example, micelles self-assemble in colloidal suspensions of surfactant molecules. Other structures such as vesicles or membranes can also form (Witten and Pincus, 2004). One set of open questions in self-assembly now con- cerns the rules of assembly when kBT is not important (Cheng et al., 2006; Hu et al., 2014). An increasingly common route to fabricating desired structures is to use

FIG.D! 18 (A) ScanningE! electronF microscope! imageG! of a bima- dynamics as a tool to force the components to assemble terial colloid with a protruding hematite cube. (B) Formation in certain directions or with certain patterns. There are of a living crystal. When illuminated, the particles assemble a variety of dynamical variables that can be controlled. into crystals. Inset shows initial homogeneous distribution. In addition to temperature, these include vibration (for (C) When the light is turned off, the crystals melt. Main large particles such as in granular media), shear, com- panel shows after 10 seconds; inset shows after 100 seconds. pression, external fields, optical tweezers, and using ac- From (Palacci et al., 2013). tive particles (including robots for larger entities for ex- ample on an architectural scale (Augugliaro et al., 2014)). and structure formation in cells. This is an attractive Particle shape is clearly an important parameter that paradigm for understanding how localized sources of en- controls the aggregate structure (Glotzer and Solomon, ergy can lead to novel and unexpected behavior. A cur- 2007). A better language than we currently have, in- rent challenge is to use the ideas and theories developed cluding a better understanding of the statistical mechan- in the biological realm to branch into other active mate- ics and thermodynamics for these far-from-equilibrium, rials. This will require developing complementary model driven, or glassy systems, is necessary to describe such systems: synthetic active particles, living organisms, re- evolution. For example, one would like to generalize the constituted systems and active granular systems. concept of minimizing the free energy to bring the col- Making new microscopic and/or macroscopic units, lective system to a specific state. which transduce energy, remains a challenge. It will re- Bringing ideas and tools from biology, where assem- quire establishing new connections between chemists who bly of specific materials and functions is constantly oc- can make new motors and biologists who can isolate, pu- curring, into soft matter is a continuing goal. Likewise, rify and modify biological force-producing units. New programmable self-assembly would allow the construc- characterization techniques will also be needed. As high- tion of more nuanced materials. For example, large two- lighted in the section below on instrumentation, new rhe- dimensional lattices have been constructed by manipu- ological capabilities would be an important advance. lating the interactions at different length scales by using biological molecules such as DNA (Reif et al., 2001; See- man, 2003; Wang et al., 2015). XI. SELF ASSEMBLY

Providing functionality to materials is an outstanding XII. PATTERNS AND STRUCTURE FORMATION challenge that incorporates the chemistry of the compo- nents and the structure of their assembly. Thus, a com- Patterns emerge from an apparently featureless back- mon goal of many areas of material science is to assemble ground. They are everywhere and many have their prove- matter into useful structures. It is particularly attrac- nance in soft matter. They play an important role in tive if the sub-units of the material (atoms, molecules, many subject areas outlined in this report: liquid crys- colloids, etc.) were to self assemble into a desired ar- tals, active matter (see Fig. 5), fluid mechanics (see rangement with as little external control exerted as pos- Figs. 1, 13 and 14), foams, packings and colloids. They sible. It has long been debated as to what constitutes are the outcome of design using origami and kirigami. self-assembly since there is always some degree of exter- Certain patterns are familiar and appear recurrently nal control exerted in setting up an experiment. One in different guises. For example, although not identical seems to know self-assembly when one sees it. in all cases, branching appears in river networks, viscous An atomic crystal may be said to self-assemble as its fingering, blood vessels and of course in tree branches constitutive atoms arrange themselves into the lowest and dielectric breakdown, (see Fig. 3). Another pervasive free-energy state. But, even there, considerable care and example is based on the Turing reaction-diffusion mech- expense must be lavished to create pure, nearly-perfect, anism (Turing, 1952) which have been demonstrated in crystals of e.g., silicon. Soft-matter can also self-organize chemical reactions (Ouyang and Swinney, 1991). An ex- into ordered arrays, albeit not as perfect as the crys- ample is shown in Fig. 19A. It has been suggested that J. Phys. Chem. B 2009, 113, 3785–3798 3785

Modulated Phases: Review and Recent Results†

David Andelman* Raymond and BeVerly Sackler School of Physics and Astronomy, Tel AViV UniVersity, Ramat AViV, Tel AViV 69978, Israel

Ronald E. Rosensweig‡ 34 Gloucester Road, Summit, New Jersey 07901 ReceiVed: September 1, 2008; ReVised Manuscript ReceiVed: September 20, 2008

19

A face tension and destabilizing viscous stresses (Saffman and Taylor, 1958). A related example occurs when the

non-linear dynamics of the system are controlled by the We consider aspects of patternings that occur in a wide array of physical systems due to interacting combinationspresence of a nearby singularity such as occurs when one liquid drop fissions into two or more smaller pieces (Eg- of dipolar, interfacial, charge exchange, entropic, and geometric influences. We review well-establishedgers, 1997; Shi et al., 1994). In that case, the Rayleigh- Plateau instability sets the scale for subsequent droplet phenomena as a basis for discussion of more recent developments. The materials of interest range from bulkformation. In other cases, such as the diffusion-limited aggregation in Fig. 4, the structure is apparently fractal with no characteristic length (Witten and Sander, 1981). inorganic solids and polymer organic melts to fluid colloids. Often, there are unifying principles behind theThe characterization of a pattern requires determin- ing its symmetries and conservation laws and identify- various modulated structures, such as the competition between surface or line tension and dipolar interactioning what order parameter describes its essential behav- ior. Sometimes the order is not obvious and new metrics in thermally reversible systems. Generally, their properties can be understood by free-energy minimization.need to be invented to describe them as in the case of hyper-uniform density fluctuations (Torquato, 2016) or in the wavy edges of a tree leaf (Sharon et al., 2007). Of- ten the patterns are caused because the system is driven B far from equilibrium, as in turbulent flow, so that new I. Introduction structures appear (Zocchi et al., 1990). Unlike amorphous materials, crystals have well-defined defects (Chaikin and Lubensky, 2000; Kleman and Lavrentovich, 2003). These defects, which often have a A diverse number of physical, chemical, and biological topological character, control deformation as the material responds to strain and determine the patterns that can systems exhibit some type of modulation in their structural be formed. This is seen dramatically in lamellar layers of block co-polymer or liquid-crystal films. The behavior of 1 the defects is intimately tied to the presence and location properties. Examples of such structures in two-dimensional of surfaces in finite-sized materials. In liquid crystals, the way in which the director field interacts with the bound- (2D) systems are elongated stripes and compact droplet-like ing surface of the liquid often requires that the defects must exist in the interior. Such defects play a special role domains, as can be seen in Figure 1. In the figure, domains in FIG. 19 Labyrinthine patterns. (A) Turing pattern. Sta- in the dynamics. tionary chemical patterns formed in a gel layer with striped Far from being a closed field, new patterns are fre- solid magnetic systems (garnet films) and in thin layers of patterns with wavelength 0.52mm. From (Vigil et al., 1992). quently being discovered. The example of proportionate (B) Modulated phase. A thin layer of ferrofluid in a mag- growth shown in Fig. 20 is an example of a pattern that netic field with a period ≈ 2mm. Adapted from Fig. 4 of had not been observed in the physical realm prior to re- ferrofluids (to be discussed in detail below) are shown side- (Rosensweig et al., 1983) as in (Andelman and Rosensweig, cent experiments on miscible viscous fingering (Bischof- 2008). berger et al., 2014). Other examples are surely waiting to by-side and exhibit striking similarity. In three-dimensional (3D) be discovered. Likewise, well-known mechanisms to cre- ate patterns, such as the reaction-diffusion mechanism of systems, the domain morphology can be more complex and such Turing patterns are responsible for many biological Turing, constantly find new applications in nature. In- shapes and forms such as pigmentation in the coats of an- deed biology is a fecund place to search for a variety of includes sheets, tubes, rods, and droplets embedded in a three- imals (for example zebras or leopards) and patterns on structures – both new and old. In order to take advan- sea shells. Similar looking labyrinthine patterns are seen tage of many of these opportunities and to search for new in modulated phases (Andelman and Rosensweig, 2008). forms of structure, it is important to understand the or- dimensional matrix. The similarity between the resulting patterns One example is in thin layers of ferrofluid (suspensions ganizational principles that appear in dissipative systems Figure 1. Domains in magnetic solidsof a magnetic and solid fluids. in a fluid solvent) (a) in Stripes a magnetic field andand (b) that govern non-equilibrium assembly processes. in certain systems of different origins is quite surprising and shown in Fig. 19B One avenue for productive exploration is the use of pat- bubble phase in a ferromagnetic garnetIn some film cases of the initiation 13 µm of patterns thickness is caused by growntern formation as a means of imparting functionality to a may allude to a common unifying mechanism. In other cases an instability as some control parameter is varied past material. Again, one can take lessons from biology. For on the 〈111〉 face of gadolinium galliuma threshold garnet. value. When Visualization there is a most unstable wave- is madeexample, inspired by embryogenesis, one could attempt additional mechanisms arise, as will be discussed. length at the threshold, a single length scale character- to use active growth patterns to localize buckling and using polarized optical microscopyizes the (Faraday onset of patterns effect). throughout the Period: system. Thisd*folding so as to create and control new structures. One occurs, for example, in viscous fingering where the finger could∼ build shape-morphing structures such as PDMS It is convenient to view these systems as systems which, due 10 µm. Adapted from ref 2. (c) Ferrofluidwidth is set by confined a competition between between stabilizing two sur- glassgels or ultrathin sheets that are able to encapsulate liq- to a competition between different interactions, achieve ther- plates exhibiting labyrinthine instability in a magnetic field.3,4 The period is d* 2mm.(d)Bubblephaseofaferrofluidconfinedinacellhaving modynamic equilibrium in a state in which the appropriate order ∼ parameter shows a spatial modulation. Examples are abundant1 agapthatincreasesfromlefttoright.Themeanbubblesizeis 1 mm. Adapted from ref 5. ∼ and include modulation of the magnetization field of ferromag- netic slabs6-8 and ferrofluids,3 the polarization field in electric 9,10 dipolar systems and certain liquid crystalline phases, the phenomena when these systems are subjected to an external superconducting order parameter in the intermediate phase of field (electric, magnetic), and we describe their morphology, 11 type I superconductors, as well as the relative composition in structure, phase separation, various instabilities, and related 12-14 block copolymer systems. phenomena. We present in this paper some of the interesting phenomena associated with modulated phases. We start by considering a II. Domains in Two-Dimensional Ferromagnetic Layers simple example explaining the underlying mechanism of Ferromagnetism is an important physical phenomenon as- wavelength selection in a quasi two-dimensional dipolar system. sociated with elements like nickel, iron, and cobalt, as well as We then address domains in magnetic garnet films and related a large number of metallic alloys that show spontaneous dipolar organic films at the water/air interface (Langmuir magnetization M in the absence of an external applied magnetic monolayers). Two other examples of systems of current field. scientific interest having many applications are subsequently In bulk magnetic systems, the uniform magnetization does discussed, mesophases in block copolymers and magnetic fluids not persist throughout the system but breaks up into spatial (ferrofluids). Competing interactions create interesting new domains, each having a specific and distinct magnetic orienta- tion. The domain size and its structure depend on competing † Part of the “PGG (Pierre-Gilles de Gennes) Memorial Issue”. * Corresponding author. E-mail: [email protected]. interactions inside of the magnet: the direct exchange interaction, ‡ E-mail: [email protected]. the demagnetization fields, and the crystal anisotropy.

10.1021/jp807770n CCC: $40.75  2009 American Chemical Society Published on Web 12/05/2008 20

A ular media simulators can vary particle shape, gravity and friction coefficient; for polymers they can vary chain length and cross-linking; for liquid crystals they can vary

molecular architecture. Because soft-matter is classical, B one does not have to resort to quantum algorithms to do an accurate simulation. The sophistication of both com- puter hardware and modeling techniques continues to in- crease rapidly so that simulations have an ever-increasing role to play in discovering new soft-matter phenomena.

Many groups now use simulations as a standard tool rather than as a research focus. The codes for molecular dynamics simulations are widely distributed and used. However, as investigators move towards more parallel processing, there is less unanimity. Some groups use FIG. 20 Proportionate growth in miscible viscous fingering. (A) The blue (darker) fluid displaces a more viscous fluid with graphic processor units (GPUs) and others use standard which it is miscible. The interface between the two fluids at an computers with e.g.,, the LAMMPS code (Frantzdale early time can be enlarged so that it lies on top of the interface et al., 2010). In simulating liquids, there are again a va- at a later time. Even the small features are nearly identical riety of techniques to simulate, for example, the Navier- after the images have been scaled. (B) This form of pattern Stokes equations. These are rapidly getting more pow- is common in the biological world (e.g., mammalian growth), erful with the ability to handle more complex problems but few physical examples exist. From (Bischofberger et al., with multiple fluids. 2014). There are, of course, new challenges. In order for a simulation to be reliable, it must be fully validated. In uids droplets (see Fig. 16). order for a simulation to be able to demonstrate novel One could also try to create functionality on new physical effects, it is essential to know what features must length and time scales. One possibility would be to be included at the microscopic level such as inter-particle employ agents with feedback that could cycle between potentials and at the macroscopic level such as boundary different states of the system. These could provide re- conditions. The best way to validate a code is to com- programmable devices or the ability to create decision- pare the simulations with benchmark experiments. This making materials. must be done on a field-by-field basis; it is valuable to have standardized codes where possible and to identify and agree on canonical experimental systems (boundary XIII. SIMULATIONS AND BIG DATA conditions, external parameters) and make sure that the data generated is reproducible. In soft matter, the syn- The sheer amount of data collected in a single soft- ergy between experiment and simulation has been ex- matter experiment can be staggering. For example, cellent; strong collaborations between experimental and many experiments depend on high-resolution, high-speed simulation groups should continue to grow. video. Manipulating a data file for even one such movie in order to identify features over hundreds of thousands While some large databases exist, such as JHTDB of frames is a daunting challenge. The era of big data (Johns Hopkins Turbulence Databases) (Li et al., 2008), has arrived in soft matter. the soft-matter community does not generally possess In addition, soft-matter systems are so complex that it and query large databases that can be used to shorten is often difficult to extract from experiments the essential the time and cost of materials development. A better roots of their distinctive behavior. As the understanding infrastructure is needed to create additional databases of of a system grows, a challenge is to devise novel methods soft-matter experimental and simulational N-body sys- that can identify new objects or descriptors that cap- tems. The question then arises of what type of data ture the essence of the underlying science. Significant should be included in such databases. At one extreme, advances along this line have been made including the one could include all the microscopic positions and inter- use of machine learning, evolutionary algorithms, com- particle forces; at the other, one could include only the putational homology and network methods. macroscopic outputs such as stress versus strain. Be- Computer simulations have an advantage over experi- cause the system response can be spatially inhomoge- ment in that they can readily investigate a series of model neous, each deformation can be different and one must systems by selectively adding or subtracting degrees of know the complete history of a system and the spatial complexity in the interactions between particles or by dependence of its response. Much can be learned from selectively tuning control parameters. Thus for gran- experience with previous databases. 21

XIV. INSTRUMENTATION

As with any branch of experimental science, signifi- cant and rapid progress often depends on the creation and development of new instruments that allow mea- surements that had previously been impossible. Soft- matter science has had an excellent history of develop- ing tools that have allowed new and precise methods of probing and manipulating matter at the mesoscopic scale. These include the invention of diffusion-wave spec- troscopy (DWS) (Maret and Wolf, 1987; Pine et al., 1988) making it possible to detect minute changes in the struc- tural arrangements of a macroscopic colloidal sample on a rapid timescale. Sophisticated laser-tweezing technol- FIG. 21 Confocal imaging of a suspension of droplets. The ogy (and related methods that use the forces from a light bright spots, where the green (lighter colored) droplets con- beam) can now manipulate enormous numbers of parti- tact each other, allows contact forces between the suspension droplets to be measured. From (Bruji´c et al., 2007). cles simultaneously (Grier, 2003). The emergence and widespread use of high-speed confocal microscopes allows precise depth resolution of structure and the tracking of individual particle positions in both soft and biolog- even to detect the transient behavior adequately, it can ical matter. Other means of visualization include holo- be necessary to extend experiments over exceedingly long graphic microscopy (for reconstructing three-dimensional times. Some spectroscopic measurements have studied volumes), magnetic resonance imaging (for imaging par- temperature-controlled glasses over periods of a month ticles and fluids containing oil or water), X-ray tomog- (Leheny and Nagel, 1997); it is prohibitive to extend that raphy, high-speed video photography, photoelastic visu- time significantly. The problem is compounded when the alization (for imaging the stress-induced birefringence of range of times is coupled with an equally large range of a two-dimensional sample), and light-sheet fluorescence lengths from 0.1nm to 1m. That is we are asking (!) for microscopy (for mapping three-dimensional particle posi- techniques that allow an order of magnitude in the order tions from two-dimensional slices). Scattering techniques of magnitude of lengths and times. have also led to important advances; for example, neu- The tactic often used in a statistical approach to a tron scattering allows the examination of individual poly- many-body system is to concentrate on average quan- mer chains by deuterium labeling and neutron reflectivity tities. However, in soft-matter systems, it is often im- provides sub-nanometer resolution of soft matter at sur- portant to resolve the behavior on the individual parti- faces and interfaces. Soft X-ray scattering enables a de- cle scale as well. The data sets become enormous when termination of the structure and morphology of complex one needs to detect many particles simultaneously over soft-mater systems with exceptional spatial resolution. long times. It is still an art to handle such big data This is not an exhaustive list. sets and new computers and algorithms are needed to Challenges remain in order to make the next advance in make such analyses routine. This is not only a prob- experimental capability. While one cannot predict where lem faced by the soft-matter community, but by other the next breakthrough will occur, we can lay out some of research and development areas as well. A synergy the areas in which progress would be welcome. with the computer-science community has already had It is common to image a sample in two dimensions large benefits in soft-matter regarding simulation tech- or at the surface of a three-dimensional material. It is niques; graphic-processing-unit (GPU) computers, origi- more difficult to resolve simultaneously the dynamics and nally developed for a different end use, has made parallel- stresses acting at the particle level in the interior of a processing easier and faster. Such hardware advances sample. One of the most used tools in colloidal science are also useful for implementing image-processing algo- is the confocal microscope which allows unprecedented rithms. Likewise, open-source software has provided so- spatial resolution. While commercial machines can now lutions to many programming challenges. operate at 1000 frames per second, higher speeds are of- As has been emphasized repeatedly throughout this ten necessary to capture fast dynamics. Figure 21 shows report, the activity of a soft material often involves its a confocal-microscope image that led to a measurement interfacial properties. The motion of contact lines is but of contact forces in a static suspension. One challenge is one example that has many consequences (Guo et al., to create a technique that can make such measurements 2013; Kaz et al., 2012; Style et al., 2013). Studying in- rapidly on a moving material in an opaque medium. terfaces requires tools that are specifically sensitive to As remarked in the introduction, soft matter is often the surface to obtain information about the local struc- also slow matter. Thus to watch a sample equilibrate, or ture, local rheology and local excitations. Being able 22 to image and characterize buried interfaces is a difficult of the structure-function relation from the sub-cellular to problem despite advances in tools. The combination of the collective scale of many cells. These are just a few rheology with other probes such as visualization or NMR of the reasons that microfluidics has become an essential can enable big advances. One example would be to study and adaptable tool in soft-matter research. Microfluidics shear thickening and flow in polymeric composites. Mea- is at the beginning of its development and many direc- suring local rheological properties of a material, active tions can be taken to produce new applications. We will microrheology is a new possibility. It is based on using mention here only a few of them. active probe particles that are sensitive to the local con- Current microfluidic devices are designed for benign ditions. Another route would be to use viscosity-sensitive environments. There is a need to develop devices that op- dyes to detect local properties. erate reliably in inhospitable environments ranging from At the level of creating new materials, there is the human bodies in which the devices generate immune re- need to develop rapid prototyping techniques. The low sponse, to the harsh conditions presented by chemical cost and availability of 3D printers has made it possi- synthesis applications involving organic solvents. ble to create particles at a small (but still not micro- One obvious direction is to scale down the size of the scopic) size. This capability opens up new types of ex- devices. Practical applications arise from the flow of flu- periments. However, it is a slow process to manufacture ids inside nanotubes composed of boron nitrides, zinc particles in bulk quantities. Thus, one advance would be oxide, and carbon nanotubes. Applications with great to create such particles more rapidly with much smaller societal need are water purification (e.g., desalination, re- feature sizes. Super-Resolution Additive Manufacturing moval of heavy metals, etc.) and efficient catalysis for fu- (SRAM) would be useful to many areas of science and els. However reducing the size introduces new challenges. technology and would be a unique contribution of soft- The large surface to volume ratios raise the fundamen- matter science. tal question of how to make entities move when surface When manufacturing new materials, one often con- forces dominate. At the nanoscale, pressure gradients fronts problems associated with the fundamental proper- are too small to drive flow. One well developed alterna- ties of the raw materials. For example, flow instabilities, tive is electrical-driven transport. Other means include rheology, and interfacial properties may limit the feature using acoustic, optic, and magnetic fields. The general size and resolution that can be achieved. Understanding question emerges of how soft matter interacts with these this collective behavior is one way in which soft-matter applied fields. At the nanoscale, there is a need for un- research continues to play an important role in the de- derstanding interactions at the molecular level and there velopment of the next stage of material fabrication. needs to be a bridge from the continuum description to Other areas of instrumentation development can come molecular physical chemistry. from exploiting developments in consumer electronics in- In order to make efficient use of these methods, stu- dustry. As already mentioned, GPUs, magnetic reso- dents and other researchers in soft matter must get nance imaging and X-ray tomography have entered the trained in the techniques. The neutron scattering com- toolkit of soft-matter experimentalists. One new area munity has been effective at establishing courses for those of cross fertilization could include the use of ultrasound using neutron sources. This can serve as a guide for how probes to study density and flow behavior deep within a to train the next generation of soft-matter scientists. The suspension (Han et al., 2016). community needs to include soft-matter topics in intro- Microfluidics: Microfluidics is one enabling technology ductory curricula in order to make full use of the burgeon- that deserves particular attention as a bridge between ba- ing opportunities in this area of science and technology. sic and applied soft-matter science (Squires and Quake, 2005). It relies on powerful and cheap photolithography techniques to cast devices that control the flow of mate- XV. OUTREACH rial. Microfluidics is an elegant, effective and precise way to organize small amounts of matter and to control com- It is the responsibility of investigators to communicate position, temperature and external fields in space and their results to other researchers who are working in the time. It is a flexible method for material synthesis and same area. However, this is only one aspect of scien- allows a large number of experiments to be performed tific communication. In addition to the research commu- while consuming only small amounts of raw material. It nity, there is a much broader audience that should be is particularly useful for exploring length scales between addressed that includes students, industry, government, 10µm and 100µm. Because these length scales are small, and the general public of interested individuals who are there is rapid transport of heat and matter and hence not necessarily expert in science (much less in a specific the ability to induce rapid changes in the material such area). For all of these communities, soft matter has a as mixing, dissolution and phase transitions. Microflu- special role to play because it deals not only with far- idics is also a powerful tool for understanding biology. It reaching scientific ideas but also with topics that are ac- allows the design of ex-vivo systems to test hypotheses cessible. 23

Soft-matter is a good place to introduce students to research. Because students can see, and sometimes even feel, the phenomenon, they can come up to speed rela- tively rapidly. Soft matter encourages students to query “why or how does this phenomenon behave the way it does?” Moreover, the research often provides tangible output – stunning visual images of structure or dynam- ics. These characteristics have been effective at attract- ing diverse students into science broadly and more specif- ically into the field of soft matter. The importance of the science – to advancing applications, to addressing broad cross-cutting scientific issues, to training students and creating a diverse workforce – should be made clear to those paying for the research. Because soft matter confronts issues found in our daily lives, the problems are approachable and not hidden in- side a “black box”. The interested layperson can get a genuine opportunity to understand a science problem, understand the way the experiment is set up, and ul- timately to understand the solution. Few other areas of science can accomplish this. In many cases, while the public can often spout some of the jargon that goes along FIG. 22 Turbulent Orb by Ned Kahn displayed in Battery with the science, it is not always clear that they have a Park City, New York in 1990. The sculpture consists a large, real understanding of what that jargon means or how it rotatable, spherical glass vessel, filled with a deep-blue fluid was arrived at and how the ideas have been or could be that reveals the currents within. From (Kahn, 1990). tested. Soft-matter research can be an antidote to such non-scientific understanding. It is often a comfort to hide behind obscurity and com- hibits devoted to the intriguing behavior of soft matter. plexity; soft matter cannot take that road! While soft- It can be displayed in curated public spaces as was the matter can be accessible to individuals with all levels of sculpture in Fig. 22. training, it is all the more essential to convey its impor- Another area in which soft-matter research can make tance. The apparent simplicity of the systems, can lead a difference is in working with science teachers in neigh- to the impression that the research is trivial, easy, or borhood schools. Here there is the opportunity to design unimportant. One challenge is therefore to convey, cor- kits for teachers and to bring teachers into the laboratory rectly and without undue hype, the importance of the for summer immersion in research. Such activities can basic science and the importance of the applications. It introduce students to exciting concepts in various areas, is hard to argue with a good technology. including extreme mechanics, robotics, and microfluidics. An effective way to accomplish this is to show the use- Likewise, bringing visitors on tours of research laborato- fulness of the research to industrial colleagues by organiz- ries has been a low-key but very effective way to com- ing short courses on different topics. What is particularly municate. Grade-school students lucky enough to have needed is to create stronger ties with different industrial had such an opportunity have sometimes advanced on to partners. While these exist between a few individual re- research internships early in their education. search groups and industrial affiliates, it is less true for More recent opportunities for outreach are enabled by a majority of soft-matter groups. It would be one im- social-media outlets and blogs on the internet. Fascinat- portant step forward to create a stronger bond between ing websites devoted to soft-matter phenomena, for ex- these two vital areas. ample, the site “fuckyeahfluiddynamics”, often incorpo- There are a number of specific ways that soft matter rate well-researched topics into interesting videos. There can achieve the goal of explaining science to a general are also special opportunities to form public networks: audience. Some of these are a continuation of what has one idea would be to create “big-data” sets that could be been so successful in the past. It should continue to ex- accessed by the general public for “crowd analysis”. ploit the visual nature of the research with captivating imagery and video. Because soft-matter phenomena ap- pear on all scales, many phenomena can be displayed XVI. CONCLUDING REMARKS with hands-on demonstrations. These can be in the form of museum exhibits where soft matter has long had a Soft matter addresses issues that are broadly relevant presence. Many science museums now have excellent ex- to many areas of science. It is, for example, arguable that 24 understanding far-from-equilibrium behavior – a topic Session chair: Granular and foams central to soft matter science – is the deepest and most Seth Fraden (Brandeis University) important question facing physics today. Soft matter has Session chair: Microfluidics had an exciting history of discovery and is also a fertile Dan Goldman (Georgia Inst. of Technology) ground for future research breakthroughs. New tools and David Grier (New York University) new ideas continue to emerge from the laboratories de- Anette (Peko) Hosoi (Massachusetts Inst. of Technology) voted to its study. This report has touched on only a few Session chair: Fluids of the directions where this research may lead. Because William Irvine (University of Chicago) it is broadly relevant to so many disciplines, to the way Heinrich Jaeger (University of Chicago) scientists can interact with the public, and to the train- Session chair: Outreach ing of a diverse next generation of students, soft-matter Ilona Kretzschmar (City College of New York) experiment occupies a very special position in the cur- Session chair: Instrumentation rent scientific endeavor. The outlook is indeed bright for Dan Lathrop (University of Maryland) continued productive and exciting activity. Wolfgang Losert (University of Maryland) Mahesh Mahanthappa (University of Minnesota) Tom Mason (University of California, Los Angeles) ACKNOWLEDGMENTS Session chair: Packing Gareth McKinley (Massachusetts Inst. of Technology) The author thanks Paul Sokol of the NSF for the op- Sidney Nagel (University of Chicago) portunity and encouragement to organize and write up Teri Odom (Northwestern University) the results of the workshop. All the workshop attendees Corey O’Hern (Yale University) (listed in the Appendix) provided their wisdom about the Session chair: Simulations and big data topics covered. I am particularly grateful to the Scien- Chinedum Osuji (Yale University) tific Advisory Committee and to the session chairs for Mark Robbins (Johns Hopkins University) organizing and outlining the discussions. The workshop, Tom Russell (University of Massachusetts Amherst) on which this report is based, was supported by the NSF Session chair: Polymers and rheology DMR-1560936 and NSF MRSEC DMR-1420709. Todd Squires (University of California, Santa Barbara) Kathleen Stebe (University of Pennsylvania) Session chair: Pattern formation Appendix A: Participants in the workshop held in Arlington Jeff Urbach (Georgetown University) VA on January 30-31, 2016 on which this report is based Eric Weeks (Emory University) 1. Scientific Advisory Committee Session chair: Glasses and jamming David Weitz (Harvard University) Paul Chaikin (New York University) Session chair: Self assembly Douglas Durian (University of Pennsylvania) Arjun Yodh (University of Pennsylvania) Sharon Glotzer (University of Michigan) Anette (Peko) Hosoi (Massachusetts Inst. of Technology) REFERENCES Heinrich Jaeger (University of Chicago) Howard Stone (Princeton University) Amend, J. R., E. Brown, N. Rodenberg, H. M. Jaeger, and H. Lipson (2012), IEEE Transactions on Robotics 28 (2), 341. 2. Attendees Andelman, D., and R. E. Rosensweig (2008), The Journal of Physical Chemistry B 113 (12), 3785. Robert Behringer (Duke University) Anderson, A., B. Golding, J. Graebner, S. Hunklinger, Dan Blair (Georgetown University) J. J¨ackle, W. Phillips, R. Pohl, M. Schickfus, D. Weaire, and W. A. Phillips (2012), Amorphous solids: low- Paul Chaikin (New York University) temperature properties, Vol. 24 (Springer Science & Busi- Session chair: Colloids ness Media). Xiang Cheng (University of Minnesota) Augugliaro, F., S. Lupashin, M. Hamer, C. Male, M. Hehn, Noel Clark (University of Colorado Boulder) M. W. Mueller, J. S. Willmann, F. Gramazio, M. Kohler, Session chair: Liquid crystals and R. D’Andrea (2014), IEEE Control Systems 34 (4), 46. Itai Cohen (Cornell University) Belousov, B. P. (1959), Compilation of Abstracts on Radia- Session chair: Mechanical meta-materials tion Medicine 147 (145), 1. Ben-Jacob, E., R. Godbey, N. D. Goldenfeld, J. Koplik, Karen Daniels (North Carolina State University) H. Levine, T. Mueller, and L. Sander (1985), Physical Zvonimir Dogic (Brandeis University) Review Letters 55 (12), 1315. Session chair: Active matter and biomaterials Berthier, L., and G. Biroli (2011), Reviews of Modern Physics Doug Durian (University of Pennsylvania) 83 (2), 587. 25

Berthier, L., G. Biroli, J.-P. Bouchaud, L. Cipelletti, and Applications 313 (1), 1. W. van Saarloos (2011), Dynamical heterogeneities in Geng, J., E. Longhi, R. Behringer, and D. Howell (2001), glasses, colloids, and granular media, Vol. 150 (OUP Ox- Physical Review E 64 (6), 060301. ford). Glotzer, S. C., and M. J. Solomon (2007), Nature Materials Berthier, L., P. Charbonneau, Y. Jin, G. Parisi, B. Seoane, 6 (8), 557. and F. Zamponi (2016), Proceedings of the National Goodrich, C. P., A. J. Liu, and S. R. Nagel (2015), Physical Academy of Sciences 113 (30), 8397. Review Letters 114 (22), 225501. Bi, D., J. Zhang, B. Chakraborty, and R. Behringer (2011), Greaves, G. N., A. Greer, R. Lakes, and T. Rouxel (2011), Nature 480 (7377), 355. Nature Materials 10 (11), 823. Bischofberger, I., R. Ramachandran, and S. R. Nagel (2014), Grier, D. G. (2003), Nature 424 (6950), 810. Nature Communications 5, 5265. Guo, S., M. Gao, X. Xiong, Y. J. Wang, X. Wang, P. Sheng, Bouchaud, J.-P., P. Claudin, D. Levine, and M. Otto (2001), and P. Tong (2013), Physical Review Letters 111 (2), The European Physical Journal E 4 (4), 451. 026101. Brown, E., and H. M. Jaeger (2012), Journal of Rheology Halsey, T., and A. Mehta (2003), Challenges in granular (1978-present) 56 (4), 875. physics (World Scientific). Bruji´c,J., C. Song, P. Wang, C. Briscoe, G. Marty, and H. A. Han, E., I. R. Peters, and H. M. Jaeger (2016), arXiv preprint Makse (2007), Physical Review Letters 98 (24), 248001. arXiv:1604.00380. Callies, M., and D. Qu´er´e(2005), Soft Matter 1 (1), 55. Henkes, S., Y. Fily, and M. C. Marchetti (2011), Physical Caswell, T. A., Z. Zhang, M. L. Gardel, and S. R. Nagel Review E 84 (4), 040301. (2013), Physical Review E 87 (1), 012303. Hosoi, A., and J. W. Bush (2001), Journal of Fluid Mechanics Cates, M. E., and M. R. Evans (2000), Soft and fragile 442, 217. matter: nonequilibrium dynamics, metastability and flow Hu, H., M. Gopinadhan, and C. O. Osuji (2014), Soft Matter (PBK) (CRC Press). 10 (22), 3867. Cavagna, A., A. Cimarelli, I. Giardina, G. Parisi, R. Santa- Hunter, G. L., and E. R. Weeks (2012), Reports on Progress gati, F. Stefanini, and M. Viale (2010), Proceedings of the in Physics 75 (6), 066501. National Academy of Sciences 107 (26), 11865. Jaeger, H. M., S. R. Nagel, and R. P. Behringer (1996), Cavagna, A., and I. Giardina (2014), Annu. Rev. Condens. Reviews of Modern Physics 68 (4), 1259. Matter Phys. 5 (1), 183. Josserand, C., A. V. Tkachenko, D. M. Mueth, and H. M. Chaikin, P. M., and T. C. Lubensky (2000), Principles of Jaeger (2000), Physical Review Letters 85 (17), 3632. condensed matter physics, Vol. 1 (Cambridge Univ Press). Kahn, N. (1990), “Turbulent orb,” http://nedkahn.com/ Chen, B. G.-g., B. Liu, A. A. Evans, J. Paulose, I. Cohen, portfolio/turbulent-orb/. V. Vitelli, and C. Santangelo (2016), Physical Review Let- Kane, C., and T. Lubensky (2014), Nature Physics 10 (1), ters 116 (13), 135501. 39. Chen, K., W. G. Ellenbroek, Z. Zhang, D. T. Chen, P. J. Kaz, D. M., R. McGorty, M. Mani, M. P. Brenner, and V. N. Yunker, S. Henkes, C. Brito, O. Dauchot, W. Van Saarloos, Manoharan (2012), Nature Materials 11 (2), 138. A. J. Liu, et al. (2010), Physical Review Letters 105 (2), Kleckner, D., and W. T. Irvine (2013), Nature Physics 9 (4), 025501. 253. Chen, K., M. Manning, P. J. Yunker, W. G. Ellenbroek, Kleman, M., and O. D. Lavrentovich (2003), Soft Matter Z. Zhang, A. J. Liu, and A. G. Yodh (2011), Physical Physics: An Introduction (Springer, New York). Review Letters 107 (10), 108301. Leheny, R., and S. Nagel (1997), EPL (Europhysics Letters) Cheng, J. Y., C. A. Ross, H. I. Smith, and E. L. Thomas 39 (4), 447. (2006), Advanced Materials 18 (19), 2505. Li, Y., E. Perlman, M. Wan, Y. Yang, C. Meneveau, R. Burns, Cheng, X., J. H. McCoy, J. N. Israelachvili, and I. Cohen S. Chen, A. Szalay, and G. Eyink (2008), Journal of Tur- (2011), Science 333 (6047), 1276. bulence 9, N31. Chuang, I., B. Yurke, R. Durrer, and N. Turok (1991), Sci- Liu, A. J., and S. R. Nagel (2010), Annual Reviews of Con- ence 251 (4999), 1336. densed Matter Physics 1 (1), 347. Debenedetti, P. G., and F. H. Stillinger (2001), Nature Majmudar, T. S., and R. P. Behringer (2005), Nature 410 (6825), 259. 435 (7045), 1079. Donev, A., F. H. Stillinger, P. Chaikin, and S. Torquato Marchetti, M., J. Joanny, S. Ramaswamy, T. Liverpool, (2004), Physical Review Letters 92 (25), 255506. J. Prost, M. Rao, and R. A. Simha (2013), Reviews of Driscoll, M. M., B. G.-g. Chen, T. H. Beuman, S. Ulrich, S. R. Modern Physics 85 (3), 1143. Nagel, and V. Vitelli (2016), Proceedings of the National Maret, G., and P. Wolf (1987), Zeitschrift f¨urPhysik B Con- Academy of Sciences 113 (39), 10813. densed Matter 65 (4), 409. Duran, J. (2012), Sands, powders, and grains: an introduction Miskin, M. Z., and H. M. Jaeger (2013), Nature Materials to the physics of granular materials (Springer Science & 12 (4), 326. Business Media). Mullin, T., S. Deschanel, K. Bertoldi, and M. Boyce (2007), Ediger, M. D., C. Angell, and S. R. Nagel (1996), The Journal Physical Review Letters 99 (8), 084301. of Physical Chemistry 100 (31), 13200. Nakagawa, M., S. Altobelli, A. Caprihan, E. Fukushima, and Eggers, J. (1997), Reviews of Modern Physics 69 (3), 865. E.-K. Jeong (1993), Experiments in fluids 16 (1), 54. Florijn, B., C. Coulais, and M. van Hecke (2014), Physical Narayan, V., S. Ramaswamy, and N. Menon (2007), Science Review Letters 113 (17), 175503. 317 (5834), 105. Frantzdale, B., S. J. Plimpton, and M. S. Shephard (2010), Nash, L. M., D. Kleckner, A. Read, V. Vitelli, A. M. Turner, Engineering with Computers 26 (2), 205. and W. T. Irvine (2015), Proceedings of the National Frenkel, D. (2002), Physica A: Statistical Mechanics and its Academy of Sciences 112 (47), 14495. 26

O’Hern, C. S., and M. D. Shattuck (2013), Nature Materials Wilen, and E. R. Dufresne (2013), Physical Review Letters 12 (4), 287. 110 (6), 066103. Onoda, G. Y., and E. G. Liniger (1990), Physical Review Tkalec, U., M. Ravnik, S. Copar,ˇ S. Zumer,ˇ and I. Muˇseviˇc Letters 64 (22), 2727. (2011), Science 333 (6038), 62. Orellana, C. S., and H. M. Jaeger (2013), Soft Matter 9 (35), Torquato, S. (2016), Physical Review E 94 (2), 022122. 8519. Torquato, S., and F. H. Stillinger (2010), Reviews of Modern Ouyang, Q., and H. L. Swinney (1991), Nature 352, 610. Physics 82 (3), 2633. P¨ahtz, T., H. J. Herrmann, and T. Shinbrot (2010), Nature Turing, A. M. (1952), Philosophical Transactions of the Royal Physics 6 (5), 364. Society of London B: Biological Sciences 237 (641), 37. Palacci, J., S. Sacanna, A. P. Steinberg, D. J. Pine, and P. M. Umbanhowar, P. B., F. Melo, H. L. Swinney, et al. (1996), Chaikin (2013), Science 339 (6122), 936. Nature 382 (6594), 793. Paulose, J., B. G.-g. Chen, and V. Vitelli (2015), Nature Vigil, R. D., Q. Ouyang, and H. L. Swinney (1992), Physica Physics 11, 153156. A: Statistical Mechanics and its Applications 188 (1), 17. Paulsen, J. D., V. D´emery, C. D. Santangelo, T. P. Russell, Wang, Y., Y. Wang, X. Zheng, E. Ducrot, J. S. Yodh, B. Davidovitch, and N. Menon (2015), Nature Materials. M. Weck, and D. J. Pine (2015), Nature Communications Pine, D., D. Weitz, P. Chaikin, and E. Herbolzheimer (1988), 6, 7253. Physical Review Letters 60 (12), 1134. Warner, M., and E. M. Terentjev (2003), Liquid crystal elas- Reif, J. H., T. H. LaBean, and N. C. Seeman (2001), in tomers, Vol. 120 (OUP Oxford). Robotics and Automation, 2001. Proceedings 2001 ICRA. Weeks, E. R., J. C. Crocker, A. C. Levitt, A. Schofield, and IEEE International Conference on, Vol. 1 (IEEE) pp. 966– D. A. Weitz (2000), Science 287 (5453), 627. 971. Wei, Y., C. M. Cejas, R. Barrois, R. Dreyfus, and D. J. Ribeiro, A. A., and V. Ortiz (2016), Chemical Reviews Durian (2014), Physical Review Applied 2 (4), 044004. 116 (11), 6488. Weitz, D., and M. Oliveria (1984), Physical Review Letters Rocks, J. W., N. Pashine, I. Bischofberger, C. P. Goodrich, 52 (16), 1433. A. J. Liu, and S. R. Nagel (2016), arXiv preprint Wereley, N. (2013), Magnetorheology: advances and applica- arXiv:1607.08562. tions, Vol. 6 (Royal Society of Chemistry). Rosensweig, R. E., M. Zahn, and R. Shumovich (1983), Jour- Witten, T., and L. M. Sander (1981), Physical Review Letters nal of Magnetism and Magnetic Materials 39 (1), 127. 47 (19), 1400. Roth, L. K., and H. M. Jaeger (2016), Soft Matter 12 (4), Witten, T. A., and P. A. Pincus (2004), Structured fluids: 1107. polymers, colloids, surfactants (Oxford University Press on Saffman, P. G., and G. Taylor (1958), in Proceedings of the Demand). Royal Society of London A: Mathematical, Physical and Wittmer, J., P. Claudin, M. Cates, J. Bouchaud, et al. (1996), Engineering Sciences, Vol. 245 (The Royal Society) pp. Nature 382 (6589), 336. 312–329. Yokoyama, H. (2012), “Liquid crystal institute, kent state Sanchez, T., D. T. Chen, S. J. DeCamp, M. Heymann, and university annual report,” http://www.lcinet.kent.edu/ Z. Dogic (2012), Nature 491 (7424), 431. facts/annual_reports/. Schoenholz, S. S., E. D. Cubuk, D. M. Sussman, E. Kaxiras, Zhabotinsky, A. M. (1964), in Proc. Acad. Sci. USSR, Vol. and A. J. Liu (2016), Nature Physics 12, 469471. 157, pp. 392–395. Seeman, N. C. (2003), Nature 421, 427. Zhang, J., S. Childress, A. Libchaber, and M. Shelley (2000), Sharon, E., B. Roman, and H. L. Swinney (2007), Physical Nature 408 (6814), 835. Review E 75 (4), 046211. Zhang, Z., N. Xu, D. T. Chen, P. Yunker, A. M. Alsayed, Shi, X., M. P. Brenner, and S. R. Nagel (1994), Science 265, K. B. Aptowicz, P. Habdas, A. J. Liu, S. R. Nagel, and 219. A. G. Yodh (2009), Nature 459 (7244), 230. Silbert, L. E., A. J. Liu, and S. R. Nagel (2009), Physical Zhao, R., Q. Zhang, H. Tjugito, and X. Cheng (2015), Pro- Review E 79 (2), 021308. ceedings of the National Academy of Sciences 112 (2), 342. Silverberg, J. L., J.-H. Na, A. A. Evans, B. Liu, T. C. Hull, Zheng, Z., F. Wang, Y. Han, et al. (2011), Physical Review C. D. Santangelo, R. J. Lang, R. C. Hayward, and I. Cohen Letters 107 (6), 065702. (2015), Nature Materials 14 (4), 389. Zocchi, G., E. Moses, and A. Libchaber (1990), Physica A: Squires, T. M., and S. R. Quake (2005), Reviews of Modern Statistical Mechanics and its Applications 166 (3), 387. Physics 77 (3), 977. Zou, L.-N., X. Cheng, M. L. Rivers, H. M. Jaeger, and S. R. Style, R. W., R. Boltyanskiy, Y. Che, J. Wettlaufer, L. A. Nagel (2009), Science 326 (5951), 408.