Electronic , Optical , Mechanical and Lattice Dynamical Properties of MgBi 2 O 6 : A First-Principles Study

Electronic structure, optical, mechanical, and lattice dynamical properties of the tetragonal MgBi2O6 are studied using a first-principles method. The band gap of MgBi2O6 calculated from the PBE0 hybrid functional method is about 1.62 eV and agrees well with the experimental value. The calculations on elastic constants show that MgBi2O6 exhibits mechanical stability and strong elastic anisotropy. The detailed analysis of calculated optical parameters and effective masses clearly indicate that MgBi2O6 has strong optical response in the visible light region and high separation efficiency of photoinduced electrons and holes.


Introduction
Currently, Bi-based oxides have received considerable attention due to their particular physical properties and wide applications in different fields like multiferroics [1,2], superconductivity [3,4], and photocatalysis [5,6].Generally, Bi exists as the trivalent state (Bi 3+ ) in most of the Bi-based oxides like Bi 2 O 3 [7], BiVO 4 [8], Bi 2 WO 6 [9], Bi 2 Sn 2 O 7 [10], BiFeO 3 [11], and BiMnO 3 [12].However, some Bi-containing oxides with the unusual pentavalent state (Bi 5+ ) have also attracted research interest.For example, NaBiO 3 has been found to show the absorption of visible light and can be used as a prominent material for photooxidation of organics [13].A recent work by Gong et al. [14] shows that AgBiO 3 can self-produce significant amounts of reactive oxygen species without light illumination or any other additional oxidant and has an excellent oxidizing reactivity.BaBiO 3 , as one kind of Bi-based oxides containing Bi 3+ and Bi 5+ mixed valent states, has been found to show the potential use for the absorber of all-oxide photovoltaics [15] and can be an active photocatalyst under visible-light irradiation [16].
MgBi 2 O 6 adopting the trirutile-type structure is also a Bi 5+ -containing compound.Kumada et al. [17] first successfully prepared MgBi 2 O 6 by the low-temperature hydrothermal method and characterized its crystal structure in detail.In 2003, Mizoguchi et al. [18] investigated the optical and electrical properties of MgBi 2 O 6 and found that MgBi 2 O 6 is a degenerate n-type semiconductor with the band gap of about 1.8 eV, and it is possible to produce an optical band gap that extends into the visible region of the spectrum by tuning the Bi 5+ 6s-O 2p interaction in order to produce a disperse conduction band.The band structure makes MgBi 2 O 6 a good candidate of visible light-sensitive photocatalysts for decomposition of organic species.Comparing with other pentavalent bismuthates such as LiBiO 3 , NaBiO 3 , KBiO 3 , ZnBi 2 O 6 , SrBi 2 O6, AgBiO 3 , BaBi 2 O 6 , and PbBi 2 O 6 , MgBi 2 O 6 does not possess the highest photocatalytic activity, but there is no adsorption observed in the decomposition of methylene blue [19].It may be a good idea to adjust the wide band gap of some traditional photocatalysts such TiO 2 by building complex compounds with MgBi2O 6 .In our recent work, the photocatalytic activity of MgBi 2 O 6 has been found to be significantly enhanced via constructing AgBr/MgBi 2 O 6 heterostructured composites [20].Theoretically, the band gap of MgBi 2 O 6 was calculated to be about 1.10 eV using the Heyd-Scuseria-Ernzerhof (HSE) functional method within the framework of the density functional theory (DFT), and was found to be widely tuned by applying external strain [21].On the basis of theoretical calculations, however, the deeper understanding of the physical properties of MgBi 2 O 6 is still in lacking.In this work, we investigated the electronic structure, optical, mechanical, and lattice dynamical properties MgBi 2 O 6 using first-principles calculations.

Computational Details
All calculations were carried out using the Vienna ab-initio simulation package (VASP) [22,23], which is an implement of DFT.The projector augmented wave (PAW) method [24,25] was used to describe the ion-electron interactions.The generalized gradient approximation (GGA) parameterized by Perdew, Burke and Ernzerhof (PBE) [26] was applied for the exchange-correlation function.The cutoff energy for the plane wave basis set was fixed at 500 eV.The Brillouin zone is sampled by a Monkhorst-Pack type k-point mesh with density of 2π × 0.03 Å −1 .Full relaxation of MgBi 2 O 6 unit cell was performed until the changes in total energy and force on each atom are less than 10 −5 eV and 10 −3 eV/Å, respectively.In this work, Mg-3s, Mg-2p, Bi-6s, Bi-6p, O-2s, and O-2p states are taken as valence electrons.Based on the DFT calculation, a more accurate screened coulomb hybrid functional (HF) developed by Heyd, Scuseria, and Ernzerhof (PBE0) [27,28] (containing 30% of the exact exchange, and 70% of the PBE exchange, and 100% of the PBE correlation energy in this work) was used to calculate the electronic structure.In the calculations of optical properties, a dense k-point density of 2π × 0.015 Å −1 including the Γ point was used to ensure the accurate precision.Moreover, 88 extra empty bands are included in the calculations to hold the excited electronic states.Phonon calculations were performed within the framework of density functional perturbation theory to investigate the lattice dynamical properties of MgBi 2 O 6 .The phonon dispersion curve and phonon density of states (PDOS) were calculated by the PHONOPY code [29] on the basis of force constants obtained from VASP code.In the phonon calculations, a 2 × 2 × 1 supercell was used.

Optimized Crystal Structure of MgBi 2 O 6
At ground state, MgBi 2 O 6 crystallizes in a tetragonal trirutile-type structure with P4 2 /mnm space group.It can be seen from the crystal structure depicted in Figure 1 that BiO 6 and MgO 6 octahedrons connect with the sharing edges and stack along the [001] direction by 2:1 ratio.In the [110] and 110 directions, BiO 6 and MgO 6 octahedrons are connected with each other by sharing O vertices.The optimized lattice constants listed in Table 1 are in good agreement with the experimental values [17].The BiO 6 octahedrons are slightly distorted and the bond lengths of Bi-O(8j) are not completely identical.

Electronic Properties of MgBi2O6
To study the electronic property of tetragonal MgBi2O6, the electronic band structure along high symmetry directions in the Brillouin zone (BZ) and density of states are calculated.The band structure calculated by conventional DFT (dotted line in Figure 2) shows that there is a bit of overlap between the conduction band and valence band, indicating the metallic feature of MgBi2O6.This contradicts the experimental findings.It is well known that DFT calculations do not take into account the effects of electron excitation and thus tend to underestimate the electronic band gap.To acquire the more accurate results, the calculations based on the HF PBE0 method are further performed.The HF corrected band structure (solid line in Figure 2) shows that MgBi2O6 is a direct semiconductor with a band gap of about 1.62 eV, which is well consistent with the experimental measured data (1.6~1.8 eV) [18][19][20] and larger than the previously calculated value (1.10 eV) based on the HF method [21].The valence band near the Fermi level (EF) is flat, while the dispersion of the conduction band close to the EF is relative strong.To study the electronic property of tetragonal MgBi 2 O 6 , the electronic band structure along high symmetry directions in the Brillouin zone (BZ) and density of states are calculated.The band structure calculated by conventional DFT (dotted line in Figure 2) shows that there is a bit of overlap between the conduction band and valence band, indicating the metallic feature of MgBi 2 O 6 .This contradicts the experimental findings.It is well known that DFT calculations do not take into account the effects of electron excitation and thus tend to underestimate the electronic band gap.To acquire the more accurate results, the calculations based on the HF PBE0 method are further performed.The HF corrected band structure (solid line in Figure 2) shows that MgBi 2 O 6 is a direct semiconductor with a band gap of about 1.62 eV, which is well consistent with the experimental measured data (1.6~1.8 eV) [18][19][20] and larger than the previously calculated value (1.10 eV) based on the HF method [21].The valence band near the Fermi level (E F ) is flat, while the dispersion of the conduction band close to the E F is relative strong.
where En(k) represents the band energy, k is the wave vector, and ℏ is the reduced Planck constant.For MgBi2O6, the effective electron mass (me * ) at the conduction band and hole mass (mh * ) at the valence band at the Γ point are calculated and listed in Table 2.The calculated tensors of m * clearly show that me * is fairly isotropic and mh * shows a relatively strong anisotropy.mh * along the [001] direction is larger than those along the [100] and [010] directions.The values of me * in Table 2 are also comparable to the previously calculated value (0.277) for MgBi2O6 [21].Moreover, the values of mh * are distinctly larger than those of me * , which indicates that the mobility of holes at the valence band is obviously slower than that of electrons at the conduction band.The big difference in mobility between electrons and holes is undoubtedly beneficial to the separation of charge carriers and the reduction of the recombination rate of electron-hole pairs.In addition, the values of mh * /me * of MgBi2O6 are in the range of 6.2~10.9,being larger than the corresponding value (2.1) of anatase TiO2 [30] widely investigated as an important photocatalytic material.Thus, semiconducting MgBi2O6 can be considered as a photocatalyst with highly efficient separation of photoinduced electrons and holes.
Table 2. Calculated effective electron mass (me * ) and hole mass (mh * ) along the three principal directions at the Γ point.All values are in units of free electron mass (m0).The effective mass (m * ) of carriers is an important parameter, because it determines the transfer and separation efficiency of electrons and holes, and directly influences the photophysical properties of semiconductors.m * can be calculated as follows:

Conduction Band
where E n (k) represents the band energy, k is the wave vector, and is the reduced Planck constant.For MgBi 2 O 6 , the effective electron mass (m e * ) at the conduction band and hole mass (m h * ) at the valence band at the Γ point are calculated and listed in Table 2.The calculated tensors of m * clearly show that m e obviously slower than that of electrons at the conduction band.The big difference in mobility between electrons and holes is undoubtedly beneficial to the separation of charge carriers and the reduction of the recombination rate of electron-hole pairs.In addition, the values of m h * /m e * of MgBi 2 O 6 are in the range of 6.2~10.9,being larger than the corresponding value (2.1) of anatase TiO 2 [30] widely investigated as an important photocatalytic material.Thus, semiconducting MgBi 2 O 6 can be considered as a photocatalyst with highly efficient separation of photoinduced electrons and holes.

Optical Properties of MgBi2O6
In order to understand the optical performance of semiconducting MgBi2O6, we calculated its dielectric functions using the hybrid functional PBE0 method, and further computed its optical properties such as complex refractive index n, extinction coefficient k, reflectivity R, absorption coefficients α, and electron energy-loss function L. The calculated dielectric functions are depicted in Figure 5.It can be clearly seen that the static dielectric function is found to be 4.1 at 0 eV.With the increase of the polarization intensity and energy, the dielectric function, the real part ε1 is also gradually increased and reaches a maximum value of 5.4 when the energy value is about 2.8 eV.When the photon energy reaches 1.6 eV, which is the width of electronic band gap, the electrons in the valence bands begin to excite and transit to conduction bands.As the carrier concentration increases, the degree of polarization decreases and the dielectric function decreases slightly.With the further increment of photon energy, the dielectric function values begins to fluctuate, which corresponds to the change of carrier concentration in the crystal.The imaginary part ε2 reflects the transition between occupied and non-occupied electrons and can be used to characterize the light absorption behavior of the crystal.The first two peaks of the imaginary part of the dielectric function

Optical Properties of MgBi 2 O 6
In order to understand the optical performance of semiconducting MgBi 2 O 6 , we calculated its dielectric functions using the hybrid functional PBE0 method, and further computed its optical properties such as complex refractive index n, extinction coefficient k, reflectivity R, absorption coefficients α, and electron energy-loss function L. The calculated dielectric functions are depicted in Figure 5.It can be clearly seen that the static dielectric function is found to be 4.1 at 0 eV.With the increase of the polarization intensity and energy, the dielectric function, the real part ε 1 is also gradually increased and reaches a maximum value of 5.4 when the energy value is about 2.8 eV.When the photon energy reaches 1.6 eV, which is the width of electronic band gap, the electrons in the valence bands begin to excite and transit to conduction bands.As the carrier concentration increases, the degree of polarization decreases and the dielectric function decreases slightly.With the further increment of photon energy, the dielectric function values begins to fluctuate, which corresponds to the change of carrier concentration in the crystal.The imaginary part ε 2 reflects the transition between occupied and non-occupied electrons and can be used to characterize the light absorption behavior of the crystal.The first two peaks of the imaginary part of the dielectric function ε 2 are, respectively, at about 2.0 and eV, which probably results from the electron transition from the top of the valence bands to the bottom of the conduction bands.ε2 are, respectively, at about 2.0 and 3.0 eV, which probably results from the electron transition from the top of the valence bands to the bottom of the conduction bands.The calculated refractive index n, extinction coefficient k, absorption coefficient α, reflectivity R, and electron energy loss function L of MgBi2O6 are shown in Figure 6.As presented in Figure 6a, the higher value of n is in the energy range of 1.5-5.8eV and the corresponding wavelength range is from 214 to 826 nm, indicating that MgBi2O6 has a strong refractive effect in both ultraviolet (UV) and visible light region.The value of k rapidly starts to increase from 2.0 eV, also showing a response in the visible and UV light region.The calculated α in Figure 6c shows that the light absorption edge is about 1.6 eV and is comparable with the band gap calculated by the PBE0 method.The value of optical absorption edge also agrees with literatural values [18,20].With the increase of energy, α also gradually increases and a series of absorption peaks appear in the energy range from 1.6 to 27 eV.Combining with the calculated dielectric function and density of states, we can find that the first two absorption peaks at 4.5 and 6.0 eV are probably related to the electron migration of the O-2p, Bi-6s, Bi-5d, and Mg-2p states.The reflectivity R and the electron energy loss function L can be used to represent the resonant frequency of the incident light and the resonant frequency of the plasma.As shown in Figure 6d, the average value of R is about 13.8%, indicating that MgBi2O6 can be used as a light absorbing material.The calculated L presented in Figure 6e completely locates in the continuous energy range of 0-45 eV, indicating that the characteristics of plasma oscillation in MgBi2O6 are not obvious.This is strongly different from the behavior appearing in Bi2Sn2O7 [31].The calculated refractive index n, extinction coefficient k, absorption coefficient α, reflectivity R, and electron energy loss function L of MgBi 2 O 6 are shown in Figure 6.As presented in Figure 6a, the higher value of n is in the energy range of 1.5-5.8eV and the corresponding wavelength range is from 214 to 826 nm, indicating that MgBi 2 O 6 has a strong refractive effect in both ultraviolet (UV) and visible light region.The value of k rapidly starts to increase from 2.0 eV, also showing a response in the visible and UV light region.The calculated α in Figure 6c shows that the light absorption edge is about 1.6 eV and is comparable with the band gap calculated by the PBE0 method.The value of optical absorption edge also agrees with literatural values [18,20].With the increase of energy, α also gradually increases and a series of absorption peaks appear in the energy range from 1.6 to 27 eV.Combining with the calculated dielectric function and density of states, we can find that the first two absorption peaks at 4.5 and 6.0 eV are probably related to the electron migration of the O-2p, Bi-6s, Bi-5d, and Mg-2p states.The reflectivity R and the electron energy loss function L can be used to represent the resonant frequency of the incident light and the resonant frequency of the plasma.As shown in Figure 6d, the average value of R is about 13.8%, indicating that MgBi 2 O 6 can be used as a light absorbing material.The calculated L presented in Figure 6e completely locates in the continuous energy range of 0-45 eV, indicating that the characteristics of plasma oscillation in MgBi 2 O 6 are not obvious.This is strongly different from the behavior appearing in Bi 2 Sn 2 O 7 [31].

Mechanical Properties of MgBi2O6
To investigate the mechanical properties of tetragonal MgBi2O6, the six independent elastic constants Cij are calculated by applying finite distortions of the lattice.The results are listed in Table 3.For the tetragonal system, the mechanical stability criterion [32] is given by the following relationships:

Mechanical Properties of MgBi 2 O 6
To investigate the mechanical properties of tetragonal MgBi 2 O 6 , the six independent elastic constants C ij are calculated by applying finite distortions of the lattice.The results are listed in Table 3.For the tetragonal system, the mechanical stability criterion [32] is given by the following relationships: The calculated values of C ij in Table 3 satisfy the stability criterion mentioned above, indicating that the tetragonal MgBi 2 O 6 has mechanical stability.As listed in Table 3, C 11 is significantly smaller C 33 , indicating that the chemical bonding strength in the (100) and (010) directions is significantly weaker than the bonding strength in the (001) direction.In addition, C 44 is obviously smaller than C 66 , which demonstrates that it is easier for shear deformation to occur along the (001) direction in comparison with the (010) direction.The shear elastic anisotropy of the material can be estimated by the relation A = 2C 66 /(C 11 − C 12 ).Typically, if A has a value of 1, meaning that the material is isotropic.The more the value of A deviates from 1, the elastic anisotropy would be more prominent.For MgBi 2 O 6 , the calculated A value is 4.6, indicating that MgBi 2 O 6 is highly anisotropic.
Young's modulus E and Poisson's ratio υ can be estimated from the bulk and shear moduli The calculated B, G, E, and υ are listed in Table 3.The B/G value is about 1.81.According to Pugh's criteria of brittleness and ductility [33], MgBi 2 O 6 exhibits some toughness, but it is not obvious, which is consistent with the Poisson's ratio υ = 0.27.In addition, the elastic anisotropy of MgBi 2 O 6 can be directly determined by the direction-dependent Young's modulus.The Young's modulus in a specific direction can be expressed by the elastic compliances (s ): where A is the matrix associated with the change of axes: Using the reduction of s ijkl for the tetragonal crystal class [34], the reduced Young's modulus with orientation can be expressed as follows: The calculated directional dependence of Young's modulus depicted in Figure 7

Lattice Dynamical Properties of MgBi2O6
The phonon dispersion curves along high symmetry directions and the phonon density of states are shown in Figure 8.There are 18 atoms in MgBi2O6 cells, so there is a total of 54 vibration modes in the phonon spectrum, including three acoustic modes and 51 optical modes.The calculated phonon spectrum shows no imaginary frequency, indicating the dynamical stability of MgBi2O6.Along F-Q and Q-Z paths in the Brillouin zone, the vibration modes are double degenerate.The acoustic modes phonons reflect the vibration of the centroid of the original cell and occupy the 0-3 THz frequency region.Among the three acoustic modes passing through the Γ point, the longitudinal mode has higher frequency in comparison with the other two transverse acoustic modes.There is no gap at the frequency range of about 3 THz between the longitudinal acoustic and transverse optical modes.Thus, phonons can transition from the acoustic mode to the optical mode without any momentum transfer [35].Further combining with the calculated PDOS shown in Figure 8

Lattice Dynamical Properties of MgBi 2 O 6
The phonon dispersion curves along high symmetry directions and the phonon density of states are shown in Figure 8.There are 18 atoms in MgBi 2 O 6 cells, so there is a total of 54 vibration modes in the phonon spectrum, including three acoustic modes and 51 optical modes.The calculated phonon spectrum shows no imaginary frequency, indicating the dynamical stability of MgBi 2 O 6 .Along F-Q and Q-Z paths in the Brillouin zone, the vibration modes are double degenerate.The acoustic modes phonons reflect the vibration of the centroid of the original cell and occupy the 0-3 THz frequency region.Among the three acoustic modes passing through the Γ point, the longitudinal mode has higher frequency in comparison with the other two transverse acoustic modes.There is no gap at the frequency range of about 3 THz between the longitudinal acoustic and transverse optical modes.Thus, phonons can transition from the acoustic mode to the optical mode without any momentum transfer [35].Further combining with the calculated PDOS shown in Figure 8

Conclusions
In this work, we have investigated the structural, electronic, optical, mechanical, and lattice dynamical properties of the trirutile-type MgBi2O6 in detail using the first-principles calculations.The calculated band gap of MgBi2O6 from the Heyd-Scuseria-Ernzerhof hybrid functional PBE0 electrical is about 1.62 eV and consistent with the experimental data (1.6~1.8 eV).The calculated effective masses show that the mobility of holes at conduction band is obviously slower than that of electrons at valence band, indicating high separation efficiency of electrons and holes in MgBi2O6.The calculated results of optical parameters clearly show that MgBi2O6 has strong light response in the visible light region and can be used as a light absorbing material.The calculated elastic constants and phonon dispersion clearly show that MgBi2O6 is mechanically and dynamically stable.Moreover, MgBi2O6 exhibits significantly elastic anisotropy.

Conclusions
In this work, we have investigated the structural, electronic, optical, mechanical, and lattice dynamical properties of the trirutile-type MgBi 2 O 6 in detail using the first-principles calculations.The calculated band gap of MgBi 2 O 6 from the Heyd-Scuseria-Ernzerhof hybrid functional PBE0 electrical is about 1.62 eV and consistent with the experimental data (1.6~1.8 eV).The calculated effective masses show that the mobility of holes at conduction band is obviously slower than that of electrons at valence band, indicating high separation efficiency of electrons and holes in MgBi 2 O 6 .The calculated results of optical parameters clearly show that MgBi 2 O 6 has strong light response in the visible light region and can be used as a light absorbing material.The calculated elastic constants and phonon dispersion clearly show that MgBi 2 O 6 is mechanically and dynamically stable.Moreover, MgBi 2 O 6 exhibits significantly elastic anisotropy.

Figure 1 .
Figure 1.Crystal structure of the trirutile-type MgBi2O6 represented as two different views.

Figure 1 .
Figure 1.Crystal structure of the trirutile-type MgBi 2 O 6 represented as two different views.

Figure 2 .
Figure 2. Calculated electronic band structure of MgBi2O6 using hybrid functional PBE0 method.Band structure from conventional DFT calculation (dotted line) is also presented for comparison.

Figure 2 .
Figure 2. Calculated electronic band structure of MgBi 2 O 6 using hybrid functional PBE0 method.Band structure from conventional DFT calculation (dotted line) is also presented for comparison.


) and hole mass (m h * ) along the three principal directions at the Γ point.All values are in units of free electron mass (m 0 ).The calculated electronic total and partial density of states (DOS) of MgBi 2 O 6 is shown in Figure 3.It can be clearly seen that the DOS near the top of the valence band is derived from the Bi-5d, Mg-2p and O-2p states.The bottom of the conduction band is mainly composed of O-2p, Bi-6s, and Mg-3s states.The strong hybridizations between Bi-5d and O-2p states in BiO 6 octahedrons and bonding interactions between Mg-2p or Mg-3s and O-2p states in MgO 6 octahedrons should be directly responsible for the structural stability of MgBi 2 O 6 .Appl.Sci.2019, 9, x FOR PEER REVIEW 5 of 13The calculated electronic total and partial density of states (DOS) of MgBi2O6 is shown in Figure3.It can be clearly seen that the DOS near the top of the valence band is derived from the Bi-5d, Mg-2p and O-2p states.The bottom of the conduction band is mainly composed of O-2p, Bi-6s, and Mg-3s states.The strong hybridizations between Bi-5d and O-2p states in BiO6 octahedrons and bonding interactions between Mg-2p or Mg-3s and O-2p states in MgO6 octahedrons should be directly responsible for the structural stability of MgBi2O6.

Figure 3 .
Figure 3. Calculated electronic density of states of MgBi2O6 using PBE0 method.

Figure 4
shows the calculated ELF for (110) plane of MgBi2O6 passing through Bi, Mg, and O atoms.The ELF values between Mg and O atoms are extremely low, indicating the ionic bonding in MgO6 octahedrons.This can be well understood on the basis of the previously obtained recognition on MgO, which is usually considered as one kind of classical ionically bonded compounds.The maximal value of ELF between Bi and O atoms reaches about 0.61, which clearly identifies the partially Bi-O covalent bonding interaction in BiO6 octahedrons.

Figure 3 .
Figure 3. Calculated electronic density of states of MgBi 2 O 6 using PBE0 method.To get a deeper understanding of the bonding nature of Bi, Mg, and O atoms, the electron localization function (ELF) was further calculated.According to the original definition, the ELF values are scaled in the range from 0 to 1.The high ELF means strong covalent bonding interaction between atoms and the very low ELF close to 0 corresponds to the ionic bonding.Figure 4 shows the calculated ELF for (110) plane of MgBi 2 O 6 passing through Bi, Mg, and O atoms.The ELF values between Mg and O atoms are extremely low, indicating the ionic bonding in MgO 6 octahedrons.This can be well understood on the basis of the previously obtained recognition on MgO, which is usually considered as one kind of classical ionically bonded compounds.The maximal value of ELF between Bi and O atoms reaches about 0.61, which clearly identifies the partially Bi-O covalent bonding interaction in BiO 6 octahedrons.

Figure 4
shows the calculated ELF for (110) plane of MgBi 2 O 6 passing through Bi, Mg, and O atoms.The ELF values between Mg and O atoms are extremely low, indicating the ionic bonding in MgO 6 octahedrons.This can be well understood on the basis of the previously obtained recognition on MgO, which is usually considered as one kind of classical ionically bonded compounds.The maximal value of ELF between Bi and O atoms reaches about 0.61, which clearly identifies the partially Bi-O covalent bonding interaction in BiO 6 octahedrons.
, we can find that the coupled vibrations of Bi, Mg, and O atoms including the O-Bi-O and O-Mg-O bending vibrations in BiO6 and MgO6 octahedrons within the frequency range from 3 to 9.8 THz and the Bi-O and Mg-O stretching vibrations in the frequency range from 12 to 16.4 THz are obvious.The phonon branches above 16.4THz are completely ascribed to the vibrations of O atoms.
, we can find that the coupled vibrations of Bi, Mg, and O atoms including the O-Bi-O and O-Mg-O bending vibrations in BiO 6 and MgO 6 octahedrons within the frequency range from 3 to 9.8 THz and the Bi-O and Mg-O stretching vibrations in the frequency range from 12 to 16.4 THz are obvious.The phonon branches above 16.4THz are completely ascribed to the vibrations of O atoms.

Figure 8 .
Figure 8. Calculated phonon dispersion curves and phonon density of states of MgBi2O6.

Figure 8 .
Figure 8. Calculated phonon dispersion curves and phonon density of states of MgBi 2 O 6 .

Table 1 .
[17]mized structural parameters of MgBi2O6, Bi-O and Mg-O interatomic distances.The experimental values[17]are listed in brackets for comparison.

Table 1 .
Optimized structural parameters of MgBi 2 O 6 , Bi-O and Mg-O interatomic distances.The experimental values [17] are listed in brackets for comparison.

Table 2 .
Calculated effective electron mass (m e *