Mechanical, elastic, anisotropy, and electronic properties of monoclinic phase of m-SixGe3−xN4
Ma Zhen-Yang, Yan Fang, Wang Su-Xin, Jia Qiong-Qiong, Yu Xin-Hai, Shi Chun-Lei
Tianjin Key Laboratory for Civil Aircraft Airworthiness and Maintenance, Civil Aviation University of China, Tianjin 300300, China

 

† Corresponding author. E-mail: zyma@cauc.edu.cn

Project supported by the National Natural Science Foundation of China (Grant No. 61601468), the Fundamental Research Funds for the Central Universities, China (Grant No. 3122014C024), and the Fund for Scholars of Civil Aviation University of China (Grant No. 2013QD06X).

Abstract

The structural, mechanical, elastic anisotropic, and electronic properties of the monoclinic phase of m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4 are systematically investigated in this work. The calculated results of lattice parameters, elastic constants and elastic moduli of m-Si3N4 and m-Ge3N4 are in good agreement with previous theoretical results. Using the Voigt–Reuss–Hill method, elastic properties such as bulk modulus B and shear modulus G are investigated. The calculated ratio of B/G and Poisson’s ratio v show that only m-SiGe2N4 should belong to a ductile material in nature. In addition, m-SiGe2N4 possesses the largest anisotropic shear modulus, Young’s modulus, Poisson’s ratio, and percentage of elastic anisotropies for bulk modulus and shear modulus , and universal anisotropic index among m-SixGe3−xN4 (x = 0, 1, 2, 3.) The results of electronic band gap reveal that m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4 are all direct and wide band gap semiconducting materials.

1. Introduction

Nitrides consisting of III A and IV A group elements can be used as high temperature materials, catalysts, light-emitting diodes and refractory ceramics because of their excellent mechanical and thermodynamic properties.[14] Silicon nitride is an important ceramic material. It is a superhard material and atomic crystal which has good lubricity and abrasion resistance. It resists oxidation at high temperature, and it can resist cold and hot shocks, it can be heated to over 1000 K in the air and will not be broken by rapid cooling and rapid heating. Because silicon nitride ceramic has such excellent characteristics, it is often used in bearings, gas turbine blades, mechanical seals, permanent die, and other mechanical components. If the silicon nitride ceramic with high temperature resistance and high heat transfer resistance is used to make the heating surfaces of the engine parts, it can not only improve the quality of diesel engine, save fuel, but also enhance the thermal efficiency. In fact, germanium is adopted because Ge atoms have higher carrier mobility than Si atoms. The semiconductor field effect transistor has superior performances. The germanium atom radius is bigger, so the critical pressure of the phase transition of Ge3N4 is less than that of Si3N4. Germanium nitride also has the advantages of corrosion resistance, high hardness, and its band gap can be adjusted. Many new structures of silicon nitride and germanium nitride have been reported in many researches.[58] It has been known for several decades that silicon nitride exists in two stable phases (namely α- and β-Si3N4) with very similar structures,[5] they both have hexagonal symmetry. The phase boundaries of the transitions are investigated using first-principles calculations within the quasi-harmonic approximation by Yu and Chen.[6] Using the plane-wave pseudopotential method, Yu and Chen[7] have investigated the structural properties, mechanical properties and thermodynamic properties of β-, γ-, wII-, and post-spinel Si3N4 and also the phase transition characters among these polymorphs. Wang et al. have conducted a comparative structural and electronic study of a new class of material, olivine-Si3N4 and olivine-Ge3N4 by using the first principles approach.[8] They found that olivine-Si3N4 has a large bulk modulus (262.9 GPa), and olivine-Si3N4 and olivine-Ge3N4 are both direct bang gap semiconductors with the band gaps of 2.99 eV and 1.55 eV, respectively.

The structural, electronic, and optical properties of the new cubic spinel nitrides c-Si3N4, c-SiGe2N4, c-Si2GeN4, and c-Ge3N4 were studied by a first-principles method.[9] The c-SiGe2N4 is a direct band gap semiconductor with a band gap of 1.85 eV, and the c-Si2GeN4 is an indirect band gap semiconductor with a band gap of 2.56 eV. In addition, the bulk moduli of c-SiGe2N4 and c-Si2GeN4 are 277.1 GPa and 258.3 GPa, respectively. Using first-principles calculations, the structural, electronic, elastic, and optical properties of cubic spinel SiGe2N4 under pressure were investigated by Moakafi et al.[10] Bouhemadou et al.[11] have investigated the structural, elastic, electronic, optical and thermal properties of c-SiGe2N4 by using the ultrasoft pseudopotential density functional method in the generalized gradient approximation under pressure. The structural and mechanical properties of the cubic spinel c-CSi2N4 were studied by first-principles total energy calculations based on the density-functional theory.[12] The Vickers hardness of c-CSi2N4 is 52.07 GPa, so it is a potential superhard material. Cui et al. explored three Si3N4 metastable phases, i.e., tetragonal t-Si3N4, monoclinic m-Si3N4 and orthorhombic o-Si3N4, by analyzing the crystal structure obtained through using the particle swarm optimization (CALYPSO) code.[13] The m-Si3N4 belongs to the space group of Cm in a monoclinic system and it is mechanically and dynamically stable under ambient pressure based on their elastic constants and phonon dispersions calculations. Cang et al. predicted the lattice structures, densities of states, phonon dispersion curves of m-, t-, and o-Ge3N4.[14] The stabilities of these Ge3N4 structures were verified by the proving of negative formation enthalpy, the satisfying of Born’s stability criteria and no imaginary frequency in the phonon dispersion curves in a pressure range of 0 GPa–20 GPa. Fan et al. investigated the electronic structures, elastic, anisotropic, and electronic properties of m-Si3N4, t-Si3N4, and o-Si3N4 under pressure,[15] and it is interesting that the shift of band gap from direct to indirect occurs for m-Si3N4 when the pressure rises up to 50 GPa. Chen et al.[16] predicted the structural, electronic, optical and thermodynamic properties of m-, t-, and o-M3N4 (M = Si, Ge, Sn).

According to the previous work, we predict two new double nitrides m-Si2GeN4 and m-SiGe2N4 in the Cm space group. In the present work, the structural, elastic anisotropic, and electronic properties of m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4 are investigated using density functional theory calculations. The Debye temperature, the ductile or brittle characters and elastic anisotropy of m-SixGe3−xN4 (x = 0, 1, 2, 3.) are discussed in detail.

2. Computational details and theory

The structural optimizations and property predictions of the monoclinic structures of m-SixGe3−xN4 (x = 0, 1, 2, 3) were performed based on density functional theory (DFT)[17,18] as implemented in the Cambridge Serial Total Energy Package (CASTEP) code.[19] The exchange–correlation potential was treated with the generalized gradient approximation (GGA) in the form of Perdew–Burke–Ernzerhof (PBE),[20] PBE for solid (PBEsol),[21] and the local density approximation (LDA) in the form of Ceperley and Alder data parameterized by Perdew and Zunger (CA-PZ).[22,23] The Broyden–Fletcher–Goldfarb–Shanno (BFGS)[24] minimization scheme was used in geometry optimization. Vanderbilt type ultrasoft pseudopotentials[25] were employed to describe the interactions between ionic core and valence electrons. The considered valence atomic configurations had been set to be Si-3s23p2, Ge-4s24p2, and N-2s22p3. In the structure calculation, a plane-wave basis set with energy cut-off 500 eV was used for all cases. The k-points of 6×14×7, 6×13×7, 6×13×7, and 5×13×7 were adopted for the structures of m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4 by using Monkhorst–Pack mesh,[26] respectively. The system reached the ground state via self-consistent calculation when the total energy, the maximum displacement of ion, the maximum ionic Hellmann–Feynman force within and the maximum stress were less than 5×10−6 eV/atom, 5×10−4 Å, 0.01 eV/Å, and 0.02 GPa, respectively. In order to calculate electronic structures, the HSE06 hybrid functional[27] was used.

3. Results and discussion

The crystal structures of m-Si3N4 and the double nitrides m-Si2GeN4 and m-SiGe2N4 are shown in Fig. 1. In a conventional cell, m-Si3N4 (m-Ge3N4) has 8 N atoms and 6 Si (Ge) atoms. For m-Si2GeN4 and m-SiGe2N4, the Ge atoms replace Si atoms with the smallest energy. The equilibrium lattice constants are determined by the minimum total energy with respect to the variation of the cell volume through utilizing both GGA and LDA. The optimized lattice parameters of m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4 at zero pressure are listed in Table 1.

Fig. 1. (color online) Crystal structures of (a) m-Si3N4, (b) m-Si2GeN4, and (c) m-Ge3N4 under ambient pressure. The silicon, germanium, and nitrogen atoms are represented as black, red, and blue spheres, respectively.
Table 1.

Calculated lattice parameters (a, b, c in unit Å and β in unit degree) of m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4 compared with previous theoretical results.

.

Along with previous theoretical results for comparison, it can be seen that the calculated lattice parameters of m-Si3N4 and m-Ge3N4 are in excellent agreement with previous results.[13,15,16] The variations of the lattice parameters (a, b, c) with increasing the composition of germanium are shown in Fig. 2. On the whole, the lattice parameters (a, b, c) increase with increasing the composition of germanium. The augment tendency of lattice parameter b is nearly linear while the lattice parameter a and c increases nonlinearly with increasing the composition of germanium. The curves of the PBE and PBEsol methods almost coincide as shown in Fig. 2. From m-Si3N4 to m-Ge3N4, the lattice parameters a, b, and c increase by 6.34% (5.10%), 5.87% (4.61%), and 6.90% (5.56%) obtained from the PBE (CA-PZ) method, respectively.

Fig. 2. (color online) Calculated lattice parameters ((a) a, (b) b, and (c) c) of m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4 by PBE, PBEsol, and CA-PZ methods.

The elastic constants are closely related to various fundamental mechanical properties of a crystal, which can provide, particularly, the information about the stability and stiffness of materials. For a monoclinic crystal, there are thirteen independent elastic constants, namely C11, C22, C33, C44, C55, C66, C12, C13, C15, C23, C25, C35, and C46. The calculated elastic constants of m-SixGe3−xN4 (x = 0, 1, 2, 3) are listed in Table 2. The elastic constants of m-Si3N4 and m-Ge3N4 are well consistent with pervious theoretical values.[13,14] For a stable crystal of monoclinic structure, it must obey the following mechanical stability according to the Born’s criteria:[28]

Table 2.

Calculated elastic constants (in unit GPa) of m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4 compared with other calculated results.

.

It can be easily verified that the elastic constants of m-SixGe3−xN4 (x = 0, 1, 2, 3) satisfy their stability requirements, which indicates that they are mechanical stable. And the dynamic stabilities of m-Si3N4 and m-Ge3N4 have been verified by the phonon spectrum calculations in Refs. [13] and [14]. So it is revealed that m-Si3N4 and m-Ge3N4 are mechanically and dynamically stable via the elastic constants and phonon calculations. The elastic constants C11, C22, C33 represent the stiffness of the material when applying stresses along the x, y, z directions,[29] respectively. The obtained result obviously shows that all the m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4 follow this ordering relationship: , which means that the resistance against deformation along the z axis is stronger than along the x axis but weaker than along the y axis. Simultaneously, with increasing the components of Ge atoms, the elastic constants Cii decrease, and thus the compression property increases.

Furthermore, according to the Voigt–Reuss–Hill (VRH) average scheme,[30,31] the bulk modulus B and shear modulus G can be estimated using calculated elastic constants via the following equations:

The Young’s modulus E and Poisson’s ratio v are obtained from the following equations: and , respectively. Comparing with the previous results,[13,14,16] the calculated values of bulk modulus B, shear modulus G, and Young’s modulus E are presented in Table 3. The bulk modulus B represents the resistance to fracture while the shear modulus G represents the resistance to plastic deformation.[32] The values of shear modulus are much larger than those of bulk modulus for m-SixGe3−xN4 (x = 0, 1, 2, 3), which means that the resistance to shape change is better than that to volume change. Generally, Young’s modulus E provides a measure of material stiffness. The larger the value of E, the stiffer the material will be.[33] The values of B, G, and E decrease with the increase of the Ge atoms except m-SiGe2N4, which suggests that the resistance to deformation becomes weaker and the stiffness decreases.

Table 3.

Calculated values of elastic modulus (in unit GPa), Poisson’s ratio v, percent compressibility of bulk modulus and shear modulus factors , , and anisotropic index for m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4, compared with other calculated results.

.

Poisson’s ratio v can be used as an indicator for ductility or brittleness. According to the Frantsevich’s rule,[34] the critical value that distinguishes between the ductile and brittle nature of material is 0.26. If , the material will behave in a ductile manner, otherwise the material demonstrates brittleness in nature. Additionally, the ratio of bulk modules to shear modulus (B/G) can also judge whether a material would be a ductile or brittle manner according to the Pugh’s criterion.[35] The critical value of this ratio which distinguishes between the ductility and brittle is 1.75. For ductile materials, the ratio of B/G is larger than 1.75. Here the calculation Poisson’s ratio v and ratio of B/G of m-SixGe3−xN4 (x = 0, 1, 2, 3) are listed in Table 3. Only the Poisson’s ratio v (B/G) of m-SiGe2N4 is larger than 0.26 (1.75), meaning that m-SiGe2N4 is ductile and m-Si3N4, m-Si2GeN4, and m-Ge3N4 are brittle in nature. The brittle materials are sensitive to thermal shocks, as the material cannot efficiently dissipate thermal stresses via plastic deformations. Thus, a brittle material can suffer limited thermal shocks before its strength drops dramatically. Hence, the capability of resisting thermal shocks for m-SiGe2N4 is stronger than for m-Si3N4, m-Si2GeN4, and m-Ge3N4.

The Debye temperature ( is one of fundamental parameters for solid materials, which is related to many physical properties, such as thermal expansion and specific heat. It is the Debye temperature that reflects the temperature of a crystal’s normal mode of vibration. The Debye temperature can be estimated from averaged sound velocity ( ), compressional sound velocity ( and shear sound velocity ( through the following formulas:[36,37] where h is the Planck’s constant, is the Boltzmann constant, is the Avogadro’s constant, M is the molecular weight, n is the number of atoms in the molecule, and ρ is the density of crystal. The calculated results of density, the Debye temperature and the sound velocities for m-SixGe3−xN4 (x = 0, 1, 2, 3) are listed in Table 4 at 0 GPa and 0 K. Following the sequence , the densities increase whereas the Debye temperature and the sound velocities decrease. According to Refs. [38]–[40], a larger suggests a higher normal vibration, which is associated with a better thermal conductivity. Meanwhile, the Debye temperature can characterize the strength of covalent bond for solid. The larger the , the stronger the covalent bond is.[41] Therefore, thermal conductivities and covalent bonds of the m-SixGe3−xN4 (x = 0, 1, 2, 3) become weaker following the above sequence.

Table 4.

Values of density (in units g/cm3), compressional, shear, and average sound velocities (in m/s) as well as Debye temperature (in K) for m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4..

.

The elastic anisotropy analysis of a crystal is of great significance for understanding the mechanisms of material microcrack and durability in engineering science.[42] There are a variety of methods to quantity the anisotropy of a crystal structure, such as the universal anisotropic index ,[43] the percent compressibility of bulk and shear modulus factors ( and ).[44] According to the following equations, these mechanical anisotropic indexes mentioned above can be evaluated from

For isotropic crystals, the percentages of elastic anisotropy for bulk modulus and shear modulus , and universal anisotropic index must be zero. The deviation from zero measures the anisotropy degree of crystal. The calculated anisotropy factors ( , , and are presented in Table 3. It is obvious that the m-SixGe3 −xN4 (x = 0, 1, 2, 3) are anisotropic. The value of is much larger than for each crystal, indicating that the anisotropy of shear modulus is more conspicuous than that of bulk modulus. The values of the universal anisotropic index are 0.968 for m-Si3N4, 0.905 for m-Si2GeN4, 2.557 for m-SiGe2N4, and 1.513 for m-Ge3N4, respectively. The anisotropy factors ( , , and of m-SiGe2N4 are maxima, which indicates the anisotropy of m-SiGe2N4 is the largest among the m-SixGe3−xN4 (x = 0, 1, 2, 3).

Simultaneously, the elastic anisotropic m-SixGe3 −xN4 (x = 0, 1, 2, 3) in Poisson’s ratio, shear modulus and Young’s modulus were investigated in the present work. The three-dimensional (3D) directional dependence and two-dimensional (2D) representation of anisotropy are calculated by utilizing the program of Elastic Anisotropy Measure (ELAM).[45] The vector direction (θ, ϕ) in the spherical coordinate system is used to describe the directional dependence of anisotropy, where θ (ϕ) represents the angle between the vector and the x-axis (z-axis) positive direction and they are expressed in radians. For an isotropic material, the 3D directional dependence would exhibit a spherical shape, while the deviation degree from the spherical shape reflects the content of anisotropy.[46]

The calculated 3D directional dependence of Poisson’s ratio and 2D representation in the (001), (010), and (100) planes for m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4 is shown in Figs. 3 and 4, respectively. It is evident that the 3D figures of the Poisson’s ratio have a large degree of deviation in shape from the sphere, which implies that the Poisson’s ratios of m-SixGe3−xN4 (x = 0, 1, 2, 3) are greatly anisotropic. The maximal (minimal) values of the Poisson’s ratio for m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4 are 0.55 (0.02), 0.55 (0.03), 0.77 (0.02), and 0.67 (0.02) at 0 GPa. Following the sequence , the maximal and minimal Poisson’s ratio directions are (3.14, 0.19), (3.14, 0.17), (0, 2.95), (3.14, 0.17) and (3.14, 0.89), (2.28, 0), (0, 2.38), (3.68, 1.35), respectively. It can be found that the maximal and minimal Poisson’s ratio directions for m-Si3N4 and m-SiGe2N4 are both in the (010) plane. From Fig. 4, it can be seen that four compounds present the different degrees of anisotropy in Poisson’s ratio in different planes. Following the above sequence, the ratios of are 0.44/0.02 = 22, 0.43/0.04 = 10.75, 0.63/0.02 = 31.5, and 0.52/0.04 = 13 in the (001) plane, 0.55/0.02 = 27.5, 0.55/0.03 = 18.33, 0.77/0.02 = 35, and 0.67/0.02 = 33.5 in the (010) plane, and 0.52/0.02 = 26, 0.53/0.04 = 13.25, 0.48/0.03 = 16.33, and 0.64/0.04 = 16 in the (100) plane, respectively. Obviously, the anisotropy of the Poisson’s ratio for m-SiGe2N4 is the largest one among m-SixGe3−xN4 (x = 0, 1, 2, 3).

Fig. 3. (color online) 3D directional dependence of Poisson’s ratio for (a) m-Si3N4, (b) m-Si2GeN4, (c) m-SiGe2N4, and (d) m-Ge3N4 at ambient pressure.
Fig. 4. (color online) 2D representations of Poisson’s ratio in the (001), (010), and (100) planes for m-SixGe3−xN4 (x = 0, 1, 2, 3). The dash–dotted and solid lines represent the maximum and minimum values in the (001), (010), and (100) planes, respectively. The black, blue, green, and red lines denote the Poisson’s ratios of m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4, respectively.

The shear moduli for all possible directions of shear strain are also calculated. The calculated 3D directional dependence of shear modulus and 2D representations in the (001), (010), and (100) planes for the m-SixGe3−xN4 (x = 0, 1, 2, 3) are plotted in Figs. 5 and 6, respectively. The 3D figures reveal that the anisotropy of shear modulus is weak since the deviation degree from the sphere shape is small. The values of ratio are as follows: for m-Si3N4, for m-Si2GeN4, for m-SiGe2N4, for m-Ge3N4 at 0 GPa. From Fig. 6, it is clear that the anisotropies of shear modulus in the (010) plane for the m-SixGe3 −xN4 (x = 0, 1, 2, 3) are slightly larger than those in other planes. From m-Si3N4 to m-Ge3N4, the values of are 150.31/63.24 = 2.38, 129.87/57.71 = 2.25, 118.83/29.98 = 3.96, and 100.23/36.25 = 2.76 in the (001) plane, 139.05/63.24 = 2.2, 117.46/57.71 = 2.04, 120.27/29.98 = 4.01, and 97.68/36.25 = 2.69 in the (010) plane, and 156.54/63.24 = 2.48, 139.54/57.71 = 2.42, 139.86/29.98 = 4.6, and 119.49/36.25 = 3.29 in the (100) plane, respectively. Consequently, m-SiGe2N4 exhibits larger anisotropy among m-SixGe3−xN4 (x = 0, 1, 2, 3).

Fig. 5. (color online) 3D directional dependence of shear modulus for (a) m-Si3N4, (b) m-Si2GeN4, (c) m-SiGe2N4, and (d) m-Ge3N4 (d) at ambient pressure.
Fig. 6. (color online) 2D representations of shear modulus in the (001), (010) and (100) planes for m-SixGe3−xN4 (x = 0, 1, 2, 3). The dash dotted and solid lines represent the maximum and minimum values in the (001), (010), and (100) planes, respectively. The black, blue, green and red lines represent the Poisson’s ratios of m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4, respectively.

The 3D directional dependencies of Young’s modulus for m-SixGe3−xN4 (x = 0, 1, 2, 3.) are plotted in Fig. 7. It can be seen that the 3D figures of the m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4 have more or less deviations in shape from sphere along the x, y, z axes, which indicates that the Young’s moduli of these compounds possess some kind of degree anisotropy. The values of ratio are as follows: for m-Si3N4, for m-Si2GeN4, for m-SiGe2N4, and for m-Ge3N4 at 0 GPa. The calculated Young’s moduli along different directions as well as in different planes for the m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4 are shown in Fig. 8. At first glance, the anisotropy of Young’s modulus in (001) plane for m-SixGe3−xN4 ( , 1, 2, 3) is similar to that in (100) plane but slightly larger than that in (010) planes. The values of for m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4 are 440.27/193.73 = 2.28, 404.44/181.34 = 2.23, 418.92/106.91 = 3.92, and 334.56/123.52 = 2.71 in the (001) plane, 344.61/190.17 = 1.81, 304.96/178.67 = 1.71, 331.71/135.61 = 2.45, and 267.23/120 = 2.23 in the (010) plane, and 440.27/227.36 = 1.94, 404.44/194.14 = 2.08, 418.92/153.44 = 2.73, and 334.56/139.31 = 2.4 in the (100) plane, respectively. It can be noted that the maximum values of Young’s modulus for m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4 in the (001) plane are the same as those in the (100) plane. A comparison of the values of leads to a conclusion that the anisotropy of Young’s modulus for m-SiGe2N4 is stronger among those of m-SixGe3 −xN4 (x = 0, 1, 2, 3), which is consistent with the pervious conclusion.

Fig. 7. (color online) 3D directional dependence of Young’s modulus for (a) m-Si3N4, (b) m-Si2GeN4, (c) m-SiGe2N4, and (d) m-Ge3N4 at ambient pressure.
Fig. 8. (color online) 2D representations of Young’s modulus in the (001), (010), and (100) planes for m-SixGe3−xN4 (x = 0, 1, 2, 3). The black, red, blue, and magenta lines represent the Poisson’s ratios of m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4, respectively.

Physical properties can be measured in different directions. Various crystal properties, including the elastic modulus, hardness, fracture resistance, thermal expansion coefficient, thermal conductivity, resistivity, and electric displacement vector, have different levels of anisotropy. If a physical property is closely related to orientation, it is anisotropic. Why does this reverse phenomenon appear? This is because of the different atomic arrangements and different atomic positions, and different atoms have different lengths in different directions, resulting in the fact that the physical quantities along the different directions possess different values. For m-SixGe3−xN4 (x = 0, 1, 2, 3), they have a monoclinic symmetry that belongs to the Cm space group, the symmetry of the monoclinic system is very low. Young’s modulus in Figs. 8(a) and 8(c) are symmetric with respect to the xy plane and yz plane, respectively, but in Fig. 8(b), the symmetry axis are not along the xz plane. And both the Poisson’s ratio, and the shear modulus have the same situation as Young’s modulus. That is to say, the symmetries of Young’s modulus, Poisson’s ratio, and shear modulus in the xz plane are lower than in the xy plane and yz plane.

The electronic band structure reflects the law of the electron motion in a solid and determines the fundamental physical and chemical properties of material. The band structures of m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4 along some high symmetry directions of the Brillouin zone are calculated with HSE06 functional at ambient pressure, and the results are shown in Fig. 9. The horizontal red dash line designates zero energy representing the Fermi lever ( ). At first glance, it appears that the four band structures are very similar. The major differences are in the width of band gap and conduction band minimum. From Fig. 9, it is clear that the valance band maximum (VBM) and the conduction band minimum (CBM) of m-SixGe3−xN4 (x = 0, 1, 2, 3) are both located at G point. Therefore, m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4 are direct semiconductors with band gaps of 5.08 eV, 4.76 eV, 4.81 eV, and 3.34 eV, respectively. It should be noted that no experimental dataare available to verify the band gaps of these compounds. All the previous results are obtained via theoretical methods, which need to be verified experimentally in future.

Fig. 9. (color online) Electronic band structures of (a) m-Si3N4, (b) m-Si2GeN4, (c) m-SiGe2N4, and (d) m-Ge3N4 with the HSE06 hybrid functional.

To further investigate the natures of the electronic band structure for m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4, the partial densities of states (PDOSs) are calculated as shown in Fig. 10. It is obvious that PDOSs of m-SixGe3−xN4 (x = 0, 1, 2, 3) look similar on a gross scale since they have the same crystal structure and Si and Ge are isoelectronic elements. The densities of electronic states are separated into three regions by the gaps. The first region is in the energy range between −20 eV and −15 eV, the second region is between −10 eV and Fermi energy ( , and the last region ranges from 5 eV to 11 eV in the conduction band. For m-Si3N4 (m-Ge3N4), the first region is mainly for N-s states with admixture from Si (Ge)-s/p states, the second region mainly involves N-p states with significant contributions from Si (Ge)-s states (between 10 eV to −7.5 eV) and Si (Ge)-p states (from −7.5 eV to , and the third region mainly contains Si (Ge)-p states and a mixture of Si (Ge)-s states and N-s/p states. For m-Si2GeN4 and m-SiGe2N4, the first region mainly contains N-s states with a mixture of Si- and Ge-s/p states, the second region mainly involes N-p states with significant contributions from Si- and Ge-s/p states, and the third region mainly covers Si/Ge-p states and a mixture of Si (Ge)-s states and N-s/p states.

Fig. 10. (color online) Comparisons of total density of states and partial density of states for m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4. The black/dash and red/solid curves represent s states and p states of silicon, germanium, and nitrogen atoms, respectively.
4. Conclusions

The structural, electronic, elastic, and anisotropic properties of the monoclinic phases of m-Si3N4, m-Si2GeN4, m-SiGe2N4, and m-Ge3N4 are investigated systemically in this work. The calculated lattice constants, elastic constants and elastic moduli of m-Si3N4 and m-Ge3N4 are in good consistence with previous theoretical results. The elastic constants of m-SixGe3 −xN4 (x = 1, 2) satisfy the stability criteria, which indicates that they are mechanically stable. According to the analyses of the obtained values of Poisson’s ratio v and B/G, it is found that only m-SiGe2N4 is ductile in nature. The densities of m-SixGe3−xN4 (x = 0, 1, 2, 3) increase whereas the Debye temperature and the sound velocities decrease with increasing the composition of germanium. The detailed analyses of the anisotropy manifest that the m-SixGe3−xN4 (x = 0, 1, 2, 3) exhibit lager anisotropy and the anisotropy of m-SiGe2N4 is the largest among the m-SixGe3−xN4 (x = 0, 1, 2, 3). The calculated band structures show that they are all direct semiconductors with band gaps of 5.08 eV, 4.76 eV, 4.81 eV, and 3.34 eV, respectively.

Acknowledgment

The authors thank Q Y Fan (School of Microelectronics, Xidian University) for being allowed to use the CASTEP code in Materials Studio.

Reference
[1] Jack K H 1976 J. Mater. Sci. 11 1135
[2] Clarke D R Lange F F Schnittgrund G D 1982 J. Am. Ceram. Soc. 65 c51
[3] Riley F L 2000 J. Am. Ceram. Soc. 83 245
[4] Xu B Dong J J McMillan P F Shebanova O Salamat A 2011 Phys. Rev. 84 014113
[5] Ching W Y Mo S D Ouyang L Z Rulis P Tanaka I Yoshiya M 2002 J. Am. Ceram. Soc. 85 75
[6] Yu B H Chen D 2013 J. Alloys Compd. 581 747
[7] Yu B H Chen D 2012 Physica 407 4660
[8] Wang H Chen Y Kaneta Y Iwata S 2006 J. Phys.: Condens. Mater. 18 10663
[9] Ching W Y Mo S D Ouyang L 2001 Phys. Rev. 63 245110
[10] Moakafi M Khenata R Bouhemadou A Benkhettou N Rached D Reshak A Elastic H 2009 Phys. Lett. 373 2393
[11] Bouhemadou A Al-Douri Y Khenata R Haddadi K 2009 Eur. Phys. J. 71 185
[12] Zhang X Y Chen Z W Du H J Ma M Z He J L Tian Y J Liu R P 2008 J. Appl. Phys. 103 083533
[13] Cui L Hu M Wang Q Q Xu B Yu D L Liu Z Y He J L 2015 J. Solid State Chem. 228 20
[14] Cang Y P Chen D Yang F Yang H M 2016 Chem. J. Chin. Univer. 37 674
[15] Fan Q Y Chai C C Wei Q Zhou P K Yang Y T 2016 AIP Adv. 6 085207
[16] Chen D Cheng K Qi B Y 2017 Chin. Phys. 26 046303
[17] Hohenberg P Kohn W 1964 Phys. Rev. 136 B864
[18] Kohn W Sham L J 1965 Phy. Rev. 137 1697
[19] Clark S J Segall M D Pickard C J Hasnip P J Probert M I J Refson K Payne M C 2005 Z. Krist.-Cryst. Mater. 220 567
[20] Perdew J P Burke K Ernzerhof M 1996 Phys. Rev. Lett. 77 3865
[21] Perdew J P Ruzsinszky A Csonka G I Vydrov O A Scuseria G E Constantin L A Zhou X L Burke K 2008 Phys. Rev. Lett. 100 136406
[22] Ceperley D M Alder B J 1980 Phys. Rev. Lett. 45 566
[23] Perdew J P Mcmullen E R Zunger A 1981 Phys. Rev. 23 2785
[24] Pfrommer B G Côté M Louie S G Cohen M L 1997 J. Comput. Phys. 131 233
[25] Vanderbilt D 1990 Phys. Rev. 41 7892
[26] Monkhorst H J Pack J D 1976 Phys. Rev. 13 5188
[27] Krukau A V Vydrov O A Izmaylov A F Scuseria G E 2006 J. Chem. Phys. 125 224106
[28] Nye J F 1985 Physical properties of crystals New York Oxford university press
[29] Hao X P Cui H L 2014 J. Korean. Phys. Soc. 65 45
[30] Wu Z J Zhao E J Xiang H P Hao X F Liu X J Meng J 2007 Phys. Rev. 76 054115
[31] Reuss A 1929 Z. Angew. Math. Mech. 9 49
[32] Vaitheeswaran G Kanchana V Kumar R S Cornelius A L Nicol M F Savane A Delin A Johansson B 2007 Phys. Rev. 76 014107
[33] Liu Y Hu W C Li D J Zeng X Q Xu C S 2013 Phys. Scr. 88 045302
[34] Duan Y H Sun Y Peng M J Zhou S G 2014 J. Alloys Compd. 595 14
[35] Pugh S F 1954 Philos. Mag. 45 823
[36] Anderson O L 1963 J. Phys. Chem. Solids 24 909
[37] Fan Q Y Chai C C Wei Q Yan H Y Zhao Y B Yang Y T Yu X H Liu Y Xing M J Zhang J Q Yao R H 2015 J. Appl. Phys. 118 185704
[38] Duan Y H Huang B Sun Y Peng M J Zhou S G 2014 J. Alloys Compd. 590 50
[39] Hu W C Liu Y Li D J Zeng X Q Xu C S 2014 Comput. Mater Sci. 83 27
[40] Li J Zhang M Luo X 2013 J. Alloys Compd. 556 214
[41] Pan Y Zheng W T Guan W M Zhang K H Fan X F 2013 J. Solid State Chem. 207 29
[42] Duan Y H Sun Y Peng M J Zhou S G 2014 J. Alloys Compd. 595 14
[43] Ranganathan S I Ostoja-Starzewski M 2008 Phys. Rev. Lett. 101 055504
[44] Fan Q Y Wei Q Yan H Y Zhang M G Zhang Z X Zhang J Q Zhang D Y 2014 Comput. Mater. Sci. 85 80
[45] Marmier A Lethbridge Z A D Walton R I Smith C W Parker S C Evans K E 2010 Comput. Phys. Commun. 181 2102
[46] Hu W C Liu Y Li D J Zeng X Q Xu C S 2014 Comput. Mater. Sci. 83 27