An Atlas of Extreme Properties in Cubic Symmetric Metamaterials

Abstract

Current research on three-dimensional metamaterial has largely focused on conventional strut, plate, and shell-based lattice designs. Although these designs offer several advantages, they possess inherent limitations that can restrict their performance in certain applications, motivating the exploration of alternative structural topologies. Here, we present a large-scale, symmetry guided framework for the generation and analysis of architected metamaterials based on all 36 cubic space groups. Using a voxel-based representation, we construct a database of approximately 1.95 million periodic unit cells spanning a broad range of relative densities and topological complexity. This dataset reveals a rich elastic property landscape shaped by crystallographic symmetry, including rare pentamode designs with high bulk to shear ratios such as $K/G \approx 166$ , isotropic-auxetic architectures with Poisson's ratio as low as $ν=-0.76$, and structures achieving up to 93% of the Hashin-Shtrikman upper bound on stiffness. Complementing the dataset, we develop a three-dimensional convolutional neural network surrogate model trained and evaluated on the full database to predict strain-energy density values under uniaxial, shear, and hydrostatic loading. Together, this work establishes a comprehensive atlas of cubic symmetric metamaterials and provides a pre-trained model for the accelerated discovery and design of 3D architected materials with extreme mechanical properties.

Figures (6)

• Figure 1: Dataset composition and elastic-property landscape. (a) Frequency of samples per crystal symmetry group across the full dataset ($N=1{,}945{,}333$). (b) Structure statistics, displaying the proportion of anisotropic, isotropic, auxetic and optimal samples within the dataset, with sample unit cells belonging to those groups displayed alongside. (c) Distribution of stiffness in $(C_{11}, C_{12}, C_{44})$ space, colored by space group, with representative unit cell topologies highlighted in the property manifold.
• Figure 2: The elastic properties of various classes of extreme structures. (a) Relative density plotted against normalized Young's modulus. (b) Poisson's ratio plotted against normalized Young's modulus. (c) Normalized Young's modulus plotted against normalized bulk modulus. (d) Normalized Young's modulus plotted against normalized shear modulus.
• Figure 3: Distributions of elastic moduli and Zener anisotropy across representative cubic space groups within the low relative density regime, balanced by random down sampling ($\rho \in [0.05,0.2]$). For clarity, only a subset of representative cubic space groups are shown. The selected space groups span multiple Bravais lattices and symmetry classes and exhibit distinct property distributions. Space groups 196, 227, 199, 230, 195, and 221 are labeled using Hermann–Mauguin notation. The red dashed lines indicate the global median across all cubic space groups visualized. (a) Zener anisotropy ratio $Z$, (b) Poisson’s ratio $\nu$, (c) normalized shear modulus $G_{norm}$, and (d) normalized Young’s modulus in the $\langle 100 \rangle$ directions $E_{norm}$ are shown for select space groups.
• Figure 4: FEA analysis of an isotropic-auxetic sample, highly anisotropic, isotropic-optimal, and pentamode sample. The stiffness tensor plot is presented to illustrate the degree of directional dependence or isotropy in the elastic response of these structures, where $\overline{E}(d)$ is the directional effective Young's modulus and $\overline{E}_{max}$ is the maximum directional Young's modulus value. (a) Left-hand side view of an isotropic-auxetic sample colored by transverse displacement in the y-axis field under uniaxial tensile loading along the x-direction (b) Warped left-hand side view of an isotropic-auxetic sample, demonstrating lateral expansion under uniaxial tensile loading along the x-direction. (c) Isotropic-auxetic stiffness tensor plot (d) Isotropic-auxetic Von Mises stress field under uniaxial tensile loading. (e) Left-hand side view of a highly anisotropic sample colored by transverse displacement in the y-axis field under uniaxial tensile loading along the x-direction (f) Warped left-hand side view of a highly anisotropic sample, demonstrating no lateral expansion under uniaxial tensile loading along the x-direction. (g) Highly anisotropic stiffness tensor plot (h) Highly anisotropic Von Mises stress field under uniaxial tensile loading. (i) Isotropic-optimal Von Mises stress field under uniaxial tensile loading. (j) Isotropic-optimal Von Mises stress field under shear loading. (k) Isotropic-optimal stiffness tensor plot (l) Isotropic-optimal Von Mises stress field under hydrostatic loading. (m) Pentamode Von Mises stress field under hydrostatic loading. (n) Pentamode Von Mises stress field under shear loading. (o) Pentamode stiffness tensor plot
• Figure 5: Overview and performance of the CNN pipeline used for stiffness prediction. (a) CNN pipeline. (b) Accuracy of predictions on uniaxial, shear and hydrostatic strain energy densities. (c) Explainability analysis on an isotropic-auxetic structure, with the gradients of the saliency map normalized from 0 to 1.