InPBi single crystals grown by molecular beam epitaxy.

InPBi was predicted to be the most robust infrared optoelectronic material but also the most difficult to synthesize within In-VBi (V = P, As and Sb) 25 years ago. We report the first successful growth of InPBi single crystals with Bi concentration far beyond the doping level by gas source molecular beam epitaxy. The InPBi thin films reveal excellent surface, structural and optical qualities making it a promising new III-V compound family member for heterostructures. The Bi concentration is found to be 2.4 ± 0.4% with 94 ± 5% Bi atoms at substitutional sites. Optical absorption indicates a band gap of 1.23 eV at room temperature while photoluminescence shows unexpectedly strong and broad light emission at 1.4-2.7 μm which can't be explained by the existing theory.

Compound semiconductors n class="Chemical">play a pivotal role in modern optoelectronics whose energy gap and lattice constant can be tailored to meet various device requirements. While semiconductors containing n class="Chemical">N, P, As and Sb in the column V of the periodic table have been extensively investigated, bismuth (Bi) containing alloys are less explored. Isoelectronic impurities in host materials, such as N and Bi in GaP and InP at the doping level, reveal rich spectroscopic features near the band gap of the host material1234. When the impurity concentration increases from the doping level to the dilute alloy or even the concentrated alloy level, interaction between the impurity energy level and the band edge of the host material results in strong band gap bowing effect in dilute nitrides56 and also in dilute GaAsBi7. This property, which can be well explained by the band anti-crossing model8, has been employed not only to demonstrate uncooled high-speed telecom lasers on GaAs9 and the record high conversion efficiency of triple junction solar cells both using dilute nitrides10, but also to trigger great intention to search for relatively unexplored dilute bismides for band engineering. Berding et al. theoretically analyzed InSbBi, InAsBi and InPBi as potential candidates for mid- and far-infrared (IR) optoelectronics applications in 1988 and pointed out that InSbBi was the easiest and InPBi was the most difficult to mix, but once formed, InPBi was the hardest among the three materials11, making it an attractive narrow-gap material candidate. Since InBi crystallizes in the tetragonal PbO-type structure rather than the cubic zinc-blende structure common for semiconducting III–V compounds, the mixing of InBi with other materials of zinc-blende structure may lead to problems12. Early work by Oe et al.13 and Noreika et al.14 showed a very narrow growth window for InSbBi due to the weak bonding of In-Bi and formation of metallic In or Bi droplets. InAsBi samples grown at low temperatures show binary InBi inclusions with a tetragonal PbO phase12. To our knowledge, InPBi single crystal has never been reported since the first theoretical prediction made 25 years ago11.

Here we report the first successful growth of InP1 − xn class="Chemical">Bix thin films with n class="Chemical">Bi concentration of x = 2.4 ± 0.4% by gas source molecular beam epitaxy (MBE). The InPBi crystal not only reveals excellent material quality, but also exhibits strong and broad room temperature photoluminescence (PL) from deep levels at 1.4–2.7 μm which provokes new theory to account for this phenomenon.

Results

Surface and structural properties

A set of samn class="Chemical">ples with thickness of 420 nm were grown by gn class="Chemical">as source MBE and their growth parameters are listed in Table 1. Details of the growth procedure can be found in Methods. Sample a, b and c have the same flux of all three elements but different growth temperatures. Atomic force microscopy (AFM) results show that droplets appear with a density in the order of 107 cm−2 as shown in Fig. 1(a) when the growth temperature is 364°C. The big ones are about 1200 nm in diameter and 300 nm in height, while the small ones are about 700 nm in diameter and 180 nm in height. Energy dispersive X-ray (EDX) measurements using a scanning electron microscope (SEM) shows the relative percentages of In and Bi atoms to be 89.4% and 7.3%, respectively, in the droplets (Fig. 1(d)), while the Bi concentration is very low between the droplets (Fig. 1(e)). A smooth surface could be obtained at 324°C with a low surface roughness value of only 0.93 nm from a 5 × 5 μm2 scan, see Fig. 1(b). If the growth temperature is further lowered to 275°C, a whisker-like surface feature with a high roughness value of 17 nm is observed as shown in Fig. 1(c). The dependence of the morphology on the growth temperature will be discussed later.

Table 1

Growth parameters and RMS roughness values measured from AFM of the three InPBi samples grown at different growth temperatures

Sample No.	Tg (°C)	TIn (°C)	TBi (°C)	PH3 (Torr)	RMS (nm)
a	364	940	465	350	84.6
b	324	940	465	350	0.93
c	275	940	465	350	17.1

Figure 1

(a), (b) and (c) show surface morphology of 5 × 5 μm2 InPBi measured by AFM using the tapping mode when grown at 364, 324 and 275°C, respectively. (d) and (e) show the EDX measurement results in an SEM from the InPBi sample grown at 364°C.

Figure 2(a) shows high-resolution X-ray diffraction (XRD) (004) ω-2θ rocking curves for the three samn class="Chemical">ples. The right narrow peak with a full width at half maximum (FWHM) of around 20 arcsec corresponds to the InP substrate and buffer, while the left peak corresponds to the n class="Chemical">InPBi layer. A weak and wide peak can be seen in InPBi grown at 364°C. A narrow FWHM of 46 arcsec is obtained from the sample grown at 324°C. Interference fringes can be clearly seen indicating a sharp InPBi/InP interface. The InPBi thickness is estimated to be about 420 nm from the simulations on the fringe separation. The strain relaxation is determined to be less than 5.0% from the two symmetric (004) scans and two asymmetric (115) scans and the lattice mismatch with respect to the InP substrate is 0.24% assuming the same Poisson ratio as InP. In order to confirm the single crystallization of the InPBi layer, the reciprocal space maps (RSM) of the 430 nm thick InPBi on both (004) and (115) planes were performed and plotted in logarithmic scale in Fig. 2(b) and 2(c). Two main peaks along with several satellite peaks can be clearly seen, in which the one with strong intensity belongs to the substrate, and the one with weak intensity belongs to the InPBi layer. The FWHM of the InP substrate peak and the InPBi is only 12.5 and 42.0 arcsec in (004) map, respectively, indicating excellent single crystal quality for a relatively thin epi-layer. The mismatch with respect to InP is 0.24% with strain relaxation of 5.4% which is the same as the value obtained using four ω-2θ scans. As the growth temperature is further lowered to 275°C, the InPBi peak becomes broadened. This layer has a mismatch of 0.25% with respect to InP and the strain relaxation increases to 21%. These results reveal that growth temperature is a critical parameter to obtain smooth InPBi single crystal.

Figure 2

(a) High resolution XRD (004) ω-2θ rocking curves from the InPBi samples. (b) and (c) are the reciprocal space maps of the 430 nm thick InPBi on both (004) and (115) planes, respectively.

After the growth temperature optimized, several samn class="Chemical">ples with similar growth parameters of samn class="Chemical">ple b were grown. The Bi concentration was determined by Rutherford Backscattering Spectroscopy (RBS), which is a powerful method to determine the atom concentration especially for heavy atoms since the He2+ beam will be easily scattered after impinging on them. A typical result is shown in Fig. 3. The Bi concentration is obtained by simulating the random spectrum. The black dotted line, the red dashed line and the blue line are the randomly oriented, the simulated and the aligned spectra, respectively. The surface spikes of the In signal at the channel number 455 of the [100] aligned curve is an artifact of the channeling analysis. Since the surface atoms are more exposed to the He2+ beam, they don't channel. From the RBS results, the Bi concentration is determined to be 2.4 ± 0.4% with 94 ± 5% Bi atoms at substitutional sites.

Figure 3

Randomly oriented (black dotted line), simulated (red dashed line) and aligned (blue line) RBS spectra for the InPBi layer grown at 324°C.

Figure 4(a) shows a typical cross-section Z-contrast scanning transmission electron microscope (STEM) image. A clear n class="Chemical">InPBi/InP interface and a lightly wavy surface are observed. The InPBi layer thickness is measured to be 423 ± 5 nm, close to the thickness determined by the XRD interference fringes. From the high resolution TEM results (not shown here), the lattice arrangement in both InP buffer side and InPBi side is very regular. The result indicates the epitaxial relation and very high interface quality. Figures 4(b)–(d) show the EDX-mapping results in STEM for Bi, In and P, respectively. Bi incorporation in the InPBi layer is clearly revealed in quite uniform distribution within the probed resolution. We note, however, information on Bi clusters or non-uniformity within the 10 nm regime can't be resolved from such a mapping, as the measured data represent the average value for each element within the e-beam probe which is typically of 10 nm. Nevertheless, we could not identify any detectable hint of Bi-related clusters from high-resolution TEM lattice imaging and high-resolution Z-contrast imaging, in contrast to our previous findings in GaAsBi samples15. For In and P mapping, the concentration of both elements is found to be reduced in InPBi compared with in InP. This is more clearly visible by the steps in Fig. 3(e) where depth profiles of the average Bi, In and P counts from EDX-mapping are shown. The detected drop of In concentration in InPBi is somewhat unexpected and can have important consequences for defect formations as will be discussed below.

Figure 4

(a) A typical STEM image and EDX-mapping results for (b) Bi, (c) In and (d) P. Figure (e) shows the depth profiles of the average Bi, In and P counts from the EDX-mapping (b)–(d).

Optical properties

Figure 5(a) shows n class="Chemical">PL spectra of sample b at 8–300 K and the InP reference samples at 300 K, respectively. The InPBi sample at 300 K reveals strong and broad PL signals at 1.4–2.7 μm, centered at 1.9 μm and its FWHM is as wide as around 660 nm (200 meV). The peak intensity is about 8 times stronger while the integrated PL intensity is almost 80 times higher compared with the InP reference sample. Careful analysis of the PL spectra reveals two peaks centered at 1.78 μm (0.70 eV) and 2.04 μm (0.61 eV) at 300 K, respectively, using the Gaussian fitting, and three peaks at T < 200 K. This strong PL signal is not observed from the InP reference sample, where only an interband transition peak at 0.92 μm (1.34 eV) can be seen. Figure 5(b) shows the normalized transmission spectra of the InPBi sample and InP reference sample of similar thickness at 300 K. The small bump observed at 0.92 eV is due to water absorption. The transmittance drops rapidly at about 1.3 eV mainly due to the large absorption coefficient near the band gap of InP16. Compared with the InP reference sample, the transmittance of InPBi sample remains essentially the same for the photon energy below 1.0 eV, and decreases faster for the photon energy above 1.0 eV, indicating stronger optical absorption in this energy regime and a smaller band gap of InPBi. Using the plot of square of absorption coefficient vs energy (not shown here), the band gap of InPBi is deduced to be 1.23 eV, giving an estimate of band gap reduction of 46 meV/%Bi. Comparing absorption with PL spectra, we notice that there is no energy overlaps of the interested parts of the signals. The absorption spectra show near band gap features at above 1.0 eV while the broad PL energy falls within 0.4–1.0 eV.

Figure 5

(a) PL spectra of the InPBi sample measured at 8-300 K and the InP reference sample at 300 K. The dashed blue lines are the Gaussian peak fitting. (b) The normalized transmission spectra of the InPBi sample and the InP reference sample at 300 K.

Discussion

To overcome the large miscin class="Chemical">bility gap and the weak bonding of III-n class="Chemical">Bi compared with those of conventional bonding of III–Vs, both the growth temperature and the V/Bi flux ratio should be kept low. This is the general strategy for almost all epitaxially grown dilute bismides. In the case of InPBi, use of a low P flux has a risk of forming In droplets. For the InPBi grown at 364°C, large metallic droplets are observed by AFM. However, the InP thin film grown under the same conditions reveals a smooth surface. This indicates that the P flux is high enough to ensure a P-rich surface reconstruction and thus a smooth surface. When Bi is added, large droplets are formed and the EDX measurement in SEM shows considerable amount of In in addition to Bi found in the droplets. This suggests that In is subject to alloying with Bi. It should be noted that both In and Bi will be in liquid form when the growth temperature is higher than their melting points (156.6°C and 271°C for In and Bi, respectively) and there indeed exist several In-Bi phases in the phase diagram17. The weak bonding nature of Bi atoms renders themselves to accumulate on the InP surface. Formation of In-Bi alloys consumes both In and Bi atoms and hinders the InP growth.

When the growth temperature is reduced to 324°C men class="Chemical">asured by the thermocoun class="Chemical">ple, high quality InPBi thin film is demonstrated. Although the real growth temperature is unknown, we believe it will be quite close to the Bi melting point. The probability of Bi atoms sticking on surface and incorporation into lattice are enhanced when decreasing the growth temperature. When the impingent Bi atoms are balanced by their incorporation to lattice sites, single crystal InPBi film is achieved and no excess Bi atoms will form Bi droplets on surface. Bismuth is known to be as an excellent surfactant, which ensures good surface and interface quality.

When the growth temperature is further lowered, excess n class="Chemical">Bi and n class="Disease">P atoms will accumulate on InP. Meanwhile, PH3 decomposition rate decreases with reducing growth temperature in gas source MBE and the resulting InPBi growth dynamics becomes rather complicated. The reference InP sample grown at 275°C with a slightly lower PH3 pressure of 320 Torr shows a similar surface feature as that of InPBi shown in Fig. 1(c). Further study is necessary to account for the observed surface morphology.

From the (004) XRD results shown in Fig. 2(a), there is some n class="Chemical">Bi incorporated in the sample grown at 364°C. As the rocking-curve signal is broad and weak, it is difficult to estimate the Bi content. The sample grown at 324°C reveals nice pendellösung fringes in the (004) scan indicating good interface. If we assume the Vegard's law is valid, the extrapolated lattice constant of zinc-blende InBi will be 6.52 ± 0.11 Å which is very close to the previously reported value of 6.5 Å18. The measured lattice mismatch is 0.24% for InPBi grown at 324°C.

The n class="Chemical">PL measurement reveals unexpectedly strong and broad signals at 0.4–1.0 eV at 300 K which are absent in the reference InP thin film grown at the same growth conditions. The absorption spectra show no difference between InPBi and InP for energy below 1.0 eV but enhanced absorption in InPBi for energy above 1.0 eV. This implies that physical origins can be different in the two energy regimes. Previous spectroscopic studies of InP:Bi grown by liquid phase epitaxy with the Bi concentration in the doping level of 10–100 ppm show many spectral features within 100 meV below the InP band gap at low temperatures, but no signals in the mid-gap at room temperature34. In these materials, Bi can be regarded as dopant and modelled by a resonant defect level within the energy band structure of the host matrix. In our case, the Bi concentration is much higher, and Bi interaction as well as cluster formation must be considered. There have been several theoretical investigations on the band structures of GaAsBi1920212223 and GaPBi23 to account for the band gap reduction of dilute bismides, but the physical explanation of whether it is due to the band anti-crossing or just the alloying and disorder is still under debate23. It is theoretically predicted that the Bi-related defect resonant level lies 80 meV below the valence band of GaAs22 and 122 meV above the valence band of GaP23. For disordered GaP1 − xBix including clustering effect with x up to 4.88%, broad band tail states exist at about 0.15–0.5 eV above the valence band23. If we assume similar band tail states also exist in InPBi, this may explain the near band gap features observed in the absorption measurement. However, it gives low energy bound of 0.84 eV (the band gap of InP is 1.34 eV at 300 K), still larger than the observed 0.70 eV and 0.61 eV in PL at 300 K. On the other hand, if the PL signals were due to the band anti-crossing, this would give exceptionally large band gap bowing of 267 meV/%Bi. Therefore, current existing theory can't explain the PL signals at the mid band gap range. Reihlen et al. reported large energy differences of 70–200 meV between PL and absorption at 300 K in InPxSb1 − x for 0.1 < x < 0.424. They attributed this difference to the combination of deep centers and intrinsic band tail recombination which were indistinguishable from the spectroscopy. The large band tail was believed to result from the compositional clustering. The compositional inhomogeneity usually induces localized states and causes S-shape in temperature dependent PL peak energy, often found in dilute nitrides25. However, in the temperature dependent PL spectra shown in Fig. 4(a), we don't observe any noticeable S-shape curves.

We tentatively attribute the strong and broad n class="Chemical">PL signals at 300 K to deep level luminescence. The deep levels can be n class="Chemical">BiIn anti-sites, P and/or In vacancies or the above combinations. The EDX-mapping results shown in Fig. 3 unambiguously confirm the decreases of both In and P atoms in InPBi compared with in InP. This gives some evidence that In and P vacancies may exist. Indeed, the scanning tunneling microscope (STM) on cross-sectional InPBi confirms the existence of P vacancies (not shown here). Bi atoms are generally believed to occupy P sites in dilute bismide. However, this will introduce very high local strain in InPBi. Our first principle calculations show that the overall formation energy is smaller for Bi taking In lattice positions instead of P positions. Also Bi can be of 3+ valence due to the relativistic effect of the two 6 s electrons. These suggest that BiIn anti-site is physically favorable and its concentration can be quite high. More experimental works, such as deep level transient spectroscopy (DLTS) and time resolved PL, are necessary to clarify the nature of those deep levels and carrier dynamics. Likewise, new theoretical model or simulations are stimulated to give physical explanations of the observed PL spectra.

In summary, new n class="Chemical">InPBi single crystal hn class="Chemical">as been successfully grown by gas source MBE for the first time. High surface and structural quality has been demonstrated and confirmed by AFM, XRD and TEM measurement. RBS measurement shows that the Bi concentration is 2.4 ± 0.4%, with 94 ± 5% Bi atoms at substitutional sites. Room temperature PL measurement reveals a very broad and strong spectrum from 1.4 to 2.7 μm while the optical absorption shows a band gap reduction of 110 meV with respect to InP. The observed PL spectrum is likely from deep levels and can't be explained by the existing theory.

Methods

Molecular Beam Epitaxy

Samn class="Chemical">ples were grown on (100) semi-insulating InP substrates by V90 gas source MBE. P2 was cracked from phosphine at 1000°C and its flux was controlled by regulating the gas pressure in the gas line. Elementary In and Bi sources were used and their fluxes were controlled by adjusting the effusion cell temperatures. Both substrate and cell temperature were measured by thermocouples.

Prior to growth, n class="Chemical">oxide desorption from the substrate surface wn class="Chemical">as carried out by heating the substrate with a P2 flux impingement on the surface. The complete desorption of the oxide layer was determined using the reflection high energy electron diffraction (RHEED), by a sharp transition to (2 × 4) pattern at about 524°C. A 75 nm un-doped InP buffer layer was grown at 474°C with the PH3 pressure of 630 Torr to improve the surface quality. 420 nm thick InP1 − xBix thin films were then grown at 275–364°C at 0.86 μm/h with a low P flux to facilitate the Bi incorporation into the lattice. InP reference samples were grown with otherwise the same growth conditions as those for InP1 − xBix for comparison.

Measurements

The surface morphology wn class="Chemical">as measured by Bruker Icon atomic force microscopy (AFM) using the tapping mode. The scan area was 5 × 5 μm2. The electron acceleration voltage for the EDX measurement in SEM was set at 10 kV to obtain good spatial resolution. Conventional X-ray diffraction ω-2θ scan was performed using a Philips X-pert high-resolution diffractometer to measure structural properties. The Cu Kα1 wavelength was selected by a Ge (220) fourfold monochrometer. The STEM and the EDX-mapping measurements were performed on a FEI Tecnai G2 F20 microscope. Cross-section STEM specimens were prepared by focus ions beam. In the RBS measurement, a 2.275 MeV4He2+ source was used to impinge on the sample surface. The PL and transmission spectra were measured using a Nicolet Megna 860 Fourier transform infrared (FTIR) spectrometer, in which a liquid-nitrogen cooled InSb detector and a CaF2 beam splitter were used. A diode-pumped solid state laser (λ = 514.5 nm) was used as an excitation source for PL measurement.

Author Contributions

S.M.W. conceived the idea behind this work, proposed the preliminary exn class="Chemical">planation for n class="Chemical">PL results and revised the manuscript. K.W. developed and grew all the samples, wrote the first version, performed the AFM and Raman measurement. Y.G. performed the PL and transmission measurement, revised the manuscript. C.Z.K. assisted in the sample growth. H.F.Z., L.Y.Z. and W.W.P. assisted in various material characterization measurements. M.J.W. performed the TEM measurement. P.F.L. performed the first principle calculations. Q.G. and all the other authors were involved in the discussion.