Structural, Electronic, Mechanical, and Thermodynamic Properties of Na Deintercalation from Olivine NaMnPO4: First-Principles Study

The impact of Na atom deintercalation on olivine NaMnPO4 was investigated in a first-principle study for prospective use as cathode materials in Na-ion batteries. Within the generalized gradient approximation functional with Hubbard (U) correction, we used the plane-wave pseudopotential approach. The calculated equilibrium lattice constants are within 5% of the experimental data. The difference in equilibrium cell volumes for all deintercalated phases was only 6%, showing that NaMPO4 is structurally more stable. The predicted voltage window was found to be between 3.997 and 3.848 V. The Na1MnPO4 and MnPO4 structures are likely to be semiconductors, but the Na0.75MnPO4, Na0.5MnPO4, and Na0.25MnPO4 structures are likely to be metallic. Furthermore, all independent elastic constants for NaxMPO4 structures were shown to meet the mechanical stability requirement of the orthorhombic lattice system.


Introduction
Sodium-ion batteries (SIBs) are considered a feasible alternative to lithium-ion batteries (LIBs) as they have the potential to be much less expensive, safer, and ecologically friendly [1]. Based on the promising results reported in previous studies [1][2][3], several discoveries revealed that there is a significant opportunity for materials design and development in low-temperature Na-ion battery chemistries. During charging and discharging, Na + ions shuttle between the positive and negative electrodes in Na-ion batteries, similarly to Li-ion batteries, with an electrolyte acting as a transport medium for those ions. A critical stage in the deployment of large-scale SIB applications is the search for suitable electrode materials that meet long-cycle stability criteria and can store and transport huge amounts of energy quickly. Several phases that allow Li-ion (de) intercalation have been explored as Na + hosts, including layered metal oxides Na x MO 2 [4][5][6] and polyanionic compounds [7,8], due to their similar intercalation chemistries. However, using standard Na equivalents, performance comparable to the high-rate capability and cycling stability of LIB electrodes is difficult to achieve [9]. The main reason for this is that Na + has a much larger ionic radius (r = 0.102 nm) than Li+ (r = 0.076 nm), resulting in slower ion mobility and larger volume changes in its host structure [10]. In layered oxides produced through cubic closed-packed oxide arrays, Na + transport is slow, and Na ion deintercalation causes complex phase transitions, resulting in rapid electrode degradation [6].
Lithium transition metal phosphates of olivine type, such as LiFePO 4 , have increased cycle and thermal stability as cathode materials compared to systems based on lithium oxide [11][12][13]. Furthermore, although the volumetric discrepancy between the olivine Materials 2022, 15, 5280 2 of 14 structured LiFePO 4 and the de-lithiated FePO 4 is 6.9%, the mismatch between the Na analogs is 17.58% [14]. During cycling, such a large volume shift has a detrimental impact on structural stability. In this setting, novel Na host materials that promote viable Na + (de)intercalation while requiring minimal volume change are crucial. The expense of electrode materials designed for SIBs is another important concern. Despite the fact that Na ions are a low-cost charge carrier, many reported cathode materials still rely on redox centers that are scarce and/or require hazardous transition metal components such as cobalt [4,15], nickel [15], vanadium [16] and chromium [6]. As a result, the cost and sustainability advantages of SIBs are considerably reduced. Large-scale applications necessitate the development of electrode materials based on Earth's abundance of transition metals, with Fe and Mn being the most common transition metals in the Earth's crust. Inspired by the LiMnPO 4 structure, NaMnPO 4 has recently attracted the attention of several researchers as a cathode material for SIBs. NaMnPO 4 has two structural variations compared to LiMnPO 4 namely, maricite and olivine type [17,18]. The phosphate group framework is the same in both configurations (space group Pnma) ( Figure 1). The difference arises from the fact that in the olivine structure (Figure 1a), Na + and Mn 2+ ions prefer to occupy two octahedral positions, but in the maricite-type structure, the opposite is true (Figure 1b). In the olivine structure, the Na + containing octahedral sites share edges and form zigzag chains along the b-axis [17].
Lithium transition metal phosphates of olivine type, such as LiFePO4, have increased cycle and thermal stability as cathode materials compared to systems based on lithium oxide [11][12][13]. Furthermore, although the volumetric discrepancy between the olivine structured LiFePO4 and the de-lithiated FePO4 is 6.9%, the mismatch between the Na analogs is 17.58% [14]. During cycling, such a large volume shift has a detrimental impact on structural stability. In this setting, novel Na host materials that promote viable Na + (de)-intercalation while requiring minimal volume change are crucial. The expense of electrode materials designed for SIBs is another important concern. Despite the fact that Na ions are a low-cost charge carrier, many reported cathode materials still rely on redox centers that are scarce and/or require hazardous transition metal components such as cobalt [4,15], nickel [15], vanadium [16] and chromium [6]. As a result, the cost and sustainability advantages of SIBs are considerably reduced. Large-scale applications necessitate the development of electrode materials based on Earth's abundance of transition metals, with Fe and Mn being the most common transition metals in the Earth's crust. Inspired by the LiMnPO4 structure, NaMnPO4 has recently attracted the attention of several researchers as a cathode material for SIBs. NaMnPO4 has two structural variations compared to LiMnPO4 namely, maricite and olivine type [17,18]. The phosphate group framework is the same in both configurations (space group Pnma) ( Figure 1). The difference arises from the fact that in the olivine structure (Figure 1a), Na + and Mn 2+ ions prefer to occupy two octahedral positions, but in the maricite-type structure, the opposite is true (Figure 1b). In the olivine structure, the Na + containing octahedral sites share edges and form zigzag chains along the b-axis [17]. Compared to olivine, maricite has been found to be the most thermodynamically stable phase. The maricite phase is 0.016 eV/formula unit more stable than the olivine phase [19,20]. This material is unlikely to be used in the construction of sodium ion batteries due to the edge-sharing MnO6 octahedrons sharing edges and the lack of cationic channels for Na diffusion. In this work, the structural, thermodynamic, electrical, and mechanical properties of olivine NaMnPO4 structures are examined to determine their stability and competency as cathode materials. Formation energies, electrical densities of states (DOS), and elastic constants will be computed to model the stability trend at 0 K. Compared to olivine, maricite has been found to be the most thermodynamically stable phase. The maricite phase is 0.016 eV/formula unit more stable than the olivine phase [19,20]. This material is unlikely to be used in the construction of sodium ion batteries due to the edge-sharing MnO 6 octahedrons sharing edges and the lack of cationic channels for Na diffusion. In this work, the structural, thermodynamic, electrical, and mechanical properties of olivine NaMnPO 4 structures are examined to determine their stability and competency as cathode materials. Formation energies, electrical densities of states (DOS), and elastic constants will be computed to model the stability trend at 0 K.

Computational Method
The Vienna ab initio simulation package (VASP) [21] program was used for all the calculation, in the projector augmented wave (PAW) pseudopotential, which is based on density functional theory (DFT) [22]. A spin-polarized generalized gradient approximation (GGA) was used to solve the Kohn-Sham equations with the Perdew-Burke-Ernzerhof functional (PBE) exchange-correlation functional, with plane-wave pseudopotential [23]. The orbitals Na (3s), Mn (3d, 4s), O (2s, 2p) and P (3p) are studied for their valence states. The Hubbard correction term U (i.e., GGA + U) was also utilized because of the high Coulomb repulsion at site among Mn 3d electrons. We use an effective U value of 5.0 eV for Mn in this study, which was chosen based on previous research [24] and validated by our test calculations. The results presented here are based on the ferromagnetic configuration because the ground-state energies of the molecule determined from other spin configurations are so close. The Monkhorst-Pack scheme was used, with a 5 × 5 × 7 k-point mesh. The cutoff energy for the plane wave function was 560 eV. When the final force on all ions is less than 0.01 eV/Å, the atomic locations and lattice parameters are considered totally relaxed. The elastic properties were calculated by means of Taylor expansion of total energy using a strain of 0.005. Lastly, the densities of states were calculated using Gaussian smearing with a width of 0.05 eV for both spin up and down orientations.

Structural Properties
The calculated structural properties of NaMnPO 4 are listed in Table 1. The calculated lattice parameters, a, b, and c are within 2% of the deviation from the experimental results, and the volume deviation was found to be within 5% of the deviation from the reported results. Differences between the calculated and experimental values are found to be small, indicating that the applied theoretical method is moderate. The Na extraction process from the NaMnPO 4 lattice is investigated to simulate the discharge process. During every extraction stage, the Na atom at the original binding site is removed assuming a topotactic approach, which is followed by the relaxation of the resulting structure. The deintercalated Na x MnPO 4 structures are shown in Figure 2. The lattice parameters and changes in cell volume of Na x MnPO 4 systems during Na extraction are listed in Table 2. During Na extraction from Na x MnPO 4 (x = 1, 0.75, 05, 0.25, 0), the lattice parameter and changes in cell volume changes are less than 5%, indicating that structural stability is well maintained throughout the various stages. However, after extraction of the fourth Na atom (x = 0), where the structure is fully deintercalated, the cell volume changed by 5.9%, which was found to be higher than the rest. Note that this volume change is not large enough to cause irreversible destruction of the structure.
The results of a fully deintercalated structure are found to have a lower percentage of deviation than those previously reported by Fang et al. [14]. The study reported on the high-performance olivine NaFePO 4 microsphere cathode synthesized by the aqueous electrochemical displacement method, with a volumetric mismatch between olivine structured LiFePO 4 and de-lithiated FePO 4 of 6.9%, while Na analogs were found to have a larger mismatch of 17.58% [10].  The results of a fully deintercalated structure are found to have a lower percentage of deviation than those previously reported by Fang et al. [14]. The study reported on the high-performance olivine NaFePO4 microsphere cathode synthesized by the aqueous electrochemical displacement method, with a volumetric mismatch between olivine structured LiFePO4 and de-lithiated FePO4 of 6.9%, while Na analogs were found to have a larger mismatch of 17.58% [10].

Electronic Properties
We calculated the densities of states (DOS) of magnetic NaxMPO4 systems in the magnetic spin-polarized state to elucidate their electronic conductivity. Total and orbital partial DOS are separated, generating a band gap near the Fermi level separating the conduction band and the conduction band, as shown in Figure 3a. To comprehend the electronic conductivity of materials, concepts of Fermi level and band gaps are required. The Fermi

Electronic Properties
We calculated the densities of states (DOS) of magnetic Na x MPO 4 systems in the magnetic spin-polarized state to elucidate their electronic conductivity. Total and orbital partial DOS are separated, generating a band gap near the Fermi level separating the conduction band and the conduction band, as shown in Figure 3a. To comprehend the electronic conductivity of materials, concepts of Fermi level and band gaps are required. The Fermi energy is employed as the energy scale's zero, and the energy band gaps for both spin-up and spin-down states were observed. The states near the Fermi level are mostly Mn 3d and O 2p, with just minor contributions from the Na and P states. Due to the smaller band gap, it can be concluded that the spin-up states are primarily insulators (according to theoretical data [26]), whereas the spin-down states are semiconductors. The partial DOS of deintercalated Na x MPO 4 structures are presented in Figure 3a-e. For the initial structure Na 1 MnPO 4 , the system was found to be a magnetic insulator with a direct band gap of 3.363 eV, as shown in Figure 3a. The valence band (spin-up) maximum is located near (0.00 0.00 0.00), at −0.203 eV with respect to the Fermi level. The conduction band (spin-up) minimum is located near (0.00 0.00 0.00), at 3.161 eV with respect to the Fermi level. The center of the gap is located at 1.479 eV with respect to the Fermi level. energy is employed as the energy scale's zero, and the energy band gaps for both spin-up and spin-down states were observed. The states near the Fermi level are mostly Mn 3d and O 2p, with just minor contributions from the Na and P states. Due to the smaller band gap, it can be concluded that the spin-up states are primarily insulators (according to theoretical data [26]), whereas the spin-down states are semiconductors. The partial DOS of deintercalated NaxMPO4 structures are presented in Figure 3a-e. For the initial structure Na1MnPO4, the system was found to be a magnetic insulator with a direct band gap of 3.363 eV, as shown in Figure 3a. The valence band (spin-up) maximum is located near (0.00 0.00 0.00), at −0.203 eV with respect to the Fermi level. The conduction band (spinup) minimum is located near (0.00 0.00 0.00), at 3.161 eV with respect to the Fermi level. The center of the gap is located at 1.479 eV with respect to the Fermi level.  Figure 3e, suggesting that some states of the valence band jump the Fermi level barrier to the conduction band. The energy band gap value and the location of the Fermi level suggest that Na0.75MnPO4 in Figure 3b is semi-metallic, resulting in good electrical conductivity in Na-ion batteries. This  Figure 3b is semi-metallic, resulting in good electrical conductivity in Na-ion batteries. This metallicity increases as the number of Na atoms extracted increases. However, interestingly, a fully deintercalated MnPO 4 system was found to be a magnetic semiconductor with a direct band gap of 0.194 eV. The maximum valence band (spin-up) is located near (0.00 0.00 0.00), at −0.100 eV with respect to the Fermi level, whereas the conduction band (spin-up) minimum is located near (0.00 0.00 0.00) at 0.094 eV with respect to the Fermi level. The center of the gap is located at −0.003263 eV with respect to the Fermi level. The Fermi energy is used as the zero of the energy scales. Regarding spin-down, it was found that Mn 3d had moved from ≈5.1 to ≈3.7 eV. Generally, it was observed that the partial density of states reveals that the Mn 3d states contribute significantly to both the conduction band and the valence band, whereas the O 2p states contribute more to the valence band. In addition, there are minimal contributions from the Na and P states. It was also noted that oxygen contributed significantly to the upper states of Na 4 MnPO 4 . Thus, most of the charge compensation for the sodium removal was obtained by oxygen oxidation. Further discussion on which reaction occurs after sodium removal would require considering the kinetics and environment, such as temperature and oxygen pressure. Additionally, it was observed that during each deintercalation stage, the magnetic moments of Mn altered, decreasing from 4.5 µ B for Na 4 MnPO 4 to 3.85 µ B for MnPO 4 . Hence, the total magnetic moments for the cell were found to be decreasing as the number of Na atoms removed increase; these total cell magnetic moments ranged from 20 µ B for Na 4 MnPO 4 to 16 µ B for MnPO 4 . We used the Bader charge analysis in this study to reveal the change in the electronic population at ionic centers throughout the Na deintercalation process. It was discovered that Mn and O engage in the redox reaction to a lesser extent, contributing an average of 1.757 e and 1.403 e, respectively.

Voltages/Redox Potential and Formation Energy
The energy change of the individual Na atom at every deintercalation stage of the NaMnPO 4 structure is calculated using analytical expression: where ∆E is the change in energy of the process that y(Na) atoms are extracted between the NaMnPO 4 layers, E Na x−y MnPO 4 , E(Na) and E(Na x MnPO 4 ) are the energy of Na x−y MnPO 4 , Na metal and Na x MnPO 4 respectively. Accordingly, the cathode voltage/redox potential is calculated as follows: The cathode voltage/redox potential is U, the absolute value of the electron charge is e, and the number of Na atoms extracted is y. Equation (2) was used to calculate the average voltage/average potential at each deintercalation stage, which is shown in Figure 4. The voltage or potential required to remove the Na ions from Na x MnPO 4 ranged between 3.997 and 3.848 V. Moreover, the Na x Mn 2 O 4 potentials are consistent with previously calculated data for Mn-based isostructures olivine LiMnPO 4 [26] and maricite NaMnPO 4 [19].  To understand the formation of solid solutions of Na removal NaxMnPO4, the formation energy per unit of formula was calculated using the expression:  To understand the formation of solid solutions of Na removal Na x MnPO 4 , the formation energy per unit of formula was calculated using the expression: where E(Na x MPO 4 ) is the energy of the partially deintercalated material, and E(NaMnPO 4 ) and E(MnPO 4 ) are the energies of the pristine and totally deintercalated structures, respectively. For Na 1 MnPO 4 , Na 0 . 75 MnPO4, Na 0 . 5 MnPO 4 , Na0.25MnPO4 and MnPO 4 , the predicted formation energies are −60.12 eV, −45.291, −30.153, and 0, respectively. With the exception of the last stage, when x = 0, we find that the formation energy of the intercalation phases is negative for all materials. The negative value indicates that these materials are predicted to have a solid solution.

Elasticity
Fundamental solid-state features, such as the equation of state, interatomic potentials, lattice constants, and phonon spectra, are all connected to elastic properties. They include crucial information about the strength of a material against an externally imposed strain and are used as stability criteria in the study of structural stability modifications of mechanical stability [27,28]. Elastic constants of a material, as defined by the bulk modulus (B), shear modulus (G), Young's modulus (E), Poisson's ratio (v) and shear anisotropy factor, describe its response to the external applied strain required to maintain a given deformation and provide useful information about the material's strength (A). Born [29] was the first to calculate elastic constants. The Born stability requirements are a collection of conditions on the elastic constants (C ij ) that are linked to a crystal's internal energy shift in the second order during formation. However, the Born stability ranges were later revealed to be sensitive to the choice of locations. Table 3 shows the elastic constants determined using a Taylor expansion obtained from Equation (4) [30].
where U(V 0 , 0) is the unstrained system energy, V 0 is the equilibrium volume, τ i is the element in the stress tensor, and δ i is a Voigt index factor. The nine independent elastic constants (C ij ) of orthorhombic NaMPO 4 structures are shown in Table 3. For orthorhombic systems, the Born mechanical stability criteria are [26,27,31,32]: We noted that all stability criteria were met, which shows that the olivine structure's NaMPO 4 compounds were mechanically stable. During the extraction phases, the first to third Na extraction phases, all Born mechanical stability criteria were met. However, for the fourth Na extraction, which is a fully deintercalated structure, the C 44 > 0 criteria were not achieved since C 44 = −39.20 GPa. This phenomenon shows that NaMnPO 4 cannot be fully deintercalated and still maintain its stability. This is the first elastic constant recorded in the Na x MnPO 4 system when Na was removed.
Using the Voigt-Reuss-Hill method, the macroscopic mechanical parameters of the bulk, shear and Young's moduli are calculated from the calculated elastic constants [30].
where B, G and E are the bulk, shear and Young's moduli, respectively, while V, R and H are the Voigt, Reuss and Hill bounds, respectively, and S ij is the inverse matrix of the elastic constant's matrix C ij , which is given by Ravindran et al. [33]. There are three sets of results in Table 4, which include the bulk, shear and Young's moduli. Depending on the material's bulk and Young's moduli, as well as its shear modulus, its hardness and stiffness may be determined. Pressure-induced volume changes are also taken into account when calculating the bulk modulus; these results were calculated using the VASP-MT package, employing Equation (4). The positive bulk, shear, and Young's moduli of Na x MnPO 4 structures are relatively large, implying hardness, great resistance to volume change, deformation, and stiffness, respectively. The elastic characteristics of the maricite NaMnPO 4 polymorph were calculated to be higher than those reported here for the olivine phase by Lethole et al. [26]. Furthermore, our calculated C 11 , C 22 , and C 33 values are much larger than C 44 , C 55 , and C 66 , implying that the materials have strong directional resistance to linear compressions against uniaxial pressures but limited resistance to shear deformations.
Furthermore, we note that B H > G H , implying that the shear modulus (G H ) is the parameter that limits the mechanical stability of the Na 1 MnPO 4 structures [33]. In addition, the Pugh ductility and brittleness criterion was calculated. Pugh proposed the bulk-to-shear modulus (B/G) ratio for polycrystalline phases, assuming that the shear modulus represents plastic deformation resistance, and the bulk modulus represents fracture resistance.
Brittleness is related with a low B/G value, whereas ductility is related with a high B/G value. The value 1.75 is the critical number that distinguishes ductility from brittleness. Because B H /G H is greater than 1.75, the structures of the Na 1 MnPO 4 , Na 0 . 75 MnPO 4 and MnPO 4 structures are ductile, which implies that these materials can bend without deformation, resulting in fewer cracks during battery operation [34]. However, because their predicted B H /G H is less than 1.75, Na 0 . 5 MnPO 4 and Na 0 . 25 MnPO 4 are both brittle.

Anisotropy in Elastic Constants
Almost all known crystals have elastically anisotropic behavior, and a proper description of this anisotropic behavior has important implications in both engineering science and crystal physics, as well as in other fields. It is possible to calculate the shear anisotropic factors by calculating the degree of anisotropy in the bonding between atoms that are in various planes. The shear anisotropic factors provide a measure of the degree of anisotropy in the bonding between atoms in different planes. The shear anisotropic factor for the {100} and 110 , 010 directions is calculated as: For the ({010}, {001}), the shear planes between the ( 011 and 010 ( 110 , 001 , respectively) directions are: Table 5 shows the shear anisotropic factors obtained from our theoretical investigations. Factors A 1 , A 2 , and A 3 must all be one for an isotropic crystal, whereas any number less than or greater than unity indicates the crystal's degree of elastic anisotropy. For Na 1 MnPO 4 , the values of A 1 (1.085), A 2 (1.145), and A 3 (1.409) diverge from the unity by 8.50%, 14.50%, and 40.09%, suggesting an isotropic characteristic. During the deintercalation stages, all the materials displayed slight deviations from the unity, whereas A 1 and A 3 for MnPO 4 showed the greatest deviation from the unity, with A 1 having a negative value. Chung and Buessem [35] developed a new term for non-cubic systems: percent elastic anisotropy [35]. It is a measure of the amount of elastic anisotropy held by the crystal under discussion. The percentage anisotropy in compressibility and shear moduli is described by the equations A B = B V − B R /B V + B R and A G = G V − G R /G V + G R , respectively, for compressibility and shear moduli. B R = B V is related with isotropic elastic constants, but a value of 100% is associated with the maximum amount of anisotropy that may be achieved. Deintercalations of Na tend to raise A B and A G values, suggesting an improvement in isotropy and a reduction in microcracks and dislocations in the charge/discharge process, apart from the MnPO 4 sample. In all stages of deintercalation, the proportion of shear modulus anisotropy A G is lower than the percentage of bulk modulus anisotropy A B . Table 5. Anisotropy in the shear elastic factor (A i with i = 1,2,3) and anisotropy in compressibility and shear moduli (A B and A G in %).

Debye Temperature
We have calculated the Debye temperature (θ θ ) from the average sound velocity (V m ) using the equation [36]: where h and k B are the Planck's and Boltzmann's constants, respectively, and V a is the atomic volume. The average sound velocity in polycrystalline systems, ν m , is evaluated by the expression: where ν t and ν l are the mean longitudinal and transverse sound velocities, which can be related by the shear and bulk moduli from Navier's equations: Table 6 shows the calculated volumetric density, sound velocities, and Debye temperature as a function of Na removal. Derived from the calculated elastic constants, the calculated Debye temperature θ D for Na 1 MnPO 4 was found to be 512.7 K. Unfortunately, we were unable to uncover any experimental or theoretical data with which to compare our computed results. As seen in Table 5, Na deintercalation resulted in a monotonic drop in all sound velocities as well as the θ D for all the systems. Considering that the θ D in a solid may be used to define its covalent strength, the drop in θ D as a results of Na deintercalation indicated that the covalent strength of NaMnPO 4 decreases with each deintercalation step. Although the thermal conductivity of the material increases with increasing θ D in general, our results indicated a modest drop in θ D , indicating that the thermal conductivity of the material was maintained after Na deintercalation. It is noteworthy that the volumetric density, sound velocities, and Debye temperature increase during the fourth step of Na deintercalation, with a Debye temperature calculated to be 497.8 K. Using the MedeA VASP program, certain essential thermodynamic parameters such as specific heat at constant volume (C v ), free energy, and entropy are calculated at 0 GPa pressure and temperatures (0-1200 K) to investigate thermodynamic properties. From Figures 5-7, the influence of temperature on C v , free energy, and entropy has been calculated and shown to be significant. At low temperatures, the rate of increase in C v is rapid for all intercalation phases, regardless of the intercalation stage. During high-temperature intercalation, the value of C v of all intercalation stages approaches the classical asymptotic limit, which is also known as the Dulong-Petit limit [37]. At low temperatures, the fluctuation of free energy versus temperature, as well as the variation of C v versus temperature, shows that the increase in free energy is significant for all intercalation stages. Unfortunately, no experimental data or other theoretical conclusions are available to compare. We hope that our research will be useful in future experiments. for all intercalation phases, regardless of the intercalation stage. During high-temperature intercalation, the value of Cv of all intercalation stages approaches the classical asymptotic limit, which is also known as the Dulong-Petit limit [37]. At low temperatures, the fluctuation of free energy versus temperature, as well as the variation of Cv versus temperature, shows that the increase in free energy is significant for all intercalation stages. Unfortunately, no experimental data or other theoretical conclusions are available to compare. We hope that our research will be useful in future experiments.

Conclusions
The effect of Na atom deintercalation on the structural, electronic, mechanical, and thermodynamic properties of NaMnPO4 has been investigated by first-principle calculations. The calculated lattice constants show good agreement with the experimental values to within 5%. During the Na deintercalation stages, the lattice parameters and the volume showed a deviation of less than 6%, which is not enough to cause irreversible distraction to the system. Electronic DOS revealed that during the Na removal stages, between the 1st and 3rd stages, the material showed an increase in metallicity, while on the other hand during the 4th stage, the material showed a semiconductor behavior with a band gap of

Conclusions
The effect of Na atom deintercalation on the structural, electronic, mechanical, and thermodynamic properties of NaMnPO 4 has been investigated by first-principle calculations. The calculated lattice constants show good agreement with the experimental values to within 5%. During the Na deintercalation stages, the lattice parameters and the volume showed a deviation of less than 6%, which is not enough to cause irreversible distraction to the system. Electronic DOS revealed that during the Na removal stages, between the 1st and 3rd stages, the material showed an increase in metallicity, while on the other hand during the 4th stage, the material showed a semiconductor behavior with a band gap of 0.194 eV. The voltage window of 3.997 to 3.848 V was obtained, and the calculated formation energy values were found to be negative, which symbolizes the prediction of a solid possessed by the material. The calculated elastic constants suggested mechanical stability for NaMnPO 4 since the stability criteria were satisfied for all deintercalated systems except fully deintercalated MnPO 4 . According to the Pugh criterion of ductility and brittleness, we note that the Na 1 MnPO 4 , Na 0 . 75 MnPO 4 , and MnPO 4 structures are ductile, while Na 0 . 5 MnPO 4 and Na 0 . 25