[en] Phase stability of NaZr₂(PO₄)₃ has been studied through density functional theory calculations of elastic constants, equation of state and enthalpies. The changes in elastic constants as a function of pressure show that the ambient rhombohedral (R3-bar c) NaZr₂(PO₄)₃ becomes unstable above 8 GPa and this instability is driven by a softening of C₄₄ elastic constant through one of the Born stability criteria. High-pressure equation-of-state and enthalpy calculations further show that the ambient rhombohedral (R3-bar c) structure transforms first into another rhombohedral (R3) phase and subsequently to LiZr₂(PO₄)₃-type orthorhombic phase at pressures above 6 and 8 GPa respectively which are in agreement with a recent X-ray diffraction study. Analysis of interatomic distances show that LiZr₂(PO₄)₃-type orthorhombic structure allows for shorter Na–O and Zr–P bonds at high pressures which appears to enable strong bonding and stability. Calculated formation enthalpy and bulk modulus of the ambient phase of NaZr₂(PO₄)₃ are found to be in reasonable agreement with the respective experimental values. (author)