Preserving elastic anisotropy with tessellations of granular packings
Pith reviewed 2026-05-10 14:50 UTC · model grok-4.3
The pith
Tessellated granular packings inherit the high elastic anisotropy of their voxels, reaching levels 100 times higher than in crystals.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Homogeneously tessellated granular systems can inherit the elastic response of the constituent voxel configurations with elastic anisotropy up to 100 times that of crystalline compounds over a range of pN². Bulk packings with N grains at pressure p reach maximum anisotropy near pN² ~ 1 and become isotropic at large pN², but the tessellations prevent this loss by limiting rearrangements.
What carries the argument
Tessellations of multiple connected grain-filled voxels that constrain grain rearrangements and allow the bulk response to match the anisotropy of the isolated voxel packings.
Load-bearing premise
The connections between voxels are strong enough to stop grains from rearranging freely and to prevent new boundary effects or inter-voxel interactions from erasing the inherited anisotropy.
What would settle it
A simulation or experiment on a large tessellated sample in which the measured A_G and A_C drop to the low values typical of bulk packings or crystals rather than staying close to the high values of the separate voxels.
Figures
read the original abstract
Multiscale periodic metamaterials have been designed for numerous applications, such as impact absorption, acoustic cloaking, photonic band gaps, and mechanical logic gates. This prior work has focused on optimizing mesoscale structure for desired bulk isotropic properties. In contrast, we seek to develop materials with highly anisotropic elastic properties. To quantify elastic anisotropy, we introduce two rotationally invariant, normalized quantities that characterize the anisotropic response to shear and compression, respectively, $A_G$ and $A_C$. We find that typical crystalline solids possess average elastic anisotropy $\overline{A}_G \approx 0.15$ and $\overline{A}_C \approx 0.09$. Compared to atomic crystals, jammed granular materials can attain elastic anisotropies that are several orders of magnitude larger. Since grain rearrangements reduce anisotropy in granular materials, to preserve strong elastic anisotropy, we design tessellated granular materials that consist of multiple connected grain-filled voxels, which limit rearrangements and enable highly anisotropic elastic properties. Bulk granular packings with $N$ grains prepared at pressure $p$ have maximal anisotropy for $pN^2\sim1$ and become isotropic in the large-$pN^2$ limit. We show that homogeneously tessellated granular systems can inherit the elastic response of the constituent voxel configurations with elastic anisotropy up to $100$ times that of crystalline compounds over a range of $pN^2$. We show further methods to tune the elastic anisotropy of tessellations by designing heterogeneously patterned voxel configurations and tessellations that allow large boundary deformations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript claims that homogeneously tessellated granular packings can inherit the high elastic anisotropy of their constituent grain-filled voxels, achieving values of A_G and A_C up to 100 times larger than those of crystalline solids (where averages are ~0.15 and ~0.09) over a range of pN^2. Bulk packings exhibit maximal anisotropy near pN^2 ~ 1 and isotropize at large pN^2 due to rearrangements; tessellations are proposed to limit such rearrangements while allowing tuning via heterogeneous voxel patterns and boundary deformations.
Significance. If validated, the result offers a practical route to engineering strongly anisotropic granular metamaterials without relying on atomic-scale crystals. The rotationally invariant, normalized measures A_G and A_C are a clear methodological advance for quantifying directional elastic response. The pN^2 scaling and its apparent preservation under tessellation provide concrete insight into the role of grain mobility, with potential relevance to impact-absorbing and logic-gate metamaterials.
major comments (2)
- [Results on tessellated systems] The inheritance claim is load-bearing for the 100x anisotropy result, yet the manuscript provides no quantitative comparison of contact networks or force-chain statistics between isolated voxels and their tessellated assemblies (e.g., in the results section discussing bulk response). Without this, it remains possible that inter-voxel grain contacts open additional shear or compression modes that reduce A_G and A_C, especially near pN^2 ~ 1 where maximal anisotropy is reported.
- [Abstract and Results] The abstract states that tessellations 'inherit the elastic response' over a range of pN^2, but no explicit bounds on that range, no error estimates on the anisotropy ratios, and no demonstration that boundary-induced softening is negligible appear in the main text or figures. This directly affects the central scaling claim and must be addressed with tabulated data or supplementary plots.
minor comments (2)
- [Methods] The explicit definitions and normalization procedures for A_G and A_C should be stated with equations in the main text (rather than deferred), including any assumptions about the stiffness tensor components used in their construction.
- [Figures] Figures illustrating tessellation geometries and corresponding anisotropy values would benefit from inclusion of multiple independent realizations with error bars to demonstrate reproducibility of the reported inheritance.
Simulated Author's Rebuttal
We thank the referee for their constructive and detailed review. We address each major comment below and agree that strengthening the evidence for inheritance via contact network comparisons and providing explicit bounds, errors, and boundary checks will improve the manuscript. We will incorporate these revisions.
read point-by-point responses
-
Referee: [Results on tessellated systems] The inheritance claim is load-bearing for the 100x anisotropy result, yet the manuscript provides no quantitative comparison of contact networks or force-chain statistics between isolated voxels and their tessellated assemblies (e.g., in the results section discussing bulk response). Without this, it remains possible that inter-voxel grain contacts open additional shear or compression modes that reduce A_G and A_C, especially near pN^2 ~ 1 where maximal anisotropy is reported.
Authors: We agree that direct quantitative comparisons of contact networks and force-chain statistics between isolated voxels and tessellated assemblies would strengthen the inheritance claim. The current manuscript demonstrates inheritance primarily through the close agreement of A_G and A_C values (and their pN² scaling) between voxels and bulk tessellations, as shown in the results figures. However, we did not include explicit contact statistics. In the revised manuscript, we will add a supplementary section with comparisons of coordination number distributions, force magnitude histograms, and force-chain visualizations for isolated voxels versus tessellated systems at pN² ≈ 1. These will show that inter-voxel contacts do not introduce modes that substantially reduce anisotropy. revision: yes
-
Referee: [Abstract and Results] The abstract states that tessellations 'inherit the elastic response' over a range of pN^2, but no explicit bounds on that range, no error estimates on the anisotropy ratios, and no demonstration that boundary-induced softening is negligible appear in the main text or figures. This directly affects the central scaling claim and must be addressed with tabulated data or supplementary plots.
Authors: We agree that the abstract and results section should provide explicit bounds on the pN² range, error estimates on A_G and A_C, and evidence that boundary softening is negligible. The manuscript shows the scaling behavior in plots but does not state numerical bounds or errors explicitly. In the revision, we will update the abstract and main text to specify the range (e.g., 0.1 < pN² < 10, where anisotropy is preserved within ~10% of voxel values), add error bars (standard deviations from ensemble averages) to the relevant figures, and include a supplementary plot of elastic moduli versus tessellation size to confirm boundary effects are negligible. A table summarizing anisotropy ratios with errors will also be added. revision: yes
Circularity Check
No circularity; derivation relies on independent definitions and simulations
full rationale
The paper introduces two new rotationally invariant normalized anisotropy measures A_G and A_C that are defined from the elastic tensor without reference to the tessellation construction or target values. The central result that homogeneously tessellated granular systems inherit voxel-level anisotropy is obtained from direct numerical simulations of the bulk response under varying pN^2, not from any algebraic reduction, fitted parameter, or self-citation chain. No load-bearing step equates the output to its inputs by construction, and the comparison to crystalline solids uses externally known low anisotropy values as a benchmark.
Axiom & Free-Parameter Ledger
free parameters (1)
- pN^2 =
~1
axioms (1)
- domain assumption Grain rearrangements reduce anisotropy in granular materials
invented entities (2)
-
A_G
no independent evidence
-
A_C
no independent evidence
Reference graph
Works this paper leans on
-
[1]
N. Liu, M. Becton, L. Zhang, K. Tang, and X. Wang, Mechanical anisotropy of two-dimensional metamateri- als: A computational study, Nanoscale Advances 1, 2891 (2019)
work page 2019
-
[2]
J.-H. Lee, J. P. Singer, and E. L. Thomas, Micro- /nanostructured mechanical metamaterials, Advanced Materials 24, 4782 (2012)
work page 2012
-
[3]
M. I. N. Rosa, K. Karapiperis, K. Radi, E. Pescialli, and D. M. Kochmann, Enhanced deformability through dis- tributed buckling in stiff quasicrystalline architected ma- terials, Advanced Materials 37, 2505125 (2025)
work page 2025
- [4]
-
[5]
M. Benedetti, A. Du Plessis, R. Ritchie, M. Dallago, N. Razavi, and F. Berto, Architected cellular materials: A review on their mechanical properties towards fatigue- tolerant design and fabrication, Materials Science and Engineering: R: Reports 144, 100606 (2021)
work page 2021
-
[6]
D. Mousanezhad, S. Babaee, R. Ghosh, E. Mahdi, K. Bertoldi, and A. Vaziri, Honeycomb phononic crys- tals with self-similar hierarchy, Physical Review B 92, 104304 (2015)
work page 2015
-
[7]
Y. Chen and L. Wang, Periodic co-continuous acous- tic metamaterials with overlapping locally resonant and bragg band gaps, Applied Physics Letters 105, 191907 (2014)
work page 2014
-
[8]
H. Chen and C. T. Chan, Acoustic cloaking in three di- mensions using acoustic metamaterials, Applied Physics Letters 91, 183518 (2007)
work page 2007
-
[9]
Z. Yang, H. M. Dai, N. H. Chan, G. C. Ma, and P. Sheng, Acoustic metamaterial panels for sound attenuation in the 50–1000 hz regime, Applied Physics Letters 96, 041906 (2010)
work page 2010
-
[10]
S. Yuan, C. K. Chua, and K. Zhou, 3d-printed mechani- cal metamaterials with high energy absorption, Advanced Materials Technologies 4, 1800419 (2019)
work page 2019
-
[11]
G. Di Frisco, R. Yousefi Nooraie, M. Guagliano, and S. Bagherifard, Structural design and characterization of hybrid hierarchical lattice structures based on sheet- network triply periodic minimal surface topology, Mate- rials & Design 246, 113336 (2024)
work page 2024
- [12]
- [13]
-
[14]
G. N. Greaves, A. L. Greer, R. S. Lakes, and T. Rouxel, Poisson’s ratio and modern materials, Nature Materials 10, 823 (2011)
work page 2011
- [15]
-
[16]
Z. Gao, D. Liu, and D. Tom´ anek, Two-dimensional me- chanical metamaterials with unusual poisson ratio behav- ior, Physical Review Applied 10, 064039 (2018)
work page 2018
-
[17]
J. Seto, H. S. Gupta, P. Zaslansky, H. D. Wagner, and P. Fratzl, Tough lessons from bone: Extreme mechanical anisotropy at the mesoscale, Advanced Functional Mate- rials 18, 1905 (2008)
work page 1905
-
[18]
T. M. Keaveny, E. F. Morgan, G. L. Niebur, and O. C. Yeh, Biomechanics of trabecular bone, Annual Review of Biomedical Engineering 3, 307 (2001)
work page 2001
-
[19]
W. P. Moestopo, A. J. Mateos, R. M. Fuller, J. R. Greer, and C. M. Portela, Pushing and pulling on ropes: Hier- archical woven materials, Advanced Science 7, 2001271 (2020)
work page 2020
- [20]
-
[21]
J. B. Berger, H. N. G. Wadley, and R. M. McMeek- ing, Mechanical metamaterials at the theoretical limit of isotropic elastic stiffness, Nature 543, 533 (2017)
work page 2017
-
[22]
E. T. Filipov, T. Tachi, and G. H. Paulino, Origami tubes assembled into stiff, yet reconfigurable structures and metamaterials, Proceedings of the National Academy of Sciences 112, 12321 (2015)
work page 2015
- [23]
-
[24]
Y. Song, R. M. Panas, S. Chizari, L. A. Shaw, J. A. Jack- son, J. B. Hopkins, and A. J. Pascall, Additively manu- facturable micro-mechanical logic gates, Nature Commu- nications 10, 882 (2019)
work page 2019
- [25]
- [26]
-
[27]
A. Ion, L. Wall, R. Kovacs, and P. Baudisch, Digital me- chanical metamaterials, in Proceedings of the 2017 CHI Conference on Human Factors in Computing Systems (ACM, Denver Colorado USA, 2017) pp. 977–988
work page 2017
-
[28]
C. P. Goodrich, S. Dagois-Bohy, B. P. Tighe, M. van Hecke, A. J. Liu, and S. R. Nagel, Jamming in finite systems: Stability, anisotropy, fluctuations, and scaling, Physical Review E 90, 022138 (2014)
work page 2014
-
[29]
M. De Jong, W. Chen, T. Angsten, A. Jain, R. Notes- tine, A. Gamst, M. Sluiter, C. Krishna Ande, S. Van Der Zwaag, J. J. Plata, C. Toher, S. Curtarolo, G. Ceder, K. A. Persson, and M. Asta, Charting the complete elas- tic properties of inorganic crystalline compounds, Scien- tific Data 2, 150009 (2015)
work page 2015
-
[30]
H. M. Jaeger, S. R. Nagel, and R. P. Behringer, Granular solids, liquids, and gases, Reviews of Modern Physics 68, 1259 (1996)
work page 1996
-
[31]
A. J. Liu and S. R. Nagel, The jamming transition and the marginally jammed solid, Annual Review of Con- densed Matter Physics 1, 347 (2010)
work page 2010
-
[32]
R. P. Behringer and B. Chakraborty, The physics of jamming for granular materials: A review, Reports on Progress in Physics 82, 012601 (2018)
work page 2018
-
[33]
M. Tsamados, A. Tanguy, C. Goldenberg, and J.-L. Barrat, Local elasticity map and plasticity in a model Lennard-Jones glass, Physical Review E 80, 026112 16 (2009)
work page 2009
-
[34]
M. L. Manning and A. J. Liu, Vibrational modes iden- tify soft spots in a sheared disordered packing, Physical Review Letters 107, 108302 (2011)
work page 2011
-
[35]
E. D. Cubuk, S. S. Schoenholz, J. M. Rieser, B. D. Mal- one, J. Rottler, D. J. Durian, E. Kaxiras, and A. J. Liu, Identifying structural flow defects in disordered solids us- ing machine-learning methods, Physical Review Letters 114, 108001 (2015)
work page 2015
-
[36]
M. L. Falk and J. S. Langer, Dynamics of viscoplastic deformation in amorphous solids, Physical Review E 57, 7192 (1998)
work page 1998
-
[37]
A. Zaccone and E. Scossa-Romano, Approximate analyt- ical description of the nonaffine response of amorphous solids, Physical Review B 83, 184205 (2011)
work page 2011
- [38]
- [39]
-
[40]
C. Maloney and A. Lemaˆ ıtre, Subextensive scaling in the athermal, quasistatic limit of amorphous matter in plas- tic shear flow, Physical Review Letters93, 016001 (2004)
work page 2004
-
[41]
G. I. Taylor, The mechanism of plastic deformation of crystals. part i.—theoretical, Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character 145, 362 (1934)
work page 1934
-
[42]
N. Pashine, D. Wang, J. Zhang, S. K. Patiballa, S. Wit- thaus, M. D. Shattuck, C. S. O’Hern, and R. Kramer- Bottiglio, Tessellated granular metamaterials with tun- able elastic moduli, Extreme Mechanics Letters 63, 102055 (2023)
work page 2023
- [43]
-
[44]
D. Wang, J. D. Treado, A. Boromand, B. Norwick, M. P. Murrell, M. D. Shattuck, and C. S. O’Hern, The struc- tural, vibrational, and mechanical properties of jammed packings of deformable particles in three dimensions, Soft Matter 17, 9901 (2021)
work page 2021
-
[45]
J. Gu´ enol´ e, W. G. N¨ ohring, A. Vaid, F. Houll´ e, Z. Xie, A. Prakash, and E. Bitzek, Assessment and optimization of the fast inertial relaxation engine (FIRE) for energy minimization in atomistic simulations and its implemen- tation in lammps, Computational Materials Science 175, 109584 (2020)
work page 2020
-
[46]
K. VanderWerf, A. Boromand, M. D. Shattuck, and C. S. O’Hern, Pressure dependent shear response of jammed packings of frictionless spherical particles, Physical Re- view Letters 124, 038004 (2020)
work page 2020
-
[47]
B. Hendrix and L. Yu †, Self-consistent elastic properties for transversely isotropic polycrystals, Acta Materialia 46, 127 (1998)
work page 1998
- [48]
-
[49]
S. I. Ranganathan and M. Ostoja-Starzewski, Universal elastic anisotropy index, Physical Review Letters 101, 055504 (2008)
work page 2008
-
[50]
C. M. Kube, Elastic anisotropy of crystals, AIP Advances 6, 095209 (2016)
work page 2016
-
[51]
H. Pan, A. M. Ganose, M. Horton, M. Aykol, K. A. Pers- son, N. E. R. Zimmermann, and A. Jain, Benchmarking coordination number prediction algorithms on inorganic crystal structures, Inorganic Chemistry 60, 1590 (2021)
work page 2021
- [52]
-
[53]
T. C. Le and D. A. Winkler, Discovery and optimization of materials using evolutionary approaches, Chemical Re- views 116, 6107 (2016)
work page 2016
-
[54]
A. Zunger, Inverse design in search of materials with tar- get functionalities, Nature Reviews Chemistry 2, 0121 (2018)
work page 2018
-
[55]
C. M. Donahue, M. C. Remillieux, G. Singh, T. Ulrich, R. J. Migliori, and T. A. Saleh, Measuring the elastic tensor of a monolithic sic hollow cylinder with resonant ultrasound spectroscopy, NDT & E International101, 29 (2019)
work page 2019
-
[56]
D. Fan, Z. Mao, J. Yang, and J.-F. Lin, Determination of the full elastic tensor of single crystals using shear wave velocities by brillouin spectroscopy, American Mineralo- gist 100, 2590 (2015)
work page 2015
-
[57]
S. Forte and M. Vianello, A unified approach to in- variants of plane elasticity tensors, Meccanica 49, 2001 (2014)
work page 2001
-
[58]
N. Auffray and P. Ropars, Invariant-based reconstruction of bidimensional elasticity tensors, International Journal of Solids and Structures 87, 183 (2016)
work page 2016
-
[59]
C. Vall´ ee and G. D. Saxc´ e, Structure of the space of 2d elasticity tensors, Discrete and Continuous Dynam- ical Systems - Series S 6, 1525 (2013)
work page 2013
-
[60]
S. Abramian, B. Desmorat, R. Desmorat, B. Kolev, and M. Olive, Recovering the normal form and symmetry class of an elasticity tensor, Journal of Elasticity 142, 1 (2020)
work page 2020
-
[61]
R. Desmorat, N. Auffray, B. Desmorat, B. Kolev, and M. Olive, Generic separating sets for three-dimensional elasticity tensors, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences 475, 20190056 (2019)
work page 2019
-
[62]
J. T. Browaeys and S. Chevrot, Decomposition of the elastic tensor and geophysical applications, Geophysical Journal International 159, 667 (2004)
work page 2004
-
[63]
B. Shoshany, OGRe: An object-oriented general relativ- ity package for Mathematica, Journal of Open Source Software 6, 3416 (2021)
work page 2021
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.