Search

Rotational invariance, a core design principle in graph neural networks (GNNs) for materials property prediction, has enabled remarkable accuracy for scalar quantities such as formation energy and band gap yet becomes a fundamental liability when the target property is a tensor of rank two or higher. Anisotropic elastic constants, which describe the directional stiffness of crystals and enter engineering design as a 6×6 Voigt matrix with up to 21 independent components in triclinic systems, must transform covariantly under rotation; any architecture that collapses all orientational information into a single invariant scalar cannot, in principle, recover the required tensorial structure. This failure-mode analysis articulates four distinct mechanisms by which invariance induces systematic error—diagonal indistinguishability, off-diagonal blindness, shear underspecification, and symmetry-induced degeneracy—and demonstrates how each arises directly from the interplay of global pooling, distance-only edge features, and symmetry averaging that characterize rotationally invariant models such as CGCNN, MEGNet, and SchNet. Detection principles based on equality tests, anisotropy-ratio checks, and differential error audits provide practitioners with practical diagnostics, while mitigation principles centered on equivariant architectures, directional encoding, and stress-tensor multi-task learning chart a clear path toward reliable tensorial prediction. The implications extend beyond elasticity to any directional property in low-symmetry crystals and call for a principled redesign of materials GNNs that treats invariance as a deliberate choice rather than a default.

Machine learning interatomic potentials (MLIPs) are widely validated against density functional theory (DFT), yet this practice risks reproducing DFT’s systematic biases rather than true experimental behavior. We present a concise four-level hierarchical framework for rigorous experimental validation focused on phonon dispersion and elastic constants. Level 1 evaluates static room-temperature elastic constants; Level 2 assesses zone-center phonon frequencies (Raman/IR); Level 3 examines full phonon dispersion curves from inelastic neutron/X-ray scattering; and Level 4 tests temperature-dependent shifts and softening. For each level we provide key observables, simple comparison metrics, uncertainty-aware criteria, and typical failure modes. A practical validation protocol and minimum reporting standards (Level 1 on ≥3 materials, Level 2 on ≥1) are proposed. This framework complements DFT benchmarks by anchoring MLIPs in experimental reality, improving reliability for thermal, mechanical, and vibrational applications in materials science.