# Large-Scale Defect Clusters with Hexagonal Honeycomb-like Arrangement in Ammonothermal GaN Crystals

Lutz Kirste, Thu Nhi Tran Thi Caliste, Jan L. Weyher, Julita

Smalc-Koziorowska, Magdalena A. Zajac, Robert Kucharski, Tomasz

Sochacki, Karolina Grabianska, Malgorzata Iwinska, Carsten Detlefs, et al.

► To cite this version:

Lutz Kirste, Thu Nhi Tran Thi Caliste, Jan L. Weyher, Julita Smalc-Koziorowska, Magdalena A. Zajac, et al. Large-Scale Defect Clusters with Hexagonal Honeycomb-like Arrangement in Ammonothermal GaN Crystals. Materials, 2022, 15 (19), pp.6996. 10.3390/ma15196996. [hal-05075692](https://hal-05075692)

HAL Id: hal-05075692

https://cnrs.hal.science/hal-05075692v1

Submitted on 20 May 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Large-Scale Defect Clusters with Hexagonal Honeycomb-like Arrangement in Ammonothermal GaN Crystals

Lutz Kirste <fnref n="1" />,* ¹D, Thu Nhi Tran Thi Caliste <fnref n="2" />, Jan L. Weyher <fnref n="3" /> ¹D, Julita Smalc-Koziorowska <fnref n="3" />, Magdalena A. Zajac <fnref n="3" />,<fnref n="4" />, Robert Kucharski <fnref n="3" />, Tomasz Sochacki <fnref n="3" /> ¹D, Karolina Grabianska <fnref n="3" /> ¹D, Malgorzata Iwinska <fnref n="3" /> ¹D, Carsten Detlefs <fnref n="2" />, Andreas N. Danilewsky <fnref n="5" /> ¹D, Michal Bockowski <fnref n="3" />,<fnref n="6" /> ¹D and José Baruchel <fnref n="2" />

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Kirste, L.; Tran Thi Caliste, T.N.; Weyher, J.L.; Smalc-Koziorowska, J.; Zajac, M.A.; Kucharski, R.; Sochacki, T.; Grabianska, K.; Iwinska, M.; Detlefs, C.; et al. Large-Scale Defect Clusters with Hexagonal Honeycomb-like Arrangement in Ammonothermal GaN Crystals.

*Materials* 2022, 15, 6996. https://doi.org/10.3390/ma15196996

Academic Editor: Alexander
A. Lebedev

Received: 31 August 2022

Accepted: 3 October 2022

Published: 9 October 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** In this paper, we investigate, using X-ray Bragg diffraction imaging and defect selective etching, a new type of extended defect that occurs in ammonothermally grown gallium nitride (GaN) single crystals. This hexagonal "honeycomb" shaped defect is composed of bundles of parallel threading edge dislocations located in the corners of the hexagon. The observed size of the honeycomb ranges from 0.05 mm to 2 mm and is clearly correlated with the number of dislocations located in each of the hexagon's corners: typically ~5 to 200, respectively. These dislocations are either grouped in areas that exhibit "diameters" of 100–250 µm, or they show up as straight long chain alignments of the same size that behave like limited subgrain boundaries. The lattice distortions associated with these hexagonally arranged dislocation bundles are extensively measured on one of these honeycombs using rocking curve imaging, and the ensemble of the results is discussed with the aim of providing clues about the origin of these "honeycombs".

**Keywords:** gallium nitride (GaN); X-ray Bragg diffraction imaging; rocking curve imaging; defect selective etching; honeycomb defect; limited size subgrain boundary; geometrically necessary dislocation; Borrmann effect

## 1. Introduction

Gallium nitride (GaN) wafers with low defect density are the substrate of choice for demanding GaN-based optoelectronic and electronic device structures. These include lateral high electron mobility transistors (HEMTs) and vertical field-effect transistors (FETs) as well as laser diodes (LDs) [1–3]. A suitable method to produce high-quality GaN crystals with low defect density is the ammonothermal approach [4–6], a solvothermal growth technique from a solution where a sufficient solubility of GaN is achieved by supercritical ammonia and mineralizers under high pressure. The ammonothermal process for growth of GaN (Am-GaN) is the ammonia counterpart of the hydrothermal crystallization of quartz from supercritical water. Such an approach is characterized by a high reproducibility. Several dozen crystals can be produced simultaneously in an autoclave in one growth run. So far, the ammonothermal growth method is the technique best suited to mass production that yields GaN crystals with the lowest threading dislocation density ($TDD < 1 \times 10^4 \text{ cm}^{-2}$) [7].

Recently, Kirste et al. demonstrated the high crystalline quality of entire 1.8-inch Am-GaN substrates by Lang-technique X-ray Bragg diffraction imaging (historically called

<footnote n="1">¹ Fraunhofer Institute for Applied Solid State Physics, Tullastraße 72, D-79108 Freiburg, Germany</footnote>

<footnote n="2">² European Synchrotron Radiation Facility, 71 Avenue des Martyrs, F-38043 Grenoble, France</footnote>

<footnote n="3">³ Institute of High-Pressure Physics, Polish Academy of Sciences, Sokołowska 29/37, 01-142 Warsaw, Poland</footnote>

<footnote n="4">⁴ Faculty of New Technologies and Chemistry, Military University of Technology, ul. Gen. Sylwestra Kaliskiego 2, 00-908 Warsaw, Poland</footnote>

<footnote n="5">5 Crystallography, Institute of Geo- and Environmental Sciences, Albert-Ludwigs-Universität Freiburg, Hermann-Herder-Straße 5, D-79104 Freiburg, Germany</footnote>

<footnote n="6">⁶ Center for Integrated Research of Future Electronics, Institute of Materials and Systems for Sustainability, Nagoya University, C3-1 Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan</footnote>

<footnote>* Correspondence: lutz.kirste@iaf.fraunhofer.de</footnote>
"X-ray topography") [8]. The main contrast mechanism for the observation of the defects was, in this case, the anomalous transmission ("Borrmann effect") predicted by the dynamical theory of diffraction (see, for instance, Authier [9]). The occurrence of the Borrmann effect in the studied substrates constitutes clear evidence of the high crystalline perfection of Am-GaN grown by a native seed approach. Using this laboratory X-ray Bragg diffraction imaging technique, not only microscopic defects such as threading dislocations (TDs), but also macroscopic defects, e.g., dislocation clusters due to preparation insufficiency, traces of facet formation, growth bands, dislocation walls, and dislocation bundles, were detected [8].

This previous study pointed out the presence, within the Am-GaN crystals, of a new type of defect cluster characterized by its hexagonal shape, macroscopic (up to 2 mm) size, and locally high defect density. We show in the present publication that this "honeycomb-like" defect cluster is mainly composed of six dislocation bundles (DBs) located in the corners of the hexagon. In order to characterise these defect clusters, several Bragg diffraction imaging techniques were combined with defect selective etching. The detailed examination of this newly discovered "honeycomb defect" (HCD) in Am-GaN crystals is the subject of the present work.

## 2. Samples and Experimental Techniques

### 2.1. Investigated Ammonothermal GaN Crystals

The investigated GaN crystals were grown by basic ammonothermal crystal growth technology in Ni-based alloy high-pressure autoclaves in a temperature range of 500–600 °C and ammonia pressure of 0.3–0.4 GPa. The growth principle is based on dissolving a polycrystalline GaN feedstock in one zone of the high-pressure autoclave filled with supercritical ammonia and transporting the dissolved material to the second zone, where the solution is supersaturated and the crystal growth of GaN takes place on native seeds [5,6,10]. Owing to a negative temperature coefficient of solubility, the so-called retrograde solubility, the chemical transport of GaN is directed from the low-temperature solubility zone to the high-temperature crystallization zone [11]. The growth process is complex: the first step aims to produce, from smaller seeds, GaN crystals with sufficient sizes. These are then used as new seeds for preparing substrates, which are obtained by applying proper wafering procedures. Further details of the basic ammonothermal growth and wafering process can be found elsewhere [5,6,10,12].

More than 40 Am-GaN substrates from different growth runs were investigated in this study. They included n-type material with a high oxygen donor concentration and semi-insulating material with manganese compensation. The diameter of the investigated Am-GaN (0001) c-plane substrates was 1 inch, 1.5 inches, 1.8 inches, or 2 inches with a thickness of 300 µm–450 µm. Both surfaces of the substrates, parallel to the (0001) and (0001̄) crystallographic planes, were prepared to exhibit a roughness in the 0.1 nm range ("epi-ready state").

### 2.2. X-Ray Bragg Diffraction Imaging

The defect analysis of these low-defect density Am-GaN crystals requires characterization techniques with a large field of view and high sensitivity to crystalline distortions, able to visualize the dislocations and grain boundaries as well as to provide a measurement of the lattice distortions. X-ray Bragg diffraction imaging (historically called X-ray topography), which allows imaging defects/dislocations in bulk crystals with low dislocation density (~10⁵ cm⁻² and below), fulfils these requirements. This X-ray imaging method is sensitive to weak misorientations (typically a few microradian) and strains (in the 10⁻⁶ range), over large sample areas (~cm²-range), with a spatial resolution in the µm range [9,13].

The classical X-ray Bragg diffraction imaging technique is based on the analysis of the X-ray diffracted beam intensity recorded by a two-dimensional detector. The "local" intensity recorded on a given area of the detector is a function of the gradient of distortion
SWB-XRT was performed using an extended parallel beam ("projection topographs") for both the transmission and reflection geometries. LauePT software was used to set the orientation of the Am-GaN substrates on the goniometer and for the indexation of the Laue pattern [21]. For the exposures in transmission geometry, the samples were tilted by either plus or minus 12 degrees about the horizontal axis; this allowed access to the Laue spots (topographs) of each of the six types of {1120} and {\bar{1}100} prismatic lattice planes. For the above-mentioned wafer thicknesses, the µt value amounts to 7.1 and 20.5 for the type 1120 reflection (λ<sub>1120</sub> ≅ 66.30 pm, 18.70 keV) and for the \bar{1}100 (λ<sub>\bar{1}100</sub> ≅ 99.45 pm, 12.47 keV) reflection, respectively. A sample-to-film distance of about 100 mm was used and the X-ray beam size was limited to 3 × 3 mm<sup>2</sup> to avoid overlapping of the various topographs on the film-detector. The 0006 reflections were measured by means of the large-area back-reflection geometry. Here, the wafers were tilted by 7 degrees about the horizontal axis, the sample-to-film distance was again 100 mm, and the X-ray beam size was limited to 5 × 5 mm<sup>2</sup>. The transmission topographs were recorded on SLAVICH VRP-M high-resolution films (Geola, Vilnius, Lithuania) and were digitized afterwards. A high-resolution Excelitas PCO 4000 CCD camera with a pixel size of 2.5 × 2.5 µm<sup>2</sup> (Excelitas PCO GmbH, Kelheim, Germany) was used for back-reflection topography.

*Synchrotron monochromatic rocking curve imaging (RCI)*: The crystalline sample is set to Bragg-diffract the large, nearly parallel, monochromatic X-ray beam and a series of X-ray diffraction topographs is recorded along the global rocking curve using a pixelated two-dimensional detector. Each pixel records an individual rocking curve originating from a defined sample volume. Real-space maps are generated from the ensemble of the recorded rocking curves: by fitting to a Gaussian profile, these maps allow for a quantitative investigation of local crystalline distortions.

The experiments were carried out in transmission geometry at the beamline BM05 at the European Synchrotron Radiation Facility (ESRF), Grenoble, France. Beams of either 18 keV (λ = 68.88 pm) or 30 keV (λ = 41.33 pm) were extracted from the synchrotron radiation using a vertically diffracting Si (111) double-crystal monochromator. We combined the measurements with these different energies to study the influence of the absorption conditions on the defect images (Appendix A). All samples were mounted in such a way that the GaN \bar{1}1\bar{2}0 or \bar{1}100 reflection was accessible in the vertical scattering plane. The beam size was 5 × 5 mm². Bragg diffraction images were recorded on a detector system comprising a scintillator screen equipped with microscope folded-relay optics to a 2048 × 2048 pixel low noise Excelitas PCO edge gold 4.2 sCMOS camera (Excelitas PCO GmbH, Kelheim, Germany). The use of this system resulted in an effective pixel size of 1.2 × 1.2 µm², so the recorded images had a corresponding field of view of ~2.4 × 2.4 mm². The distance between the sample and the detector was 50 mm. The sample was rotated along the rocking curve with an angular step of 10⁻⁴ degrees over a range of 0.01°. For every step, an image was recorded, i.e., a total of 100 X-ray diffraction images were captured for each scan. The generated stacks of images contain, for each of the pixels, a "local" rocking curve. These ~4 × 10⁶ rocking curves were processed by the Rocking Curve Imaging Analysis (RCIA) software tool [22]. Each of the rocking curves, corresponding to a single pixel, was fitted independently using a Gaussian function. Maps of the integrated intensity (INT), the full-width at half-maximum (FWHM), and the peak position angle (PPOS) of the diffracted intensity are then calculated from the ensemble of these rocking curves.

The INT map is comparable to a classical X-ray topograph and provides information on the defect structure of the crystals. The FWHM and the PPOS maps provide, in addition, quantitative values related to the local effective misorientation δΘ of the used lattice planes [9]. This effective misorientation is mainly constituted by the local variation in the rotation Δφ (or tilt) of these lattice planes and the local relative variation in the lattice parameter spacing Δd/d. These maps are, therefore, quantitative measurements of the local crystal lattice distortion and indicate the local lattice quality of the investigated crystal sample [23].
Under low-absorption conditions (μt < 1), both higher diffracted intensity and higher FWHM values correspond to an increased gradient of distortion within the diffracting volume. This is not what occurs for higher absorption conditions (typically μt > 6), when Borrmann conditions are present, as is the case for some of our experiments on Am-GaN, the INT map will behave reciprocally to the FWHM map, i.e., for areas of high defect density, the FWHM map will display high values, but the INT map will exhibit low intensity. This is because of the fact that defects "disturb" the anomalous transmission; when the strain gradient around the defect is high, new wavefields are created and propagate. Most of these new wavefields are strongly absorbed and this leads to a locally reduced intensity ("Borrmann shadow") [9,24]. The PPOS map provides information on the local departure from the nominal Bragg angle. As indicated above, this departure includes contributions from lattice rotations (tilts) and from lattice plane spacing variations (strain) [23,25]. By combining the INT, FWHM, and PPOS maps, we can obtain extra quantitative information with respect to that provided by classical topography techniques.

The principle of the RCI technique and evaluation procedure are presented schematically in Figure 1. Table 1 gives an overview of the measurement conditions for all three XRT methods used for this study. The values were calculated using LauePT software [21].

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Schematic drawing of Bragg diffraction imaging using (a) rocking curve technique with a wide monochromatic incident beam; (b) acquisition of an image series (each pixel records its own local rocking curve); and (c) reconstructions of maps of the integrated intensity (INT), the full-width at half-maximum (FWHM), and the diffraction peak position (PPOS), extracted on a pixel-by-pixel basis from the local rocking curves.</figcaption></figure>

<table><caption>**Table 1.** Overview of the used measurement conditions for L-XRT, SWB-XRT, and RCI.</caption><thead><tr><th>Technique</th><th>Reflection<br>hkil</th><th>Energy<br>[keV]</th><th>Wavelength<br>λ [pm]</th><th>Bragg Angle<br>Θ<sub>B</sub> [°]</th><th>Absorption<br>μt*</th></tr></thead><tbody><tr><td rowspan="2">L-XRT</td><td>1120</td><td>8.05 (Cu-Kα<sub>1</sub>)</td><td>154.06</td><td>28.89</td><td>9.1–13.6</td></tr><tr><td>1120</td><td>18.70</td><td>66.30</td><td>12.00</td><td>7.1</td></tr><tr><td rowspan="2">SWB-XRT</td><td>11100</td><td>12.47</td><td>99.45</td><td>10.37</td><td>20.5</td></tr><tr><td>1120</td><td>18.00</td><td>68.88</td><td>12.47</td><td>7.9</td></tr></tbody></table>
<table><caption>**Table 1.** *Cont.*</caption><thead><tr><th>Technique</th><th>Reflection<br>hkl</th><th>Energy<br>[keV]</th><th>Wavelength<br>λ [pm]</th><th>Bragg Angle<br>Θ<sub>B</sub> [°]</th><th>Absorption<br>μt *</th></tr></thead><tbody><tr><td rowspan="2">RCI</td><td>11̅20</td><td>30.00</td><td>41.33</td><td>7.45</td><td>1.9</td></tr><tr><td>1̅100</td><td>30.00</td><td>41.33</td><td>4.29</td><td>1.9</td></tr></tbody></table>

* μ = linear X-ray absorption coefficient, t = thickness of GaN substrate, - thickness of the GaN substrates investigated by L-XRT using Cu-Kα₁: 300 µm-450 µm, - thickness of the GaN substrate with the intensively investigated defect designated “HCD-A” that was the subject of detailed investigations using L-XRT (Cu-Kα₁) and, in addition, SWB-XRT as well as RCI: 300 µm.

## 2.3. Defect Selective Etching

Defect selective etching (DSE) was used as a complementary tool with respect to Bragg diffraction imaging. DSE is an established characterization technique for the analysis of defects in GaN thin films or bulk crystals. The method is suitable for revealing a wide range of different defects in GaN, such as dislocations, precipitates, grain boundaries, inversion domains, striations, and others. DSE is very sensitive to the local deformation of the GaN lattice and is mainly used to reveal dislocations. A particular advantage of the method is the relatively easy identification of the different observed types of TDs in the frequently used (0001) GaN crystals or substrates. Three types of TDs occur in (0001) GaN, namely, screw type (TSD) with Burgers vector **b** = ⟨0001⟩; edge type dislocation (TED), **b** = 1/3⟨1120⟩; and mixed-type (TMD), **b** = 1/3⟨1123⟩. When etching Ga-polar (0001) GaN, the etch pits differ significantly in terms of size, shape, and inclination angle of the side walls of the pits, depending on the doping and impurity level, the composition of the etchant, the etching temperature and time, and the type of the dislocations. However, for the three types of TDs in Ga-polar (0001) GaN, the size grading is generally such that TSDs form the largest etch pits, followed by TMDs, and TEDs form the smallest etch pits [26-29].

The approach of the DSE investigations in this work was to first use L-XRT to localize the HCDs on the Am-GaN wafers and then to apply target-oriented DSE at the relevant sites. For revealing dislocations, classical orthodox etching in molten E + M etch was used (eutectic mixture (E) of KOH + NaOH with 10% of MgO (M) [27]). In the present work, etching was performed on Ga-polar surface at a temperature of ~550 °C for 20 min. The temperature was measured inside a hot steel plate. After etching, the samples were examined by differential interference contrast (DIC) microscopy (Nikon Eclipse 80i). More details of the principle of DSE of GaN and other semiconductors can be found in the work of Weyher and Kelly [30].

# 3. Experimental Results

## 3.1. Results of X-ray Bragg Diffraction Imaging

### 3.1.1. Tracking of Honeycomb Defects (HCDs) by Lang X-ray Topography (L-XRT)

Figure 2 shows L-XRT images of three 1-inch Am-GaN substrates recorded using the 11̅20 reflection. The topographs in Figure 2a,b are from semi-insulating substrates and the topograph in Figure 2c is from an n-type one. We use the usual photographic convention, where higher X-ray intensity corresponds to increased blackening.

The experimental conditions described above coupled with the high crystalline quality of the crystals led to the fact that all the topographs were recorded under anomalous transmission (Borrmann effect) conditions. This means, as indicated above, that defects become visible through a locally reduced intensity. In the topographs of the Am-GaN substrates shown on Figure 2, different types of defects are visible on a µm-scale up to a mm-scale. The defects on the µm-scale are threading dislocations (TDs) of the screw (TSD), edge (TED), or mixed (TMD) type. In addition, macroscopic defects with mm-scale such as subsurface damage and scratches, traces of facets, growth bands, and various forms of dislocation clusters can be observed. The dislocation clusters can be divided into single locally isolated clusters; chains clusters; and a particular defect that consists of six clusters with a hexagonal, honeycomb-like, arrangement. These “honeycomb defects” (HCDs)
constitute the topic of the present work. Descriptions and discussions regarding the other
defect types visible in the overview topographs can be found elsewhere [8].

<figure><img src="image_5.png" /><figcaption>**Figure 2.** L-XRT 1120 reflection topographs of three GaN substrates with HCDs (marked by red hexagons); (a) and (b) are topographs of semi-insulating Am-GaN substrates, while (c) is a topograph of an n-type Am-GaN substrate. The defect cluster denoted as HCD-A in Figure 2b is subject of detailed studies in the following and is, therefore, highlighted.</figcaption></figure>

The HCDs we identified in the three substrates of Figure 2 are indicated by red hexagons. It can be seen that the HCD clusters display different sizes and distributions. In Figure 2b,c, the HCDs are rather isolated, whereas in Figure 2a, there is an area in the right part of the substrate where the HCDs are arranged close to one another. It should be noted that, from a statistical point of view, although the HCD is found in many Am-GaN substrates, the “density” of this defect is rather low; we observed, on the more than 40 studied substrates, a density between ~0.2 and 2.4 cm⁻².

3.1.2. Diffraction Imaging Contrast of the Honeycomb Defect (HCD) in Transmission L-XRT

Figure 3a shows the magnified cluster HCA-A of the overview topograph of Figure 2b (L-XRT, 1120 reflection). The entire HCD cluster has a “diameter” of ~1.8 mm. Six areas of distinct light-dark contrast and hexagonal arrangement are observed, each forming large, many-lobed contrast rosettes. The arrangement of the six contrast rosettes is in a hexagonal shape and the outer edges of the honeycomb run along the *a*-plane {1120} prismatic facets. Figure 3b shows the same topograph as Figure 3a with the six contrast rosettes numbered (1–6) as well as indexed planes for a better understanding of the following discussions in this paper. The contrast of the six rosettes does not look all the same, but there are three pairs of contrast rosettes with similar appearance, namely, (5, 6), (2, 3), and (1, 4), respectively. Using the orientation of the diffraction vector *h* and the numbering scheme of Figure 3b, the contrast rosettes can be described as follows: rosettes 5 and 6: two large lobes of dark contrast nearly perpendicular to *h* and two small lobes of bright contrast nearly parallel to *h*; rosettes 2 and 3: two large lobes of bright contrast nearly perpendicular to *h* and two small lobes of dark contrast nearly parallel to *h*. Rosettes 2 and 3 are a kind of contrast-inverted mirror image of rosettes 5 and 6 and rosettes 1 and 4 are an intermediate appearance of rosette types 5, 6 and 2, 3, but with much smaller dark and bright contrast lobes. Figure 3c shows magnified examples of each of the three different contrast rosette pairs. The size of the single rosettes exhibits a width of ~500 µm.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** The HCD cluster shown in (a) is the enlarged L-XRT image, 1120 reflection, of the one labeled HCD-A on Figure 2b; (b) for a better understanding of the explanations in the text, the contrast rosettes were numbered (1 to 6) and the crystallographic planes are indicated; (c) details of some of the corners of the hexagon are shown.</figcaption></figure>

The six contrast rosettes of the HCD are embedded in a matrix containing isolated TDs. The TDD around the HCD shown here is ~5 × 10⁴ cm⁻². Under Borrmann contrast imaging conditions, TDs become visible as many-lobed rosette-shaped contrasts in low defect crystals. Similar images have already been observed for TDs in highly absorbing semiconductors such as Ge [31], GaAs [32,33], CdTe [34], and ZnGeP₂ [35–37]. Recently, many-lobed rosette-shaped contrasts of TDs were reported for low-defect Am-GaN for the first time [8]. The appearance of the many-lobed contrast rosettes of the HCDs is comparable to the contrast rosettes of single TDs in Am-GaN, but their dimensions are about 10 times larger. There is no doubt, however, that the defect images of large, many-lobed contrast rosettes observed here by L-XRT under high µt conditions correspond to a strong and very complex Borrmann effect of a specific defect cluster in GaN, suggesting a similarity in defect type of HCDs to TDs that was not tangible at this point.

For further discussion of the HCD contrast rosettes in transmission XRT, see Appendix A.

### 3.1.3. Diffraction Imaging Contrast of the Honeycomb Defect (HCD) in Back-Reflection SWB-XRT

To obtain a simpler image of the HCD clusters, mainly by restricting to the surface, a SWB-XRT topograph was recorded in back-reflection geometry. Figure 4a shows the
0006 reflection topograph of the defect cluster labeled HCD-A in Figure 2b. Six contrasts displaying a hexagonal arrangement are clearly visible. The contrasts located at the six corners of the hexagon show two parts: a dark inner part and a bright, radially aligned, outer part. It has been shown that, using the SWB-XRT back-reflection technique, the threading dislocations (TDs) become visible as more or less point-shaped contrasts [38]. Three types of TDs, namely, screw, edge, or mixed, can be associated with varying size, brightness (light or dark), and uniformity of the contrast. The contrasts seen in Figure 4 are similar to what was observed before for TDs in Am-GaN crystals. The images we observe on the corners of the HCD could, therefore, result from the sum of the contrasts of aligned TDs.

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Back-reflection geometry 0006 SWB-XRT topograph showing (a) the hexagon labeled HCD-A in Figure 2b; (b) an enlarged view of corners 1 and 2.</figcaption></figure>

Let us note that some of the six images display a single bright straight line (contrasts 1, 5, and 6 in Figure 4a), whereas others show a series of bright lines, suggesting a different structure of the associated defects (contrasts 2–4 in Figure 4a). The described features are easy to observe on the magnified image of the honeycomb corners 1 and 2 shown in Figure 4b. The corner 1 contrast exhibits a bright straight radial line, whereas in corner 2, the bright contrast splits into several nearly radial lines that, as will be shown in Section 3.2, correspond to several dislocation bundles.

In addition, we observe images of TDs present outside the corners of the HCD, which display a black dot–white dot contrast. The orientation of the line joining the center of the dots varies between TDs to the other, depending on their Burgers vector direction [39].

### 3.1.4. Variations in the Honeycomb Defect Images

Although the HCDs observed by L-XRT in the different Am-GaN substrates show basically the same defect pattern, they display striking differences in the details of their image contrast. Figure 5 shows eight different HCDs, again marked with a red hexagon for better visibility. The occurrence of HCDs is independent of the doping type used in the substrates: the HCDs of Figure 5b–e were observed in n-type samples, whereas the HCDs of Figure 5a,f–h were found in semi-insulating (Mn-doped) substrates. The image size and contrast vary considerably from HCD to HCD. The sizes of the observed HCD vary from 50 µm to more than 2 mm. The rosettes present in the HCD images are not associated with the size; Figure 5a displays large and strong contrast rosettes, whereas Figure 5g shows small and weak contrasts, with both HCDs having a size of more than 2 mm. The contrast of the rosettes can be very different from one HCD to the other; Figure 5h shows, for instance, a HCD with a rather weak contrast. In addition to the six rosettes that form the hexagon, we observe, in some cases, other rosettes that are marked by small blue circles
(Figure 5c–f), which occur both inside and outside the HCD. The hexagonal symmetry of the HCD cluster does not appear to be affected by these additional rosettes or by the TDs present in the crystal matrix. Figure 5e,f show, for instance, that, inside the HCD indicated by the red lines parallel to the *a*-planes, six additional rosettes form a second hexagon (marked in green), which is rotated by 30 degrees and bounded by the six {1100} *m*-planes. In addition, another rosette occurs exactly in the center of this defect cluster. Quite different types of HCDs can arise simultaneously on the same single wafer, so that no direct correlation to specific growth conditions can be found for the individual types of HCDs. As the distortion fields observed on the topographies vary between the HCDs, the features leading to the hexagonal defects may also vary within a general framework.

<figure><img src="image_8.png" /><figcaption>**Figure 5.** L-XRT of eight different HCDs (marked with a red hexagon for better visibility), which exhibit different sizes and different contrasts with light and dark areas arranged hexagonally, as well as often large many-lobed rosettes. All topographs were measured using a type 1120 reflection. Additional clusters of threading dislocations are indicated by blue circles; the green hexagons mark a second kind of honeycomb bounded by the six {1100} *m*-planes. (a) HCD with strong contrast rosettes in semi-insulating Am-GaN, (b) HCD with strong contrast rosettes in n-type Am-GaN, (c–f) HCDs with in addition to the six rosettes that form the hexagon, other rosettes marked by small blue circles, (e,f) HCDs with six additional rosettes that form a second hexagon (marked in green), (g,h) HCDs with weak contrast rosettes in semi-insulating Am-GaN.</figcaption></figure>

## 3.2. *Results of Defect Selective Etching (DSE)*

A total of seven HCDs of different sizes, from different Am-GaN substrates, and with different contrast details observed by L-XRT, as described in Section 3.1.4, were investigated by DSE. There are various types of HCDs; examples with the typical differences will be described in detail below.

Figure 6a shows the DIC optical overview image after orthodox etching of the HCD-A (Am-GaN substrate from Figure 2b). The use of DSE allows to identify the defects responsible for the hexagonal pattern of six multi-lobed contrast rosettes observed on Bragg
diffraction images; these are bundles of TDs that are concentrated in the corners of the honeycomb in very restricted locations (marked by blue circles).

<figure><img src="image_9.png" /><figcaption>**Figure 6.** DIC optical images recorded after etching the HCD-A present in the Am-GaN substrate shown in Figure 2b: (a) overview image that allows identifying the bundles of TEDs that are concentrated in the corners of the honeycomb in very restricted locations (marked by blue circles); (b) enlarged images of the etch pits of the six dislocation bundles (DBs) located in the corners of the hexagon (a highly magnified image of the etched DB 5 is shown in Appendix B, Figure A4).</figcaption></figure>

The analysis of the pit size and morphology of the pits indicates that these are TEDs. Figure 6b shows scaled-up images of the etch pits of the six dislocation bundles (DBs). The number of edge dislocations present per DBs is high, amounting to 127–199 individual TEDs (average number: 154 TEDs). There are also differences in the distribution of the TEDs in the DBs, with two different groups. The DBs labeled 1, 5, and 6 consist of TEDs that line up predominantly along single, straight long chains (~150 µm in length) and only
a few single TEDs are located next to these chains. In a similar way, the edge dislocations of DBs 2, 3, and 4 are again predominantly arranged in several, nearly parallel, but shorter chains with a length of ~50 µm. The dislocation chains are radially oriented and parallel to one of the six crystallographic ⟨1100⟩ directions. Within the chains, the distance between TEDs varies between 0.5 µm and 2 µm. DSE shows, inside and outside of the hexagon, other pits associated with further TDs, which were already mentioned when describing the L-XRT and SWB-XRT results. Besides TEDs, there is a large number of TMDs and few TSDs in the matrix. The TMDs and TSDs differ from the TEDs in the DSE analysis by forming larger etch pits. In some cases, TDs bundle together in the matrix. However, these are different from the six DBs that form the HCD cluster. The number of TDs is much lower, with around ~10–30 dislocations per bundle and, besides TEDs, TMDs also occur in these bundles. Some of these additional DBs are marked by yellow circles or ellipsoids in Figure 6a.

We pointed out in Section 3.1.4 that the individual HCDs display quite different details in terms of size, contrast strength, contrast of the dislocation images associated with the corners of the hexagon, and extra images associated with other dislocations. These differences in the XRT images led to extending the DSE analysis to other HCDs. Figure 7 shows a DIC optical image of the etched HCD shown in Figure 5c; etch pits of the six DBs of the TED type are again observed, but their number now ranges from 47 to 63 per DB, with an average number of 53. These are arranged, within the bundles, in several short chains with lengths of 10–20 µm. The L-XRT image (Figure 5c) showed another rosette in the center of the hexagon. DSE reveals that these are also pits from a large bundle of TDs (labelled as 7 in Figure 7). The pits of this central DB are larger than the ones forming the HCD, this indicating that they are TMDs. Further smaller groups of clustered TMDs are also present in the interior of the HCD hexagon. We can also observe in the interior of this HCD two large deep pits formed on screw dislocations that appear in the center of former growth hillocks (TSD-1 and TSD-2, Figure 7a,b). To complete the description of the DSE analysis of this sample, we also observe defects related to growth facets (F) and scratches (S).

<figure><img src="image_10.png" /><figcaption>**Figure 7.** (a) DIC optical image of the etched HCD shown in Figure 5c. The corners of the hexagon are formed by TEDs. Growth facets (F) and scratches (S) are also visible. (b) Enlarged images of the etch pits formed on the bundles of TEDs (1-6). Other pits of dislocation bundles (like the one labeled 7) are of mixed type and two pits of single threading screw dislocations (TSD-1 and TSD-2) are observed.</figcaption></figure>

Figure 8 shows, as a third example, the DIC optical image after etching of the HCD shown in Figure 5h, which was detected in a semi-insulating Am-GaN substrate. Here, the etch pits corresponding to TEDs in the six corners of the honeycomb are spread over a large area (with a "diameter" of ~250 µm, Figure 8b) and show almost no chain-like arrangement.
The number of TEDs per DB is between 80 and 120 (average: 101). The count of TEDs for this HCD is not so precise owing to the wide dispersion of the distribution.

<figure><img src="image_11.png" /><figcaption>**Figure 8.** (a) DIC optical image of the etched HCD shown in Figure 5h. TEDs of each DB in six corners of the honeycomb cluster are spread over an area with a “diameter” of ~250 µm, as shown for DB-3 in the enlarged image (b).</figcaption></figure>

For the seven HCD clusters analyzed by DSE, it can be stated that, in each case, the six DBs forming the hexagon are always composed of pure TEDs. On the other hand, DBs inside and outside the HCD cluster are mostly bundles of TMDs, but some TED bundles were also observed. Pure TSDs were not observed as bundles; these occurred only in isolation.

The examples of HCDs examined here using DSE show the following characteristics:

(i) single TEDs in the DBs arrange as a few (or even single) long chains, many short chains, or are distributed over an area of up to 250 µm. There appears to be a correlation between the way the TEDs are located and the images they produce on the topographs: when more TEDs in the DBs are arranged as chains, stronger contrast is recorded on L-XRT;

(ii) the number of TEDs in the DBs differs significantly between HCDs. We can extract a correlation, shown in Figure 9, between the number of TEDs in the DBs and the size of the HCD clusters: the higher the number of TEDs in the corners of the hexagon, the larger the size of the corresponding HCD cluster.
<figure><img src="image_12.png" /><figcaption>**Figure 9.** Correlation between the number of TEDs in the DBs located in the corners of the HCD and the size of the HCD clusters: the larger an HCD cluster, the higher the number of TEDs in the corresponding six DBs. The blue hexagons represent the data points and the red line is a guide for the eyes.</figcaption></figure>

### 3.3. Use of RCI for a Quantitative Determination of the HCD's Associated Lattice Distortion

The synchrotron RCI measurements were performed on the HCD-A indicated in Figure 2b, which was also investigated previously by L-XRT, SWB-XRT, and DSE. For these studies, the 1120 reflection and 1100 reflection of the prismatic GaN planes were used. In order to try to facilitate the image interpretation, we decided to mainly use a 30 keV energy incident X-ray beam (μt ~1.9, 30 keV), which allows to avoid being in Borrmann conditions. A comparison of the images recorded under high absorption (μt~7.9, 18 keV) and the medium absorption displayed here is shown in Appendix A.

From the images recorded along the rocking curve, as already indicated, maps of intensity (INT) and, with the help of Gaussian fits of the rocking curves recorded in each individual pixel, maps of full-width at half-maximum (FWHM) as well as peak position (PPOS) were extracted.

The first, simple approach to gain information on the DBs forming the HCD cluster can be achieved through the evaluation of the original individual topographs recorded along the rocking curve for the 30 keV measurement of the 1120 reflection. This enables to observe "simple" images of defects because, if an image is recorded at an angular position ΔΘ very far away from Θcentre, the central position of the diffraction curve, the main contribution to Bragg diffraction originates from the distorted regions lying around the defect, and not from the perfect-crystal matrix. This XRT variant is termed "weak beam topography" [40,41].

Figure 10 shows a series of five individual images taken along the rocking curve for the area of HCD-A (see Figure 2b). The image shown in Figure 10c was taken in the center of the rocking curve (Θcentre); images from Figure 10b,d in the steeply rising and steeply falling flanks of the rocking curve, respectively; and images from Figure 10a (left tail) and Figure 10e (right tail) far away (ΔΘ~ ±60 μradian) from the maximum Θcentre of the rocking curve. The left tail image (Figure 10a) and right tail image (Figure 10e), respectively, show the weak beam images of the six DBs of the HCD-A. These topographs were taken outside the exact Bragg conditions and, therefore, no dynamical-theory related contrast is visible. Instead, only a "direct contrast image" of the lattice distortion associated with the DBs is observed. Further analysis of the original individual INT images captured along the rocking curve shows that the areas inside and outside the HCD become visible on different topographs (in particular Figure 10b,d). This indicates that the Am-GaN crystal regions inside and outside the honeycomb are in diffraction position for different angles, probably corresponding to a small misorientation.
<figure><img src="image_13.png" /></figure>

**Figure 10.** A series of individual images recorded on the HCD-A (from Figure 2b) for several angular positions along the rocking curve: going from its left tail, through the top, and to its right tail. The detailed description of the individual images (a–e) is given in the text. (c) image taken in the center of the rocking curve (Θcentre), (b,d) images from the steeply rising and steeply falling flanks of the rocking curve, (a) image from the left tail and (e) image from the right tail far away (ΔΘ~ ±60 μradian) from the maximum Θcentre of the rocking curve, respectively.

Figures 11 and 12 show the RCI maps at 30 keV of the 1120 and 1100 reflections, respectively. The INT maps presented in Figures 11a and 12a are, in spite of the difference in absorption, qualitatively similar to the topographs of the L-XRT and SWB-XRT experiments. Indeed, the HCD shows in the INT maps six large multi-lobed rosettes, which appear to be the images of the DBs composed of 127–199 individual TEDs. The formation of these distinctive contrast rosettes is, as already indicated, associated with dynamical theory-related interferences that produce high-structural-quality crystals. These interferences also show up, in a reverse way, in the FWHM map, but it can be observed that, in addition, the DBs located in the six corners of the hexagon exhibit a higher FWHM (10–15 μradians) than the one of most of the pixels in the image (~5 μradians). The PPOS map shows, around most of the corners of the hexagon, a blue-red contrast corresponding to a variation in peak position of ~30–40 μradians. It is known that a variation in peak position appears around a single dislocation and its value (typically 5 μradians) is proportional to the Burgers vector. Tsoutsouva et al. [25] have reported that, for the case of parallel dislocations located within one small area (a few μm² range), misorientation, measured a few μm away from this area, corresponds to an effective Burgers vector (the sum of the Burgers vector of the individual dislocations). In our case, this would lead to the conclusion that the sum of the Burgers vectors of the more than one hundred dislocations present on each of the corners is less than ten times the one of an individual dislocation. However, this procedure does not correspond to our case, because the dislocations present in the corners of HCDs are either arranged in “walls” or are located within relatively large areas (~200 μm).

We retrieve, in these PPOS maps, the variation in orientation, already observed in the images taken along the rocking curve (Figure 10), outside and inside of the HCD cluster. The mutually misoriented crystal volumes are separated by the formation of a kind of subgrain boundary lying along the notional connecting lines between the six DBs of the HCD cluster. Therefore, these subgrain boundaries run along the *a*-plane {1120} prismatic planes. A clear complete hexagon cannot be observed for each diffraction vector used, but, depending on the orientation with respect to the diffraction vector, only certain subgrain boundaries of the honeycomb become visible.
<figure><img src="image_14.png" /><figcaption>**Figure 11.** RCI maps of the 1̅120 reflection at 30 keV: (a) integrated intensity map (INT); (b) full-width at half-maximum map (FWHM); (c) peak position map (PPOS).</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>**Figure 12.** RCI maps of the 1̅100 reflection at 30 keV: (a) integrated intensity map (INT); (b) full-width at half-maximum map (FWHM); (c) peak position map (PPOS).</figcaption></figure>

Figure 13a,b show the PPOS maps of the 1̅120 and the 1̅100 reflection, respectively, with superimposed indications of the extra features we wish to point out. For better understanding, a scheme of the observed features and the corresponding crystallographic orientations is shown in Figure 13c. For the PPOS map corresponding to the 1̅120 reflection (Figure 13a), the subgrain boundaries of the (2110), (2110), (1210), and (1210) facets forming the hexagon are clearly visible, whereas for the PPOS map of the 1̅100 reflection (Figure 13b), the subgrain boundaries of the (1120) and (1120) facets stand out. The visibility of the subgrain boundaries is significantly reduced when they are perpendicular to the direction of the diffraction vector. This indicates that the rotation of these planes one side and the other of the subgrain boundary is such that the rotation axis has a main component parallel to the diffraction vector. From the observations, it follows that the misorientation of the volume is not preferentially in one of the six *a*-plane directions, but it is a three-dimensional misorientation. Let us note, however, that the average misorientation of the hexagon area relative to the surrounding Am-GaN matrix is small; the measurement of the average position inside and outside the honeycomb shows a variation of ~3 µradians for the 1̅120 reflection and ~8 µradians for the 1̅100 reflection. We can also observe that, except for the large six DBs in the corners of the hexagon with the high numbers of TDs, almost no TDs occur on the boundary courses of the subgrain boundaries between the DBs.
<figure><img src="image_16.png" /><figcaption>**Figure 13.** Peak position maps for (a) the 1120 reflection and (b) the 1100 reflection at 30 keV. The subgrain boundaries along the *a*-plane {1120} prismatic planes forming the hexagon are indicated as black dashed lines and the subgrains present within the HCD are indicated as orange and violet dashed lines. These last subgrain boundaries mostly go from one TD bundle to the other and can even form smaller honeycombs within the big one. The misorientations among these subgrains are extremely small, in the very few µradian range. (c) Stand-alone subgrain boundaries shown. For better understanding, the orientations of the *a*-plane facets (black) and *m*-plane facets (blue) are shown.</figcaption></figure>

In addition to the subgrain boundaries along the *a*-plane prismatic planes forming the hexagon (black dashed lines), the HCD cluster appears to be subdivided into sub-subgrains with boundaries (indicated by orange and violet dashed lines in Figure 13a-c), which mostly go from one TD to the other. These TDs act as geometrically necessary dislocations (GNDs; see, for instance, [42,43]) and can even form, as shown in Figure 13a-c, smaller hexagons. Most of the inner sub-subgrain boundaries run along the *a*-plane prismatic planes, but a few of them run along the *m*-plane (1100) prismatic plane (marked with a dashed violet line). These "sub-subgrains" within the HCD cluster exhibit very weak misorientations, in the µradian range, smaller than that of the entire HCD cluster relative to the surrounding Am-GaN matrix.

Figure 14 shows the magnified 1120 and 1100 PPOS map images of two opposite DBs of the HCD cluster, together with the corresponding DSE images, which indicate the occurrence of what we call a “limited size subgrain boundary”. Indeed, the edge dislocations that are aligned (Figure 14f) can, if they have the same Burgers vector (or at least when one Burgers vector is predominant), act as GNDs to produce a rotation of the lattice planes on one side of the alignment with respect to the other side. This variation in orientation is visible in Figure 14e as a red-blue contrast. However, in Figure 14d, we can notice that the contrast is actually more complicated, with a succession, within a globally “light blue” matrix, of “light blue-blue-red-blue-light blue” angular positions. This is expected for a limited size subgrain boundary because the induced misorientation must
come back to the matrix orientation when far from the subgrain, implying a rotation in the opposite sense, as schematically indicated in Figure 15. This configuration is confirmed by the weak beam images of the RCI (Figure 10a,e), which show a double point on each of the corners, indicating that there are two areas (the red or blue ones in Figure 15 schematic drawing) around the limited size subgrain boundary displaying the same orientation.

<figure><img src="image_17.png" /><figcaption>**Figure 14.** Magnified $\overline{1\overline{1}20}$ and $\overline{1}\overline{1}00$ PPOS map images of two DBs (indicated as 2 and 5 in, e.g., Figures 3b and 6a) of the HCD-A cluster, together with the corresponding DSE images (a highly magnified image of the etched DB 5 (Figure 14f) is shown in Appendix B, Figure A4). The explanation for the individual figures is given in the text. (a) $\overline{1\overline{1}20$ reflection PPOS map of DB-2, (b) $\overline{1}\overline{1}00$ reflection PPOS map of DB-2, (c) DSE image of DB-2, (d) $\overline{1\overline{1}20$ reflection PPOS map of DB-5, (e) $\overline{1}\overline{1}00$ reflection PPOS map of DB-5, (f) DSE image of DB-5.</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>**Figure 15.** Lattice plane rotations around a limited size subgrain boundary; (a) RCI PPOS measurement of the lattice rotation; (b) schematic drawing: black lines—no rotation of the lattice planes, blue lines—downward rotation of the lattice planes, red lines—upward rotation of the lattice planes, and purple line—limited size subgrain boundary.</figcaption></figure>
We can, in addition, infer from Figure 14 that the aligned or partially aligned TEDs (Figure 14f,c, respectively) that constitute the DBs located in the corners of the hexagon induce qualitatively similar distortions. However, the limited subgrain boundary behavior schematically presented in Figure 15 is clearly visible only for the aligned dislocations. On the other hand, the Burgers vectors of the dislocations located in the opposite corners of the hexagon do not appear to be aligned. If they were aligned, the contrast of the DBs images would be the same (or reversed if the Burgers vector is reversed) in Figure 14a,d, which is not observed.

4. Discussion

Bragg diffraction imaging and DSE analyses revealed that all of the HCDs occurring
in Am-GaN have a series of common features. Indeed, they are mainly constituted by
six dislocation bundles, each consisting of edge-type threading dislocations (TEDs) with
dislocation lines nearly parallel to the <0001> direction. The number of TEDs present in
each corner varies between HCDs, but the higher this number, the bigger the HCD size.
This size is a function of the number of edge dislocations, but their arrangement (lines,
succession of small lines, more diffuse cluster) or the number of other type of dislocations
in their neighbourhood appear to have a reduced influence.

The origin of the formation of the HCDs during ammonothermal GaN growth is not clear at the present stage. However, information can be provided on the stage of formation within the crystal growth process chain for these defect clusters. A complex technology is used to produce large area Am-GaN crystals, suitable for substrates with (0001) orientation, from small seeds. Process steps that have to be performed include the following: (1) multiple regrowth of slender GaN seed in the lateral <1120> directions; (2) joining two crystals by tiling technology; (3) continuing GaN grow in the lateral <1120> directions as well as in the vertical [0001] direction in multiple regrowth steps to a size sufficient for substrate fabrication; and (4) multiple regrowth in the vertical [0001] direction (with lateral growth blocked to increase only the thickness in the vertical [0001] direction) [5,6]. As has recently been shown, the history of the grown crystals, or more precisely, the individual process steps of crystal growth used for preparing these Am-GaN samples, can be traced in detail on the basis of observed XRT defect patterns [8]. Macroscopic defects, such as growth bands and chains of threading dislocations, related to regrowth steps as well as chains of threading dislocation bundles related to tiling seams of the joining process of two seed crystals for crystal area enlargement act as chronological markers of the Am-GaN growth. Figure 16 shows an HCD formed on chains of threading dislocation bundles (marked by red arrows) of a tiling seam. This observation is clear evidence that, chronologically, this HCD was formed after joining two seed crystals, i.e., in a process step that served to increase the thickness in the vertical [0001] direction. It can thus be ruled out that the HCD cluster formed in a process step for lateral crystal enlargement during growth in the <1120> directions. Further research is needed in order to explain the reason for the formation of the HCDs.

In the case of extensively studied HCD-A, the subgrain boundaries form along the
{1120} prismatic planes and the total crystal volume within the hexagon is slightly mis-
oriented relative to the surrounding Am-GaN matrix. The average misorientation of the
investigated HCD-A volume is very small. At this stage, it is not possible to say if this
low misorientation is caused by the DBs of the HCD. The origin can also be connected to
isolated threading dislocations or other DBs, which do not belong to the six DBs forming the
honeycombs. In addition, the volume of the hexagon is intersected by further sub-subgrain
boundaries with even smaller misorientations.
<figure><img src="image_19.png" /><figcaption>**Figure 16.** HCD (marked with a hexagon) formed on chains of threading dislocation bundles (red arrows) of a tiling seam. This allows a chronological assignment of the development of this defect cluster in the stages of seed enlargement.</figcaption></figure>

We do not have, presently, a convincing model to explain the HCDs and further work
is in progress to determine the Burgers vectors of the observed TEDs. Notwithstanding,
from the observed misorientation of the bounded subgrain boundaries in the RCI PPOS
maps, we can conclude that the Burgers vectors *b* in opposite corners cannot be collinear.
On the other hand, no contrast is observed for subgrain boundaries lying along a side of the
honeycomb when the diffraction vector *h* is perpendicular to the boundary. If we assume
that the subgrain rotation originates from the DBs that are present in the corners of the HCD,
which act as GNDs, this non-visibility could correspond to an "effective Burgers vector"
in the middle of this subgrain boundary perpendicular to the diffraction vector *h*. This
effective *b* can be approximated by the sum of two effective *b* corresponding to two DBs in
the neighboring corners, providing further clues for possible Burgers vectors directions.
However, all of these considerations remain highly speculative and more information is
required for preparing a model. We can nevertheless say that HCDs should result from
attractive and repulsive interactions between the TED bundles sitting in various hexagon
corners. Additionally, these interactions should have a different dependence on the size of
the hexagon in order to have a minimum energy size, which is related to the number of
dislocations in each corner.

**5. Conclusions**

In this paper, we analyzed defect clusters with hexagonal honeycomb-like arrangement
in ammonothermally grown GaN substrates. These defect clusters vary in size from
0.05 mm to 2 mm. In each of the hexagons' corners, we find a large concentration of
edge-type threading dislocations forming bundles. Within these bundles, the dislocations
are either grouped in areas or they show up as straight long chain-like radial alignments of
the same size that behave like limited subgrain boundaries. No defect concentration is seen
at the edges of the hexagons. The size of the honeycomb defect clusters correlates with
the number of dislocations in the bundles. At present, no clear model for the formation
of this complex defect type in GaN can be given. Further work is in progress to collect
more statistics for these defect clusters as well as to determine by transmission electron
microscopy the direction of the Burgers vectors of the edge-type threading dislocations
in the bundles forming the geometrically necessary dislocations. Understanding and,
possibly in future, preventing the development of HCDs during ammonothermal GaN
crystal growth is of great importance, both because Am-GaN crystals are presently used to manufacture electronic or optoelectronic devices and because they are used as seed crystals for new generations of GaN bulk crystals, e.g., by ammonothermal or hydride vapor-phase epitaxy (HVPE) growth processes.

**Author Contributions:** Conceptualization: L.K.; methodology: L.K., T.N.T.T.C. and J.L.W.; investigation—Bragg diffraction imaging: L.K., T.N.T.T.C., C.D., A.N.D. and J.B.; investigation—defect selective etching: J.L.W. and M.A.Z.; investigation—crystal growth and substrate preparation: R.K., T.S., K.G. and M.B.; writing—original draft preparation, L.K. and J.B.; writing—review and editing: L.K., T.N.T.T.C., J.L.W., J.S.-K., C.D., M.I., M.B. and J.B.; supervision: L.K. All authors have read and agreed to the published version of the manuscript.

**Funding:** The research leading to these results has received funding by the German Federal Ministry of Defense under the PureGaN project (E/E810/LC056/LF060). In addition, this research was funded by the Polish National Centre for Research and Development through project TECHMATSTRATEG-III/0003/2019-00. This research also received funding from the Department of the Navy, Office of Naval Research Global under ONRG award number N62909-21-1-2063, as well as the ECSEL Joint Undertaking (JU) under grant agreement No 101007310. The JU receives support from the European Union´s Horizon 2020 research and innovation program and from the National Authorities of Italy, Germany, France, Poland, Czech Republic, and The Netherlands.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data presented in this study are available upon request from the corresponding author when reasonably indicated.

**Conflicts of Interest:** The authors declare no conflict of interest.

## Abbreviations

The following abbreviations are used in this manuscript:

<table><tr><td>Am-GaN</td><td>ammonothermal gallium nitride</td></tr><tr><td>DB</td><td>dislocation bundle</td></tr><tr><td>DSE</td><td>defect selective etching</td></tr><tr><td>FWHM map</td><td>full-width at half-maximum map</td></tr><tr><td>GaN</td><td>gallium nitride</td></tr><tr><td>GND</td><td>geometrically necessary dislocation</td></tr><tr><td>HCD</td><td>honeycomb defect</td></tr><tr><td>HVPE-GaN</td><td>hydride vapor-phase epitaxy gallium nitride</td></tr><tr><td>INT map</td><td>integrated intensity map</td></tr><tr><td>L-XRT</td><td>Lang technique X-ray topography</td></tr><tr><td>PPOS map</td><td>peak position angle map</td></tr><tr><td>RCI</td><td>synchrotron monochromatic rocking curve imaging</td></tr><tr><td>SWB-XRT</td><td>synchrotron white-beam X-ray topography</td></tr><tr><td>TD</td><td>threading dislocation</td></tr><tr><td>TED</td><td>threading edge dislocation</td></tr><tr><td>TMD</td><td>threading mixed dislocation</td></tr><tr><td>TSD</td><td>threading screw dislocation</td></tr><tr><td>XRT</td><td>X-ray topography</td></tr></table>

## Appendix A. Contrast of the DBs Located at the HCDs' Corners

### Appendix A.1. Contrast of the DBs as a Function of the Reflection Type and Orientation of the Diffraction Vector

In general in XRT, the contrast image of dislocations of a crystal changes depending on the direction as well as on the type of the diffraction vector. This topic was also investigated below for the HCD defect using SWB-XRT. For this purpose, topographs were taken containing all reflections of type {1120} (a-plane) and {1100} (m-plane) prismatic planes.
Figure A1a shows topographs of the six type 1120 reflections. The SWB-XRT topographs are similar to the L-XRT image (Figure 3). However, the six contrast rosettes of the DBs appear less filigree and somewhat poorer in details. In these six topographs of the same type 1120 reflection, it can be seen that the contrast rosettes of the DBs have qualitatively the same appearance, but they are arranged differently depending on the direction of the diffraction vector. Thus, the contrast rosettes consisting of two large lobes with dark contrast and two small lobes with light contrast are always facing the diffraction vector, and vice versa. The contrast rosettes consisting of two large lobes with light contrast and two small lobes with dark contrast are always facing away from the diffraction vector. The DBs also show rosette-like contrasts in the topographs of the type 1100 reflections, but these are solely dark (Figure A1b). The only difference between the six individual rosette contrasts in a corresponding topograph is that their expansion changes with position relative to the diffraction vector. Those dark contrast rosettes facing the diffraction vector show a larger broadening than the rosettes facing away from the diffraction vector.

<figure><img src="image_20.png" /><figcaption>**Figure A1.** Transmission SWB-XRT of the HCD-A (Figure 2b): (a) 1120 reflections family, (b) 1100 reflection family. For the 1120 reflection family, the position of the light and dark contrast lobes changes depending on the direction of the diffraction vector. For the 1100 reflection family, only dark contrasts are observed; these show differences in their broadening depending on the orientation to the diffraction vector.</figcaption></figure>

### Appendix A.2. Contrast of the DBs as a Function of the X-ray Energy

As described for the L-XRT in Section 3.1.2 (and for the SWB-XRT in Appendix A1), depending on their orientation with respect to the diffraction vector, the rosettes do not all have the same appearance. On the two 1120 reflection INT maps (18 keV and 30 keV), there are three pairs of contrast rosettes visible that have a similar appearance, while for the 1100 reflection INT map imaged at 30 keV, there are two groups of three with similar contrasts. The comparison of the two INT maps of the 1120 reflection, measured at 18 keV and 30 keV, respectively, shows the characteristic of the Borrmann effect on the contrast pattern of the defects for Am-GaN investigated here as a function of energy. At first sight, the contrast patterns on the two 1120 reflection INT maps are similar. However, closer inspection reveals that the contrasts in the INT map for the 18 keV imaging with a µt of 7.9 are much more complex and filigree owing to the occurrence of anomalous transmission (Borrmann effect), which is not present in the image at 30 keV with a µt of 1.9. This is especially the case for the lobes of the six large contrast rosettes of the HCD, but also for
the matrix of TDs inside and outside the HCD cluster, which exhibit rosette-like contrasts, albeit much smaller. For the INT maps under the 30 keV conditions, the lobes of the six contrast rosettes of the HCD appear much plainer.

<figure><img src="image_21.png" /><figcaption>**Figure A2.** RCI maps of the 1120 reflection at 18 keV: (a) integrated intensity map (INT); (b) full-width at half-maximum map (FWHM); (c) peak position map (PPOS) (compare to Figure 11).</figcaption></figure>

### Appendix A.3. Phenomenological Model for the Emergence of Contrast Differences in Symmetrically Equivalent DBs of the HCDs

The defect pattern of the DBs leads to interesting diffraction imaging contrasts, which can be summarized phenomenologically as follows: the contrasts of the DBs resemble the dynamic contrast rosettes of single TEDs, but these are significantly larger owing to the accumulation of sometimes almost 200 individual TEDs in one DB. The observation of these large dynamic contrast rosettes in strongly absorbing GaN is possible, because the actual matrix of the investigated material is of extraordinary structural perfection and, therefore, a pronounced Borrmann effect is present. Although the individual six DBs are comparable in first approximation with respect to number and distribution in the respective HCDs, the contrast rosettes show phenomenologically clearly different appearances. As described in Section 3.1.2, among the DBs, groups of the same appearance are formed and there is a dependence of the appearance on the direction and distance of the diffraction vector. As shown in Figure A3, using the 1120 reflection as an example, the appearance of individual contrast rosettes can be inverse. These extraordinarily large contrast rosettes are created by the dynamically diffracted X-ray beam, which impinges on very large symmetrically arranged DBs. The X-ray beam has to pass through two subgrain boundaries, i.e., when entering and leaving the honeycomb volume. The X-ray beam has to pass through two subgrain boundaries, i.e., when entering and when leaving the honeycomb volume. The matrix volume and the honeycomb volume, each considered separately, are of high structural perfection, except that the two volumes are minimally misoriented with respect to each other. To understand the complicated Borrmann contrasts for the HCD cluster more precisely, appropriate simulations for dynamic diffraction conditions would have to be performed. However, this is outside the scope of this work.
<figure><img src="image_22.png" /><figcaption>**Figure A3.** Phenomenological model concerning the contrast inversion of the rosettes of the DBs owing to two subgrain boundaries that have to be passed by the X-ray beam (see text for explanation).</figcaption></figure>

## Appendix B. Distances of the TEDs in the Limited Size Subgrain Boundaries

<figure><img src="image_23.png" /><figcaption>**Figure A4.** (a) DIC optical image of the etched DB-5 of the HCD-A, forming a limited size subgrain boundary (compare Figures 6 and 14f); (b) magnified partial view—these are threading edge dislocations (TEDs), which are predominantly aligned in a line in the *a*-direction. The distance in the line between individual TEDs is between ~0.5 and 1 µm. Owing to the small distances, trenches are formed during etching. The smaller the distance between the TEDs, the wider the trenches are etched. In addition to the TEDs, a single threading mixed dislocation (TMD) can also be seen in Figure A4a.</figcaption></figure>

## References

1. Meneghini, M.; De Santi, C.; Abid, I.; Buffolo, M.; Cioni, M.; Abdul Khadar, R.; Nela, L.; Zagni, N.; Chini, A.; Medjdoub, F.; et al. GaN-based power devices: Physics, reliability, and perspectives. *J. Appl. Phys.* **2021**, *130*, 181101. [CrossRef]

2. Bojarska-Cieślińska, A.; Marona, Ł.; Smalc-Koziorowska, J.; Grzanka, S.; Weyher, J.; Schiavon, D.; Perlin, P. Role of dislocations in nitride laser diodes with different indium content. *Sci. Rep.* **2021**, *11*, 21. [CrossRef]
3. Zhang, Y.; Zubair, A.; Liu, Z.; Xiao, M.; Perozek, J.A.; Ma, Y.; Palacios, T. GaN FinFETs and trigate devices for power and RF applications: Review and perspective. *Semicond. Sci. Technol.* **2021**, *36*, 054001. [CrossRef]

4. Bockowski, M.; Iwinska, M.; Amilusik, M.; Fijalkowski, M.; Lucznik, B.; Sochacki, T. Challenges and future perspectives in HVPE-GaN growth on ammonothermal GaN seeds. *Semicond. Sci. Technol.* **2016**, *31*, 093002. [CrossRef]

5. Kucharski, R.; Sochacki, T.; Lucznik, B.; Bockowski, M. Growth of bulk GaN crystals. *J. Appl. Phys.* **2020**, *128*, 050902. [CrossRef]

6. Grabianska, K.; Jaroszynski, P.; Sidor, A.; Bockowski, M.; Iwinska, M. GaN Single Crystalline Substrates by Ammonothermal and HVPE Methods for Electronic Devices. *Electronics* **2020**, *9*, 1342. [CrossRef]

7. UNIPRESS. Available online: http://www.unipress.waw.pl/growth/index.php/ammono-gan-wafers-sales (accessed on 31 August 2022).

8. Kirste, L.; Grabianska, K.; Kucharski, R.; Sochacki, T.; Lucznik, B.; Bockowski, M. Structural Analysis of Low Defect Ammonothermally Grown GaN Wafers by Borrmann Effect X-ray Topography. *Materials* **2021**, *14*, 5472. [CrossRef] [PubMed]

9. Authier, A. *Dynamical Theory of X-ray Diffraction*, 2nd ed.; Oxford University Press: Oxford, UK, 2003. [CrossRef]

10. Doradziński, R.; Dwiliński, R.; Garczyński, J.; Sierzputowski, L.P.; Kanbara, Y. Ammonothermal Growth of GaN Under Ammono-Basic Conditions. In *Technology of Gallium Nitride Crystal Growth*; Ehrentraut, D., Meissner, E., Bockowski, M., Eds.; Springer: Heidelberg, Germany, 2010; pp. 137–160. [CrossRef]

11. Hashimoto, T.; Nakamura, S. A Pathway Toward Bulk Growth of GaN by the Ammonothermal Method. In *Technology of Gallium Nitride Crystal Growth*; Ehrentraut, D., Meissner, E., Bockowski, M., Eds.; Springer: Heidelberg, Germany, 2010; pp. 161–182. [CrossRef]

12. Zajac, M.; Kucharski, R.; Grabianska, K.; Gwardys-Bak, A.; Puchalski, A.; Wasik, D.; Litwin-Staszewska, E.; Piotrzkowski, R.; Domagala, J.Z.; Bockowski, M. Basic ammonothermal growth of Gallium Nitride—State of the art, challenges, perspectives. *Prog. Cryst. Growth Charact. Mater.* **2018**, *64*, 63–74. [CrossRef]

13. Klapper, H. Defects in Non-metal Crystals. In *Characterization of Crystal Growth Defects by X-Ray Methods*; Tanner, B.K., Bowen, D.K., Eds.; Plenum Press: New York, NY, USA, 1980; pp. 133–160. [CrossRef]

14. Lübbert, D.; Baumbach, T.; Härtwig, J.; Boller, E.; Pernot, E. µm-resolved high-resolution X-ray diffraction imaging for semiconductor quality control. *Nucl. Inst. Meth. Phys. Res. B* **2000**, *160*, 521–527. [CrossRef]

15. Lang, A.R. Direct Observation of Individual Dislocations by X-Ray Diffraction. *J. Appl. Phys.* **1958**, *29*, 597–598. [CrossRef]

16. Borrmann, G. Die Absorption von Röntgenstrahlen im Fall der Interferenz. *Z. Phys.* **1950**, *127*, 297–323. [CrossRef]

17. Rozgonyi, G.A.; Miller, D.C. X-ray characterization of stresses and defects in thin films and substrates. *Thin Solid Films* **1976**, *31*, 185–216. [CrossRef]

18. Rack, A.; Weitkamp, T.; Bauer, S.; Modregger, P.; Cecilia, A.; Rolo, T.; Rack, T.; Haas, D.; Simon, R.; Hedele, R.; et al. The micro-imaging station of the TopoTomo beamline at the ANKA synchrotron light source. *Nucl. Inst. Meth. Phys. Res. B* **2009**, *267*, 1978–1988. [CrossRef]

19. Tuomi, T.; Naukkarinen, K.; Rabe, P. Use of synchrotron radiation in X-ray diffraction topography. *Phys. Status Solidi A* **1974**, *25*, 93–106. [CrossRef]

20. Miltat, H. White beam synchrotron radiation topography. In *Characterization of Crystal Growth Defects by X-ray Methods*; Tanner, B.K., Bowen, D.K., Eds.; Plenum Press: New York, NY, USA, 1980; pp. 401–420. [CrossRef]

21. Huang, X.R. LauePt, a graphical-user-interface program for simulating and analyzing white-beam X-ray diffraction Laue patterns. *J. Appl. Cryst.* **2010**, *43*, 926–928. [CrossRef]

22. Tran Thi Caliste, T.N. RCIA Software. Available online: https://gitlab.com/l_sim/scripts/rcia (accessed on 31 August 2022).

23. Tran Thi, T.N.; Morse, J.; Caliste, D.; Fernandez, B.; Eon, D.; Härtwig, J.; Barbay, C.; Mer-Calfati, C.; Tranchant, N.; Arnault, J.C.; et al. Synchrotron Bragg diffraction imaging characterization of synthetic diamond crystals for optical and electronic power device applications. *J. Appl. Cryst.* **2017**, *50*, 561–569. [CrossRef]

24. Klapper, H. X-Ray Diffraction Topography: Principles and Techniques. In *X-Ray and Neutron Dynamical Diffraction: Theory and Applications*; Authier, A., Lagomarsino, S., Tanner, B.K., Eds.; Plenum Press: New York, NY, USA, 1996; pp. 137–146. [CrossRef]

25. Tsoutsouva, M.G.; Amaral De Oliveira, V.; Baruchel, J.; Camel, D.; Marie, B.; Lafford, T. Characterisation of defects in mono-like silicon for photovoltaic applications using X-ray Bragg diffraction imaging. *J. Appl. Cryst.* **2015**, *48*, 645–654. [CrossRef]

26. Weyher, J.L.; Brown, P.D.; Rouvière, J.L.; Wosinski, T.; Zauner, A.R.A.; Grzegory, I. Recent advances in defect-selective etching of GaN. *J. Cryst. Growth* **2000**, *210*, 151–156. [CrossRef]

27. Kamler, G.; Weyher, J.L.; Grzegory, I.; Jezierska, E.; Wosinski, T. Defect-selective etching of GaN in a modified molten bases system. *J. Cryst. Growth* **2002**, *246*, 21–24. [CrossRef]

28. Weyher, J.L.; Lazar, S.; Macht, L.; Liliental-Weber, Z.; Molnar, R.J.; Müller, S.; Sivel, V.G.M.; Nowak, G.; Grzegory, I. Orthodox etching of HVPE-grown GaN. *J. Cryst. Growth* **2007**, *305*, 384–392. [CrossRef]

29. Knoke, I.Y.; Berwian, P.; Meissner, E.; Friedrich, J.; Strunk, H.P.; Müller, G. Selective etching of dislocations in GaN grown by low-pressure solution growth. *J. Cryst. Growth* **2010**, *312*, 3040–3045. [CrossRef]

30. Weyher, J.L.; Kelly, J.J. Defect-selective etching of semiconductors. In *Springer Handbook of Crystal Growth*; Dhanaraj, G., Byrappa, K., Prasad, V., Dudley, M., Eds.; Springer: Berlin/Heidelberg, Germany, 2010; pp. 1453–1476. [CrossRef]

31. Okkerse, B. Anomalous transmission of X-rays in germanium. *Philips Res. Rpts.* **1962**, *17*, 464–478.

32. Builov, A.N.; Danil’chuk, L.N. Edge dislocations with large burgers vector in sphalerite crystals. *Tech. Phys. Lett.* **2002**, *28*, 762–764. [CrossRef]
33. Möck, P. Analysis of thermal-treatment-induced dislocation bundles in GaAs wafers by means of X-ray transmission topography and complementary methods. *J. Appl. Cryst.* **2001**, *34*, 65–75. [CrossRef]

34. Yildirim, C.; Gout, E.; Pagot, A.; Tran Thi Caliste, T.N.; Baruchel, J.; Brellier, D.; Ballet, P. Bragg Diffraction Imaging of CdZnTe Single Crystals. *J. Electron. Mater.* **2020**, *49*, 4550–4556. [CrossRef]

35. Okunev, A.O.; Verozubova, G.A.; Trukhanov, E.M.; Dzjuba, I.V.; Galtier, P.R.J.; Said Hassani, S.A. Study of structural defects in ZnGeP₂ crystals by X-ray topography based on the Borrmann effect. *J. Appl. Cryst.* **2009**, *42*, 994–998. [CrossRef]

36. Okunev, A.O.; Verozubova, G.A. X-ray topography contrast from edge dislocations in ZnGeP₂ single crystals. *J. Appl. Cryst.* **2015**, *48*, 1228–1233. [CrossRef]

37. Lei, Z.; Okunev, A.; Zhu, C.; Verozubova, G.; Yang, C. Low-angle boundaries in ZnGeP₂ single crystals. *J. Appl. Cryst.* **2018**, *51*, 361–367. [CrossRef]

38. Sintonen, S.; Rudziński, M.; Suihkonen, S.; Jussila, H.; Knetzger, M.; Meissner, E.; Danilewsky, A.; Tuomi, T.O.; Lipsanen, H. Synchrotron radiation x-ray topography and defect selective etching analysis of threading dislocations in GaN. *J. Appl. Phys.* **2014**, *116*, 083504. [CrossRef]

39. Yao, Y.; Ishikawa, Y.; Sugawara, Y.; Takahashi, Y.; Hirano, K. Observation of Threading Dislocations in Ammonothermal Gallium Nitride Single Crystal Using Synchrotron X-ray Topography. *J. Electron. Mater.* **2018**, *47*, 5007–5012. [CrossRef]

40. Barrett, R.; Baruchel, J.; Härtwig, J.; Zontone, F. The present status of the ESRF diffraction topography beamline: New experimental results. *J. Phys. D Appl. Phys.* **1995**, *28*, A25O–A255. [CrossRef]

41. Baruchel, J.; Härtwig, J. X-Ray Topography. In *Encyclopedia of Condensed Matter Physics*; Bassani, F., Liedl, G.L., Wyder, P., Eds.; Elsevier: Amsterdam, The Netherlands, 2005; pp. 342–348. [CrossRef]

42. Hughes, D.A.; Hansen, N.; Bammann, D.J. Geometrically necessary boundaries, incidental dislocation boundaries and geometrically necessary dislocations. *Scr. Mater.* **2003**, *48*, 147–153. [CrossRef]

43. Cai, W.; Nix, W.D.; *Imperfections in Crystalline Solids*; Cambridge University Press: Cambridge, UK, 2016. [CrossRef]