Saharat
Chomdech
a,
Chalermchai
Himwas
*a,
Wenich
Pumee
a,
Suphakan
Kijamnajsuk
b,
Waraporn
Tanthanuch
c and
Songphol
Kanjanachuchai
*a
aSemiconductor Device Research Laboratory, Department of Electrical Engineering, Faculty of Engineering, Chulalongkorn University, 254 Phayathai Road, Bangkok 10330, Thailand. E-mail: Chalermchai.H@chula.ac.th
bNational Metal and Materials Technology Center, Thailand Science Park, 114 Paholyothin Rd., Klong 1, Klong Luang, Pathumthani 12120, Thailand
cSynchrotron Light Research Institute (Public Organization), Muang, Nakhon Ratchasima 30000, Thailand
First published on 13th August 2024
GaAsPBi is a comparatively novel quaternary III–V compound semiconductor with attractive properties that suit optoelectronic applications. However, the quaternary compound synthesized on GaAs (001) is not straightforward and requires elaborately adjusted growth parameters. The quaternary compound usually presents inhomogeneity and defects, making the compositional estimation of the quaternary compound challenging. This work proposes a method to estimate the quaternary compound using synchrotron-based extended X-ray absorption fine structure (EXAFS) around the P-K edge. This technique considers only bound atoms in the local structure, making it suitable for estimating the composition of both defective and homogeneous compounds. From the EXAFS estimations, the average composition of defective GaAsPBi (homogeneous GaAsPBi) is GaAs0.24P0.68Bi0.08 (GaAs0.66P0.25Bi0.09). The oxidation states of Bi and P in the quaternary compounds are also probed by synchrotron-based X-ray absorption near edge structure. Finally, we show that the ex situ annealing of lattice-matched GaAsPBi epitaxial layer on GaAs (001) under the flow of hydrogen-containing gas can only slightly improve its luminescence (up to 1.3 times as compared to the as-grown condition) owing to the low defect density.
One way to increase the Bi concentration without sacrificing the crystal quality is to synthesize the GaAsPBi epitaxial layer instead of the GaAsBi. However, the synthesis of GaAsPBi is not straightforward. There is a catalytic effect of the Bi atoms towards the existing P atoms, making it hard to control the alloy composition. Furthermore, orthogonal dislocations and running droplets readily occur if grown without using the reported strategic procedure.16 Another challenge when incorporating more elemental atoms in an epitaxial layer is the uncertainty in the compositional estimation. Nattermann et al. estimated the GaAsPBi alloy composition by comparing the optical results against the structural model17 when our earlier works relied on synchrotron-based X-ray photoelectron spectroscopy (XPS) for the quaternary alloy compositional estimation under a cation–anion stoichiometric assumption.16,18 XPS is a surface-sensitive technique that measures all ejected electrons from the surface, hence unsuitable for probing inhomogeneous epitaxial layers. Besides, Bi-containing materials grown at low temperatures usually host more defects/dislocations if the growth parameters are not elaborately adjusted. These defects/dislocations invalidate the carriers participating in the radiative luminescent mechanism. Thermal annealing in hydrogen-containing gas aims to permeate hydrogen atoms into the epilayer to neutralize these defects/dislocations and free up carriers, thus improving luminescent intensity.19,20 Since high-temperature device fabrication processes can easily disrupt the epitaxial layer and degrade the device's quality, another benefit of the annealing observations is to set the maximum temperature that the epitaxial can tolerate. Here, we propose a means to estimate the compositions of the GaAsPBi epitaxial layers and study the ex situ annealing of the lattice-matched GaAsPBi epitaxial layer. This work reports the study of X-ray absorption near-edge structures (XANES) and compositional estimations using extended X-ray absorption fine structures (EXAFS) for defective and lattice-matched GaAsPBi epitaxial layers by implying information obtained from surface and structural investigations. This innovative technique considers only bound atoms in the local structure, making it suitable for estimating the composition of both defective and homogeneous compounds. Subsequently, we report the effect of ex situ annealing on the structural properties, optical properties, XANES spectra, and EXAFS spectra of the lattice-matched GaAsPBi on GaAs (001).
The surfaces of samples A and B were measured using an atomic force microscope (AFM, SPA-400, Seiko Instruments) and a differential interference contrast (DIC) optical microscope (OM). The structural properties of both samples were probed by high-resolution X-ray diffraction (HR-XRD, Rigaku). The chemical formations were characterized by advanced X-ray absorption spectroscopy (XAS). XAS spectra were acquired at the SUT-NANOTEC-SLRI XAS beamline (BL5.2) of the SLRI, Thailand. The X-ray light has an energy ranging from 1240–12100 eV and a photon flux ranging from 108–1010 photons per s/100 mA. An InSb (111) double-crystal monochromator was used for energy scanning. The Bi–M5 edge XANES calibration was performed using a Mo foil in transmission mode. Spectral acquisition was conducted from −30 eV to 100 eV relative to E0 at 2580 eV, with an energy step of 0.2. For P-K edge XAS measurement, the elemental phosphorus standard was used for calibration. P-K edge XANE acquisition covered a range from −20 eV to 90 eV of E0 at 2146 eV with an energy step of 0.2. In P-K edge EXAFS spectral acquisitions, the energy scan was set up in a range of −100, −20, 30, and 10.6 k of E0 at 2146 eV with an energy step of 5, 0.3, 0.05 k, with a corresponding time step of 1, 1, and 3 s. Data preprocessing and theoretical fitting were performed consecutively using Athena-IFEFFIT and Artemis-IFEFFIT packages.21
For the ex situ annealing study of lattice-matched GaAsPBi, sample B is cut into four pieces, then annealed for ∼3 hours under the flow of hydrogen-containing gas at different annealing temperatures (TA): 450 °C (sample C), 550 °C (sample D), and 650 °C (sample E). The structural properties of samples B, C, D, and E were measured using AFM, OM, and high-resolution XRD. The optical quality of the samples was assessed by photoluminescence at 30 K. The samples were mounted on the cold finger of a close cycle He cryostat and excited by Ventus solid-state laser (λ = 532 nm) at an average power density of 50 W cm−2. Focused PL emitted from the samples was dispersed in a monochromator (Horiba iHR320) and collected by a cooled (77 K) InGaAs detector using a standard lock-in detection technique. PL intensity acquired for each sample was normalized against a standard sample to allow meaningful comparisons. Finally, we report the effect of ex situ annealing on the XAS spectra.
We use X-ray absorption spectroscopy (XAS) to characterize the structures of the quaternary alloy. Depending on the energy range extended from the absorption edge of an atom, XAS is categorized into XANES and EXAFS. The sample's thickness exceeding 300 μm, restricts the measurement of XAS spectra to the fluorescent mode. We skip characterizing the absorption edges for Ga and As atoms since the measuring signal potentially comes from the substrate, disturbing the fitting and interpretation. Furthermore, the beamline 5.2 configuration gives the highest photon flux within a limited energy range, making it suitable for study only at Bi M5- and P K-edges for this experiment.
The oxidation states of the synthesized quaternary alloy are measured using XANES at the Bi M5-edge and P K-edge and compared against their known standards. Fig. 2(a) shows Bi M5-edges for sample A (E0 = 2599 eV) and for sample B (E0 = 2601 eV). Additionally, Bi M5-edge for Bi2O3 (E0 = 2598 eV), which theoretically corresponds to Bi3+, is added in the figure for comparison. Since the absorption edges of both samples are close to that of Bi2O3, the oxidation state of Bi in samples A and B is likely 3+. Fig. 2(b) shows P K-edges for samples A and B (E0 = 2146 eV). Additionally, the P K-edges for Ca3(PO4)2 (E0 = 2154 eV), KH2PO4 (E0 = 2154 eV), and phosphorus powder (P4, E0 = 2146 eV) are included in the figure for comparison. Theoretical assignments suggest that these correspond to P5+, P5+, and 0 oxidation states, respectively. It indicates that the synthesized quaternary alloy possesses a phosphorus oxidation state of 0, as the absorption energy aligns with that of P4.
For the EXAFS, the measuring spectrum in terms of energy is converted to k-space, then to R-space by applying a Fourier transform. The R-space represents the radial distribution of neighboring atoms around the absorbing atom.
The k-space and R-space are fitted with a known model that includes parameters such as bond lengths, coordination numbers, and the Debye–Waller factor (σ2). The fitting results consider these parameters. The shifted energy (ΔE0) should be within ±10 eV for well-fitted EXAFS. The coordination number represents the number of neighboring atoms around the absorber. Bond length (R) is the distance between the absorbing atom and its neighboring atoms. σ2 value represents the Debye–Waller factor or mean square relative displacement (MSRD), which accounts for the thermal and static disorder in the distances between the absorbing atom and its neighboring atoms. Typical values of σ2 for well-ordered structured structures are on the order of 0.001–0.01, while higher values indicate more disorder or higher temperature.21,23 The amplitude reduction factor (S02) adjusts the overall amplitude of the EXAFS oscillations, which typically ranges between 0.7 and 1.0.21,24R-Factor quantifies the agreement between the fitted model and the experimental data, aiming for a low value to indicate close agreement between the experimental and theoretical EXAFS spectra.25
We characterize the local structures of samples A and B by using synchrotron-based EXAFS around the P K-edge. EXAFS around the Bi M5-edge is neglected because its energy is too close to the Bi M4-edge (<250 eV difference). Fig. 3(a) and (b) consecutively depict Fourier transformed of P K-edge EXAFS spectra (black spectra) for samples A and B. Fitting the spectra in R space allows revealing the local structures around the P atoms. For the fitting, we initially generate the theoretical EXAFS signal by the FEFF code performed on the cubic GaP crystal (prevalently experimentally observed structure). We execute the fitting using the theoretical signal in conjunction with the lattice constant value resulting from the X-ray diffractogram shown in Fig. 1(e).
Fig. 3 Experimental and simulated Fourier transformed P K-edge EXAFS of GaAsPBi epitaxial layers of (a) sample A and (b) sample B. |
For sample A, highly mismatched probed by XRD, cracks, and orthogonal dislocations evidenced in AFM attest to perform the diffractogram/lattice constant conversion under a fully relaxed assumption. The theoretical EXAFS signal assumes that the quaternary alloy is stoichiometric, i.e., P, As, and Bi atoms (anions) are interchangeable and are balanced by Ga atoms (cations). The Ga atoms in the first shell are intact, and the P atoms in the second shell are interchangeable with As and Bi atoms until the experimental and fitting spectra converge (red scatters in Fig. 3(a)). For sample A, the best fit is obtained when the second shell comprises 2.83 As atoms, 8.18 P atoms, and 0.99 Bi atoms (equivalent to GaAs0.24P0.68Bi0.08). The resulting high-P-content quaternary alloy agrees with the XRD diffractogram. This composition, which measures the local structure of the four atomic types binding together, represents an average composition of the quaternary alloy. The discrepancy of the composition acquired using EXAFS and XPS is associated with the nature of the acquired signal. The composition estimated using XPS spectroscopy relies on all the atoms that appeared on the surface. Hence, the method is suitable only for homogeneous alloys.
For sample B, the adjacent diffraction peaks of the substrate and epitaxial layer (Fig. 1(e)) and the smooth AFM surface (Fig. 1(c)) suggest that the grown layer is pseudo-morphically synthesized on the substrate. We obtain the best fit (scatters in Fig. 3(b)) by applying the identical fitting technique reported earlier when the second shell comprises 7.95 As atoms, 3.00 P atoms, and 1.05 Bi atom (equivalent to GaAs0.66P0.25Bi0.09).
To study the ex situ annealing effect of the lattice-matched GaAsPBi on GaAs (001), we cut sample B into four pieces, kept one piece as a controlled sample, and annealed the other three samples at 450 °C (sample C), 550 °C (sample D), and 650 °C (sample E). To anneal each sample, we placed it in the middle of a quartz tube, increased the TA to the set point for one hour, dwelled at the set point for one hour, and decreased the TA to room temperature for another hour. Thus, the complete annealing profile lasts three hours under hydrogen-containing gas (95% N2, 5% H2). The dissociated hydrogen atoms can neutralize dislocations and dangling bonds and free the carriers, thus improving the optical emission. Indeed, we annealed the samples using the identical chamber and annealing profiles as in the case of reported dislocations embedded GaAsPBi epilayer,26 allowing us to compare the annealing effects from both studies.
Fig. 4(a)–(f) show the AFM micrographs (OM images) of latticed-matched films annealed at 450, 550, and 650 °C, respectively. The AFM images confirm that the surfaces of all the annealed samples remain unchanged from the as-grown surface (Fig. 1(c)). Their RMS roughness is in the range of 0.36–0.62 nm. For a lower magnification field of view, the OM images of annealed samples at TA = 450 °C and 550 °C are similar to that of the as-grown sample, in agreement with the AFM results. Blistering defects emerge on the surface of sample E (TA = 650 °C) (see Fig. 4(f)). The low density of the defects hinders us from observing them by AFM. The average diameter and density of the blisters are 41.65 μm and 1.31 × 10−5 μm−2, respectively. The defects distributed randomly for the lattice-matched GaAsPBi/GaAs, differently from those that occurred on the surface with orthogonal dislocations reported elsewhere.26
A 90°-tilted SEM micrograph of sample D (TA = 550 °C, Fig. 4(g)) illustrates a smooth surface of the 280 nm-thick GaAsPBi epitaxial layer. The measured thickness is thicker than the 200 nm-thick prior determined by RHEED oscillations, owing to the different incorporating mechanisms during the growth and calibration explained earlier.18 90°-tilted SEM micrographs of samples B (as-grown) and C (TA = 450 °C), having the same structure as sample D, are not shown here. Fig. 4(h) shows a 90°-tilted SEM micrograph of sample E (TA = 450 °C) with a blistering defect on the surface. This blistering defect has a depth of 3 μm, indicating that it destroys both the epitaxial layer and GaAs substrate.
HR-XRD scans around the symmetric (004) reflections of the annealed GaAsPBi epitaxial layers are depicted in Fig. 5. For samples C (TA = 450 °C) and D (TA = 550 °C), the GaAs and the quaternary alloy diffraction peaks are identical to those of the as-grown (sample B), confirming that the structures after annealing at TA = 450 °C and 550 °C are intact. The structural stabilities agree with the surface results shown earlier. For sample E (TA = 650 °C), the 2θ diffraction peak of GaAs is stable at 66.05° but that of GaAsPBi shifts to 66.36°. The XRD peak shift confirms that the structural degradation occurred at TA = 650 °C, agreeing with the reported surface results. The XRD peak position suggests that the degraded film has comparatively lower As or Bi concentrations than other samples in the annealing series.
Fig. 5 HR-XRD diffractograms scanned around the symmetric (004) reflections of as-grown GaAsPBi epitaxial layers (samples B) and of annealed GaAsPBi epitaxial layers (samples C, D, and E). |
Fig. 6(a) shows normalized PL spectra of samples B (as-grown), C (TA = 450 °C), D (TA = 550 °C), and E (TA = 650 °C). The PL spectra of all the samples show two luminescent peaks with energies centered around 1.32 eV and 1.47 eV, emitting from the epitaxial layer and GaAs substrate, respectively. Single Gaussian luminescence from the epitaxial layer suggests that the epitaxial layer is homogeneous, unlike sample A, which possesses multiple peak luminescence originating from the in-plane inhomogeneity.26 The optical improvement of the annealing series is analyzed by comparing their normalized integrated photoluminescence (IPL) with the IPL of sample B (as-grown) (see Fig. 6(b), black bar graph). For comparison, we also summarize the optical improvement in the case of the annealed series of the inhomogeneous GaAsPBi epitaxial layers from ref. 26 (see Fig. 6(b), gray bar graph). For lattice-matched GaAsPBi/GaAs, annealing with hydrogen-containing gas can only slightly enhance the luminescence. The most luminescent improvement occurs at TA = 550 °C, where the IPL is improved 1.3 times compared to the as-grown sample. The optimal annealing temperature (TA = 550 °C) is identical to the case of the inhomogeneous quaternary alloy in that the IPL is improved 34 times from the as-grown condition that occurred along with the luminescent spectral feature modification.26
Fig. 6 (a) Normalized PL spectra of GaAsPBi epitaxial layers (samples B, C, D, and E). (b) Normalized integrated PL (IPL) of GaAsPBi epitaxial layers (samples B, C, D, and E) as a function of annealing temperature (TA) (black bar graph). The gray bar graph represents the normalized IPL of defective GaAsPBi epitaxial layers as a function of TA taken from ref. 26. |
Since the improvement comes mainly from the carriers released from the neutralized defects, the current annealing recipe can significantly improve the optical quality for the sample inherited with high defect density, for example, high orthogonal dislocation density as sample A. When the annealing temperature increases to TA = 650 °C, the IPL decreases to 0.07 times compared to the as-grown sample. The reduction of the luminescence agrees with the degraded structural quality. After the annealing process, the XANES signal for all the annealed samples is similar to that for sample B (not shown). The oxidation state of Bi (P) for all the annealed samples C, D, and E is 3+ (0) since the Bi M5-edges (P K-edges) for the samples are at E0 = 2601 eV (2146 eV).
Fig. 7(a)–(c) consecutively depict Fourier transformed of P K-edge EXAFS spectra (black spectra) for the annealed GaAsPBi samples C, D, and E. We execute the fitting using the theoretical signal with the lattice constant value resulting from the X-ray diffractogram shown in Fig. 5 – the identical fitting technique used earlier. The best fits (scatters in Fig. 7(a)–(c)) are obtained when the compositions of samples C, D, and E are equivalent to GaAs0.66P0.25Bi0.09. The EXAFS fitting parameters for samples A, B, C, D, and E are listed in Table 1.
Fig. 7 Experimental and simulated Fourier transformed P K-edge EXAFS of annealed GaAsPBi samples (a) sample C (TA = 450 °C), (b) sample D (TA = 550 °C), and (c) sample E (TA = 650 °C). |
Path | Sample A | Sample B | Sample C | ||||||
---|---|---|---|---|---|---|---|---|---|
N | R (Å) | σ 2 | N | R (Å) | σ 2 | N | R (Å) | σ 2 | |
P–Ga | 4.00 | 2.362 | 0.004 | 4.00 | 2.361 | 0.002 | 4.00 | 2.362 | 0.003 |
P–P | 8.18 | 3.901 | 0.028 | 3.00 | 3.873 | 0.005 | 1.74 | 3.786 | 0.008 |
P–As | 2.83 | 3.932 | 0.009 | 7.95 | 3.925 | 0.007 | 9.84 | 3.965 | 0.009 |
P–Bi | 0.99 | 3.706 | 0.015 | 1.05 | 3.772 | 0.013 | 1.23 | 3.885 | 0.003 |
P–Ga | 12.00 | 4.546 | 0.018 | 12.00 | 4.593 | 0.027 | 12.00 | 4.627 | 0.005 |
R factor | 0.018 | 0.045 | 0.044 | ||||||
S 0 2 | 0.929 | 0.963 | 0.963 | ||||||
ΔE | 3.759 | 3.256 | 3.860 |
Path | Sample D | Sample E | |||||||
---|---|---|---|---|---|---|---|---|---|
N | R (Å) | σ 2 | N | R (Å) | σ 2 | ||||
P–Ga | 4.00 | 2.358 | 0.0031 | 4.00 | 2.361 | 0.0033 | |||
P–P | 2.64 | 3.830 | 0.017 | 2.64 | 3.813 | 0.0022 | |||
P–As | 8.99 | 3.936 | 0.0105 | 8.84 | 3.951 | 0.0105 | |||
P–Bi | 0.30 | 4.065 | 0.0006 | 0.53 | 4.088 | 0.0002 | |||
P–Ga | 12.00 | 4.624 | 0.0129 | 12.00 | 4.616 | 0.0125 | |||
R factor | 0.049 | 0.059 | |||||||
S 0 2 | 0.963 | 0.963 | |||||||
ΔE | 2.867 | 3.690 |
The room temperature band gap of GaAs0.66P0.25Bi0.09 calculated by equations proposed by Luo et al.27 is 1.17 eV. The GaAs0.66P0.25Bi0.09 band gap at 25 K is estimated at 1.27 eV (an increase of ∼0.1 eV as the temperature decreases from 300 K to 25 K, the same as GaAs). The marked band gap position in Fig. 5(a) is within the luminescent curve but still lower than the experimental luminescent peak of 50 meV. The discrepancy is attributed either to inappropriate parameters in the fitting packages or to the equations proposed by Luo et al.,27 resulting in a slight overestimation of the Bi content (or underestimation of the P content). The local structure of the lattice-matched GaAsPBi alloy does not change even at the highest annealing temperature (TA = 650 °C), where we observe blisters on the surface and degraded luminescent quality, implying that the degradation occurred macroscopically while the local structure was intact. Here, we do not estimate the band gap of GaAs0.24P0.68Bi0.08 for two reasons. First, the P concentration (0.68) is out of the application range for the equations proposed by Luo et al. (0.40).27 Second, the sample is inhomogeneous, which makes the comparison between the optical and structural asserted band gaps invalid.
This journal is © The Royal Society of Chemistry 2024 |