Evolutionary search for superconducting phases in the lanthanum-nitrogen-hydrogen system with universal neural network potential

Abstract: Recently, Grockowiak $\textit{et al.}$ reported "hot superconductivity" in ternary or multinary compounds based on lanthanum hydride [A. D. Grockowiak $\textit{et al.}$, Front. Electron. Mater. $\textbf{2}$, 837651 (2022)]. In this paper, we explored thermodynamically stable phases and superconducting phases in the lanthanum-nitrogen-hydrogen system (La${x}$N${y}$H${1-x-y}$, $0 \leq x \leq 1$, $0 \leq y \leq 1$) at pressure of 20$\,$GP. We rapidly and accurately constructed the formation-enthalpy convex hull using an evolutionary construction scheme based on density functional theory calculations, extracting the candidates for stable and moderately metastable compounds by the universal neural network potential calculations. The convex hull diagram shows that more than fifty compounds emerge as stable and moderately metastable phases in the region of $\Delta H \leq 4.4$$\,$mRy/atom. In particular, the compounds are concentrated on the line of $x = 0.5$ connecting between LaH and LaN. We found that the superconductivity is gradually enhanced due to N doping for LaH and the superconducting critical temperature $T{\rm c}$ reaches 8.77$\,$K in La$2$NH with $y = 0.25$. In addition, we predicted that metastable La$_2$NH$_2$ shows the highest $T{\rm c}$ value, 14.41$\,$K, of all the ternary compounds predicted in this study. These results suggest that it is difficult to obtain the hot superconductivity in the La-H compounds with N at 20$\,$GPa.