Skip to main content
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian
  • Log in
  • My Cart

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses

New Research In

Physical Sciences

Featured Portals

  • Physics
  • Chemistry
  • Sustainability Science

Articles by Topic

  • Applied Mathematics
  • Applied Physical Sciences
  • Astronomy
  • Computer Sciences
  • Earth, Atmospheric, and Planetary Sciences
  • Engineering
  • Environmental Sciences
  • Mathematics
  • Statistics

Social Sciences

Featured Portals

  • Anthropology
  • Sustainability Science

Articles by Topic

  • Economic Sciences
  • Environmental Sciences
  • Political Sciences
  • Psychological and Cognitive Sciences
  • Social Sciences

Biological Sciences

Featured Portals

  • Sustainability Science

Articles by Topic

  • Agricultural Sciences
  • Anthropology
  • Applied Biological Sciences
  • Biochemistry
  • Biophysics and Computational Biology
  • Cell Biology
  • Developmental Biology
  • Ecology
  • Environmental Sciences
  • Evolution
  • Genetics
  • Immunology and Inflammation
  • Medical Sciences
  • Microbiology
  • Neuroscience
  • Pharmacology
  • Physiology
  • Plant Biology
  • Population Biology
  • Psychological and Cognitive Sciences
  • Sustainability Science
  • Systems Biology
Research Article

Structural hierarchy confers error tolerance in biological materials

Jonathan A. Michel and View ORCID ProfilePeter J. Yunker
PNAS February 19, 2019 116 (8) 2875-2880; first published February 5, 2019; https://doi.org/10.1073/pnas.1813801116
Jonathan A. Michel
aSchool of Physics, Georgia Institute of Technology, Atlanta GA 30332
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Peter J. Yunker
aSchool of Physics, Georgia Institute of Technology, Atlanta GA 30332
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Peter J. Yunker
  • For correspondence: peter.yunker@gatech.edu
  1. Edited by Andrea J. Liu, University of Pennsylvania, Philadelphia, PA, and approved January 3, 2019 (received for review August 15, 2018)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Significance

Structural hierarchy is ubiquitous in nature and an emerging trend in engineered materials. Despite their many virtues, hierarchical materials also add complexity and dramatically increase the opportunities for random errors in assembly. Nonetheless, highly hierarchical tissues have evolved many times in diverse lineages; this prevalence suggests a common source of mechanical robustness. In this work, we introduce a tractable, model hierarchical lattice with controllable attributes on each length scale. We find, contrary to intuition, that adding additional levels of structure actually reduces the relative variation in mechanical properties, despite an increase in assembly errors. This finding informs our understanding of the emergence of hierarchically structured biological materials, while also offering a practical heuristic in materials design.

Abstract

Structural hierarchy, in which materials possess distinct features on multiple length scales, is ubiquitous in nature. Diverse biological materials, such as bone, cellulose, and muscle, have as many as 10 hierarchical levels. Structural hierarchy confers many mechanical advantages, including improved toughness and economy of material. However, it also presents a problem: Each hierarchical level adds a new source of assembly errors and substantially increases the information required for proper assembly. This seems to conflict with the prevalence of naturally occurring hierarchical structures, suggesting that a common mechanical source of hierarchical robustness may exist. However, our ability to identify such a unifying phenomenon is limited by the lack of a general mechanical framework for structures exhibiting organization on disparate length scales. Here, we use simulations to substantiate a generalized model for the tensile stiffness of hierarchical filamentous networks with a nested, dilute triangular lattice structure. Following seminal work by Maxwell and others on criteria for stiff frames, we extend the concept of connectivity in network mechanics and find a similar dependence of material stiffness upon each hierarchical level. Using this model, we find that stiffness becomes less sensitive to errors in assembly with additional levels of hierarchy; although surprising, we show that this result is analytically predictable from first principles and thus potentially model independent. More broadly, this work helps account for the success of hierarchical, filamentous materials in biology and materials design and offers a heuristic for ensuring that desired material properties are achieved within the required tolerance.

  • network mechanics
  • evolution of biomaterials
  • structural hierarchy
  • biophysics
  • soft matter

Living systems organize across many distinct levels, spanning from molecular to macroscopic scales. Such hierarchical arrangements endow organisms with many beneficial material properties; they may have high strength-to-weight ratios, exhibit strain stiffening, or be robust against fracture (1⇓⇓⇓–5). A seeming drawback of this approach, however, is the enormous amount of information needed to specify the structure of highly hierarchical tissues and the increased number of opportunities for stochastic errors. Even for a self-assembled material, each hierarchical level likely increases the number of local minima in the free energy landscape, increasing the opportunity for kinetic errors in assembly (6, 7). While one may reasonably fear that this cascade of errors will undermine the reliable realization of self-assembled hierarchical materials, structural hierarchy is used effectively by organisms belonging to many diverse evolutionary lineages (2, 8, 9). Such widespread success suggests the presence of an underlying mechanism responsible for this emergent robustness. However, the number of elements necessary to describe a hierarchical structure grows geometrically with the number of hierarchical levels; thus, a 10-level structure is currently computationally inaccessible. While identification of the underlying principles responsible for hierarchical robustness would greatly aid in explaining the ubiquity of natural hierarchical structures, this objective first requires developing a mechanistic understanding of how each scale contributes to a material’s overall properties.

To gain a foothold in the study of hierarchical materials mechanics, we focus on a highly tractable model system: a triangular lattice of nodes connected by harmonic springs. Frames made of slender, elastic beams have long been of interest in technical mechanics (10⇓⇓⇓⇓⇓⇓⇓–18) and the physics of living tissue (1⇓–3, 9, 19⇓⇓⇓–23), and recent work has demonstrated that fibers can generically emerge from diverse building blocks (24). Further, the mechanics of elastic networks are easily interpretable through the Maxwell counting heuristic; briefly, to constrain every degree of freedom in the network, there must be 2d bonds per node, where d is the dimensionality of the system (12, 25). While much work has been done to characterize elastic networks constructed with a single important length scale (10, 13⇓–15, 26, 27), we lack a general characterization of elastic networks constructed with multiple disparate length scales; a priori, it is unclear how Maxwell counting applies to hierarchical structures. Are there distinct degrees of freedom associated with “large” nodes, just as there are with “small” nodes? How do constraints on large and small length scales compare? Identification of underlying mechanisms that make hierarchical structures robust first requires developing a comprehensible hierarchical model.

Here, we introduce a model system with a nested, dilute triangular lattice structure, in which distinct network connectivities can be defined on multiple scales. We examine the dependence of tensile stiffness on each of these connectivities and capture this relationship with a simple model. Using this model, we then assess the resilience of a hierarchical material’s mechanical properties in the presence of random errors in assembly.

Geometrical Characteristics of the Model System

We consider an extension of the well-studied dilute, triangular lattice in 2D. Nodes arranged in a triangular Bravais lattice are connected to nearest neighbors, and bonds are then removed at random such that some fraction, referred to as the bond portion, p, remains. The infinite triangular lattice has a connectivity of six bonds per node when p=1, while in 2D Maxwell counting dictates a minimum connectivity of four bonds per node; thus, the infinite, dilute triangular lattice should lose stiffness when p falls below 23. Lattices of finite size would require a slightly higher bond portion, due to the presence of underconstrained nodes at the boundary. This prediction has been thoroughly confirmed for ball-and-spring networks, via simulation and mean-field theoretic approaches (10, 13⇓–15).

We create hierarchical triangular lattices through an iterative process, in which the bonds of the lattice at one length scale are in turn crafted from smaller-scale triangular lattices. In principle, this process can be iterated ad infinitum; in practice, if the number of large bonds is held constant, the total number of nodes grows geometrically with the number of hierarchical levels. This places a practical limit on the number of levels that can be considered in simulations. We have numerically constructed and simulated lattices with one, two, and three levels of structural hierarchy (Fig. 1); bond portion is independently set on each hierarchical level.

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

(A) A dilute triangular lattice with one level of structure. Missing bonds are indicated with dashed lines. (B) A two-level triangular lattice, with each bond replaced by a smaller-scale, dilute triangular lattice. (C) An extension to three levels of structural hierarchy.

Hierarchical Model of Stiffness

We propose to model the stiffness of our networks by generalizing the scaling law proposed by Garboczi and Thorpe (10), for a single-scale, diluted triangular lattice. Garboczi and Thorpe (10) found that, for the ball-and-spring case, components of the elastic constant tensor should have the formK=kp−pc1−pc,p≥pc0,p<pc,[1]where pc is the minimum bond portion necessary for marginal stiffness, and k is the value of the modulus when the network is fully connected. We propose to describe large-scale bonds using an effective stiffness with the form of Eq. 1 and introduce plarge and psmall, the portion of bonds retained on the large and small scales. Because of the finite width of large-scale bonds, we do not assume that the minimum small- or large-scale bond portions needed for marginal stiffness are 23. The stiffness of a large-scale bond will be inherited from its small-scale structure, such that the overall stiffness scales asK=kplarge−pc,largepsmall−pc,small1−pc,large1−pc,small,[2]where K is tensile stiffness and k is the stiffness for a network fully connected on all scales. Now consider a general network with N distinct length scales. The stiffness ki of bonds at scale i will be inherited from the structure at scale i−1, so thatki∝pi−1−pc,i−1.[3]The overall stiffness for some general number N levels of structural hierarchy will then beK=k∏i=1Npi−pc,i1−pc,i.[4]

Simulation Procedure

Networks were simulated in 2D, with ball-and-spring interactions between pairs of connected nodes. To measure stiffness, nodes along the tops of networks were uniformly displaced along the vertical direction, after which the y coordinates of the top and bottom nodes were fixed. Next, the x coordinates of top and bottom nodes, as well as both coordinates of all other points, were relaxed using the fast inertial relaxation engine algorithm (28). Each bond is modeled as a fiber that resists stretching with a 1D stretching modulus, μ. Let each bond be a parametric curve r(s), where s ranges from 0 to the length of the bond, and let u→(s) be a field describing the displacement of a point on the bond in the strained state. The energy for a bond of length l is described by the functionalU=∫0lμ2du→ds2ds.[5]A uniform, 1D stretching modulus of unity was assigned to each bond, and nodes were relaxed until the rms residual force in the network was less than 1×10−10 in units of stretching modulus in the case of the one- and two-level networks and 1×10−9 in units of stretching modulus for three-level networks. We extract the tensile stiffness by fitting plots of elastic energy vs. strain to parabolas.

Simulation Results

We simulated one-, two-, and three-level lattices with a wide range of bond portions and measured their stiffness (Fig. 2). Large-scale bonds are 10 small-scale bonds long and contain three rows of small-scale bonds. Before comparing to our model, we must identify each critical bond portion for one-, two-, and three-level lattices.

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

(A) Simulated stiffness plotted vs. bond portion for one length scale. (B) Heat map of simulated stiffness as a function of small and large bond portion for a two-level network. Stiffness is plotted in simulation units, as indicated in the key. (C) Heat map of simulated stiffness for a slice in bond portion space for three-level networks with full connectivity on the largest scale. Stiffness is plotted in simulation units, as indicated in the key. (D) Heat map of simulated stiffness for a slice in bond portion space for three-level networks with full connectivity on the smallest scale. Stiffness is plotted in simulation units, as indicated in the key. (E) Stiffness normalized by maximum stiffness for networks with one, two, and three levels of structure.

For the one-level lattice, we recover the expected linear relationship between stiffness and bond portion (Fig. 2A). We find the critical bond portion to be 0.670 (SI Appendix).

For the two-level lattice, we find that stiffness increases as either the small or the large bond portion is increased (Fig. 2B). To determine the threshold bond portion for the large (small) scale, we find those points in bond portion space for which the small (large) bond portion is equal to 1. We then fit a plot of large (small) bond portion vs. stiffness and fit these data to an equation of the form Eq. 1. We find that the critical bond portions are 0.83 and 0.58 for the small and large scales, respectively (SI Appendix). Initially, it may be surprising to find that the critical bond portion on the large scale is less than 0.67. However, as large-scale bonds are endowed with a finer-scale structure, they acquire an effective bending stiffness, rather than being governed strictly by harmonic, central force interactions. Notably, networks with bonds possessing bending stiffness have been demonstrated to be rigid even below the classic isostatic point (14, 15). Thus, it is crucial that pc be directly measured and not assumed from Maxwell counting.

For the three-level lattice, we find that stiffness increases as any bond portion is increased (Fig. 2 C and D and SI Appendix). Following an approach similar to that used for the two-level lattice, we find the critical bond portions are 0.83, 0.72, and 0.62 for the small, medium, and large scales, respectively.

To test our model, we compare the stiffness obtained from simulations to the stiffness predicted by our model. In Fig. 2E, we plot stiffness values from our model against stiffness values from simulations for one, two, and three levels. For ease of comparison, stiffness values for an N-level network are normalized to the maximum attainable stiffness for an N-level network. Each point in the plot compares a simulation result with its model prediction counterpart. There are no free parameters in our model; we normalize stiffnesses by the maximum values for one-, two-, or three-level lattices and use the identified critical bond portions. We find remarkable agreement; linear fits between simulated and predicted quantities have r2 values and slopes, respectively, of 0.983 and 1.003 for one level, 0.989 and 1.01 for two levels, 0.988 and 0.98 for three levels, and 0.993 and 1 for all data combined.

Experimental Prospects

While the attempt to manipulate actual tissues as we have manipulated our networks may prove elusive, we envision, as an accessible proof-of-principle experiment, cutting networks described in this work from elastic sheets. Such sheets would ideally have an out-of-plane thickness that is large in comparison with the in-plane width of bonds, to frustrate buckling, as described in ref. 29. Existing work has also established tapering of bonds at junctions as a useful technique for controlling the ratio of bending to stretching stiffness (30).

Model Applicability

We next performed a series of tests and analyses to determine how broadly our model can be applied. As a first step, we varied the size and aspect ratio of large-scale bonds. Up to this point we have shown results for networks in which large-scale bonds are 10 small-scale bonds long and contain three rows of small-scale bonds; however, we have also considered two-level networks whose large-scale bonds are 20 small-scale bonds long and contain five rows of small-scale bonds. In both cases, we find similarly strong agreement. More details are provided in SI Appendix.

While these two systems captured the same phenomena, we cannot test all possible configurations and must instead carefully consider the circumstances under which our model is applicable. Large-scale bonds can be too wide or too narrow. If the width of large-scale bonds is too great in comparison with their length, then the resulting structure will effectively be a single-scale lattice with a regular array of small perforations. If large-scale bonds contain two or fewer rows of small-scale bonds, the coordination number of the network with all bonds present will be four, making the network susceptible to immediate loss of stiffness upon dilution of bonds. Moreover, stretching of networks will lead to transverse contraction, producing a lateral load on horizontal bonds. If large-scale bonds buckle, our picture of a primarily stretching-stabilized network will be invalid. We here propose methods to identify the range of aspect ratios for which our model should be expected to apply.

We performed tensile simulations to determine the change in the displacement field of a single-scale lattice under tension due to the presence of a hole of width w. Note that a two-scale lattice with bonds of length l and width w, in units of small-scale bond length, may be viewed as a single-scale network with a periodic array of equilateral triangular holes with side length s=l−w3. We argue that our model, with a separation of scales, is appropriate when the change to the displacement field is of the order of the applied strain times the small-scale bond length over a length equal to or greater than the separation between neighboring holes. This analysis, explained in detail in SI Appendix, yieldsl/w≳2+3/w.[6]Large-scale bonds may also be too narrow and thus buckle. Using the theory of slender rods as discussed in ref. 31, we find buckling will become a concern when the ratio of length to width is of the order of the square root of the applied strain (ε); thus, to avoid buckling, the ratio of length to width must satisfylw≲ε−1/2.[7]In our case, large-scale bonds have length 10 and width 3, in units of small-scale bond length, for an aspect ratio of 5.77. This is above the limit separating hierarchical networks from perforated sheets (3.73). Further, as strains used in simulations do not exceed 5×10−3, this aspect ratio is below the limit at which buckling may be expected to occur (14.1). This places us comfortably within the constraints on the applicability of our model discussed above.

Thus, our proposed formulation for stiffness of a hierarchical lattice is expected to apply to a wide range of networks. Crucially, we observe that it accurately describes 1-, 2-, and 3-level lattices (Fig. 2E), suggesting that the smallest length scale can always be replaced by a network of even smaller bonds, and the stiffness will remain the product of all excess bond portions. This general formulation facilitates investigation of highly hierarchical (e.g., 10-level) lattices.

Structural Error Tolerance

Now that we have obtained a general relationship for the stiffness of a hierarchical structure, we are primed to consider the possibility of random errors in assembly, a likely complication in any real assembly process. We focus in particular on how stochastic deviation from a targeted set of connectivities (on all relevant length scales) results in a deviation in network stiffness. We consider two distinct regimes. In the first case, we consider a target point near the isosurface along which stiffness vanishes. In the second one, we consider a target point in bond portion space far from both the limiting case of full connectivity on any length scale and the contour of vanishing stiffness.

We first consider target points in bond portion space corresponding to marginally stiff structures. Such points are of interest, as highly compliant materials have critical biological roles (1) and are attracting increasing attention for applications (11). Such materials typically must not be susceptible to critical transitions in their elastic moduli as a result of small fluctuations in their fine-scale structure.

We use a numerical technique to estimate the distribution of stiffness arising from random errors. First, we choose a nominal point in bond portion space, such that a certain target stiffness is achieved, and the excess bond portion is the same for each length scale. We next add Gaussian random noise to the bond portion on each length scale. The stiffness of the resulting “noisy” point is then estimated by means of an interpolated function computed from simulation data for one- and two-level lattices and a fitted model for three-level lattices (Fig. 2); for lattices with more than three levels, we use Eq. 4. This process is carried out for 50,000 trials.

Strikingly, the variance in stiffness is greatly reduced with each additional level of hierarchical structure (Fig. 3A). Further, the stiffness distribution for the single-level network exhibits a large peak at zero, which is absent in the stiffness distribution of the two- and three-level networks. Thus, despite having opportunities for errors at three separate stages of assembly, the three-level network is more reliably constructed than the one-level network and can more reliably avoid stochastically generating a floppy network.

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

(A) Stiffness probability distribution functions (PDFs) estimated from histogram data for networks with one, two, and three levels of structural hierarchy. Points in 1D, 2D, and 3D bond portion with the same minute nominal stiffness were chosen, and Gaussian random variables were added to each bond portion. Note the spike in the PDF for the one-level network at zero stiffness. In this case, we consider a nominal stiffness of 0.025 and Gaussian random noise for each bond portion with zero mean and a SD of 0.005. Other cases are addressed in SI Appendix. (B) Relative error in stiffness vs. levels of hierarchy is plotted for K¯=0.001, σ=0.0001. As additional levels of structural hierarchy are added, the relative error in the tensile stiffness decreases precipitously at first, and the effect saturates at a certain number of levels. Provided the assumptions leading to Eq. 9 hold, our analytical theory and numerical approach are in close agreement (r2≈1). Here, the product of excess bond portions is 0.001, and the noise has amplitude 0.0001 on each scale. (C) The number of all small-scale bonds in a network, divided by the area enclosed by the outer perimeter of the network, is shown vs. number of hierarchical levels for stiffness values of 0.02, 0.03, and 0.04, in units of stretching modulus over small-scale bond length. Networks were chosen to have the same bond portion on all three levels. In each case, increasing hierarchy leads to markedly lower density of small-scale bonds, attesting to the ability of hierarchy to confer both robustness and efficiency.

Next, we seek a general understanding of target points far from any boundary in bond portion space (see SI Appendix for detailed derivation). For N levels, there is a nominal excess bond portion, pe, for each level; we consider identically and independently distributed deviations from the nominal bond portions according to a normal distribution with zero mean and standard deviation (SD) σ. Referring to Eq. 4, we define the reduced stiffness:K¯=Kk∏i=1N1−pc,i.[8]With this definition we find the expected deviation in the stiffness of a network with N hierarchical levels to beΔK¯K¯≈NσK¯1/N,[9]where the approximation holds when σ≪pe. For a target K¯, this functional form predicts the optimal number of levels to beN*=⌊−2⁡lnK¯⌋,[10]where “⌊⌋” denotes the floor.

To test our analytical result, we again use the above numerical approach to calculate SD in stiffness for networks with 1–15 hierarchical levels. For lattices with more than 3 levels, we use Eq. 4 to estimate the stiffness of the resulting noisy point. We find very good agreement between the numerically generated data and our analytical prediction (r2≈1 for the case shown in Fig. 3B).

The above derivations were performed assuming identically and independently distributed random errors in bond portion; however, the results presented above do not strictly require such stringent conditions. Similar robustness against fluctuation in stiffness can occur for error rates that vary on different scales and for distributions of random errors in which errors in bond portion on different length scales are correlated (see SI Appendix for an in-depth treatment).

Interestingly, investigating networks with different error rate distributions allows us to identify a useful heuristic. Consider a two-level network with the same error rate in its large-scale bond portion as a one-level network. The two-level network will have a smaller variance in its stiffness than the one-level network as long as its small bond portion error rate is less than three times larger than its large bond portion error rate (note that this argument also works if the roles of the small and large scales are swapped; see SI Appendix for details).

These observations do not simply arise from the generic properties of random variables. Although a relation of the form Eq. 9 with N=1 would hold for a general product of random numbers, the dependence upon N of the divisor indicates a scaling of the relative fluctuation specific to a stiffness with the functional form of Eq. 4. To produce a supple solid of a particular stiffness, the portion of bonds that must be retained on each scale to achieve this stiffness increases with increasing hierarchy. This can ensure that the connectivity of the network on each scale exceeds the threshold connectivity by an amount large in comparison with the the typical fluctuations in connectivity.

We further find that hierarchical structures offer superior economy of material, in keeping with previous studies (16). To quantify this benefit, we define the number density, ns, of small-scale bonds as follows: Let ls be the length of a single small-scale bond, Ns denote the total number of small-scale bonds, and A denote the area enclosed by the perimeter of the network. Then, ns is given byns=NsA[11]and has units of ls−2. We regard small-scale bonds as the fundamental structural units of a network, such that ns provides a measure of the number of connections per unit area needed to achieve a target set of mechanical properties. To compare one-, two-, and three-level networks, we chose target stiffness values and identified the points in bond portion space with equal bond portions on each level corresponding to these stiffness values. We then estimated accompanying values of ns for networks with these bond portions by interpolation. For each stiffness chosen, there was a sharp decline in the number of bonds per unit area needed with increasing numbers of hierarchical levels. This result emphasizes that hierarchical materials’ greater error tolerance is a benefit inherent to their geometry, rather than deriving from redundancy. A nested, modular assembly process thus enables the creation of a structure that is both more reliable and more frugal.

Discussion and Conclusion

Contrary to expectation, the elastic moduli of hierarchical materials are more reliably controlled than the elastic moduli of materials with one relevant length scale. Thus, with respect to random errors, it is easier to make hierarchical structures than to make single–length-scale structures. This finding may have wide-ranging implications for evolutionary biology and materials design.

It may seem that evolving progressively larger, more complex bodies is accompanied, and impeded, by a growing assortment of mechanical challenges. To the contrary, this work provides evidence that adding hierarchical complexity can actually reduce stochastic variation in material properties. This effect decreases the need for coevolving error-correcting mechanisms, thus facilitating the evolution of new traits that are “good enough” for an organism to survive.

Upon successful assembly, structural hierarchy is known to endow materials with a host of desirable properties that are unattainable with single–length-scale structures. Our finding that hierarchy also reduces susceptibility to stochastic errors suggests that bottom–up production processes for synthetic hierarchical materials may produce finished products which perform as intended despite modest, but discernible errors.

While we have explicitly addressed 2D, lattice-based networks in this study, we anticipate that our findings will have broader applicability. Feng et al. (13) have demonstrated that an analogous model to Eq. 1 describes the dependence of the elastic moduli of a dilute face-centered cubic lattice in 3D on the portion of bonds retained. Moreover, connectivity has been established as a useful control parameter for predicting the mechanical properties of many disordered solids, such as jammed packings of hard (32, 33) and soft (34) particles and semiflexible biopolymer networks (21).

More generally, we have broadened the scope of Maxwell’s visionary means of characterizing frames to account for structural hierarchy, facilitating the understanding of cases in which it has proved difficult to relate materials’ emergent properties to their fine-scale structure. While recent computational advancements have enabled study of biological macromolecules over experimentally relevant time scales (35), comprehensive understanding at the level of an organism demands a coarse-graining procedure for which our model may offer a useful road map. A generalized counting heuristic may also offer a means of expediting feasible, yet cumbersome calculations in materials design. Hierarchical materials necessarily have many design attributes, but our accessible model may considerably narrow the search of parameter space needed to reach a goal.

Materials and Methods

Network Creation.

First, a large-scale lattice is created and diluted. Dilution begins with the shuffling of all bonds, after which a random minimum spanning tree is created using Kruskal’s algorithm. Bonds are then drawn at random from those bonds not used to create the spanning tree and added to the network until the desired bond portion is reached. Next, a small-scale lattice is overlaid such that the large-scale bond length is an integer multiple of the small-scale bond length, and the position of each large-scale node coincides with the position of a small-scale node. Each small-scale bond lying within a large-scale bond is retained, after which small-scale bonds are diluted to the desired bond portion. Small-scale bond dilution is carried out in such a way that a system-spanning contact network remains, no large-scale bond is severed, and all adjacent large-scale bonds remain connected. In each successive iteration, a smaller scale may then be introduced, with the previously smallest scale taking the role of the large scale. This process is described schematically in Fig. 1.

Finding the Critical Bond Portions.

Critical bond portions for a network with N levels are computed from simulation data for each level by choosing all points in bond portion space for which the network is fully connected on the other N−1 levels, and the stiffness is nonzero. For each level, we then fit a plot of stiffness vs. the bond portion of interest to an equation of the form Eq. 1. For three-level networks, data for stiffness vs. bond portion for three-level networks were fitted to a line of the formK(p)=a⋅p+k0[12]and the x intercept of this line was taken to be the critical bond portion. For more details, please see SI Appendix.

Acknowledgments

The authors thank Alberto Fernandez-Nieves and Zeb Rocklin for helpful comments. The authors acknowledge funding from Georgia Institute of Technology’s Soft Matter Incubator and the NSF through Grant IOS-1656549.

Footnotes

  • ↵1To whom correspondence should be addressed. Email: peter.yunker{at}gatech.edu.
  • Author contributions: J.A.M. and P.J.Y. designed research; J.A.M. performed research; J.A.M. analyzed data; and J.A.M. and P.J.Y. wrote the paper.

  • The authors declare no conflict of interest.

  • This article is a PNAS Direct Submission.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1813801116/-/DCSupplemental.

Published under the PNAS license.

References

  1. ↵
    1. Fratzl P,
    2. Weinkamer R
    (2007) Nature’s hierarchical materials. Prog Mater Sci 52:1263–1334.
    OpenUrlCrossRef
  2. ↵
    1. Meyers MA,
    2. McKittrick J,
    3. Chen P
    (2013) Structural biological materials: Critical mechanics-materials connections. Science 339:773–780.
    OpenUrlAbstract/FREE Full Text
  3. ↵
    1. Launey ME,
    2. Buehler M,
    3. Ritchie RO
    (2010) On the mechanistic origins of toughness in bone. Annu Rev Mater Res 40:25–53.
    OpenUrlCrossRef
  4. ↵
    1. Yao H,
    2. Gao H
    (2007) Multi-scale cohesive laws in hierarchical materials. Int J Solids Struct 44:8177–8193.
    OpenUrl
  5. ↵
    1. Sen D,
    2. Buehler MJ
    (2011) Structural hierarchies define toughness and defect-tolerance despite simple and mechanically inferior brittle building blocks. Sci Rep 1:35.
    OpenUrlPubMed
  6. ↵
    1. Hormoz S,
    2. Brenner MP
    (2011) Design principles for self-assembly with short-range interactions. Proc Natl Acad Sci USA 108:5193–5198.
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Zeravcic Z,
    2. Manoharan VN,
    3. Brenner MP
    (2014) Size limits of self-assembled colloidal structures made using specific interactions. Proc Natl Acad Sci USA 111:15918–15923.
    OpenUrlAbstract/FREE Full Text
  8. ↵
    1. Harris J,
    2. Böhm C,
    3. Wolf S
    (2017) Universal structure motifs in biominerals: A lesson from nature for the efficient design of bioinspired functional materials. Interface Focus 7:20160120.
    OpenUrlCrossRef
  9. ↵
    1. Dunlop JWC,
    2. Fratzl P
    (2010) Biological composites. Annu Rev Mater Res 40:1–24.
    OpenUrlCrossRef
  10. ↵
    1. Garboczi EJ,
    2. Thorpe MF
    (1985) Effective-medium theory of percolation on central force elastic networks. II. Further results. Phys Rev B 31:7276–7281.
    OpenUrl
  11. ↵
    1. Wegst U,
    2. Bai H,
    3. Saiz E,
    4. Tomsia A,
    5. Ritchie R
    (2015) Bioinspired structural materials. Nat Mater 14:23–26.
    OpenUrlCrossRefPubMed
  12. ↵
    1. Maxwell JC
    (1864) On the calculation of the equilibrium and stiffness of frames. Philos Mag 27:294–299.
    OpenUrl
  13. ↵
    1. Feng S,
    2. Thorpe MF,
    3. Garboczi E
    (1985) Effective-medium theory of percolation on central-froce elastic networks. Phys Rev B 31:276–280.
    OpenUrl
  14. ↵
    1. Das M,
    2. MacKintosh FC,
    3. Levine AJ
    (2007) Effective medium theory of semiflexible filamentous networks. Phys Rev Lett 99:038101.
    OpenUrlCrossRefPubMed
  15. ↵
    1. Mao X,
    2. Stenull O,
    3. Lubensky TC
    (2013) Effective-medium theory of a filamentous triangular lattice. Phys Rev E 87:042601.
    OpenUrl
  16. ↵
    1. Lakes R
    (1993) Materials with structural hierarchy. Nature 361:511–515.
    OpenUrlCrossRef
  17. ↵
    1. Wilhelm J,
    2. Frey E
    (2003) Elasticity of stiff polymer networks. Phys Rev Lett 91:108103.
    OpenUrlCrossRefPubMed
  18. ↵
    1. Feng J,
    2. Levine H,
    3. Mao X,
    4. Sander L
    (2016) Nonlinear elasticity of disordered fiber networks. Soft Matter 12:1419–1424.
    OpenUrl
  19. ↵
    1. Head DA,
    2. Levine AJ,
    3. MacKintosh FC
    (2003) Distinct regimes of elastic response and deformation modes of cross-linked cytoskeletal and semiflexible polymer networks. Phys Rev E 68:061907.
    OpenUrl
  20. ↵
    1. Lindstrom SB,
    2. Vader DA,
    3. Kulachenko A,
    4. Weitz DA
    (2010) Biopolymer network geometries: Characterization, regeneration, and elastic properties. Phys Rev E 82:051905.
    OpenUrl
  21. ↵
    1. Broedersz CP,
    2. MacKintosh FC
    (2014) Modeling semiflexible polymer networks. Rev Mod Phys 86:995–1036.
    OpenUrlCrossRef
  22. ↵
    1. Amuasi HE,
    2. Heussingler C,
    3. Vink RLC,
    4. Zippelius A
    (2015) Nonlinear and heterogeneous elasticity of multiply crosslinked biopolymer networks. New J Phys 17:083035.
    OpenUrl
  23. ↵
    1. Licup AJ, et al.
    (2015) Stress controls the mechanics of collagen networks. Proc Natl Acad Sci USA 112:9573–9578.
    OpenUrlAbstract/FREE Full Text
  24. ↵
    1. Lenz M,
    2. Witten TA
    (2017) Geometrical frustration yields fibre formation in self-assembly. Nat Phys 13:1100–1105.
    OpenUrl
  25. ↵
    1. Broedersz CP,
    2. Mao X,
    3. Lubensky TC,
    4. MacKintosh FC
    (2011) Criticality and isostaticity in fibre networks. Nat Phys 7:983–985.
    OpenUrlCrossRef
  26. ↵
    1. Zhang T,
    2. Schwarz JM,
    3. Das M
    (2014) Mechanics of anisotropic spring networks. Phys Rev E 90:062139.
    OpenUrl
  27. ↵
    1. Lubensky TC,
    2. Kane CL,
    3. Mao X,
    4. Souslov A,
    5. Sun K
    (2015) Phonons and elasticity in critically coordinated lattices. Rep Prog Phys 78:073901.
    OpenUrlCrossRefPubMed
  28. ↵
    1. Bitzek E,
    2. Koskinen P,
    3. Gaehler F,
    4. Moseler M,
    5. Gumbsch P
    (2006) Structural relaxation made simple. Phys Rev Lett 97:170201.
    OpenUrlCrossRefPubMed
  29. ↵
    1. Coulais C,
    2. Sabbadini A,
    3. Vink F,
    4. van Hecke M
    (2018) Multi-step self-guided pathways for shape-changing metamaterials. Nature 561:512–515.
    OpenUrl
  30. ↵
    1. Rocks JW, et al.
    (2017) Designing allostery-inspired response in mechanical networks. Proc Natl Acad Sci USA 114:2520–2525.
    OpenUrlAbstract/FREE Full Text
  31. ↵
    1. Landau LD,
    2. Lifshitz EM
    (2012) Theory of Elasticity (Pergamon Press, Oxford), 2nd Ed, pp 97–98.
  32. ↵
    1. Liu A,
    2. Nagel S
    (2010) The jamming transition and the marginally jammed solid. Annu Rev Condens Matter Phys 1:347–369.
    OpenUrlCrossRef
  33. ↵
    1. O’Hern C,
    2. Liu A,
    3. Nagel S
    (2003) Jamming at zero temperature and zero applied stress: The epitome of disorder. Phys Rev E Stat Nonlin Soft Matter Phys 68:011306.
    OpenUrlCrossRefPubMed
  34. ↵
    1. van Hecke M
    (2010) Jamming of soft particles: Geometry, mechanics, scaling and isostaticity. J Phys Condens Matter 22:033101.
    OpenUrlCrossRefPubMed
  35. ↵
    1. Phillips JC, et al.
    (2005) Scalable molecular dynamics with NAMD. J Comput Chem 26:1781–1802.
    OpenUrlCrossRefPubMed
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Structural hierarchy confers error tolerance in biological materials
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Structural hierarchy confers error tolerance in biological materials
Jonathan A. Michel, Peter J. Yunker
Proceedings of the National Academy of Sciences Feb 2019, 116 (8) 2875-2880; DOI: 10.1073/pnas.1813801116

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Structural hierarchy confers error tolerance in biological materials
Jonathan A. Michel, Peter J. Yunker
Proceedings of the National Academy of Sciences Feb 2019, 116 (8) 2875-2880; DOI: 10.1073/pnas.1813801116
Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley
Proceedings of the National Academy of Sciences: 116 (8)
Table of Contents

Submit

Sign up for Article Alerts

Article Classifications

  • Physical Sciences
  • Physics
  • Biological Sciences
  • Biophysics and Computational Biology

Jump to section

  • Article
    • Abstract
    • Geometrical Characteristics of the Model System
    • Hierarchical Model of Stiffness
    • Simulation Procedure
    • Simulation Results
    • Experimental Prospects
    • Model Applicability
    • Structural Error Tolerance
    • Discussion and Conclusion
    • Materials and Methods
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Surgeons hands during surgery
Inner Workings: Advances in infectious disease treatment promise to expand the pool of donor organs
Despite myriad challenges, clinicians see room for progress.
Image credit: Shutterstock/David Tadevosian.
Setting sun over a sun-baked dirt landscape
Core Concept: Popular integrated assessment climate policy models have key caveats
Better explicating the strengths and shortcomings of these models will help refine projections and improve transparency in the years ahead.
Image credit: Witsawat.S.
Double helix
Journal Club: Noncoding DNA shown to underlie function, cause limb malformations
Using CRISPR, researchers showed that a region some used to label “junk DNA” has a major role in a rare genetic disorder.
Image credit: Nathan Devery.
Steamboat Geyser eruption.
Eruption of Steamboat Geyser
Mara Reed and Michael Manga explore why Yellowstone's Steamboat Geyser resumed erupting in 2018.
Listen
Past PodcastsSubscribe
Multi-color molecular model
Enzymatic breakdown of PET plastic
A study demonstrates how two enzymes—MHETase and PETase—work synergistically to depolymerize the plastic pollutant PET.
Image credit: Aaron McGeehan (artist).

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Special Feature Articles – Most Recent
  • List of Issues

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Librarians
  • Press
  • Site Map
  • PNAS Updates

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490