Atomic stiffness for bulk modulus prediction and high-throughput screening of ultraincompressible crystals

Determining bulk moduli is central to high-throughput screening of ultraincompressible materials. However, existing approaches are either too inaccurate or too expensive for general applications, or they are limited to narrow chemistries. Here we define a microscopic quantity to measure the atomic stiffness for each element in the periodic table. Based on this quantity, we derive an analytic formula for bulk modulus prediction. By analyzing numerous crystals from first-principles calculations, this formula shows superior accuracy, efficiency, universality, and interpretability compared to previous empirical/semiempirical formulae and machine learning models. Directed by our formula predictions and verified by first-principles calculations, 47 ultraincompressible crystals rivaling diamond are identified from over one million material candidates, which extends the family of known ultraincompressible crystals. Finally, treasure maps of possible elemental combinations for ultraincompressible crystals are created from our theory. This theory and insights provide guidelines for designing and discovering ultraincompressible crystals of the future.