The investigation was divided into four parts. In the first part, we investigated the effect of different misorientations; in the second part, we focussed on the effect of the growth conditions; in the third part, we looked at the optical properties; and in the final part, we combined this to obtain smooth N-polar GaN on low misorientation sapphire while trying to maintain good layer quality.
3.1. Impact of Sapphire Misorientation
Figure 1 shows the surface of N-polar GaN layers after the simultaneous growth of ≈1.4
m in 45 min growth at 1200 °C with a V/III ratio of 1650 on sapphire wafers with various misorientations in the A or M direction. Under these conditions, flat surfaces have been obtained already for misorientation angles of 0.8° and higher.
When looking at the surfaces using AFM in
Figure 2, one can clearly see that under our conditions, the surfaces grown on sapphire with M misorientation had a generally lower root mean square roughnesses over
m
than the ones along A. There is, however, a large-scale waviness on the M misorientation samples that leads to higher roughnesses when measured over larger areas (e.g., 8–9 nm for 1°M over
m
m compared to 2–3 nm for 2°A) as reported before [
4,
5].
In addition to the different roughness, there is also a different morphology. The wafer flats (sapphire a-plane) were at 45°. Thus, the step edges should run at 45° on M misorientations (top row in
Figure 2) and −45° on A (middle row in
Figure 2). However, the steps are not very visible for M misorientation in the top row. Our explanation is that when the step-edges are parallel to the GaN m-plane (as on A misorientation), they are straight and visible. However, when the step-edges run along the a-plane, they become unstable and instead form a zig-zag of different equivalent m-planes, which causes the larger waviness. This is more clearly seen in the large-area surface maps in
Figure 2 bottom, again very straight steps from 1.5°A and a strong meandering of steps and higher RMS roughness for 1°M.
The morphology is also reflected in the XRD rocking curves FWHM in
Figure 3; the smoother a sample, the narrower the 0002 FWHM becomes. When the HHLs vanish, the FWHM reduces from 700–900″ to ≈400–550″. The reason is that the HHLs all have a slightly different twist and tilt. Hence, a coalesced sample without HHLs has larger-sized domains with similar twist and tilt and thus less broadening. However, the XRD FWHMs increase again together with the RMS roughness at 4°A misorientation, which is probably due to the very pronounced step-bunching reducing again the average domain size.
Thus, for our growth condition, we found that sapphire misorientations along the m-direction result in smoother surfaces in small scales with narrower XRD FWHM. Smooth surfaces were obtained for misorientations as low as 0.8°, which is less than what had been previously reported. Apparently, the growth conditions are most important in order to obtain smooth N-polar surfaces and with them a good crystal quality.
3.2. Impact of Growth Conditions
Therefore, we investigated the impact of growth conditions, essentially the variation of the NH flow. For this, we used 2°A sapphire together with 0.3°M or 0.2°M. There were two series with a nominal growth rate of 0.45 nm/s at 1150 °C and a second using a 1.5× higher TMGa flow at 1180 °C to obtain smoother surfaces at low misorientations. For each of the series, the ammonia flow was changed and all other parameters were kept constant.
Even to the naked eye, the surfaces of the N-polar GaN on the 0.3°M change strongly in
Figure 4. There are two trends visible. At high V/III ratios (above 500), the HHLs became smaller, and additional smaller HHLs appear. Thus, the height differences and with it the RMS roughness decrease. Below 500, the tops of the HHLs become first flatter and then rounder until at a V/III ratio of 33 just some <50 nm high pancake-like structures remain on 0.3°M at 1180 °C. The same but a little weaker is observed on 0.2°M, too.
In comparison, the GaN surfaces grown on 2°A sapphire remained flat for almost all conditions with an average RMS roughness of nm. The waviness only increased at the lowest and highest V/III ratio, which is in reverse to the trend for the 0.2°M and 0.3°M surfaces.
Figure 5 shows the same trend in numbers for the growth at 1150 °C. The roughness on 0.2°M sapphire and the growth rate on 2°A correlate nicely; both increase with the increasing V/III ratio until 500. At higher V/III ratios, the overall roughness and the growth again rate decrease. The XRD FWHM instead decreases until a V/III ratio of 500 and then increases again.
The same reversal of trends at a V/III ratio of about 500 is seen in
Figure 6 for the growth with a higher TMGa flow at 1180 °C on 0.3°M and 2°A sapphire. At this high temperature, the GaN layer grown on the 2°A sapphire substrates developed larger holes at V/III ratios below 100 and thus, XRD and electrical measurements were not reliable, and only partial data are available for 2°A. Nevertheless, the trends reverse for nearly all parameters at a V/III ratio of around 500, which is similar to the 1150 °C growth series. The growth rate first increases and then decreases with the increasing V/III ratio, the XRD rocking curve FWHM first decreases and then increases, and there is an opposite trend in RMS roughness between low and 2°A misorientation with either maximum respective minimum at V/III ratios of 500. A reversal of trends was also found in the PL data, which is discussed in the next
Section 3.3.
At first glance, it seems surprising that the XRD FWHM for samples with hillocks (
Figure 6c) and with a smooth surface on 2°A (
Figure 5c) show exactly the same trends despite a reverse trend in RMS roughness (
Figure 5a and
Figure 6a). However, for the GaN hetero-epitaxy, the initially high dislocation densities are reduced by a 3D growth phase at the onset of buffer growth, where dislocations are bent towards each other and can annihilate. Too long diffusion lengths (at very low V/III ratios) lead to a very quick coalescence and thus to a very short 3D growth phase with many remaining dislocations. Hence, a surface with a long and strong 3D growth phase is the one with the lowest number of dislocations and thus the one with the smallest XRD FWHM. For low misorientations, this results in the largest HHLs.
The initial increase of the growth rate with increasing V/III ratio can be explained by the effect of a reduced Ga desorption with increased NH flow. Since the Ga species stay a shorter time on the surface until they meet with NH, the probability of desorption decreases compared to incorporation. Shorter times on the surface with increasing NH flow also translates into shorter diffusion lengths, which increases the 3D growth that is essential for dislocation annihilation, as discussed above. Thus, the XRD FWHM reduces too with the increasing V/III ratio. However, the reversal of trends starting at V/III ratio of 300 shows that there is more at work than supersaturation alone.
One way to understand this would be a change of the surface termination. On an ideal (000
) surface, nitrogen is bound with three bonds downwards, while Ga has three upwards dangling bonds and is bound with just a single bond. However, this surface is not the most stable structure. According to DFT calculations, there are two more stable surfaces on N-polar (000
) GaN: a Ga-adlayer at Ga-rich conditions and a N-terminated 3N-H
surface at N-rich conditions [
15]. The N-polar Ga-adlayer surface has only half the energy barrier for oxygen adsorption compared to the Ga-polar Ga-adlayer surface [
16] and thus will incorporate much more oxygen than the Ga-polar (0001) Ga-adlayer surface, leading to a very high oxygen background at lower V/III ratios. On the other hand, all adsorption sites on the 3N-H
surface are Ga sites, which makes oxygen incorporation more difficult since oxygen would prefer a nitrogen site. Furthermore, it seems likely that oxygen reaching the 3N-H
surface rather reacts with the adsorbed hydrogen and should further reduce oxygen incorporation.
Therefore, we assume that at low V/III ratios, a Ga-adlayer covers the entire N-polar surface. However, when increasing the V/III ratio above 300, domains of the N-rich 3N-H
surface would appear as well and grow in size, until at V/III ratios above 2000, at which point the N-rich 3N-H
would cover the entire surface. The change of the dominating surface would affect the adsorption and desorption probabilities of oxygen and Ga. The oxygen incorporation would reduce from the high value of the Ga-adlayer at V/III ratios below 300 to the low value of the 3N-H
for V/III ratios above 2000. Indeed, the data in
Figure 5d) show a high value for low V/III ratios and then a strong reduction for the V/III ratio between 300 and 2000 to a constant low value.
The growth rate is more complex, since it overlaps with the thermodynamic effect of NH
on Ga surface lifetimes and Ga desorption as discussed above. Furthermore, no calculated data on Ga adsorption exists. However, a N-terminated surface such as the 3N-H
means that Ga can adsorb only with a single bond compared to three bonds on a Ga-adlayer surface. Moreover, the mixed N/hydrogen termination of the 3N-H
surface means that even the dangling bonds are most often saturated by hydrogen. Together, it is reasonable to assume a lower Ga adsorption probability on the 3N-H
surface. Then, the adsorption probability of Ga would start to reduce above a V/III ratio of 300 with the increasing 3N-H
coverage. This reduction compensates for the increase from shorter surface lifetimes of Ga with increased NH
flow. With the increasing V/III ratio, the growth rate reduces with increasing 3N-H
area until at V/III ratios above 2000, only the N-rich 3N-H
surface covers everything. For higher V/III ratios, the growth rate should stay constant, and indeed, the last point of
Figure 5 at the maximum possible V/III ratio of 2800 is changing less than the trend line.
The simultaneous presence of domains with two different surfaces at the same time should also shorten the diffusion length because of the high number of dangling bonds at the border of a domain. Thus, the diffusion length is at a minimum between 500 and 1000, where the coverage with either surface would be equal. The diffusion length would increase again with higher V/III ratios until the whole surface is only covered by the 3N-H
surface. Further increasing the V/III ratio should again decrease the diffusion length due to a reduced surface lifetime with more incoming N-species, but this was not possible with the given TMGa flow rates. The increasing diffusion length between the V/III ratios of 500 and 2000 due to the larger single surface domains could explain the smaller and denser HHLs on low misorientations (
Figure 4 top and middle) and the increasing XRD FWHM at very high V/III ratios (
Figure 5c) and
Figure 6c due to weaker 3D growth as the 3N-H
domain size increases.
To summarize our model: at lower V/III ratios, the surface is entirely Ga-adlayer terminated, and the oxygen incorporation is at maximum. Since the Ga diffusion length equals the time to meet a N species to incorporate, the Ga diffusion length increases with the decreasing V/III ratio. However, a longer Ga diffusion length also means that the mobile Ga rests on the surface for longer, thus increasing the probability of Ga desorption. Hence, the growth rate decreases with decreasing V/III ratios until for very low V/III ratios below 30, the growth becomes N-limited and the growth rate drops sharply (below V/III = 40 in
Figure 5b and
Figure 6b). When increasing the V/III ratio above 300, domains of the 3N-H
surface start to appear. This 3N-H
surface has a much lower probability of oxygen incorporation, and thus, the residual oxygen starts to decrease with increasing 3N-H
coverage (
Figure 5d). In addition, Ga adsorption is reduced on the 3N-H
surface, reducing the growth rate with increasing 3N-H
coverage with increasing V/III ratio. At the same time, the presence of two different surfaces strongly reduces the Ga diffusion length due to a high probability of nucleation at domain boundaries and enhances the initial 3D growth phase, leading to a minimum in XRD FWHM and a maximum in HHL sizes on low misorientations at V/III ratios around 500. Increasing the V/III ratio further reduces the Ga-adlayer surface area until only the 3N-H
surface covers everything for V/III ratios above 2000. Until this point, the oxygen incorporation and growth rate decrease, the Ga diffusion length increases again, the FWHM increases (due to less 3D growth), and the HHLs sizes decrease again.
It must be emphasized that the conditions for reaching a nitrogen-terminated surface depend not only on the TMGa and NH
flow rates but also on the reactor design. Consequently, the V/III ratio where the 3N-H
surface dominates (and thus, the oxygen incorporation has reached its minimum) differs a lot among the literature judging from the reported oxygen incorporation [
2,
8,
9,
10]. Even for our two series at 1150 °C and 1180 °C, the value of the turning point at 500 was rather a coincidence: the higher temperatures and thus higher NH
decomposition was offset by the higher growth rate for the second series.
3.3. Optical Properties of N-Polar GaN
The room-temperature photoluminescence data are complex.
Figure 7 shows three typical spectra at the low (50) (a), trend reversal (500) (b) and highest (1750) (c) V/III ratio for the 1180 °C series. Most striking is a strong shift of the band edge for the low miscut samples; see
Figure 7d. The shift of the band edge is mostly related to strain, since the different surface morphologies (smooth versus HHLs) resulted in a different strain, as confirmed in
Figure 7e: the shift of the band edge follows nicely the strain along [0001]. The slope of the dotted line is 14 eV which is between the literature values for the Ga polar GaN of 12 eV [
17] and 15.4 eV [
18]. However, at low V/III ratios (and high carrier concentrations), other effects such as band filling may additionally shift the bands to even lower values and cause a deviation from the strain-related shift.
The yellow luminescence band (YL) in
Figure 7f is the highest in relation to the band edge at the V/III ratio of 500. This was independent of the surface roughness; similar trends were observed for 2°A and 0.2°M/0.3°M despite their opposite trends in roughness. This would indicate that the YL is unrelated to HHLs. Following our assumption of a surface transition around this V/III ratio, then the YL could be related to a defect from domain boundaries at each surface. Carbon alone is unlikely to cause the YL, since the carbon background is low (<10
cm
) and is almost not changing with the V/III ratio (
Figure 5d). Furthermore, the XRD FWHM is at minimum at the highest YL (
Figure 5c and
Figure 6c), ruling out a connection with extended structural defects or dislocations. This is quite different from the behaviour of Ga-polar GaN where a high YL is usually associated with carbon-related complexes and poor crystal quality.
The YL of our samples is centred around 2.2 eV (570 nm) which according to a recent report for the molecular beam epitaxy of N-polar GaN might be related to Gallium vacancies [
19]. However, there are reports of YL of N-polar GaN related with carbon. After very high C- and O-implantation into N-polar GaN, the YL only increased with C-implantation [
20]. Moreover, a strong YL for N-polar GaN growth using Na-flux was correlated with high carbon incorporation [
21]. The same report noted that the YL was shifted to longer wavelengths for N-polar GaN [
21], while the YL shifted with heavy carbon implantation to slightly shorter values around 550 nm [
20]. Given the very low carbon incorporation for V/III ratios as low as 100 (
Figure 5d), which is typical for N-polar GaN MOVPE growth [
2,
9], for our series, a direct correlation of carbon and YL seems unlikely. The best explanation seems that there is more than one complex active in the same energy range, and Ga vacancies would agree with an increasing YL due to growth on a disordered surface with two simultaneous terminations.
There is further evidence for a different active complex.
Figure 7a–c also show a second band centred at 430 nm (hence, blue luminescence band, BL) for the 0.3°M (and 0.2°M) samples, which is most prominent at the lowest and highest V/III ratios. This BL is barely visible on the higher misorientations (and thus smoother surfaces) and thus was often hidden in the donor–acceptor-related transitions with phonon replicas close to the band edge [
2,
22]. However, a strong BL was reported on N-polar GaN with many HHLs [
23]. In another series of samples at lower growth temperatures and even lower V/III ratios, the BL increased with increasing growth temperatures (
Figure 7i).
Moreover, the BL anticorrelates with the YL and the roughness but correlates with the net donor concentration (
Figure 7h). An early publication suggested indeed a higher carbon incorporation with HHLs on N-polar GaN [
8] and speculated that the apex of the HHLs could favour carbon incorporation. In a recent report, a BL at 430 nm can be observed to increase after heavy carbon implantation and annealing of N-polar GaN (
Figure 5 in [
20]). Finally, on Ga-polar GaN, the blue luminescence is related to carbon [
24].
Since the BL increases in intensity with growth temperature and carrier concentration, we speculate that the BL might probably come from electrically active carbon incorporated on a Ga site or maybe an interstitial. However, the low-carbon background makes this questionable, and more data are needed.
3.4. Smooth N-Polar GaN with Low Misorientation
As shown in
Section 3.1, misorientations in m-directions are favourable to reduce HHLs at low misorientations. Since typical (0001) sapphire substrates are usually misoriented along this direction, we used 0.2°M and 0.3°M in our study as well.
The best sample from these series is shown in
Figure 8, which was obtained deep in the desorption-limited growth regime at the stability edge of even N-polar GaN. (Ga-polar GaN layers were completely desorbed at the conditions within a few minutes.) The large area RMS roughness is 13 nm, which is three to four times higher than for typical GaN surfaces (4–6 nm) but could be still accepted for devices. Moreover, the surface of this sample shows even some traces of atomic steps. This was the only, low misorientation N-polar GaN surface where we have ever seen step edges. Typically, the surfaces are featureless with 1 nm deep holes such as the surface shown in
Figure 9.
However, the smooth surface comes with a price. As shown, in the bottom right of
Figure 8, the XRD 0002
rocking curve FWHM strongly increases when the RMS roughness decreases. This decrease of crystal quality is due to a lack of 3D growth, i.e., the layer coalesces too fast, and the dislocations do not have enough time to bend sidewards and cannot annihilate.
Hence, we have tried a multi-step approach for the high-temperature growth of N-polar GaN. The first step was with a high V/III ratio to start with many nuclei. In the next step, the layer is smoothed by overgrowth with a very low V/III ratio. By this, reasonable flat surfaces with acceptable XRD FWHM around 700–900 arcsec have been achieved.
Figure 9 shows a typical resulting surface. The surface does not show atomic steps. Instead, there are many 1–2 nm deep holes. Such a morphology has been reported before for N-polar GaN [
25]. The density of the holes is in the
–
cm
range and thus might be related to dislocations.
There remains another challenge: While somewhat flat and with not too bad crystal quality, these layers had high background carrier concentrations. Thus, we tried regrowth on these surfaces with the optimised V/III ratio of 1000. This must be performed in a third step because half of the sample thickness was desorbed when heating to 1180 °C for the regrowth.
The third layer has lower carrier densities close to
cm
. However, small HHLs reappear on its surface, as in
Figure 9 bottom right. The density of these small HHLs is 1–2
cm
with a height of 20–40 nm after 30 min of regrowth. Due to the above-discussed strong GaN desorption before the start of the regrowth, we have not systematically studied their shapes as a function of growth time.
Therefore, it is possible to obtain a smooth N-polar GaN with decent structural quality on low misorientation. However, it seems very challenging to also obtain a low background carrier concentration.