First-Principles Studies for Electronic Structure and Optical Properties of p-Type Calcium Doped α-Ga2O3

Gallium oxide (Ga2O3) is a promising wide-band-gap semiconductor material for UV optical detectors and high-power transistor applications. The fabrication of p-type Ga2O3 is a key problem that hinders its potential for realistic power applications. In this paper, pure α-Ga2O3 and Ca-doped α-Ga2O3 band structure, the density of states, charge density distribution, and optical properties were determined by a first-principles generalized gradient approximation plane-wave pseudopotential method based on density functional theory. It was found that calcium (Ca) doping decreases the bandgap by introducing deep acceptor energy levels as the intermediate band above the valence band maximum. This intermediate valence band mainly consists of Ca 3p and O 2p orbitals and is adequately high in energy to provide an opportunity for p-type conductivity. Moreover, Ca doping enhances the absorptivity and reflectivity become low in the visible region. Aside, transparency decreases compared to the pure material. The optical properties were studied and clarified by electrons-photons interband transitions along with the complex dielectric function’s imaginary function.

Band structure engineering would open new avenues in device applications [10]. Doping with particular elements in practical applications can boost the material's electrical and optical properties and maximize the devices' efficiency. Mostly, α-Ga 2 O 3 is privileged with good n-type doping conductivity due to its deep donor oxygen defect nature [11]. Silicon (Si) and Tin (Sn) are highly studied n-type shallow donors with a small ionization energy [11,12]. In contrast, high-quality p-type α-Ga 2 O 3 is challenging due to the lack of hole conductivity, low solubility, and activation rate of dopants which largely impede the formation of the p-n junction bipolar semiconductor devices [13][14][15][16]. Despite several experimental reports on α-Ga 2 O 3 , its theoretical analysis is relatively small and started to receive attention recently. The theoretically predicted Gr-I and Gr-II metals shows p-type nature by doping on Ga 2 O 3 as reported by Tang et al. [17]. In 2019, Dong et al. investigated the magnesium (Mg) element substitute at the gallium site that had proven to be a good p-type conductivity in α-Ga 2 O 3 through first-principles studies [18]. Tang et al. reported that the same alkaline earth metal group, calcium (Ca), plays a deep acceptor element and studied the electronic structures and optical properties of p-type Ca doped on a β-Ga 2 O 3 first-principle calculation [17]. There is a lack of in-depth analysis and systematic study of new strong acceptor doping studies for the bipolar p-n junction of electronic and optical α-Ga 2 O 3 devices, which is of high demand for further exploration of potential α-Ga 2 O 3 films.
In this research, Gr-IIA alkaline metal Ca doping α-Ga 2 O 3 had reported the electronic structure and optical properties by the first-principles studies based on density functional theory (DFT). The work aims to find Ca potential by substituting the Ga site and supporting this theoretical investigation for future experimental work. Ca attributes of small ionization energy quickly release holes and show p-type conductivity nature. The doping of Ca contributes to a slight reduction in the optical bandgap due to the addition of the impurity band into the valance band (VB). The optical absorption spectra have shown a redshift of the absorption edge towards the visible-infrared region.

Computational Details
In this study, all the calculations were performed by using the Cambridge Serial Total Energy Package (CASTEP) code in Materials Studio (MS) 6.1 software based on the DFT [19]. The ultra-soft pseudopotential approach is used for general energy. The electronic wave function is unfolded with plane wave base groups, and ion potential is replaced by ultra-soft pseudopotential [20]. Ultrasoft pseudopotentials technology interacts with ions and electrons with total energy convergence 5 × 10 −6 eV/atom [21]. Cut-off plane-wave energy is set at 380 eV for total-energy, band structure, and spectra calculations. It has been shown that this functional exchange-correlation gives rise to widegap semiconductor atomic geometries with strong ionic bonds in excellent accordance with measured values. The exchange-correlation potential is defined by the generalized gradient approximation (GGA) in the Perdew-Burke-Ernzerhof (PBE) function and the local density approximation (LDA) [22]. The LDA utilization in this simulation work was used to compare the parameters, cell volume, and cell angle from the feature of the GGA-PBE. The electronic and optical characterization is then determined based on the structural geometry optimization to make the calculations effective and trackable [18]. The optimization of all calculated structures was done until the Hellmann Feynman force became less than 0.01 eV/Å (1 Å = 0.1 nm). Broyden-Fletcher-Goldfarb Shanno (BFGS) is a reliable approach to optimize the structure by full relaxation of the lattice constant and internal coordinates [23]. For maximum stress and displacement, the convergence criterion was 0.02 GPa and 5 × 10 −4 Å, respectively. A primitive cell composed of 30 atoms (12 Ga atoms and 18 O atoms) was used to study the electronic structure and optical properties of the Ca doped α-Ga 2 O 3 . A 3 × 3 × 2 k-point mesh of Monkhorst-Pack methods was used to integrate the Brillouin zone [24]. Ga, O, and Ca valence shell electrons configuration were [Ar] 3d 10 4s 2 4p 1 , [He] 2s 2 2p 4 , and [Ar] 4s 2 , respectively.

Geometrical Structures
Firstly, geometry optimization was implemented using the BFGS method [25]. The most stable structure was achieved after the fully relaxed cell parameter and volume by reducing the energy state [20]. The investigated optimized crystal structure of intrinsic α-Ga 2 O 3 and Ca-doped α-Ga 2 O 3 is shown in Figure 1a,b. The lattice parameters a, b, c, and unit cell volume V are increased after Ca doping due to the larger atomic radius of Ca than that of Ga. Meanwhile, the bandgap E g decreases with Ca doping. Ga2O3 and Ca-doped α-Ga2O3 is shown in Figure 1a,b. The lattice parameters a, b, c, and unit cell volume V are increased after Ca doping due to the larger atomic radius of Ca than that of Ga. Meanwhile, the bandgap Eg decreases with Ca doping.  Figure 1 shows the typical corundum structure of α-Ga2O3 with R3 c symmetry, and the crystalline cell comprises six Ga2O3 formula units. In the unit cell of α-Ga2O3, oxygen anions are almost closely hexagonally packed, and gallium atoms occupy octahedral twothird sites. Every gallium octahedron shares one face and three edges with the other three octahedra sites. The gallium octahedral sites are distorted by the crystal lattice [7]. We investigated and explored the electronic and optical properties of Ca doped α-Ga2O3 by substitution of Ga atoms lowest energy preferred sites. Table 1 lists all the optimized lattice parameters, and cell volume of pure and Ca doped α-Ga2O3 material as compared to experimental results (a = b = 4.983, c = 13.433) [26]. The experimental value for α-Ga2O3 has a lattice parameter and cell volume closer to that of the LDA calculation. In addition, the error is less than 1.5%, ensuring the reliability of our results. It has been shown that the calculated LDA is more precise and accurate than GGA-PBE according to the experimental reports on the lattice parameter of hexagonally structured α-Ga2O3. Therefore, in this study, the electronic structure and optical properties of Ca doped α-Ga2O3 were calculated using the DFT-LDA approach. The bond lengths for pure and Ca-doped α-Ga2O3 are shown in Table 2. It is shown that the atomic radius of the Ga atom is smaller than that of the Ca atom according to the periodic table. Therefore, the Ga-O bond length is shorter than the Ca-O and O-O bonds in Ca-doped α-Ga2O3. This situation is the same as the pure α-Ga2O3, where the Ga-O bond length is shorter than O-O bonds. After Ca doping, the overall bond length is larger than before doping due to the concentration of Ca 2+ ions.  Figure 1 shows the typical corundum structure of α-Ga 2 O 3 with R3c symmetry, and the crystalline cell comprises six Ga 2 O 3 formula units. In the unit cell of α-Ga 2 O 3, oxygen anions are almost closely hexagonally packed, and gallium atoms occupy octahedral twothird sites. Every gallium octahedron shares one face and three edges with the other three octahedra sites. The gallium octahedral sites are distorted by the crystal lattice [7]. We investigated and explored the electronic and optical properties of Ca doped α-Ga 2 O 3 by substitution of Ga atoms lowest energy preferred sites. Table 1 lists all the optimized lattice parameters, and cell volume of pure and Ca doped α-Ga 2 O 3 material as compared to experimental results (a = b = 4.983, c = 13.433) [26]. The experimental value for α-Ga 2 O 3 has a lattice parameter and cell volume closer to that of the LDA calculation. In addition, the error is less than 1.5%, ensuring the reliability of our results. It has been shown that the calculated LDA is more precise and accurate than GGA-PBE according to the experimental reports on the lattice parameter of hexagonally structured α-Ga 2 O 3 . Therefore, in this study, the electronic structure and optical properties of Ca doped α-Ga 2 O 3 were calculated using the DFT-LDA approach. Table 1. The lattice parameters (a, b, and c) and unit cell volume (V) of basic α-Ga 2 O 3 and Ca doped α-Ga 2 O 3 crystal structure. The values were compared with theoretical and experimental results.

Phase a (Å) b (Å) c (Å) V (Å 3 )
Pure The bond lengths for pure and Ca-doped α-Ga 2 O 3 are shown in Table 2. It is shown that the atomic radius of the Ga atom is smaller than that of the Ca atom according to the periodic table. Therefore, the Ga-O bond length is shorter than the Ca-O and O-O bonds in Ca-doped α-Ga 2 O 3 . This situation is the same as the pure α-Ga 2 O 3, where the Ga-O bond length is shorter than O-O bonds. After Ca doping, the overall bond length is larger than before doping due to the concentration of Ca 2+ ions. The ionic and covalent bonding-based charge density contours of intrinsic α-Ga 2 O 3 and Ca-doped α-Ga 2 O 3 are shown in Figure 2a,b, respectively. The charge density distributions can be used to assess crystal bonding characteristics. In α-Ga 2 O 3 , the ionic or covalent structures are a contentious subject and have vital facts to understand the charge transport property [23,24]. The nuclei are bound by the charge density, which is shared between them. It has been shown that certain covalent bonding characters can exist between Ga and O atoms [24]. The electrons are located around the O atoms in the hexagonal phases, and there is no region of bond localization, indicating a predominantly ionic bond.

Type of Bond Length
Pure α-Ga2O3

Electron Charge Density
The ionic and covalent bonding-based charge density contours of intrinsic α-Ga2O3 and Ca-doped α-Ga2O3 are shown in Figure 2a,b, respectively. The charge density distributions can be used to assess crystal bonding characteristics. In α-Ga2O3, the ionic or covalent structures are a contentious subject and have vital facts to understand the charge transport property [23,24]. The nuclei are bound by the charge density, which is shared between them. It has been shown that certain covalent bonding characters can exist between Ga and O atoms [24]. The electrons are located around the O atoms in the hexagonal phases, and there is no region of bond localization, indicating a predominantly ionic bond.
In Figure 2b, a major reflection is observed in the distribution of charge density between the doping Ca atom and native atoms. By incorporating Ca, the oxygen atom's nearest electrons population reduce compared to distant oxygen atoms, as presented in   In Figure 2b, a major reflection is observed in the distribution of charge density between the doping Ca atom and native atoms. By incorporating Ca, the oxygen atom's nearest electrons population reduce compared to distant oxygen atoms, as presented in

Band Structure and Density of States
The first-principal calculations of the electronic band structure of undoped α-Ga 2 O 3 are shown in Figure 3a. The density of states is shown in Figure 4a,c at different energy windows. The bandgap energy is usually determined between the valence band maximum (VBM) and the conduction band minimum (CBM), placed at the Brillion region's G point. The Fermi level is situated at 0 eV of the energy scale. The large disparity of theoretically calculated DFT-LDA bandgap (2.950 eV) and experimental band gap (5.3 eV) is due to the ground state DFT results of strong Coulomb correlation and exchange-correlation potential between the excited electrons in the underestimated LDA bandgap [27]. pant has been introduced above the valence band and intersected with the Fermi level. The Ca-α-Ga2O3 has a direct bandgap with the most upper part of the valence band and the lowest part of the G point's conduction band. It is noticed that the DFT-LDA calculated bandgap of Ca-doped α-Ga2O3 is 2.733 eV which is 0.217 eV less than that of pure α-Ga2O3. In Ca-doped systems, the Ca 3p state contributes to the VBM, and the CBM is contributed to by Ga-4s and Ca-3p states. Ca creates dopant states near the topmost of VBs and shift the valence band edge at higher energy and reduce the bandgap.   (Ga1−xCax)2O3 alloys bandgap energy with the incorporation of Ca impurity is schematically shown in Figure 5. In α-Ca2O3, 3.33% of Ca concentrations were produced by replacing one Ga atom with Ca, respectively. Ca substitutes Ga octahedral sites and fully occupied band energy higher than the original valence band of the Ga2O3. This intermediate band originates from the coupling of Ca 3p, which lies below the O 2p band that The electronic band structure of pure α-Ga 2 O 3 , as shown in Figure 3a, is represented by the crystal structure made with the 12 Ga atoms and 18 O atoms of primitive cells. Figure 4a shows the uppermost valence band is a major contribution of O 2p states with a width from about −0.12 to −7.43 eV and a minor contribution of Ga 3d, 4p, 4s, and O 2s orbitals which are located between −11.5 and −18.5 eV. CBM consists of Ga 4s orbital states, and lies between 2.32 to 8.46 eV [26]. The valence band's edge is reasonably flat, and the effective mass at VBM is relatively heavy. Characteristics of CBM is almost equally dispersed everywhere.
The electronic band structure of Ca-doped α-Ga 2 O 3 is shown in Figure 3b, and its corresponding DOS is shown in Figure 4b

Theoretical Description and Optical Properties
The optical properties of pure α-Ga2O3 and Ca doped α-Ga2O3 can be studied based

Theoretical Description and Optical Properties
The optical properties of pure α-Ga 2 O 3 and Ca doped α-Ga 2 O 3 can be studied based on the dielectric function, absorption coefficient, and reflectivity. The optical transitions are related to inter-band and intra-band transitions in the crystal structure of solid materials. The inter-band transitions mainly contributed to the semiconductor rather than intraband transitions to semiconducting crystal structure. The dielectric function ε(ω) has a frequency dependence linear electromagnetic interaction response of incident photons with electrons and is separated into the real part ε 1 (ω) and imaginary part iε 2 (ω). It is difficult to achieve accurate optical constants due to the underestimation of the bandgap. Hence, the optical property results were corrected by scissors operators according to experimental results. Scissors operators = 5.3 eV−2.950 eV = 2.35 eV, where 2.950 is theoritical DFT-LDA calculated bandgap and 5.3 eV is experimental bandgap [28]. The scissors operators used to shift all conduction levels to match the measured value of the bandgap. All-optical properties can be described from the complex dielectric constant in Equation (1) [29]: The imaginary dielectric part ε 2 (ω) semiconductors are an essential parameter for determining optical properties for different designs of optoelectronic devices. The dielectric function imaginary part ε 2 (ω) is associated with the dielectric loss of energy or absorption of light described by the absorption coefficient in the material while the real part ε 1 (ω) is associated with the stored energy within the material (degree of polarization). The complex dielectric function imaginary part ε 2 (ω) is determined by summing the transitions between occupied and unoccupied electronic states, as reported in Equation (2) [22]: where m is the mass of free electrons, e is the electron charge, M is the dipole matrix, and ω is the frequency of incident photons. i and j are the initial and final states, respectively, f i denotes the Fermi distribution function and i-th state with wave function vector k. The dielectric function real part ε 1 (ω) can be derived from the imaginary part iε 2 (ω) by the Kramers-Kronig dispersion Equation ( where P is the principal value of the integral. The absorption coefficient can be derived from the dielectric function as shown in Equation (4) [31]:

Dielectric Function
The dielectric constant is directly proportional to the crystal's polarizability, representing the deformability of the electronic distribution, and relates the shape of the valence band charge density. Oxygen ions are supposed to contribute to polarizability in such a highly ionic substance as Ga2O3. As was predicted in the current study, a difference of dielectric properties between Ca and α-Ga2O3 is projected. A considerably larger dielectric constant of α-Ga2O3 compared to Ca doped α-Ga2O3 can be drawn to a lightly dense ion package and a large number of Ga atoms coordinated in the α-Ga2O3 lattice [7].
From Figure 6, the calculated imaginary dielectric function of intrinsic α-Ga2O3 and Ca-α-Ga2O3 were plotted and compared. The main peaks are 10.12 eV of pure and Cadoped α-Ga2O3, respectively, corresponding to the electrons' transitions from the upper state O-2p valence band to the lower state Ga-4s conduction band. The Ca-doped α-Ga2O3 peak is smaller than that of intrinsic α-Ga2O3 and is enhanced towards the lower energy site. Moreover, this Ca doping absorption peak at 2.9 eV indicates that the absorption spectra increase from the UV to the visible region, and the average optical transmittance (3) where P is the principal value of the integral. The absorption coefficient can be derived from the dielectric function as shown in Equation (4) [31]:

Dielectric Function
The dielectric constant is directly proportional to the crystal's polarizability, representing the deformability of the electronic distribution, and relates the shape of the valence band charge density. Oxygen ions are supposed to contribute to polarizability in such a highly ionic substance as Ga 2 O 3 . As was predicted in the current study, a difference of dielectric properties between Ca and α-Ga 2 O 3 is projected. A considerably larger dielectric constant of α-Ga 2 O 3 compared to Ca doped α-Ga 2 O 3 can be drawn to a lightly dense ion package and a large number of Ga atoms coordinated in the α-Ga 2 O 3 lattice [7].
From Figure 6, the calculated imaginary dielectric function of intrinsic α-Ga 2 O 3 and Ca-α-Ga 2 O 3 were plotted and compared. The main peaks are 10.12 eV of pure and Cadoped α-Ga 2 O 3 , respectively, corresponding to the electrons' transitions from the upper state O-2p valence band to the lower state Ga-4s conduction band. The Ca-doped α-Ga 2 O 3 peak is smaller than that of intrinsic α-Ga 2 O 3 and is enhanced towards the lower energy site. Moreover, this Ca doping absorption peak at 2.9 eV indicates that the absorption spectra increase from the UV to the visible region, and the average optical transmittance decreases.
such a highly ionic substance as Ga2O3. As was predicted in the current study, a difference of dielectric properties between Ca and α-Ga2O3 is projected. A considerably larger dielectric constant of α-Ga2O3 compared to Ca doped α-Ga2O3 can be drawn to a lightly dense ion package and a large number of Ga atoms coordinated in the α-Ga2O3 lattice [7].
From Figure 6, the calculated imaginary dielectric function of intrinsic α-Ga2O3 and Ca-α-Ga2O3 were plotted and compared. The main peaks are 10.12 eV of pure and Cadoped α-Ga2O3, respectively, corresponding to the electrons' transitions from the upper state O-2p valence band to the lower state Ga-4s conduction band. The Ca-doped α-Ga2O3 peak is smaller than that of intrinsic α-Ga2O3 and is enhanced towards the lower energy site. Moreover, this Ca doping absorption peak at 2.9 eV indicates that the absorption spectra increase from the UV to the visible region, and the average optical transmittance decreases. Figure 6. The imaginary part of the complex dielectric function for intrinsic α-Ga2O3 and Ca doped α-Ga2O3.

Absorption and Reflectivity
The absorption spectra of intrinsic α-Ga2O3 and Ca-doped α-Ga2O3 in the range of energy (0 to 40 eV) and wavelength (0-1000 nm) are shown in Figure 7a,b. The presence of a high energy peak at 11.5 eV into the deep ultraviolet region is shown in both cases of α-Ga2O3 and Ca-doped α-Ga2O3. These peaks are associated with the inter-band transition from the valance band O 2p occupied states to the conduction band Ga 4s unoccupied states. From Figure 7a, it is worth noticing that the absorptivity of α-Ga2O3 is in the deep UV spectrum due to its semi-insulating nature. Meanwhile, Ca-doped α-Ga2O3 appears as

Absorption and Reflectivity
The absorption spectra of intrinsic α-Ga 2 O 3 and Ca-doped α-Ga 2 O 3 in the range of energy (0 to 40 eV) and wavelength (0-1000 nm) are shown in Figure 7a,b. The presence of a high energy peak at 11.5 eV into the deep ultraviolet region is shown in both cases of α-Ga 2 O 3 and Ca-doped α-Ga 2 O 3 . These peaks are associated with the inter-band transition from the valance band O 2p occupied states to the conduction band Ga 4s unoccupied states. From Figure 7a, it is worth noticing that the absorptivity of α-Ga 2 O 3 is in the deep UV spectrum due to its semi-insulating nature. Meanwhile, Ca-doped α-Ga 2 O 3 appears as a sharp new small peak at 2.9 eV which originated from the inter-band transitions of C 3p from the VBM to these CBM that promote ε 2 and enhanced absorption in the visible region. Figure 7a,b show that the absorption spectra have a redshift with the Ca-doping. This shifting is mainly because of the hole carriers generated after the Ca doping so that the Fermi level enters the valence band, which is realized in the density of the states. The same phenomenon was observed in the absorption spectrum wavelength window in both cases of undoped and Ca doped α-Ga 2 O 3 ( Figure 7b). As shown in Figure 7b, the absorption spectra indicate that the transparency decreases with Ca doping due to enhancement of absorbance in the visible region. Figure 7c,d show the reflection spectra of intrinsic α-Ga 2 O 3 and Ca-doped α-Ga 2 O 3 in the energy and wavelength window regime. For pure α-Ga 2 O 3 , a sharp main intense peak appears at 16 eV in the energy range and 70.45 nm in the wavelength range, and it decreases from the deep UV-visible to infrared region. Meanwhile, for the Ca doped α-Ga 2 O 3 system, two major peaks appear at the same position as the pure phase, but reflectivity is lower than that of the pure one. These characters drag material potential towards optoelectronics and bipolar device application. Figure 7c,d show the reflection spectra of intrinsic α-Ga2O3 and Ca-doped α-Ga2O3 in the energy and wavelength window regime. For pure α-Ga2O3, a sharp main intense peak appears at 16 eV in the energy range and 70.45 nm in the wavelength range, and it decreases from the deep UV-visible to infrared region. Meanwhile, for the Ca doped α-Ga2O3 system, two major peaks appear at the same position as the pure phase, but reflectivity is lower than that of the pure one. These characters drag material potential towards optoelectronics and bipolar device application. Figure 7. The absorption spectrum of (a) intrinsic α-Ga2O3 and Ca-doped α-Ga2O3 in energy, (b) intrinsic α-Ga2O3 and Cadoped α-Ga2O3 in wavelength, the reflectivity of (c) intrinsic α-Ga2O3 and Ca-doped α-Ga2O3 in energy, and (d) intrinsic α-Ga2O3 and Ca-doped α-Ga2O3 in wavelength window regime.

Conclusions
The substitutional Ca doping on the α-Ga 2 O 3 electronic, structural, and optical properties was studied through a first-principles calculation based on the DFT. The Ca doping predominantly affects the band structure and optical properties of pure α-Ga 2 O 3 . The DFT-LDA calculated bandgap of α-Ga 2 O 3 is decreased to 2.950 eV after Ca doping and the forbidden bandwidth appeared at the top of the valance band. The results obtained show that Ca-doped α-Ga 2 O 3 is an obvious variation of energy band structures and the charge density distribution. The acceptor energy states are composed of Ca 3p and O 2p orbital states. The Ca-doped α-Ga 2 O 3 dynamic dielectric function is observed. The most important absorption peak shifted in the visible regime; it realizes p-type α-Ga 2 O 3 but not a strong p-type candidate for optoelectronic device application. However, the lack of experimental results about the optical properties of Ca doped α-Ga 2 O 3 requires further experimental works to compare with our simulation results.