First Principles Study on the Thermodynamic and Elastic Mechanical Stability of Mg 2 X (X = Si,Ge) Intermetallics with (anti) Vacancy Point Defects

: In this paper, based on the density functional theory, through thermodynamic and mechanical stability criteria, the crystal cell model of intermetallic compounds with vacancy and anti-site point defects is constructed and the lattice constant, formation heat, binding energy, elastic constant, and elastic modulus of Mg 2 X (X = Si, Ge) intermetallics with or without point defects are calculated. The results show that the di ﬀ erence in the atomic radius leads to the instability and distortion of crystal cells with point defects; Mg 2 X are easier to form vacancy defects than anti-site defects on the X (X = Si, Ge) lattice site, and form anti-site defects on the Mg lattice site. Generally, the point defect is more likely to appear at the Mg position than at the Si or Ge position. Among the four kinds of point defects, the anti-site defect x Mg is the easiest to form. The structure of intermetallics without defects is more stable than that with defects, and the structure of the intermetallics with point defects at the Mg position is more stable than that at the Si / Ge position. The anti-site and vacancy defects will reduce the material’s resistance to volume deformation shear strain, and positive elastic deformation, and increase the mechanical instability of the elastic deformation of the material. Compared with the anti-site point defect, the void point defect can lead to the mechanical instability of the transverse deformation of the material and improve the plasticity of the material. The research in this paper is helpful for the analysis of the mechanical stability of the elastic deformation of Mg 2 X (X = Si, Ge) intermetallics under the service condition that it is easy to produce vacancy and anti-site defects.


Computational Method
In this paper, based on the density functional theory (DFT) [26,27], we use the CASTEP [28] module in Materials Studio 8.0. The exchange-correlation potential was in the form of Perdew-Burke-Ernzerh (PBE) of the generalized gradient approximation (GGA) [29,30]. Considering computational accuracy and efficiency [31], the final cutoff energy of plane wave E cut was set to be 380 eV. The Brillouin zone integral is in the form of 8 × 4 × 8 monkhorst-pack [32], and the super soft pseudopotential of the reciprocal space is selected as the potential function [33].
The parameter conditions are set as follows: the total energy of the system converges within 5.0 × 10 −6 ev/atom, the force on each atom is less than 0.01 eV/Å, the stress deviation is lower than 0.02 GPa, the tolerance deviation is less than 5 × 10 −4 Å.; the maximum convergence tolerance of energy and self-consistent field (SCF) convergence threshold are 5.0 × 10 −6 eV/atom and 5.0 × 10 −7 eV/atom, respectively. Both Mg 2 Si and Mg 2 Ge belong to the face-centered cubic (FCC) structure of anti-fluorite (CaF 2 ).
In this paper, we use 1 × 1 × 2 and 2 × 2 × 2 supercell models, so that we can more easily verify whether the calculated data is reasonable, and can compare the influence of different supercell sizes on the calculation results. The 1 × 1 × 2 Mg 16 X 8 supercell crystals are shown in Figure 1. There are two vacancy defects V mg , V x , two anti-site defects Mg X , X Mg in the crystals, and the small red box indicates the location of the point defect. compounds. References [22][23][24] used the first principle to calculate the interface structure of Mg2Si and related compounds. Kessair [25] studied the structural characteristics of seven different types of Mg2Si and found that the thermoelectric potential of cubic CaF2 structure is very high, which may be the preferred material in thermoelectric applications. Vacancy point defects play a decisive role in the physical properties of thermoelectric materials and are used as a new way to improve thermoelectric properties. However, there are few studies on the structures of Mg2X (X = Si, Ge) with vacancy and anti-site defects. Therefore, in this paper, the equilibrium lattice constant, formation heat, binding energy, and elastic constant of Mg2X (X = Si, Ge) system with vacancy defect structure are calculated by Density Functional Theory (DFT). The effect of vacancy defect on the structure and properties of Mg2X (X = Si, Ge) intermetallics is explored from the microscopic point of view.

Computational Method
In this paper, based on the density functional theory (DFT) [26,27], we use the CASTEP [28] module in Materials Studio 8.0. The exchange-correlation potential was in the form of Perdew-Burke-Ernzerh (PBE) of the generalized gradient approximation (GGA) [29,30]. Considering computational accuracy and efficiency [31], the final cutoff energy of plane wave Ecut was set to be 380 eV. The Brillouin zone integral is in the form of 8 × 4 × 8 monkhorst-pack [32], and the super soft pseudopotential of the reciprocal space is selected as the potential function [33].
The parameter conditions are set as follows: the total energy of the system converges within 5.0 × 10 −6 ev/atom, the force on each atom is less than 0.01 eV/Å, the stress deviation is lower than 0.02 GPa, the tolerance deviation is less than 5 × 10 −4 Å.; the maximum convergence tolerance of energy and self-consistent field (SCF) convergence threshold are 5.0 × 10 −6 eV/atom and 5.0 × 10 −7 eV/atom, respectively. Both Mg2Si and Mg2Ge belong to the face-centered cubic (FCC) structure of anti-fluorite (CaF2).
In this paper, we use 1 × 1 × 2 and 2 × 2 × 2 supercell models, so that we can more easily verify whether the calculated data is reasonable, and can compare the influence of different supercell sizes on the calculation results. The 1 × 1 × 2 Mg16X8 supercell crystals are shown in Figure 1. There are two vacancy defects Vmg, Vx, two anti-site defects MgX, XMg in the crystals, and the small red box indicates the location of the point defect.   The 2 × 2 × 2 supercell model with 96 atoms is shown in Figure 2. The small red circle indicates the position of point defect, and the point defect concentration is reduced to a quarter of the former. We have made some comparative studies with the related calculation results of 1 × 1 × 2 supercells in this work. We have made some comparative studies with the related calculation results of 1 × 1 × 2 supercells in this work.

Equilibrium Lattice Stability of Intermetallics with or without Point Defects
The equilibrium lattice constants of Mg2X (X = Si, Ge) intermetallic compounds have been calculated and shown in Table 1.
The lattice constants of Mg16Si8 and Mg16Ge8 with point defects are listed in Table.1, and Figure  2 is the comparison of lattice constants of Mg2X (X = Si, Ge) with point defects. It can be seen from Table.1 that the equilibrium lattice constants calculated in this paper are very close to the experimental values [34,35], which shows that the method adopted is reliable. It also can be seen from Table 1 that when the cell contains vacancy or anti vacancy point defects, lattice distortion occurs in the cell, resulting in cell expansion and a lattice constant increase. At the same time, as shown in Figure 3, the lattice constant of Mg16Ge8 with vacancy and anti-site defects is larger than that of Mg16Si8, which is due to the following relationship between the relative atomic radius: RSi (0.118 nm) < RGe (0.123 nm) < RMg (0.160 nm).
The calculated results of lattice parameters of 1 × 1 × 2 supercells and 2 × 2 × 2 supercells (as shown in Table 2) are compared, and the results are consistent. Table 1. Calculated equilibrium lattice constants of 1 × 1 × 2 Mg2X (X = Si, Ge) cells with or without two types of point defects.

Equilibrium Lattice Stability of Intermetallics with or without Point Defects
The equilibrium lattice constants of Mg 2 X (X = Si, Ge) intermetallic compounds have been calculated and shown in Table 1.  Table 1, and Figure 2 is the comparison of lattice constants of Mg 2 X (X = Si, Ge) with point defects. It can be seen from Table 1 that the equilibrium lattice constants calculated in this paper are very close to the experimental values [34,35], which shows that the method adopted is reliable. It also can be seen from Table 1 that when the cell contains vacancy or anti vacancy point defects, lattice distortion occurs in the cell, resulting in cell expansion and a lattice constant increase. At the same time, as shown in Figure 3, the lattice constant of Mg 16 Ge 8 with vacancy and anti-site defects is larger than that of Mg 16 Si 8 , which is due to the following relationship between the relative atomic radius: R Si (0.118 nm) < R Ge (0.123 nm) < R Mg (0.160 nm).

Effects of Vacancy and Anti-Site Points Defects on Thermodynamic Stability
The smaller the bond energy is, the less energy is needed to break the bond, the easier it is to be broken, the more unstable it is, so the harder it is to generate. The larger the bond energy is, the more energy is needed to destroy the bond. The less likely the bond is to be destroyed, the more stable the bond is, and the easier the crystal formed by the bond is. The calculation formula of the formation heat of Mg16X8 (X = Si, Ge) intermetallics with point defects is [36]:  The calculated results of lattice parameters of 1 × 1 × 2 supercells and 2 × 2 × 2 supercells (as shown in Table 2) are compared, and the results are consistent. Table 2. Calculated equilibrium lattice constants of 2 × 2 × 2 Mg 2 X (X = Si, Ge) cells with or without two types of point defects.

Effects of Vacancy and Anti-Site Points Defects on Thermodynamic Stability
The smaller the bond energy is, the less energy is needed to break the bond, the easier it is to be broken, the more unstable it is, so the harder it is to generate. The larger the bond energy is, the more energy is needed to destroy the bond. The less likely the bond is to be destroyed, the more stable the bond is, and the easier the crystal formed by the bond is. The calculation formula of the formation heat of Mg 16 X 8 (X = Si, Ge) intermetallics with point defects is [36]: where, E AB tot represents the average total energy of each crystal cell of Mg 2 X intermetallic compound, N indicates the number of atoms in the crystal cell, and E A solid , E B solid represent the ground state energy of each atom in the solid cell, E Mg solid = −973.95 eV, E Si solid = −107.26 eV, and E Ge solid = −107.29 eV. The binding energy E coh of a crystal is the energy needed to decompose the crystal into a single atom, which is usually negative. Therefore, the larger the absolute value of the binding energy, the higher the stability of the intermetallics. The calculation formula of binding energy is [37]: where, E AorB atom represents the energy of each atom in the crystal cell in the free state, E Mg atom = −972.48 eV, E Si atom = −101.85 eV and E Ge atom = −102.71 eV. The calculated formation heat and binding energy of Mg 2 X (X = Si, Ge) cells with or without point defects are listed in Table 3. The more negative ∆H f of the crystal is, the easier it is to form. That is to say, the absolute ∆H f of the crystal is larger, it is not easy to be damaged, and it is easier to form. Table 3. Computed formation heat (∆H f ) and binding energy (E coh ) of Mg 16 X 8 and Mg 64 X 32 (X=Si, Ge) with or without (anti) vacancy point defects. As shown in Table 3 that the absolute value of ∆H f of Mg 2 X (X = Si, Ge) with point defects is less than that without point defects, which means that the structure with point defects is easy to be destroyed, and it is easy to be in an unstable state before formation, that is to say; it is difficult to form compounds (For example, it can forms a stable structure with point defects at 0 K and 0 GPa, and the formed structure may become unstable under the influence of external factors, which is not contradictory to the relationship between binding energy and crystal stability in the following section of this paper).
It can also be seen from Table 3 that on the X (X = Si, Ge) lattice site, the absolute value of ∆H f of intermetallic compounds with vacancy defects is larger than that with anti-site defects, indicating that the metal bond energy with vacancy defects on the X lattice site is larger than that with anti-site defects, and it is not easy to be damaged, and it is easier to form crystals, that is to say, intermetallic compounds are easier to form vacancy defects than anti-site on the X (X = Si, Ge) lattice site. However, Mg 2 X is easier to form anti-site defects on the Mg lattice site.
It is also found that for 1 × 1 × 2 Mg 2 Si intermetallics, the ∆H f (−0.092 eV/atom) of vacancy point defects at Mg position is smaller than that at Si position (−0.084 eV/atom), the absolute value is larger, and it is not easy to be destroyed, so it is easier to form; The formation energy (−0.108 eV/atom) of the anti-site defect at Mg position is smaller than that at Si position (−0.072 eV/atom), the absolute value is larger, it is not easy to be damaged, and it is easier to form, which indicates that the point defect is more likely to appear at the Mg position than at the Si position; The same analysis shows that for Mg 2 Ge intermetallics, point defects are more likely to form at the Mg position than at the Ge  Table 3, we can also find that this conclusion is also applicable to Mg 64 X 32 (X = Si, Ge) intermetallics with 2 × 2 × 2 supercells.
It is also found that no matter 1 ×1× 2 supercells or 2 × 2 × 2 supercells, the absolute value of the binding energy of the intermetallics without point defects is higher than that of the intermetallics with point defects, which indicates that the energy required for atom bonding is higher and the structure is more stable when there is no defect. For the same type of defects, the absolute value of the binding energy at the Mg position is larger than that at the Si/Ge position, so the structure of the intermetallics forming point defects at the Mg position is more stable than that at the Si/Ge position.
The defect formation energies of vacancy and anti-site can be calculated by the following expression: where, E total (per f ect) and E total (de f ect) represent the total energy of perfect supercell and defect supercell of Mg 2 X (X = Si, Ge), respectively; n i represent the number of atoms added or removed in Mg 2 X (X = Si, Ge) intermetallics, µ i represent the corresponding chemical potential of the atom, i.e., the defect formation energies of vacancy and anti-site in Mg 2 X (X = Si, Ge) are shown in Table 4. As shown in Table 4, for Mg 2 X (X = Si, Ge), the ∆E vac and ∆E anti at the Si or Ge position is bigger than that at the Mg position, which indicates that the point defect is more likely to appear at Mg position than at Si or Ge position. Among the four kinds of point defects, the defect forming energy of the anti-site defect ∆E anti Mg is the smallest, which indicates that the anti-site defect x Mg is the easiest to form in Mg 2 X. Compared with other point defects, the ∆E anti Mg is negative, which means it can form spontaneously in 1 × 1 × 2 and 2 × 2 × 2 supercells. The defect formation of vacancy ∆E vac and inversion ∆E anti (except Mg Si ) in Mg 2 X (X = Si, Ge) can decrease with the increase of supercell size, which shows that only the anti-site defect Mg Si is easy to form in 1 × 1 × 2 supercells, while other point defects tend to form in 2 × 2 × 2 supercell.

Elastic Constants of Stability Criteria
Elastic constants are used to measure the ability of materials to resist external force and deformation and can be used to predict the mechanical stability of materials. To determine the elastic stability of a crystal, three independent elastic constants [38] are needed for the cubic structure: C 11 , C 12 , C 44 . For tetragonal structures, six independent elastic constants [39] are needed: C 11 , C 12 , C 13 , C 33 , C 44 , C 66 . In this paper, Mg 2 X (X = Si, Ge) cells with point defects belong to tetragonal structure, and Mg 2 X (X = Si, Ge) cells without point defects belong to cubic structure [40,41].
The calculated elastic constants of Mg 2 X (X = Si, Ge) cells with (without) point defects are listed in Table 5. It can be seen that our calculation results are close to others' theoretical calculation results [42], indicating that the selection of the parameters is reasonable. In the case of cubic crystals, the conditions of elastic stability are well known [43]: The elastic Born stability criteria for the tetragonal system are [43]: Based on the computed elastic constants, the effects of point defects on the mechanical stability of Mg 2 X (X = Si, Ge) intermetallics can be evaluated. Obviously, the elastic constants C ij calculated in Table 3 meet the mechanical stability criteria. Therefore, the Mg 2 X (X = Si, Ge) intermetallic compounds with and without vacancy or anti-vacancy point defects are all mechanically stable structures.

Elastic Modulus of Mechanical Stability
The yield limit σ s and strength limit σ b obtained in the tensile test reflect the bearing capacity of the material to the force, while the elongation δ or section shrinkage ψ reflect the plastic deformation capacity of the material. In order to express the difficulty of material resisting deformation in the elastic range, the meaning of material elastic modulus E is usually reflected by the rigidity of parts in practical engineering, because once the parts are designed according to the stress, the rigidity is judged by the deformation generated by the load in the service process in the elastic deformation range. Generally, the load-causing unit strain is the rigidity of the part. For example, for tension and compression parts, the rigidity is EA 0 , and A 0 is the cross-section area of the part. It can be seen that in order to improve the stiffness EA 0 of the parts, which is to reduce the elastic deformation of the parts, the materials with high elastic modulus can be selected, and the cross-sectional area A 0 of the bearing can be appropriately increased. The importance of stiffness is that it determines the stability of parts in service, especially for slender rod-shaped parts and thin-walled parts. Therefore, for the theoretical analysis and design calculation of components, the elastic modulus E is an important mechanical performance index that is often used. Elastic modulus E is defined as the ratio of stress to corresponding strain when the ideal material has small deformation. The nature of modulus depends on the nature of deformation. Elastic modulus E (unit: GPa) includes the Young's modulus E (tensile stress/tensile strain, also called positive elastic modulus, tensile modulus), shear modulus G (shear stress/shear strain angle), compression modulus K (compressive stress/shrinkage strain), bulk modulus B (bulk pressure/bulk strain)-the reciprocal of the modulus is called the flexibility expressed in J. When it is not easy to cause confusion, the elastic modulus of general metal materials refers to the young's modulus, that is, the positive elastic modulus.
From the macroscopic point of view, the modulus of elasticity is a measure of the ability of an object to resist elastic deformation. From the microscopic point of view, it is a reflection of the bonding strength among atoms, ions, or molecules. All the factors that affect the bonding strength can affect the elastic modulus of the material, such as bonding mode, crystal structure, chemical composition, microstructure, temperature, etc. The young's modulus E of metal materials will fluctuate by 5% or more due to different alloy compositions, heat treatment state and cold plastic deformation. However, generally speaking, the elastic modulus of metal material is a mechanical property index that is not sensitive to the structure. Alloying, heat treatment (fiber structure), cold plastic deformation, and other factors have little influence on the elastic modulus E, and the external factors such as temperature and loading rate have little influence on it, so the elastic modulus E is regarded as a constant in general engineering application.
Poisson's ratio ν refers to the ratio of the absolute value of the transverse positive strain and the positive axial strain when the material is under tension or compression in a single direction, also known as the transverse deformation coefficient, which is the elastic constant reflecting the transverse deformation of the material. When the material extends (or shortens) along the load direction, it will also shorten (or lengthen) in the direction perpendicular to the load. The negative value of the ratio of the strain ε l in the vertical direction to the strain ε in the load direction is called Poisson's ratio of the material, i.e., ν = −ε l /ε. Poisson's ratio ν of materials is generally determined by test method. For traditional materials, ν is generally constant in the elastic working range, but beyond the elastic range, ν increases with the increase of stress until v = 0.5. For the material with large Poisson's ratio, the amount of transverse deformation is larger than that of the longitudinal deformation after the material is stressed and before plastic deformation, otherwise, the amount of transverse deformation is smaller than that of the longitudinal deformation.
In Table 6, we compare the results of calculation when there are vacancy point defects in different supercells of Mg 2 X (X = Si, Ge) intermetallic compounds; it is found that for Mg 2 Si type intermetallics, the B, G, and E of Mg 2 X with Mg vacancy (V Mg ) are higher than those of Mg 2 X with X vacancy (V X ). The above results show that the vacancy defects in Mg position, compared with those in the Si position and Ge position, will lead to the stronger resistance to volume deformation, shear strain and elastic deformation of crystal materials, the greater the rigidity, the less deformation, the stronger the brittleness and the worse the plasticity.
For intermetallic compounds with anti-site point defects, the B, G and E of Mg 2 X (X = Si, Ge) with the X Mg anti-site are higher than those of Mg 2 X with the Mg X anti-site, indicating that when the Si or Ge atom replace the Mg atom to form anti-vacancy defects, compared with Mg atom replacing Si or Ge atom to form anti vacancy defect, the material has a strong resistance to volume deformation, shear strain and elastic deformation, but poor plasticity. Besides, for Mg 2 X (X = Si, Ge) with same point defect structure, the B, G, and E of 1 × 1 × 2 supercells are smaller than those of 2 × 2 × 2 supercells, which indicates that the Mg 2 X (X = Si, Ge) with point defects are easier to form in 2 × 2 × 2 supercells. Therefore, the point defects generated on the Mg position of 2 × 2 × 2 supercells can enhance the mechanical stability of the material, while the plasticity is poor.  Table 6 also shows that the Poisson's ratio ν (0.168 and 0.164) of Mg 16 X 8 without point defects is the smallest, and the Poisson's ratio ν with defects is larger (greater than 0.18), indicating that the formation of defects makes the transverse deformation of the material larger and the transverse elastic mechanical instability increased. It was also found that the Poisson's ratio of Mg 2 X (X = Si, Ge) with X vacancy is higher than that of Mg 2 X with the Mg X anti-site. Similarly, the Poisson's ratio of Mg 2 X with Mg vacancy was slightly higher than that of Mg 2 X with X Mg anti-site. Thus, the above results show that when there are defects in the Mg, Si, and Ge position, the void point defects increase the instability of transverse deformation more than the anti-void point defects. Note that for Mg 2 X (X = Si, Ge) with the same point defect structure (except Mg Si ), the Poisson's ratio ν of the 2 × 2 × 2 supercells is smaller than that of the 1 × 1 × 2 supercells, indicating that the point defects in larger supercells can only increase the instability of transverse deformation slightly, compared with 1 × 1 × 2 supercells.
According to Pugh's report [44], G/B can predict the ductility or brittleness of materials, and materials with G/B value higher than the critical value of 0.57 belong to brittle materials; it can be seen from Table 6 that both Mg 2 X (X = Si, Ge) with and without point defects are brittle materials. The maximum G/B value of Mg 2 X without point defects indicates that the formation of defects will improve the ductility of materials; the G/B value of vacancy point defects is smaller than that of anti-site defects, indicating that vacancy point defects can improve the plasticity of materials more than anti-site defects. Compared with 2 × 2 × 2 supercells, the G/B of Mg 2 X with the same point defect structure (except Mg Si ) in 1× 1 × 2 supercells can enhance the plasticity of materials.
In theory, only two of the three elastic constants E, G and ν of isotropic materials are independent, because they have the following relationship: G = E/[2 (1 + ν)]. This is consistent with the essence of the definition of Poisson's ratio. For anisotropic materials, the elastic modulus in x, y and z directions, shear modulus and Poisson's ratio in xy, xz, and yz planes are required

Elastic Anisotropy
Anisotropy exists widely in all kinds of materials. In short, crystals have different properties in different directions. It is easy to confuse anisotropy and inhomogeneity, which are two completely different concepts. Anisotropy means that the properties of materials are related to directions, while inhomogeneity means that the properties of materials are related to parts. In the case of single crystals, at any point in the interior, the structures and properties are the same, but it has different performance in different directions.
Crystal is anisotropic, with different properties in different directions, and has strict symmetry. Here, the symmetry is a very remarkable property, which is widely used in mathematics and physics, and also commonly used in the analysis of material properties. Polycrystalline materials have preferred orientation and certain anisotropy.
All kinds of materials are elastic, and the elastic properties of most materials are anisotropic [45,46]. For example, in cubic crystals, the [111] direction is usually more difficult to compress than the [100] direction. The three-dimensional surface Young"s modulus E of intermetallics is calculated as follows [47] 1 where S 11 , S 12 , and S 44 represent the flexibility coefficient, l 1 , l 2 , and l 3 indicate direction cosines relative to the x, y, and z directions, respectively. Figures 4 and 5 show the three-dimensional Young's modulus E surface of each intermetallic compound Mg 2 X (X = Si, Ge) at 0 GPa with 1 × 1 × 2 supercells and 2 × 2 × 2 supercells, respectively. The near sphere shows that these intermetallic compounds are elastic and anisotropic [48]. Compared with other point defects of Mg 2 X, the Mg 2 X with X vacancy exhibits a relatively obvious deviation from spherical shape, which indicates that the Mg 2 X with X vacancy has a relatively strong anisotropy. The degree of spherical deviation of Mg 2 X decreases with the increase of supercell size, implying that Mg 2 X with point defects show a relatively weak anisotropy in the 2 × 2 × 2 supercells. Anisotropy exists widely in all kinds of materials. In short, crystals have different properties in different directions. It is easy to confuse anisotropy and inhomogeneity, which are two completely different concepts. Anisotropy means that the properties of materials are related to directions, while inhomogeneity means that the properties of materials are related to parts. In the case of single crystals, at any point in the interior, the structures and properties are the same, but it has different performance in different directions.
Crystal is anisotropic, with different properties in different directions, and has strict symmetry. Here, the symmetry is a very remarkable property, which is widely used in mathematics and physics, and also commonly used in the analysis of material properties. Polycrystalline materials have preferred orientation and certain anisotropy.
All kinds of materials are elastic, and the elastic properties of most materials are anisotropic [45,46] where 11 , 12 , and 44 represent the flexibility coefficient, 1 , 2 , and 3 indicate direction cosines relative to the x, y, and z directions, respectively. Figure 4 and Figure 5 show the three-dimensional Young's modulus surface of each intermetallic compound Mg2X (X = Si, Ge) at 0 GPa with 1 × 1 × 2 supercells and 2 × 2 × 2 supercells, respectively. The near sphere shows that these intermetallic compounds are elastic and anisotropic [48]. Compared with other point defects of Mg2X, the Mg2X with X vacancy exhibits a relatively obvious deviation from spherical shape, which indicates that the Mg2X with X vacancy has a relatively strong anisotropy. The degree of spherical deviation of Mg2X decreases with the increase of supercell size, implying that Mg2X with point defects show a relatively weak anisotropy in the 2 × 2 × 2 supercells.

Conclusions
In this paper, based on DFT, the influence of point defects on the structure and elastic mechanical stability of Mg2X (X = Si, Ge) intermetallic is revealed from the microscopic point of view: (1) The difference of atomic radius, the lattice constant of Mg2Ge type intermetallics with (without) point defects is larger than that of Mg2Si type, and the existence of point defects leads to the instability of crystal cells and the expansion of lattice distortion.
(2) The predicted ∆Hf and Ecoh are consistently negative for all Mg2X (X = Si, Ge) with point defects, implying the thermodynamic stability of Mg2X (X = Si, Ge) with point defects. The absolute value of ∆ of Mg2X (X = Si, Ge) with point defects is smaller than that of the cells without point defects, which indicates that the structure with point defects is relatively difficult to form; intermetallics are easier to form vacancy defects than anti-site on the X (X = Si, Ge) lattice site and form anti-site defects on the Mg lattice site; the ∆Hf of the vacancy or anti-site defects of Mg2X at the Si or Ge position is bigger than that at the Mg position, which indicates that the point defect is more likely to appear at the Mg position than at the Si or Ge position; the absolute value of the binding energy of the intermetallics without point defects is higher than that of the intermetallics with point defects, which indicates that the energy required for atom bonding is higher and the structure is more stable when there is no defect; the structure of the intermetallics forming point defects at the Mg position is more stable than that at the Si/Ge position.
(3) When the point defects exist, the bulk modulus , shear modulus and Young's modulus of intermetallic compounds become smaller, which indicates that the void and anti-void point defects will reduce the volume deformation resistance, shear strain resistance, and positive elastic deformation resistance of materials, and increase the mechanical instability of elastic deformation of materials. When there is a vacancy or anti-site defect in the Mg position, the resistance to volume deformation, shear strain and elastic deformation of the crystal is stronger than that in the Si or Ge position. Therefore, the point defects in the Mg position enhance the mechanical stability of the material, while the plasticity is poor. The maximum / value of Mg2X without point defects indicates that the formation of defects will improve the ductility of materials; the / value of vacancy point defects is smaller than that of anti-site defects, indicating that vacancy point defects can improve the plasticity of materials more than anti-site defects.
(4) The Poisson's ratio of Mg2X without point defects is the smallest, which indicates that the formation of defects makes the transverse deformation of the material larger, and the transverse elastic instability increased. Compared with 2 × 2 × 2 supercells, the G/B of Mg2X with the same point defect structure (except ) in 1 × 1 × 2 supercells can enhance the plasticity of materials.

Conclusions
In this paper, based on DFT, the influence of point defects on the structure and elastic mechanical stability of Mg 2 X (X = Si, Ge) intermetallic is revealed from the microscopic point of view: (1) The difference of atomic radius, the lattice constant of Mg 2 Ge type intermetallics with (without) point defects is larger than that of Mg 2 Si type, and the existence of point defects leads to the instability of crystal cells and the expansion of lattice distortion.
(2) The predicted ∆H f and E coh are consistently negative for all Mg 2 X (X = Si, Ge) with point defects, implying the thermodynamic stability of Mg 2 X (X = Si, Ge) with point defects. The absolute value of ∆H f of Mg 2 X (X = Si, Ge) with point defects is smaller than that of the cells without point defects, which indicates that the structure with point defects is relatively difficult to form; intermetallics are easier to form vacancy defects than anti-site on the X (X = Si, Ge) lattice site and form anti-site defects on the Mg lattice site; the ∆H f of the vacancy or anti-site defects of Mg 2 X at the Si or Ge position is bigger than that at the Mg position, which indicates that the point defect is more likely to appear at the Mg position than at the Si or Ge position; the absolute value of the binding energy of the intermetallics without point defects is higher than that of the intermetallics with point defects, which indicates that the energy required for atom bonding is higher and the structure is more stable when there is no defect; the structure of the intermetallics forming point defects at the Mg position is more stable than that at the Si/Ge position.
(3) When the point defects exist, the bulk modulus B, shear modulus G and Young's modulus E of intermetallic compounds become smaller, which indicates that the void and anti-void point defects will reduce the volume deformation resistance, shear strain resistance, and positive elastic deformation resistance of materials, and increase the mechanical instability of elastic deformation of materials. When there is a vacancy or anti-site defect in the Mg position, the resistance to volume deformation, shear strain and elastic deformation of the crystal is stronger than that in the Si or Ge position. Therefore, the point defects in the Mg position enhance the mechanical stability of the material, while the plasticity is poor. The maximum G/B value of Mg 2 X without point defects indicates that the formation of defects will improve the ductility of materials; the G/ B value of vacancy point defects is smaller than that of anti-site defects, indicating that vacancy point defects can improve the plasticity of materials more than anti-site defects.
(4) The Poisson's ratio ν of Mg 2 X without point defects is the smallest, which indicates that the formation of defects makes the transverse deformation of the material larger, and the transverse elastic instability increased. Compared with 2 × 2 × 2 supercells, the G/B of Mg 2 X with the same point defect structure (except Mg Si ) in 1 × 1 × 2 supercells can enhance the plasticity of materials.