Structure prediction of stable sodium germanides at 0 and 10 GPa

In this work we used ab-initio random structure searching (AIRSS) to carry out a systematic search for crystalline Na-Ge materials at both 0 and 10 GPa. The high-throughput structural relaxations were accelerated using a machine-learned interatomic potential (MLIP) fit to density-functional theory (DFT) reference data, allowing ∼1.5 million structures to be relaxed. At ambient conditions we predict three new Zintl phases, Na_3Ge_2, Na_2Ge and Na_9Ge_4, to be stable and a number of Ge-rich layered structures to lie in close proximity to the convex hull. The known Na_δGe_34 clathrate and Na_4Ge_13 host-guest structures are found to be relatively stabilized at higher temperature by vibrational contributions to the free energy. Overall, the low energy phases exhibit exceptional structural diversity, with the expected mixture of covalent and ionic bonding confirmed using the electron-localisation function (ELF). The local Ge structural motifs present at each composition were determined using Smooth Overlap of Atomic Positions (SOAP) descriptors and the Ge-K edge was simulated for representatives of each motif, providing a direct link to experimental x-ray absorption spectroscopy (XAS). Two Ge-rich phases are predicted to be stable at 10 GPa; NaGe_3 and NaGe_2 have simple kagome and simple hexagonal Ge lattices respectively with Na contained in the pores. NaGe_3 is isostructural with the MgB_3 and MgSi_3 family of kagome superconductors and remains dynamically stable at 0 GPa. Removing the Na from NaGe_2 results in the hexagonal lonsdalite Ge allotrope, which has a direct band gap.