Density Functional Theory Study of Electronic Structure and Optical Properties of Ln 3+ -Doped γ -Bi 2 MoO 6 (Ln=Gd, Ho, Yb)

: Based on density functional theory (DFT), theoretical models of three kinds of lanthanide rare earth metal ion-doped γ -Bi 2 MoO 6 were constructed (Ln-BMO (Ln=Gd, Ho, Yb)). The geometric structure, electronic structure, and optical properties of the model were calculated, and the inﬂuence of doped Ln 3+ ions on the structures and properties of the system was analyzed. The results revealed that the substitution of smaller ionic radius Ln 3+ ions for Bi 3+ ions caused a contraction of the lattice parameters. At the same time, the contribution of the [Ln]4d near valence band and conduction band reduced the bandwidth of γ -Bi 2 MoO 6 , forming the Ln-O ionic bond with different strengths to obtain higher charge conductivity and charge-separation ability. Secondly, Ln 3+ ions have a strongly ionic charge, which leads to the appearance of optical absorption bands in the infrared region and part of the visible region. This reduces the reﬂection in the visible region, improves the utilization rate, delays the loss of electron energy, and promotes phase matching in the visible region. And the Gd 3+ - doped system has better photocatalytic activity than the other Ln 3+ -doped system. This research provides theoretical insights into doped lanthanide rare earth ions and also provides strategies for the modiﬁcation of γ -Bi 2 MoO 6 nanomaterials.


Introduction
Semiconductor crystal materials possess the characteristics of periodicity and symmetry in the atomic structure, enabling interaction and transformation of various forms of energy within the crystal, including power, heat, light, electricity, etc. [1][2][3][4].It is the exceptional versatility and remarkable stability of semiconductor crystals that make indispensable core materials in advanced fields such as laser technology, atomic energy technology, and optoelectronic technology.
In the field of photocatalytic degradation, scholars have investigated the crystal structure of the bismuth molybdate material and found that among the three crystal structures, γ-Bi 2 MoO 6 with an orthogonal structure has a suitable band gap width (2.5-2.8 eV).This structure has a strong response to visible light, making it a photocatalytic material which is driven by visible light.However, due to the large electron effective reciprocity mass, electron-hole-pair recombination, and small specific surface area, the photocatalytic activity of bismuth molybdate crystals is limited [32,33].Therefore, it is necessary to modify the intrinsic material by doping, recombination, modification, and loading to improve its properties.Doping is one of the effective methods to modify intrinsic materials.
Meng et al. [34].successfully synthesized the Fe 3+ -doped photocatalytic material Fe-BMO, which exhibited excellent photocatalytic activity for photocatalytic nitrogen fixation.According to the results of experiments and DFT, the introduction of the iron medium reduces the work function of the intrinsic material and promotes the migration of surface electrons, while accelerating the redox between Fe 3+ /Fe 2+ to facilitate charge collection.Yu et al. [35] fluorinated Bi 2 MoO 6 materials that were prepared by the combination of solvothermal and calcination methods, which showed excellent degradation efficiency for dyes, antibiotics, phenols, and other pollutants.The study demonstrated that doped F - ions promoted the ability of the intrinsic material surface and oxygen vacancies to capture electrons, and increased the output of active functional groups.Secondly, the introduction of F -ions made the valence band of the intrinsic material have a stronger positive potential and improved the oxidation ability of holes.Li et al. [36] successfully prepared Lu 3+ -doped Bi 2 MoO 6 nanorods.According to the combination of experiment and DFT, it was found that by doping different amounts of Lu 3+ ions, the intrinsic material gradually transformed from nano-sheet to a rod-like morphology.At the same time, the lattice surface selectively exposed the (0 0 1) surface, and Lu 3+ ions also promoted the generation of oxygen vacancy in the intrinsic lattice to promote charge separation and light response.
The electronic and optical properties of compounds can be improved by doping rare earth metal ions, transition metal ions, and non-metal ions.Alemi et al. [37] prepared three different lanthanide rare earth ions Ln 3+ (Ln=Gd, Ho, Yb)-doped γ-Bi 2 MoO 6 by thermal polymerization, and found that different rare earth metal ions doped to form flakes with different shape structures; all three types of doping promoted the photocatalytic activity of the intrinsic material.It was mainly attributed to the enlarged specific surface area and reduced bandgap, resulting in a larger reaction contact area and better electron separation efficiency.At the same time, in the photocatalytic experiment, through the comparison of different concentrations of dopant ions, it was found that the photocatalytic performance of the product Gd 0.5 Bi 1.5 MoO 6 was improved the most, which was attributed to the half-filled electronic configuration of Gd 3+ ions.
From the experimental point of view, this paper comprehensively analyzes the influence and improvement of three kinds of the doping of lanthanide rare earth metal ions on the photocatalytic performance of γ-Bi 2 MoO 6 , but does not involve the corresponding theoretical part.Considering the significance of lanthanide rare earth metal ions for improving the photocatalytic activity of γ-Bi 2 MoO 6 , this paper constructs theoretical models for three lanthanide rare earth metal ion-doped γ-Bi 2 MoO 6 (Ln-BMO (Ln=Gd, Ho, Yb)), and calculates and analyzes the geometric structure, electronic structure, and optical properties of Ln-BMO by density functional theory (DFT).Compared with the experimental data, this study not only provides theoretical support for the experiment, but also provides the basis and guidance for rare earth metal ion-doped γ-Bi 2 MoO 6 .Crystals 2023, 13, x FOR PEER REVIEW 3 data, this study not only provides theoretical support for the experiment, but also vides the basis and guidance for rare earth metal ion-doped γ-Bi2MoO6.

Theoretical Model
As described in the introduction, γ-Bi2MoO6 can be represented by the genera mula (Bi2O2)(MonO3n+1), where n = 1, in which the structure is formed by the interlaci MoO6 octahedron and Bi-O layer (Bi2O2).As for the initial file of γ-Bi2MoO6 crystal, t file (mp-23064) is downloaded from the materials project resource library [38][39][40].crystal structure belongs to the orthorhombic structure, and the space group

Computational Details
Vienna Ab-initio Simulation Package (VASP) is a first-principle program based o combination of plane wave basis sets and ultrasoft pseudopotential.It is widely us material simulation calculations due to its high precision and convenient use [41] DFT calculation method of this work was carried out using the VASP software comb with the projected augmented wave (PAW) method, and the related functional descri is calculated using the generalized gradient approximation (GGA-PBE) [42].

Computational Details
Vienna Ab-initio Simulation Package (VASP) is a first-principle program based on the combination of plane wave basis sets and ultrasoft pseudopotential.It is widely used in material simulation calculations due to its high precision and convenient use [41].The DFT calculation method of this work was carried out using the VASP software combined with the projected augmented wave (PAW) method, and the related functional description is calculated using the generalized gradient approximation (GGA-PBE) [42].

SCF Structure Converges
Firstly, the crystal structures of BMO and Ln-BMO were optimized and converged.In previous work, the lattice parameters of pure BMO and Ln-BMO with different proportions were calculated using X-ray diffraction spectral data.The results were compared to relevant information after DFT structure convergence, as shown in Table 1.In this table, EXP represents experimental data and DFT represents structure convergence data.From Table 1, it is found that the lattice parameters of the intrinsic materials tend to shrink with the doping of lanthanide rare earth metal ions, and also shrink with the increase in doping ion concentration.Through the doping of three lanthanide rare earth ions, lanthanide rare earth metal ions with smaller radius replace Bi 3+ ions in the intrinsic material, resulting in the reduction of lattice parameters, which is consistent with the experimental analysis.At the same time, the average relative errors between theoretical calculation and experimental data are also given in the table, and the error is controlled within 3%.The possible causes of the error involve experimental measurement error and convergence error.It shows that the model and convergence accuracy of theoretical calculation are reliable.

Electronic Structure Information
To better illustrate how the doping of lanthanide rare earth metal ions affects the electronic structure of intrinsic materials, static calculations were carried out through the converged structural model, and the band structure, electron state density, and differential charge density of BMO and Ln-BMO were obtained.The detailed calculation results and analysis are as follows.

Band Structure
The band structures are obtained through theoretical calculation and simulation, and are shown in Figure 2a-d.Firstly, it should be noted that the maximum conduction band (CBM) and minimum valence band (VBM) of the band structure of BMO and Ln-BMO are in the same position in K-space, indicating that the materials before and after doping belong to the direct band gap semiconductor.It has better responsive activity than other semiconductor materials [46].Compared to BMO, the electronic states of Ln-BMO in the VB and CB are denser, indicating that more electron states are generated in the valence and conduction bands.The introduction of the Ln 3+ ions is the main reason for the appearance of these new electronic states.Secondly, the band gap value (Eg) of the band structure can reflect the carrier migration efficiency.In the figure, the undoped BMO band gap is Eg = 2.435 eV, while the calculated band gaps of Gd-BMO, Ho-BMO, and Yb-BMO doped with lanthanide rare earth metal ions are 2.407 eV, 2.368 eV, and 2.268 eV, respectively.Table 2 lists the comparison of the band gap between experimental and theoretical calculations.The experimental and theoretical trends are unified, and the errors are controlled within 4%.In addition, the carrier mobility is inversely proportional to the effective mass of the electron [47].The equation (y = A + Bx + Cx 2 ) is used for the quadratic fitting of the valence band top, and then the relevant electron effective mass can be obtained by fitting coefficient C. The relevant results are shown in Table S1.The effective mass of electron of the intrinsic material is effectively reduced by Ln 3+ doping.Among them, 0.5 Gd-BMO has the lowest electron effective mass (m = 0.624 m e ), indicating that its carrier migration efficiency is the best.These results show that the bandwidth decreases by the doping of lanthanide rare earth ions, which promotes the visible light response and promotes the carrier mobility of the intrinsic material.within 4%.In addition, the carrier mobility is inversely proportional to the effective mass of the electron [47].The equation ( =  +  +  2 ) is used for the quadratic fitting of the valence band top, and then the relevant electron effective mass can be obtained by fitting coefficient C. The relevant results are shown in Table S1.The effective mass of electron of the intrinsic material is effectively reduced by Ln 3+ doping.Among them, 0.5 Gd-BMO has the lowest electron effective mass (m = 0.624 me), indicating that its carrier migration efficiency is the best.These results show that the bandwidth decreases by the doping of lanthanide rare earth ions, which promotes the visible light response and promotes the carrier mobility of the intrinsic material.To determine the contribution of Ln atoms and other atomic orbitals to the electronic structure of BMO and Ln-BMO, the total density of state (TDOS) and the partial density of state (PDOS) are calculated, and the calculated results are shown in Figure 2e-h As for the doped materials with three lanthanide rare earth metal ions, compared with BMO, it is found that the introduction of Ln 3+ ions does not significantly effect the orbital contribution of Bi, Mo, and O atoms, indicating that a small amount of Ln 3+ doping

Density of Electronic States
To determine the contribution of Ln atoms and other atomic orbitals to the electronic structure of BMO and Ln-BMO, the total density of state (TDOS) and the partial density of state (PDOS) are calculated, and the calculated results are shown in Figure 2e-h As for the doped materials with three lanthanide rare earth metal ions, compared with BMO, it is found that the introduction of Ln 3+ ions does not significantly effect the orbital contribution of Bi, Mo, and O atoms, indicating that a small amount of Ln 3+ doping (0.25 Ln-BMO and 0.5 Ln-BMO) cannot change the crystal structure of γ-Bi 2 MoO 6 , which is consistent with the experimental analysis.Furthermore, the orbital contributions of [Gd] 4d, [Ho] 4d, and [Yb] 4d are found at −5.0 eV to 0.0 eV and 0.0 eV to 8.0 eV, and their contributions near CBM and VBM lead to the change in BMO band gap.It can be inferred that the Ln 3+ doping generates new interactions, which makes it easier for the intrinsic structure to excite electrons into the conduction band, and obtain higher charge conductivity and effective charge separation ability, thus proving that the doping of Ln 3+ ion can effectively improve the photocatalytic ability [49].

Differential Charge Density
The differential charge density of Ln 3+ ions in Ln-BMO is calculated by DFT, which reveals the influence of lanthanide rare earth metal ions doping on the electron transport performance of γ-Bi 2 MoO 6 (Figure 3a-c,g-i).The yellow part represents electron depletion and the blue part represents electron accumulation.Three lanthanide rare earth ions show similar differential charge densities, indicating that the doping of the three ions has similar properties.The differential charge density shows that electrons accumulate near the three types of Ln atoms, while the depletion of electrons occurs in O near Ln atoms, and the nearby Bi atoms and Mo atoms are also disturbed.In particular, the interaction between the Ln atoms forms an Ln-Ln ionic bond in 0.5 Ln-BMO, which further facilitates electron transfer.In addition, the influence range of different Ln atoms on the electron density of the intrinsic material is estimated.This shows that the influence of Gd 3+ doping on the electron density of the system reaches 5.103 Å and 6.175 Å, which are higher than the values for Ho 3+ and Yb 3+ .The result shows that the photogenerated electrons of BMO are transferred to doped-ion sites, which inhibits the recombination of electron-hole pairs, and the influence range of 0.5 Gd-BMO on the electron density is stronger, which is consistent with the experimental analysis [50].The 3D differential charge cross section of (j) (0 0 1) 0.5 Gd-BMO; (k) (0 0 1) 0.5 Ho-BMO; (l) (0 0 1) 0.5 Yb-BMO.

Optical Structure Information
As a photocatalytic material that relies on visible light response, it is necessary to study its optical properties.The dielectric constant, as a bridge between the process of inter-band transition and the electronic structure of a solid, reflects the band structure and other spectral information of a solid [52].Thus, we analyze the influence of the doping of At the same time, the 3D differential charge cross-section is drawn based on differential charge density.The (0 0 1) (8.100 Å away from the origin), (0 1 0) (4.900 Å away from the origin), and (1 0 0) (4.320 Å away from the origin) cross-section of convergence structure is selected.As shown in Figure 3d-f,j-l and Figure S1, the blue part represents electron depletion and the red part represents electron accumulation.The darker the color, the higher the peak value, which represents more electron transfer.The figure reveals that the doping of Ln 3+ ions forms ionic bonds with different strengths of Ln-O, and the strength of the ionic bond is an important factor affecting the change in charge compensation [51].According to the peak value, the peak intensity of 0.25 Gd-BMO reaches + 0.03926 e/bohr 3 (0 0 1), + 0.03819 e/bohr 3 (0 1 0), and + 0.03824 e/bohr 3 (1 0 0), which is higher than the other two doped systems and shows isotropy under different cross sections.Similarly, the ionic bond strength of 0.5 Gd-BMO is also higher than the two other doping systems, showing strong electron transfer ability.The results show that Ln-O ionic bonds with different strengths are formed and charge compensation is generated in BMO system, showing the improvement of photo response and decolorization ability, and the Gd 3+ doping effect is more obvious.

Optical Structure Information
As a photocatalytic material that relies on visible light response, it is necessary to study its optical properties.The dielectric constant, as a bridge between the process of inter-band transition and the electronic structure of a solid, reflects the band structure and other spectral information of a solid [52].Thus, we analyze the influence of the doping of lanthanide rare earth metal ions on the optical properties of BMO and Ln-BMO materials.For the calculation of optical properties, we calculate the complex dielectric function, complex refractive index, reflectivity, absorption coefficient, and electron energy loss function of the material using Equations ( 1)-( 9).Among them, the Kramers-Kronig transformation dispersion relation of the imaginary part is used to solve the solution of the real part of the complex permittivity [53].The detailed formula is explained in the Supplementary Materials. (5) Crystals 2023, 13, 1158 8 of 15

Complex Dielectric Constant ε(ω)
The complex dielectric constant is mainly manifested in the binding and dissipative ability of charge, which can indirectly reflect the light utilization ability of semiconductor photocatalytic materials.According to Equations ( 3) and ( 4), the complex dielectric constant is shown in Figures 4a-c and S2a-c.Firstly, we find the real part and imaginary parts of the complex dielectric constant of BMO and Ln-BMO show similar curves in the whole photon energy range.In the real part curve, BMO shows obvious peaks at 3.082 eV, 5.674 eV, and 7.713 eV, while the real part of the dielectric constant of Ln-BMO migrates to high energy.For the imaginary part curve, BMO has an obvious peak at 4.167 eV and 6.442 eV.Interestingly, the doping of Ln 3+ ions induces the peak value at the lower energy to shift positively, while the opposite trend is observed for peak values at higher energies.These phenomena may be attributed to the interaction between [Ln] 4d and the intrinsic material.
For the photon energy range of the visible region (1.64-3.19eV), the real part of the complex dielectric constant can approximately describe the energy stored in a material upon exposure to photon radiation.In the infrared region, the real part of Ln-BMO shows high peaks, which are not present in the intrinsic materials.And for 0.5 Ln-BMO, the peak value of the infrared region has a different degree of change.These peaks correspond to electron transitions between the band gaps, which are attributed to the long-lived electronic energy levels of [Ln] 4d and the interaction between [O] 2p and [Ln] 4d [54,55].
When the photon energy is 0, the actual value is expressed as the static dielectric constant, in which Gd 3+ and Ho 3+ ion doping significantly improve it, and the higher the static dielectric constant is, the better the refraction coefficient is.In the visible region, the real part of the permittivity of Ln-BMO decreases slightly compared with BMO.
The imaginary part reflects the energy dissipated from the material to the outside upon exposure to photon radiation.Compared with the intrinsic material, Ln-BMO is enhanced in the infrared and visible regions, indicating that the carriers of the doped materials with Ln 3+ ions are more active and the optical activity is enhanced [56].The direct band gap of the corresponding system can be determined by the imaginary part of the dielectric constant.According to the calculation, the direct optical band gap of the BMO, Gd-BMO, Ho-BMO, and Yb-BMO systems is 2.454 eV, 2.436 eV, 2.417 eV, and 2.402 eV, respectively.The optical band gap is generally larger than the electronic band gap [57].At the same time, it is found that the imaginary part of the dielectric constant of 0.5 Gd-BMO increases more strongly in the visible region, indicating that the appropriate amount of Gd 3+ ions doped can improve the photo-response intensity of the material and obtain better photocatalytic activity.

Complex Refractive Index n(ω) and Birefringence ∆n
The complex refractive index and birefringence index of materials are important optical constants for analyzing the absorption spectrum.The parameters calculated according to Equations ( 5) and ( 6) are shown in Figures 4d,e and S2d-f.The real part determines the propagation speed of light in the material, while the extinction coefficient of the imaginary part represents the attenuation of light propagation in the material, and also indirectly represents the absorption capacity of the material.The extinction coefficient is the main analysis part of the complex refractive index.Among them, the attenuation of Ln-BMO is enhanced in the infrared and visible regions, which further indicates that the doping of Ln 3+ ions leads to the enhancement of the light absorption capacity in the intrinsic material.In particular, 0.5 Gd-BMO has the largest enhancement amplitude in the visible region, indicating that Gd 3+ ion doping promotes the absorption of visible light and light response.respectively.The optical band gap is generally larger than the electronic band gap the same time, it is found that the imaginary part of the dielectric constant of 0.5 G increases more strongly in the visible region, indicating that the appropriate am Gd 3+ ions doped can improve the photo-response intensity of the material and obt ter photocatalytic activity.

Complex Refractive Index 𝑛(𝜔) and Birefringence ∆𝑛
The complex refractive index and birefringence index of materials are import tical constants for analyzing the absorption spectrum.The parameters calculated ing to Equations ( 5) and ( 6) are shown in Figures 4d,e and S2d-f.The real part dete the propagation speed of light in the material, while the extinction coefficient of th inary part represents the attenuation of light propagation in the material, and als rectly represents the absorption capacity of the material.The extinction coefficien main analysis part of the complex refractive index.Among them, the attenuation BMO is enhanced in the infrared and visible regions, which further indicates that th ing of Ln 3+ ions leads to the enhancement of the light absorption capacity in the in material.In particular, 0.5 Gd-BMO has the largest enhancement amplitude in the region, indicating that Gd 3+ ion doping promotes the absorption of visible light an response.
As for the real part of the complex refractive index, it can be analyzed by solv birefringence.The birefringence can be obtained by the difference of the real par complex refractive index of the three spindles of the crystal.For BMO crystals, th fringence can be obtained from ∆ =  ∥ −  ⊥ ; the results are shown in Figure 4 the visible region, the birefringence of BMO is relatively stable around 0.1, bu proaches 0 near 2.854 eV, which may affect the phase matching in this region [58].BMO, this point is shifted to 3.086 eV, and the Gd-BMO is shifted more obviousl shift phenomenon promotes the phase matching of the intrinsic material in the region.As for the real part of the complex refractive index, it can be analyzed by solving the birefringence.The birefringence can be obtained by the difference of the real part of the complex refractive index of the three spindles of the crystal.For BMO crystals, the birefringence can be obtained from ∆n = n z − n ⊥z ; the results are shown in Figure 4f.Near the visible region, the birefringence of BMO is relatively stable around 0.1, but it approaches 0 near 2.854 eV, which may affect the phase matching in this region [58].For Ln-BMO, this point is shifted to 3.086 eV, and the Gd-BMO is shifted more obviously.This shift phenomenon promotes the phase matching of the intrinsic material in the visible region.

Reflectivity R(ω) and Absorption Coefficient α(ω)
The reflectivity and absorption coefficients of materials are crucial parameters in the calculation of optical properties.The percentage of the radiant energy reflected by the material in the total radiant energy is called reflectivity, and the amount of radiant energy absorbed by the material is called the absorption coefficient.The reflectivity and absorption coefficient of intrinsic material and doped material are calculated according to Equations ( 7) and (8).Figures 5a and 6a shows the calculated reflectance.The intrinsic material and the doped material exhibit similar reflectivity curves, indicating that the doping of a small amount of Ln 3+ ions (0.25 Ln-BMO and 0.5 Ln-BMO) does not change the optical structure of the intrinsic material, which is consistent with the previous analysis.The reflectivity of all calculated materials reaches its maximum at 7.850 eV, reaching approximately 40%.And, characteristic peaks different from those of the intrinsic materials appear at 28.471 eV, 30.219 eV, and 32.943 eV, which correspond to the doping of Gd 3+ , Ho 3+ , and Yb 3+ , respectively.With the increase in doping concentration, the corresponding characteristic peak is more obvious.Secondly, element doping reduces the reflectivity of intrinsic materials in non-characteristic bands, which confirms the enhancement of optical absorption.Figures 5c and 6c shows the calculated absorption coefficient.The valence and conduction bands of BMO and Ln-BMO are continuous, and the corresponding absorption curves are also continuous [59].The increase in photon energy promotes light absorption in the material, indicating the direct transition process of electrons.At the same time, the Ln-BMO observes characteristic peaks similar to the corresponding energy value of the reflectivity, which promotes the intrinsic material to have stronger light absorption in the higher photon energy band.
For the photon energy range of visible light, its reflectivity and absorption coefficient are shown in Figures 5b,d, and 6b,d.Firstly, the reflectivity of Ln-BMO is lower than BMO in the visible light region, but higher than BMO in the energy range less than 0.6 eV.And with the increase in doping concentration, its reflectivity in the visible region is lower, indicating that the doping of lanthanide rare earth metal ions can reduce the reflection in the visible light region and increase the utilization rate of the material for visible light.As for the calculated absorption coefficient, it can be observed that the doping of Ln 3+ ions presents a new absorption peak that cannot appear in the intrinsic material in the infrared region and visible region (less than 2.0 eV), which is attributed to the introduction of the Ln 3+ ions.And Ln 3+ ions have a smaller ionic radius and a stronger ionic charge that leads to the extra electronic state in the band gap, resulting in the appearance of the low-energy absorption band [54].Especially in the visible region, 0.25 Gd-BMO and 0.5 Gd-BMO exhibit a strong absorption coefficient and low reflectance, which further indicates that Gd-BMO has a strong light response and transmittance, and explains the reason that Gd-BMO has strong light catalytic ability, which is consistent with the previous analysis.

Energy Loss Function L(ω)
The electron energy loss function is mainly manifested in the primary process of energy loss caused by inelastic scattering, which cannot involve the characteristic excitation of the ground state process.According to Equation ( 9), the energy loss function of intrinsic material and doped material is calculated with the energy of the incident photon.Figures 5e,f and 6e,f display the loss functions of the calculation systems.In the figure, the main peaks are in 24.090 eV, 24.042 eV, 23.885 eV, and 23.794 eV, corresponding to BMO, Yb-BMO, Ho-BMO, and Gd-BMO, respectively, which are represented as plasma peaks.These peaks illustrate the collective excitation of electron charge density in the crystal structure [60].And Ln 3+ doping causes the plasma peak value to shift to the lower energy level, which may be attributed to the new spin-orbit effect of Ln 3+ and leads to the change in crystal field effect.Secondly, similar to the reflectance and absorption coefficient, the characteristic peaks of corresponding elements are also observed.No significant electron energy loss is observed in the infrared and visible energy ranges.Although the electron energy loss function increases with the doping of Ln 3+ , such an enhancement has no obvious effect on the optical response, and the energy loss is the result of ionization and electric dipole oscillation by external fields [61].
As described in the introduction, γ-Bi 2 MoO 6 can be represented by the general formula (Bi 2 O 2 )(Mo n O 3n+1 ), where n = 1, in which the structure is formed by the interlacing of MoO 6 octahedron and Bi-O layer (Bi 2 O 2 ).As for the initial file of γ-Bi 2 MoO 6 crystal, the cif file (mp-23064) is downloaded from the materials project resource library [38-40].This crystal structure belongs to the orthorhombic structure, and the space group number is 29, corresponding to Pca21.The lattice parameters of the original file are a = 5.527 Å, b = 5.528 Å, c = 16.298Å, α = β = γ = 90 • , Vo = 497.973Å 3 .The crystal structure contains two non-equivalent Bi 3+ and one Mo 6+ which bond with six non-equivalent O 2− , respectively, forming the Bi 2 O 2 tetrahedron and MoO 6 octahedra [38].The 2 × 1 × 1 supercell is used to construct the model.The doping model of three rare earth metal ions is constructed by replacing Bi atoms near the center of the supercell system, which is named Ln-BMO (Ln=Gd, Ho, Yb), and the undoped system is named BMO.The structural models are shown in Figure 1.
. For BMO, the Fermi energy from −20.0 eV to −15.0 eV is mainly contributed by O atoms; Mo atoms also contribute in the interval, and mainly come from [O] 1s and [Mo] 4d orbitals, while only [Bi] 6s orbital contributes near 10.0 eV.Secondly, from the energy of −6.0 eV to 0.0 eV in PDOS, Bi, Mo, and O can be found to participate in contribution, which can be composed of [O] 2p-[Mo] 4d and [O] 2p-[Bi] 6s orbital interactions [48].The similar pattern accounts for the orbital contribution from 0.0 eV to 8.0 eV, which explains the formation of a conductive state between the [Bi] 6p, [Mo] 4d, and [O] 2p orbits.
. For BMO, the Fermi energy from −20.0 eV to −15.0 eV is mainly contributed by O atoms; Mo atoms also contribute in the interval, and mainly come from [O] 1s and [Mo] 4d orbitals, while only [Bi] 6s orbital contributes near 10.0 eV.Secondly, from the energy of −6.0 eV to 0.0 eV in PDOS, Bi, Mo, and O can be found to participate in contribution, which can be composed of [O] 2p-[Mo] 4d and [O] 2p-[Bi] 6s orbital interactions [48].The similar pattern accounts for the orbital contribution from 0.0 eV to 8.0 eV, which explains the formation of a conductive state between the [Bi] 6p, [Mo] 4d, and [O] 2p orbits.

Figure 5 .
Figure 5.The reflectivity of BMO and 0.25 Ln-BMO: (a) full spectrum; (b) near the visible The absorption coefficient of BMO and Ln-BMO: (c) full spectrum; (d) near the visible reg energy loss function of BMO and Ln-BMO: (e) full spectrum; (f) near the visible region.

Figure 5 .
Figure 5.The reflectivity of BMO and 0.25 Ln-BMO: (a) full spectrum; (b) near the visible region.The absorption coefficient of BMO and Ln-BMO: (c) full spectrum; (d) near the visible region.The energy loss function of BMO and Ln-BMO: (e) full spectrum; (f) near the visible region.

Figure 5 .Figure 6 .
Figure 5.The reflectivity of BMO and 0.25 Ln-BMO: (a) full spectrum; (b) near the visible re The absorption coefficient of BMO and Ln-BMO: (c) full spectrum; (d) near the visible region energy loss function of BMO and Ln-BMO: (e) full spectrum; (f) near the visible region.

Figure 6 .
Figure 6.The reflectivity of 0.25 Ln-BMO and 0.5 Ln-BMO: (a) full spectrum; (b) near the visible region.The absorption coefficient of 0.25 Ln-BMO and 0.5 Ln-BMO: (c) full spectrum; (d) near the visible region.The energy loss function of 0.25 Ln-BMO and 0.5 Ln-BMO: (e) full spectrum; (f) near the visible region.

Table 1 .
This experimental result and DFT structure convergence result of lattice parameters of BMO and Ln-BMO.

Table 2 .
Experimental result and DFT structure convergence result of band gap of BMO and Ln-BMO.

Table 2 .
Experimental result and DFT structure convergence result of band gap of BMO and Ln-BMO.