pith. sign in

arxiv: 2604.24604 · v1 · submitted 2026-04-27 · ❄️ cond-mat.mtrl-sci

Step- and terrace-resolved crystal truncation rod scattering from vicinal surfaces under coherent heteroepitaxy

Pith reviewed 2026-05-08 03:05 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords crystal truncation rod scatteringvicinal surfacesheteroepitaxycoherent strainelastic deformationtriclinic distortionterrace orderingInGaN/GaN
0
0 comments X

The pith

An elasticity-based model for coherently strained films on vicinal surfaces adds a shear-induced triclinic deformation that selectively alters non-specular crystal truncation rods.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a unified theory of crystal truncation rod scattering from vicinal surfaces that includes a coherently strained heteroepitaxial film. The formalism combines film-induced interference fringes, complete elastic lattice distortion, terrace ordering, surface reconstruction, and real-time growth evolution. Direct comparison of the Nagai model with an elasticity-based model finds nearly identical lattice tilt in both, but the elasticity approach introduces an extra triclinic deformation caused by shear strain. This added deformation barely changes specular rods yet strongly modifies non-specular rods, turning the latter into a direct probe of the film's full elastic state. The same framework preserves its original sensitivity to step and terrace details, as illustrated by representative calculations for InGaN/GaN that support quantitative analysis of both static and dynamic measurements.

Core claim

We develop a general theory of crystal truncation rod (CTR) scattering from vicinal surfaces with a coherently strained heteroepitaxial film. The formalism incorporates film-induced interference fringes, full elastic lattice distortion, terrace ordering, surface reconstruction, and real-time growth evolution within a unified description. Comparison between Nagai model and elasticity-based model shows that the lattice tilt is nearly identical in the two approaches, whereas the elasticity-based model predicts an additional triclinic deformation arising from shear strain. This deformation has little effect on specular CTRs but strongly modifies non-specular rods, making them a sensitive probe.

What carries the argument

The elasticity-based model of full lattice distortion in a coherently strained heteroepitaxial film on a vicinal surface, which introduces a triclinic deformation from shear strain in addition to the lattice tilt shared with the Nagai model.

Load-bearing premise

The heteroepitaxial film remains fully coherent with no plastic relaxation or defects, allowing direct comparison of Nagai and elasticity models without additional relaxation mechanisms.

What would settle it

Non-specular CTR intensity profiles measured on a confirmed-coherent InGaN/GaN film that match Nagai-model predictions but deviate from the elasticity-based predictions of triclinic deformation effects.

Figures

Figures reproduced from arXiv: 2604.24604 by Erqi Xu, Guangxu Ju, Jiale Wang, Jiaqing Yue, Junlin Wu, Qihui Lin, Zihao Xu.

Figure 1
Figure 1. Figure 1: FIG. 1. Sketch of (a) the Nagai model and (b) the elasticity-b view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Lattice tilt angle view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Typical CTRs from a vicinal surface with a het view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Calculated reflectivities for the (00 view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Schematic illustration of the unit-cell configurati view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Calculated reflectivities for the (00 view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Calculated CTR reflectivities for a vicinal GaN(0001 view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Calculated reflectivities for the (00 view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Calculated reflectivities for the (00 view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Color map of the logarithm of the CTR reflectivity for view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Simulated CTR intensity during growth in the ( view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. Simulated time evolution of CTR intensity at fixed view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13. Complex-plane trajectory of (a) view at source ↗
Figure 14
Figure 14. Figure 14: FIG. 14. Complex-plane trajectories of the bulk structure fa view at source ↗
Figure 15
Figure 15. Figure 15: FIG. 15. Calculated amplitude and phase of the bulk and epilay view at source ↗
Figure 16
Figure 16. Figure 16: FIG. 16. Complex-plane trajectories of view at source ↗
read the original abstract

We develop a general theory of crystal truncation rod (CTR) scattering from vicinal surfaces with a coherently strained heteroepitaxial film. The formalism incorporates film-induced interference fringes, full elastic lattice distortion, terrace ordering, surface reconstruction, and real-time growth evolution within a unified description. Comparison between Nagai model and elasticity-based model shows that the lattice tilt is nearly identical in the two approaches, whereas the elasticitybased model predicts an additional triclinic deformation arising from shear strain. This deformation has little effect on specular CTRs but strongly modifies non-specular rods, making them a sensitive probe of the full elastic state of coherent epitaxial films. We further show that the characteristic sensitivity of vicinal CTRs to terrace ordering, surface reconstruction, and terrace-resolved compositional modification remains robust in the presence of a coherent film. Representative calculations for InGaN/GaN demonstrate that the framework enables quantitative interpretation of both static and real-time CTR measurements and provides access to step- and terrace-resolved structural and kinetic information during heteroepitaxial growth.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 3 minor

Summary. The paper develops a general theory of crystal truncation rod (CTR) scattering from vicinal surfaces with a coherently strained heteroepitaxial film. The formalism unifies film-induced interference fringes, full elastic lattice distortion, terrace ordering, surface reconstruction, and real-time growth evolution. Comparison between the Nagai model and an elasticity-based model shows nearly identical lattice tilts, but the elasticity model introduces an additional triclinic deformation from shear strain. This has little effect on specular CTRs but strongly modifies non-specular rods, making them a sensitive probe of the full elastic state. Representative calculations for InGaN/GaN demonstrate utility for quantitative interpretation of static and real-time CTR measurements, providing step- and terrace-resolved structural and kinetic information.

Significance. If the derivations hold, this work is significant for advancing quantitative CTR analysis in heteroepitaxy on vicinal surfaces by providing a unified framework that incorporates multiple effects. The key insight that non-specular rods probe the full elastic state (including triclinic shear) is valuable for coherent films. Strengths include explicit derivations of the distorted lattice and CTR intensity formulas, incorporation of established formalisms for terraces and reconstruction, and representative InGaN/GaN calculations that support reproducibility.

minor comments (3)
  1. The abstract is information-dense; consider separating the description of the unified formalism from the model comparison results for improved readability.
  2. In the InGaN/GaN example calculations, figure captions or the text should explicitly list the input parameters (e.g., film thickness, misfit strain, terrace width) to facilitate reproduction of the shown CTR profiles.
  3. Ensure consistent notation for strain components and lattice distortions when contrasting the Nagai and elasticity models; a side-by-side table of key parameters would aid clarity.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive assessment of our manuscript, recognition of its significance for quantitative CTR analysis in heteroepitaxy, and recommendation for minor revision. We are pleased that the unified framework, the insight regarding non-specular rods as probes of the full elastic state, and the representative InGaN/GaN calculations were viewed favorably.

Circularity Check

0 steps flagged

Derivation chain is self-contained; no circular reductions identified

full rationale

The paper develops an explicit formalism for CTR scattering from vicinal surfaces with coherent heteroepitaxial films, incorporating elastic lattice distortion, terrace ordering, and growth evolution. The central comparison between the Nagai tilt model and the elasticity-based model derives the additional triclinic shear deformation directly from continuum elasticity under the stated full-coherence assumption, with explicit statements that this deformation has negligible effect on specular CTRs but modifies non-specular rods. No equations reduce by construction to fitted parameters or prior self-citations; the InGaN/GaN calculations are presented as representative applications of the derived intensities. The coherence premise is declared as the modeling scope rather than smuggled in. This satisfies the criteria for an independent derivation without load-bearing self-references or ansatz smuggling.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on standard elasticity theory applied to coherent films and conventional CTR scattering principles; no new entities or fitted parameters are introduced in the abstract.

axioms (1)
  • domain assumption Elasticity theory fully describes the lattice distortion in coherently strained heteroepitaxial films without relaxation.
    Invoked when comparing Nagai and elasticity-based models and when claiming sensitivity of non-specular rods to the full elastic state.

pith-pipeline@v0.9.0 · 5504 in / 1272 out tokens · 67539 ms · 2026-05-08T03:05:07.848709+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

60 extracted references · 1 canonical work pages

  1. [1]

    Compared with an exactly ori- ented surface, the vicinal geometry separates rods that would otherwise overlap and provides enhanced sensitiv- ity to step-resolved surface structure. Relative to our previous treatment of a vicinal substrate without an epi- taxial layer[ 15], the key new feature here is the emergence of interference fringes around each Brag...

  2. [2]

    ( 14) and ensures proper normalization of the scattering amplitude

    arises from the summation over the step periodicity ( N3) in Eq. ( 14) and ensures proper normalization of the scattering amplitude. For a given Bragg peak H0K0L0, the corresponding CTR satisfies K = K0 + L − L0 M , (18) reflecting its orientation along the vicinal surface nor- mal. In the limit M → ∞ , corresponding to an exactly oriented surface, the expr...

  3. [3]

    from the substrate region ( m ≤ 0) into the film region ( m = 1 to J M), thereby preserving the phase continuity across the substrate-film interface while intro- ducing distinct phase factors in the epilayer containing lattice distortion. This yields repi = rf Fepi M JM∑ m=1 Y m epi = rf Fepi M Yepi(Y JM epi − 1) Yepi − 1 , (20) where J denotes the film thic...

  4. [4]

    The surface reconstruction is taken to be 3H(T1)[ 9]

    For all calculations presented in this work, the growth temperature is fixed at 1076 K. The surface reconstruction is taken to be 3H(T1)[ 9]. A photon energy of 25 . 78 keV ( λ = 0. 4809 ˚ A) is used, con- sistent with recent experiments.[ 9] Atomic form factors are taken from Ref. 40, including resonant corrections at this energy.[ 41] Absorption lengths ...

  5. [5]

    Specifically, the modified reconstructions on the α and β terraces are denoted as α ∗ and β ∗ , respectively

    Building upon the α/β terrace-resolved CTR formalism developed above, we in- troduce modified reconstruction domains that account for terrace-dependent indium incorporation. Specifically, the modified reconstructions on the α and β terraces are denoted as α ∗ and β ∗ , respectively. Within a double-step spacing of M unit cells in the y direction, the modified...

  6. [6]

    Mw0 gives the starting row index of region w; Mw gives the number of rows in that region

    for each reconstruction w = α , α ∗ , β , and β ∗ . Mw0 gives the starting row index of region w; Mw gives the number of rows in that region. w M w0 Mw α (J − 1)M N − Nα ∗ α ∗ (J − 1)M + N − Nα ∗ Nα ∗ β (J − 1)M + N M − N − Nβ ∗ β ∗ (J − 1)M + M − Nβ ∗ Nβ ∗ near the step edges. Figure 9 shows the calculated CTR profiles for a vicinal GaN(0001) surface in t...

  7. [7]

    These minima are determined by the bulk GaN crystal structure (see Appendix B 3)

    The phase of the substrate contribution changes by π not only when crossing the Bragg peak, but also at intermediate positions where the scattering amplitude passes through minima. These minima are determined by the bulk GaN crystal structure (see Appendix B 3). In addition, due to terrace-dependent height offsets, the effective film thick- ness deviates fro...

  8. [8]

    V A or in situ laser reflectometry)

    provides a direct route to determine the 16 indium composition xIn without full curve fitting, pro- vided that the growth rate is independently known (e.g., from Sec. V A or in situ laser reflectometry). The key quantity is the oscillation period ∆ d, which can be di- rectly obtained from the temporal evolution of the CTR intensity. The oscillation period ∆...

  9. [9]

    All asymptotic orders below are taken with respect to the small off-cut angle θ

    Details of the elastic-based model The indices 1 , 2, 3 and x, y, z are used interchangeably. All asymptotic orders below are taken with respect to the small off-cut angle θ. The normal strains ǫii are treated as finite parameters, while the shear components scale as ǫ23, ǫ ′ 23 = O(θ). The two coordinate systems mentioned in Sec. II B are related by   ˆe...

  10. [10]

    Thus, the solution is uniquely determined

    and ( 9), ǫx′x′ = ǫm1; ǫy′y′ = ǫm2; ǫx′y′ = 0; σx′z′ = 0; σy′z′ = 0; σz′z′ = 0. Thus, the solution is uniquely determined. For later use, we summarize several relations: ǫ23 = σ23/ 2c44 = − sin θ0 epi cos θ0 epiσ ′ 22/ 2c44; (A6) ǫ′ 23 = sin θ0 epi cos θ0 epi(ǫ22 − ǫ33) + cos 2θ0 epiǫ23; (A7) ǫ′ 11 = ǫ11; (A8) ǫ′ 22 = (cos θ0 epi)2ǫ22 − sin 2θ0 epiǫ23 + (...

  11. [11]

    Here we assume that the off-cut angle θ ∼ 1/M is small

    Small-quantity expansion The purpose of this subsection is twofold: first, to show that the Nagai and elasticity-based models give the same lattice tilt to leading order; and second, to identify which components of the deformation matrix control the differ- ence in CTR scattering. Here we assume that the off-cut angle θ ∼ 1/M is small. In practice, θ usually...

  12. [12]

    + cos γǫ′ 23 cos γ (1 + ǫ′

  13. [13]

    (A12) Under the small-quantity approximation, γ = θ0 epi − θ + ǫ′ 23 1 + ǫ′ 22 + O(θ3)

    − sin γǫ′ 23 = tan ( θ0 epi − θ ) . (A12) Under the small-quantity approximation, γ = θ0 epi − θ + ǫ′ 23 1 + ǫ′ 22 + O(θ3). (A13) The angle δ is the angle between the substrate b vec- tor and the corresponding bepi vector in the strained epi- layer. It describes the lattice tilt of the epilayer. Geo- metrically, δ = γ + arctan[ǫ23/ (1 + ǫ22)]. (A14) 18 Th...

  14. [14]

    Substituting Eq

    has been used. Substituting Eq. ( A17) and (A18) then yields δN = δe + O(θ2). (A23) In Sec. III B, we noted that the difference in CTR in- tensity between two crystal models is negligible in the specular case but more significant in the non-specular case. The underlying reason is that the scalar product q ·(∆ ′ e − ∆ ′ N ) is close to zero when qy = 2 πK/b ...

  15. [15]

    In the expression for the total reflectivity amplitude rtot, the only term that depends on J is Y JM epi , which appears in both repi and rrec

    Phase accumulation and thickness oscillation period We first analyze the dependence of the CTR phase on the film thickness J. In the expression for the total reflectivity amplitude rtot, the only term that depends on J is Y JM epi , which appears in both repi and rrec. Since |Yepi| ≈ 1, its magnitude deviates only weakly from unity, and the oscillatory behav...

  16. [16]

    V B, we noted that the trajectories of constant phase lie close to the extrema of the reflectivity

    Complex-plane interpretation of equiphase trajectories and reflectivity extrema In Sec. V B, we noted that the trajectories of constant phase lie close to the extrema of the reflectivity. To clar- ify the origin of this behavior, we consider a simplified case in which surface reconstruction is neglected. This behavior can be understood geometrically in the c...

  17. [17]

    Since Ybulk, Y epi ≈ 1, the factor Y / (Y − 1) amplifies the difference in magnitude between the bulk and epilayer contributions, while preserving their similar phase

    49(rad), respectively, reflecting the similar atomic con- figurations of the GaN substrate and the InGaN epilayer. Since Ybulk, Y epi ≈ 1, the factor Y / (Y − 1) amplifies the difference in magnitude between the bulk and epilayer contributions, while preserving their similar phase. As a result, rbulk lies close to the extremal direction of the rtot circle, le...

  18. [18]

    This behavior is particularly evident near L = 1 and L = 3 in Fig

    Phase jumps and extrema switching near substrate minima We observe a rapid phase shift along the equiphase trajectories near minima of the total reflectivity R(L). This behavior is particularly evident near L = 1 and L = 3 in Fig. 10 for the (00L) CTR with fα = 0 and fα = 1, where the trajectories that initially correspond to reflectivity minima switch to m...

  19. [19]

    Specifically, for 2 < L < 2

    As established in Appendix B 2, the relative loca- tion of equiphase trajectories with respect to reflectivity extrema is determined not only by the phase difference between rbulk and zepi, but also by their relative magni- tudes. Specifically, for 2 < L < 2. 87, the phases of rbulk and zepi are nearly identical, while |rbulk| > |zepi|. In this regime, the b...

  20. [20]

    Fong and C

    D. Fong and C. Thompson, In situ syn- chrotron x-ray studies of ferroelectric thin film, Annu. Rev. Mater. Sci. 36, 431 (2006)

  21. [21]

    A. S. Disa, F. J. Walker, and C. H. Ahn, High- resolution crystal truncation rod scattering: Ap- plication to ultrathin layers and buried interfaces, Adv. Mat. Interf. 7, 1901772 (2020)

  22. [22]

    I. K. Robinson, Crystal truncation rods and surface roughness, Phys. Rev. B 33, 3830 (1986)

  23. [24]

    J. Wang, Q. Lin, E. Xu, B. Shen, and G. Ju, Determination of off-cut angle and direction on vicinal surfaces via x-ray crystal truncation rods, physica status solidi (b) 262, 2500177 (2025)

  24. [25]

    Thompson, C

    C. Thompson, C. M. Foster, J. A. Eastman, and G. B. Stephenson, Observation of the polarization of do- mains in ferroelectric thin films using x-ray interference, Appl. Phys. Lett. 71, 3516 (1997)

  25. [26]

    Fenter and C

    P. Fenter and C. Park, Termination interference along crystal truncation rods of layered crystals, J. Appl. Crystallogr. 37, 977 (2004)

  26. [27]

    T. A. Petach, A. Mehta, M. F. Toney, and D. Goldhaber- Gordon, Crystal truncation rods from miscut surfaces, Phys. Rev. B 95, 184104 (2017)

  27. [28]

    G. Ju, D. Xu, C. Thompson, M. J. Highland, J. A. Eastman, W. Walkosz, P. Zapol, and G. B. Stephen- son, In situ microbeam surface x-ray scattering re- veals alternating step kinetics during crystal growth, Nat. Commun. 12, 1721 (2021)

  28. [29]

    G. Ju, Y. Honda, M. Tabuchi, Y. Takeda, and H. Amano, In situ x-ray investigation of changing barrier growth temperatures on ingan single quantum wells in metal- organic vapor phase epitaxy, Journal of Applied Physics 115, http://dx.doi.org/10.1063/1.4867640 (2014)

  29. [30]

    I. K. Robinson, R. T. Tung, and R. Feidenhans, X- ray interference method for studying interface structures , Phys. Rev. B 38, 3632 (1988)

  30. [31]

    V. M. Kaganer, B. Jenichen, R. Shayduk, and W. Braun, Structure of Fe 3Si/GaAs (001) epitaxial films from x-ray crystal truncation rods, Phys. Rev. B 77, 125325 (2008)

  31. [32]

    J. M. LeBeau, R. Engel-Herbert, B. Jalan, J. Cagnon, P. Moetakef, S. Stemmer, and G. B. Stephenson, Stoichiometry optimization of ho- moepitaxial oxide thin films using x-ray diffraction, Appl. Phys. Lett. 95, 142905 (2009)

  32. [33]

    G. Ju, S. Fuchi, M. Tabuchi, and Y. Takeda, In situ x- ray measurements of movpe growth of In xGa1− xN single quantum wells, J. Cryst. Growth 370, 36 (2013)

  33. [34]

    G. Ju, D. Xu, C. Thompson, M. J. Highland, J. A. East- man, W. Walkosz, P. Zapol, and G. B. Stephenson, Crys- tal truncation rods from miscut surfaces with alternating terminations, Phys. Rev. B 103, 125402 (2021)

  34. [35]

    J. Y. Tsao, Ch 6: Surface Morphology, in Materials Fundamentals of Molecular Beam Epitaxy (Academic Press, Inc., San Diego, CA, 1993) pp. 201–257

  35. [36]

    Sarzynski, M

    M. Sarzynski, M. Leszczynski, M. Krysko, J. Z. Doma- gala, R. Czernecki, and T. Suski, Influence of GaN sub- strate off-cut on properties of InGaN and AlGaN layers, Cryst. Res. Technol. 47, 321 (2012)

  36. [37]

    G. A. Held and J. D. Brock, Calculation of x-ray-diffraction profiles from vicinal surfaces, Phys. Rev. B 51, 7262 (1995)

  37. [38]

    T. P. Trainor, P. J. Eng, and I. K. Robin- son, Calculation of crystal truncation rod structure factors for arbitrary rational surface terminations, J. Appl. Crystallogr. 35, 696 (2002)

  38. [39]

    Wollschlager and C

    J. Wollschlager and C. Tegenkamp, Diffraction from disordered vicinal surfaces with alternating terraces, Phys. Rev. B 75, 245439 (2007)

  39. [40]

    G. Ju, D. Xu, C. Thompson, M. J. Highland, J. A. Eastman, W. Walkosz, P. Zapol, and G. B. Stephenson, Burton-Cabrera-Frank theory for surface with alternat- ing step types, Phys. Rev. B 5, 054312 (2022)

  40. [41]

    Nagai, Structure of vapor-deposited Ga xIn1− xAs crys- tals, J

    H. Nagai, Structure of vapor-deposited Ga xIn1− xAs crys- tals, J. Appl. Phys. 45, 3789 (1974)

  41. [42]

    A. T. Macrander and K. E. Strege, X-ray double-crystal characterization of highly perfect InGaAs/InP grown by vapor-phase epitaxy, J. Appl. Phys. 59, 442 (1986)

  42. [43]

    D. A. Neumann, H. Zabel, and H. Morko¸ c, Terracing in strained-layer superlattices, J. Appl. Phys. 64, 3024 (1988)

  43. [44]

    Kry´ sko, J

    M. Kry´ sko, J. Z. Domagala, R. Czernecki, M. Leszczyn- ski, P. Perlin, T. Suski, S. Grzanka, G. Tar- gowski, I. Grzegory, M. Bo´ ckowski, and S. Porowski, Tilt of InGaN layers on miscut GaN substrates, Phys. Stat. Sol. Rap. Res. Lett. 4, 142 (2010)

  44. [45]

    Suzuki, T

    H. Suzuki, T. Sasaki, M. Takahasi, Y. Ohshita, N. Kojima, I. Kamiya, A. Fukuyama, T. Ikari, and M. Yamaguchi, Effect of substrate orientation on strain relaxation mechanisms of InGaAs layer grown on vicinal GaAs substrates measured by in situ x-ray diffraction, Japanese Journal of Applied Physics 56, 08MA06 (2017)

  45. [46]

    Krysko, J

    M. Krysko, J. Z. Domagala, R. Czernecki, and M. Leszczynski, Triclinic deformation of InGaN lay- ers grown on vicinal surface of GaN (00.1) substrates, J. of Appl. Phys. 114, 113512 (2013)

  46. [47]

    M. E. Liao and M. S. Goorsky, A refined model for epi- taxial tilt of elastically strained epilayers grown on misc ut substrates, J. Appl. Phys. 129, 025105 (2021)

  47. [48]

    Jiang, R.-V

    F. Jiang, R.-V. Wang, A. Munkholm, S. K. Streif- fer, G. B. Stephenson, P. H. Fuoss, K. Latifi, and C. Thompson, Indium adsorption on GaN un- der metal-organic chemical vapor deposition conditions, Appl. Phys. Lett. 89, 161915 (2006)

  48. [49]

    A. E. Romanov, T. J. Baker, S. Nakamura, J. S. Speck, and ERATO/JST UCSB Group, Strain-induced polarization in wurtzite III-nitride semipolar layers, Journal of Applied Physics 100, 023522 (2006)

  49. [50]

    A. F. Wright, Elastic properties of zinc- blende and wurtzite AlN, GaN, and InN, J. Appl. Phys. 82, 2833 (1997)

  50. [51]

    Richard, M

    M.-I. Richard, M. J. Highland, T. T. Fister, A. Munkholm, J. Mei, S. K. Streiffer, C. Thompson, P. H. Fuoss, and G. B. Stephenson, In situ synchrotron x-ray studies of strain and composition evolution dur- ing metal-organic chemical vapor deposition of InGaN, Appl. Phys. Lett. 96, 051911 (2010)

  51. [52]

    S. R. Andrews and R. A. Cowley, Scattering of x-rays from crystal surfaces, J. Phys. C 18, 6427 (1985) . 23

  52. [53]

    Munkholm, S

    A. Munkholm, S. Brennan, and E. Carr, A com- parison of surface roughness as measured by atomic force microscopy and x-ray scattering, J. Appl. Phys. 82, 2944 (1997)

  53. [54]

    Inatomi and Y

    Y. Inatomi and Y. Kangawa, Theoretical study of adatom stability on polar gan surfaces during mbe and movpe, Applied Surface Science 502, 144205 (2020)

  54. [55]

    R. G. van Silfhout, J. F. van der Veen, C. Norris, and J. E. Macdonald, Structure of the Ge(111)- c(2 × 8) surface as determined from scattered X-ray intensities along crystal truncation rods, Faraday Discuss. Chem. Soc. 89, 169 (1990)

  55. [56]

    Gu´ enard, G

    P. Gu´ enard, G. Renaud, A. Barbier, and M. Gautier- Soyer, Determination of the α -Ai2 O3 (0001) Surface Re- laxation and Termination by Measurements of Crystal Truncation Rods, Surf.Rev.Lett. 05, 321 (1998)

  56. [57]

    P. J. Eng, T. P. Trainor, G. E. Brown Jr., G. A. Way- chunas, M. Newville, S. R. Sutton, and M. L. Rivers, Structure of the Hydrated α -Al2 O3 (0001) Surface, Science 288, 1029 (2000)

  57. [58]

    M. H. Xie, S. M. Seutter, W. K. Zhu, L. X. Zheng, H. Wu, and S. Y. Tong, Anisotropic step-flow growth and is- land growth of GaN(0001) by molecular beam epitaxy, Phys. Rev. Lett. 82, 2749 (1999)

  58. [59]

    Waasmaier and A

    D. Waasmaier and A. Kirfel, New analytical scattering-factor functions for free atoms and ions, Acta Crystallogr. Sec. A 51, 416 (1995)

  59. [60]

    B. L. Henke, E. M. Gullikson, and J. C. Davis, X-ray interactions: photoabsorption, scattering, trans- mission, and reflection at E=50-30000 eV, Z=1-92, At. Data Nucl. Data Tables 54, 181 (1993) , updates at https://henke.lbl.gov/optical_constants/asf.html

  60. [61]

    T. K. Uˇ zdavinys, S. Marcinkeviˇ cius, M. Mensi, L. La- hourcade, J.-F. Carlin, D. Martin, R. Butte, and N. Grandjean, Impact of surface morphology on the properties of light emission in ingan epilayers, Applied Physics Express 11, 051004 (2018)