research papers\(\def\hfill{\hskip 5em}\def\hfil{\hskip 3em}\def\eqno#1{\hfil {#1}}\)

Journal logoFOUNDATIONS
ADVANCES
ISSN: 2053-2733

Theoretical study of the properties of X-ray diffraction moiré fringes. II. Illustration of angularly integrated moiré images

CROSSMARK_Color_square_no_text.svg

aSakai 5-13-2-A322, Musashino, Tokyo 180-0022, Japan
*Correspondence e-mail: j-yoshimura@voice.ocn.ne.jp

Edited by K. Tsuda, Tohoku University, Japan (Received 16 January 2019; accepted 4 April 2019; online 26 June 2019)

Using a theory of X-ray diffraction moiré fringes developed in a previous paper, labelled Part I [Yoshimura (2015[Yoshimura, J. (2015). Acta Cryst. A71, 368-381.]). Acta Cryst. A71, 368–381], the X-ray moiré images of a silicon bicrystal having a weak curvature strain and an interspacing gap, assumed to be integrated for an incident-wave angular width, are simulation-computed over a wide range of crystal thicknesses and incident-wave angular width, likely under practical experimental conditions. Along with the simulated moiré images, the graphs of characteristic quantities on the moiré images are presented for a full understanding of them. The treated moiré images are all of rotation moiré. Mo Kα1 radiation and the 220 reflection were assumed in the simulation. The results of this simulation show that fringe patterns, which are significantly modified from simple straight fringes of rotation moiré, appear in some ranges of crystal thicknesses and incident-wave angular width, due to a combined effect of Pendellösung oscillation and an added phase difference from the interspacing gap, under the presence of a curvature strain. The moiré fringes which slope to the perpendicular direction to the diffraction vector in spite of the assumed condition of rotation moiré, and fringe patterns where low-contrast bands are produced with a sharp bend of fringes arising along the bands are examples of the modified fringe pattern. This simulation study provides a wide theoretical survey of the type of bicrystal moiré image produced under a particular condition.

1. Introduction

In a previous paper (Yoshimura, 2015[Yoshimura, J. (2015). Acta Cryst. A71, 368-381.]), hereafter referred to as Paper I, a theory of X-ray diffraction moiré fringes from a bicrystal specimen was given, and the properties of the moiré image derived by this theory were explained by showing examples of plane-wave moiré images computed by the theory. However, the practically observed moiré images are integrated images for the angular spread of the incident X-ray wave. Therefore, the illustration of plane-wave images alone is incomplete for the study of moiré images, although they are important as the basics in moiré-fringe study. By showing the simulation of integrated images, the theory would become better understood, thus making it useful for practical problems. Therefore, in this paper, similar to the previous presentation of plane-wave moiré images, a series of integrated moiré images from a bicrystal specimen having a weak curvature strain and an interspacing air gap are simulation-computed and surveyed according to this moiré-fringe theory.

A newly occurring problem in treating integrated images is the effect of gap phase difference (hereafter gap phase) which is involved in the total interference phase of the wavefield [see equations (8)[link] and (9b)[link] shown later]. When we treat plane-wave images, the gap phase does not become a significant problem, except when the front component crystal of the bicrystal (hereafter front crystal) is strained (see Fig. I-12; equations and figures in Paper I are indicated by the Roman numeral I attached to the equation or figure numbers). However, in integrated images the gap phase plays an essential role in the formation of the interference fringe pattern. As explained in Paper I (see Fig. I-2 and the succeeding equations), the gap phase is produced when the diffracted waves from the front crystal are propagated through the interspacing gap, and is added to the total phase difference of the wavefield. This gap phase varies with the glancing angle of the incident wave [see equation (9b)[link] later], thereby causing a variation in the fringe position through a variation in the total phase difference of interference. By the integration of such small fringe-position variations, the resulting moiré-fringe pattern is significantly modified and the fringe contrast is decreased. To expand on the gap phase, its effect was first studied as gap interference fringes in X-ray topographs in theory and experiment (Authier et al., 1968[Authier, A., Milne, A. D. & Sauvage, M. (1968). Phys. Status Solidi, 26, 469-484.]; Hart & Milne, 1970[Hart, M. & Milne, A. D. (1970). Acta Cryst. A26, 223-229.]). The expression of the gap phase in plane-wave X-ray diffraction and a rocking-curve measurement from a gapped bicrystal of silicon have been reported by Yoshimura (1991[Yoshimura, J. (1991). Phys. Status Solidi A, 125, 429-440.]).

2. Theoretical

The integrated moiré images observed and discussed in this paper were computed as an integral of the plane-wave image intensity (diffracted-wave image) [see equation (I-20)] as follows:

[\eqalignno{{R_g}({{{\bf{r}}_{b'}}} )& = (1/ \Delta \Theta _{\rm inc})\textstyle\int\limits_{\Delta {\theta _i}}^{\Delta {\theta _f}} \big[ {I_{og'}}({\bf{r}}_{b'},u ) + I_{gg'}({\bf{r}}_{b'},u) &\cr &\quad+ {A_g}({\bf{r}}_{b'},u)\cos {\Psi _g}({\bf{r}}_{b'},u)&\cr &\quad + {B_g}({\bf{r}}_{b'},u)\sin {\Psi _g}({\bf{r}}_{b'},u)\big]\,{\rm d}\Delta \theta .&(1)}]

This integration may be understood to be a convolution integral of the plane-wave intensity function and a rectangular function of the peak height Iinc = 1 and angular width [\Delta {\Theta _{\rm inc}}] of the incident wave. The variable of integration [\Delta \theta] represents the deviation angle from the exact Bragg position when the X-ray wave is incident on the front crystal of the bicrystal. When we denote the upper and lower limit of integration by [\Delta {\theta _f}] and [\Delta {\theta _i}], respectively, the angular width of the integration, [\Delta {\Theta _{\rm inc}}], and the mid position of the integration width, [\Delta {\theta _{\rm m}}], are given, respectively, as follows:

[\Delta {\Theta _{\rm inc}} = \Delta {\theta _f} - \Delta {\theta _i}\eqno(2)]

[\Delta {\theta _{\rm m}} = (\Delta {\theta _i} + \Delta {\theta _f})/ 2. \eqno(3)]

For the meaning and expressions of [{I_{og'}}({{{\bf{r}}_{b'}},u} )], [{I_{gg'}}({{{\bf{r}}_{b'}},u} )], [{A_g}({{{\bf{r}}_{b'}},u} )] and [{B_g}({{{\bf{r}}_{b'}},u} )] in the integrand in equation (1)[link], see equations (I-22a,b) and (I-23a,b). [However, it must be stated here that the present expression of equation (I-22a) has an error: `[\cosh ({2K{\alpha _{21,r}}{t_1}} )]' in the expression must be corrected to `[\cos ({2K{\alpha _{21,r}}{t_1}} )]'. Although the error is present in the text, the computations in the paper were all made correctly using the correct expression.] In the calculation of the functions [{I_{og'}}({{{\bf{r}}_{b'}},u} )] etc. above, the following three deviation parameters are relevant [see equations (I-45), (I-46a,b)]:

[u=K\Delta\theta\sin2\theta_{\rm B}+\textstyle{1\over 2}K\chi_o(1-\gamma_g/\gamma_o)-2\pi (\Delta{\bf g}_1\cdot \hat{\bf K}_g),\eqno(4)]

[u_o=u-2\pi(\Delta{\bf g}\cdot \hat{\bf K}_g)\eqno(5a)]

[u_g=u-2\pi(\gamma_g/\gamma_o)(\Delta {\bf g}\cdot\hat{\bf K}_o).\eqno(5b)]

Here, u is the deviation parameter corresponding to the deviation angle [\Delta \theta] above; uo and ug are the deviation parameters when the transmitted and diffracted waves (O and G waves, respectively) emerging from the front crystal are incident on the rear component crystal (hereafter rear crystal) to excite the transmitted and diffracted waves once more; K is the wavenumber in vacuum, and [{\theta _{\rm{B}}}] is the Bragg angle; [{\chi _o}] is the 0th Fourier component of the dielectric susceptibility of the crystal; [{\gamma _o}] and [{\gamma _g}] are the direction cosines of the transmitted and diffracted beams, respectively, with respect to the normal to the entrance surface of the front crystal; [\Delta {{\bf{g}}_1}] is a local change in the reciprocal-lattice vector from the [{\bf{g}}] vector for the perfect region in the front crystal; [{\hat {\bf K}}_g] = [(\sin \theta _{\rm B},0,\cos \theta _{\rm B})] and [{\hat {\bf K}}_o] = [( - \sin \theta _{\rm{B}},0,\cos \theta _{\rm{B}})] denote unit vectors along the direction of the diffracted and transmitted waves, respectively. The argument [{{\bf{r}}_{b'}}] in the functions [{I_{og'}}({{{\bf{r}}_{b'}},u} )] etc. is a vector referring to a position on the exit surface of the rear crystal. [\Delta {\bf{g}}] in equations (5a[link]), (5b)[link] is defined as

[\Delta {\bf{g}} = \Delta {{\bf{g}}_2} - \Delta {{\bf{g}}_1}.\eqno(6)]

[\Delta {{\bf{g}}_2}] is a change in the reciprocal-lattice vector in the rear crystal, relative to the reciprocal-lattice vector [{\bf{g}}] in the front crystal. This [\Delta {\bf{g}}] is the reciprocal-vector difference between the front and rear crystals, which is relevant in the production of moiré-fringe patterns. It is expressed as

[\Delta {\bf{g}} = (1/ d)\left ( - {\Delta d}/d,\Delta \rho, \Delta \omega \right) \eqno(7)]

in the orthogonal coordinate system with the yz plane placed on the diffracting lattice plane (see Figs. I-2, I-3); [({\Delta d}/ d)] = [({\Delta d}/ d)_2 - ({\Delta d}/d)_1], [\Delta \rho] = [\Delta {\rho _2} - \Delta {\rho _1}], [\Delta \omega] = [\Delta {\omega _2} - \Delta {\omega _1}]; [({\Delta d}/ d )_i], [\Delta {\rho _i}] and [\Delta {\omega _i}] (i = 1, 2) denote, respectively, the lattice-spacing difference, and the rotation of the lattice plane about the z and y axes in the front crystal (i = 1) and rear crystal (i = 2), relative to an unstrained region in the front crystal; d is the lattice spacing.

When we assume symmetric Laue geometry where the diffracting lattice plane is perpendicular to the crystal surfaces, the term representing the phase difference of the interference in the integrand function in equation (1)[link] is expressed as

[{\Psi _g}({{{\bf{r}}_{b'}},u} ) = \left\{ ({2\pi } / d) [- ({\Delta d}/d)x + \Delta \rho y] - {\phi _{{\rm{gap}}}} \right\} \eqno(8)]

[see equations (I-34)–(I-37) for details]. Here, [ -({\Delta d}/ d)/ d] and [{\Delta \rho }/ d] in the first term on the right side are the x and y components, respectively, of [\Delta {\bf{g}}], as shown in equation (7)[link]; the moiré-fringe pattern is intrinsically drawn by these two quantities. The second term represents the gap phase as mentioned in Section 1[link], and is given by [see equations (I-11), (I-21) and (I-34)]:

[{\phi _{{\rm{gap}}}} = {t_{{\rm{gap}}}}\left[u - \textstyle{1\over 2}K{\chi _o}(1 - {\gamma _g}/ {\gamma _o}) \right]/{\gamma _g}\eqno(9a)]

[=t_{\rm gap}\left[K\Delta \theta\sin2\theta_{\rm B}-2\pi(\Delta{\bf g}_1\cdot \hat{\bf K}_g)\right]/\gamma_g.\eqno(9b)]

Here, tgap is the width of the gap. If the gap is filled with some high-density material such as silicon oxide, instead of air, an added phase owing to the passage through it should be taken into account. Then, the gap phase above is modified to

[\eqalignno{\phi_{\rm gap}&=t_{\rm gap}\left[K\Delta \theta\sin2\theta_{\rm B}-2\pi(\Delta{\bf g}_1\cdot \hat{\bf K}_g)\right]/\gamma_g &\cr &\quad \times (1 + \textstyle{1\over 2}{\chi _{s,o}}). &(10)}]

Here, [{\chi _{s,o}}] is the Fourier component of the dielectric susceptibility of the filled material. However, while [\Delta \theta] remains a minute quantity, this correction is an even more minute quantity. Therefore, the correction is considered to be practically unnecessary.

The preparations for computing moiré images are ready now. However, to discuss the characteristics of the interference patterns of the computed moiré images, we also need the wavefield expression in which the two phase terms in equation (1)[link] are arranged to form one term. For this reason, we prepare the following expression by expanding and rearranging equations (1)[link] and (8)[link]:

[\eqalignno{R_g({\bf{r}}_{b'} )& = R_{og'}({\bf{r}}_{b'} ) + R_{gg'}({\bf{r}}_{b'})&\cr &\quad + \left\{[RAA_g({\bf{r}}_{b'})]^2 + [RBB_g({\bf{r}}_{B'})]^2 \right\}^{1/2}&\cr &\quad\times \cos \left\{ ({2\pi }/ d)\left [ - ({\Delta d}/ d)x + \Delta \rho y \right] - \Phi _{{\rm P},g}({\bf{r}}_{b'}) \right\} &(11)}]

with

[R_{og'}({\bf{r}}_{b'}) = (1/\Delta \Theta _{\rm inc})\textstyle\int\limits_{\Delta {\theta _i}}^{\Delta {\theta _f}} I_{og'}({\bf{r}}_{b'},u )\,{\rm d}\Delta \theta \eqno(12a)]

[R_{gg'}({\bf{r}}_{b'}) = (1/\Delta \Theta _{\rm inc})\textstyle\int\limits_{\Delta {\theta _i}}^{\Delta {\theta _f}} I_{gg'}({\bf{r}}_{b'},u )\,{\rm d}\Delta \theta \eqno(12b)]

[\Phi _{{\rm P},g} ({\bf{r}}_{b'}) = {\rm{atan}}\left[RBB_g ({\bf{r}}_{b'})/ RAA_g ({\bf{r}}_{b'}) \right]\eqno(13)]

[\eqalignno{RAA_g({{\bf{r}}_{b'}})& = (1/ \Delta \Theta _{\rm inc}) \textstyle\int\limits_{\Delta {\theta _i}}^{\Delta {\theta _f}} \Big [{A_g}({{\bf{r}}_{b'}},u)\cos (\phi _{\rm gap})&\cr &\quad - B_g({{\bf{r}}_{b'}},u)\sin (\phi _{\rm gap}) \Big] \,{\rm d}\Delta \theta & (14a)}]

[\eqalignno{RBB_g({{\bf{r}}_{b'}})& = (1/ \Delta \Theta _{\rm inc}) \textstyle\int\limits_{\Delta {\theta _i}}^{\Delta {\theta _f}} \Big [{A_g}({{\bf{r}}_{b'}},u)\sin (\phi _{\rm gap})&\cr &\quad + B_g({{\bf{r}}_{b'}},u)\cos (\phi _{\rm gap}) \Big] \,{\rm d}\Delta \theta. &(14b)}]

As shown in the expression in equation (11)[link], [{\Phi _{{\rm{P}},g}}({{{\bf{r}}_{b'}}} )] in equation (13)[link] represents the contribution of Pendellösung oscillation (hereafter PL oscillation) to the total interference phase of the wavefield. In this sense it may be called the Pendellösung-connected term (hereafter the PL phase), although it is also connected with the gap phase as seen from equations (13)[link] and (14a)[link], (14b)[link]. In combination with the effect of the gap phase, the phase [{\Phi _{{\rm{P}},g}}({{{\bf{r}}_{b'}}} )] becomes a source of oscillations that are more complicated than the simple PL oscillation. From the expression in equation (11)[link], the fringe contrast in moiré images is given as follows:

[\eqalignno{{V_g}({\bf{r}}_{b'} )& = \left[R_{g,\max }({\bf{r}}_{b'}) - R_{g,\min }({\bf{r}}_{b'})\right]/\left[R_{g,\max }({\bf{r}}_{b'}) + R_{g,\min }({\bf{r}}_{b'})\right]&\cr &=\left\{ [RAA_g({\bf{r}}_{b'})]^2 + [RBB_g({\bf{r}}_{b'})]^2\right\}^{1/2}&\cr &\quad/\left[R_{og'}({\bf{r}}_{b'}) + R_{gg'}({\bf{r}}_{b'})\right]. &(15)}]

As to the numerical conditions for computation, the Si 220 reflection with Mo Kα1 radiation (0.70926 Å) was assumed throughout this paper; [{\theta _{\rm{B}}}] = 10.64°. The linear absorption coefficient is μ = 1.462 mm−1. [All these values are taken from Pinsker (1978[Pinsker, Z. G. (1978). Dynamical Scattering of X-rays in Crystals. Berlin: Springer-Verlag.]).] As already mentioned above, the diffracting (220) lattice plane was assumed to lie perpendicular to the crystal surfaces in the symmetric Laue geometry ([{\gamma _o}] = [{\gamma _g}]). The crystal surfaces in the bicrystal specimen are all parallel to one another (see Fig. I-2). Except for the case of Fig. 7, the differences in the spacing and orientation of the lattice plane from the reference region were assumed as follows: the front crystal is strain free, namely [\Delta {{\bf{g}}_1}] = 0 [[({\Delta d}/ d)_1] = [\Delta {\rho _1}] = [\Delta {\omega _1}] = 0] and s1 = 0 [si (i = 1, 2) denotes the strength of the curvature deformation of the crystals around the y axis]; the rear crystal, on the other hand, has a non-zero reciprocal-vector difference [\Delta {{\bf{g}}_2}] [\ne] 0, i.e. [({\Delta d}/ d)_2] = 0, but [\Delta {\rho _2}] = d/0.045 (rad); in addition, the crystal has a curvature deformation of s2 = 0.045′′ mm−1 (radius of curvature r = 4600 m) and thereby a rotation of the lattice plane around the y axis [\Delta\omega_2] = s2·(x-xo) is induced over the whole crystal. Here, xo denotes a position on the crystal surface corresponding to the centre of curvature, and was assumed to be xo = 9.0 (mm). The deviation parameters in equations (4)[link] and (5a)[link], (5b)[link] become

[u = K\Delta \theta \sin 2 \theta _{\rm{B}} \eqno(16)]

[{u_o} = {u_g} = u - (2\pi/d)\Delta {\omega _2}\cos \theta _{\rm{B}}\eqno(17)]

based on the assumptions above. All these conditions are the same as those in the computation of the plane-wave moiré images in Paper I. A rotation-moiré pattern should appear in the diffracted images, similar to the case of the previous plane-wave images. Hereafter, the thicknesses of the component crystals, the gap width, and the x position in the crystals and moiré images are given in units of mm (the y position is not referred to in this paper). The computation of moiré images was made with Visual Basic.NET version 2003, and the intensity and contrast in the output images were adjusted in image processing so as to be best suited for observation. Computation of the graphs of the characteristic quantities [R_{og'}({\bf{r}}_{b'} ) + R_{gg'}({\bf{r}}_{b'} )], [{V_g}({\bf{r}}_{b'})] etc. was made using SigmaPlot version 11.0.

3. Results and discussion

The commonly assumed numerical conditions for the computation were as stated above. In the following we observe and discuss the moiré images obtained by changing the thicknesses of the component crystals t1 and t2 of the front and rear crystals, respectively, and the incident-wave angular width [\Delta {\Theta _{\rm inc}}]. The gap width tgap was also set at different values in the computation.

3.1. Moiré images of bicrystals of t1 = t2 = 0.8 (μt1,2 = 1.170)

3.1.1. Images with tgap = 0

First, we show some integrated moiré images when tgap = 0 in Figs. 1[link](a)–1[link](c). These are representative of images obtained when [\Delta {\Theta _{\rm inc}}] is increased gradually, starting from a minute value. To compare unambiguously with previous plane-wave images, the mid deviation angle in the integration was set to be [\Delta\theta_{\rm m}=0.32''] for which angle the plane-wave images Fig. I-8 etc. were computed; it is slightly off the peak position of the diffracted-intensity curve. In accordance with the assumption of [({\Delta d}/ d) = 0] and [\Delta \rho = d/ 0.045], moiré fringes with a fringe spacing Λ = 0.45 mm in the y direction are produced. Needless to say, these images are topographs of the diffracted wave (G image). The images are shown in the way that they are viewed from the emerging-beam side, and the x- and y-coordinate axes are taken as shown in Fig. 1[link](c). The image contrast is reproduced in such a way that white contrast indicates higher intensity, which is the opposite to convention. The diffraction vector [{\bf{g}}], which is not shown in the images, is directed from the left to the right along the x axis (see e.g. Figs. I-1 and I-3).

[Figure 1]
Figure 1
Computer-simulated integrated moiré images of X-ray diffraction for a silicon bicrystal specimen with component-crystal thicknesses of t1 = t2 = 0.8 mm and gap width of tgap = 0. The 220 reflection with Mo Kα1 radiation was assumed. Angular widths of the incident X-rays were set to be [\Delta\Theta_{\rm inc}] = 0.04′′, 0.50′′ and 10.0′′ for images (a), (b) and (c), respectively. The mid deviation angle was set to be [\Delta\theta_{\rm m}] = 0.32′′. The graduations at the bottom of the image (c) are given in mm. For further information see text.

Fig. 1[link](a) shows the image computed with angular width [\Delta\Theta_{\rm inc}=0.04''], which is very close to the plane-wave condition. In fact, this moiré image is very similar to the plane-wave image in Fig. I-8(a). When [\Delta {\Theta _{\rm inc}}] was increased, the integrated image remained nearly the same as the plane-wave image up to about [\Delta {\Theta _{\rm inc}}=0.12''], and a somewhat clear difference between the integrated and plane-wave images began to be found from about [\Delta {\Theta _{\rm inc}}=0.16'']. When [\Delta {\Theta _{\rm inc}}] was increased further, an oscillation in the fringe lines and a change in the image intensity distribution with [\Delta {\Theta _{\rm inc}}] continued to be observed up to about [\Delta {\Theta _{\rm inc}}\simeq 2''], although the oscillation was of very small amplitude and the intensity distribution change was very gentle. Fig. 1[link](b) shows an example of a moiré image in such a middle region of [\Delta {\Theta _{\rm inc}}] before reaching [\Delta {\Theta _{\rm inc}}\simeq 2'']. The fringe lines still show a slight oscillation in this image [see Fig. 2[link](b) later]. With the angular width of [\Delta {\Theta _{\rm inc}}\ge 2''], the fringe pattern had almost settled to the one shown in Fig. 1[link](c). Even when [\Delta {\Theta _{\rm inc}}] was increased further, the fringe pattern did not change significantly from the pattern in Fig. 1[link](c), although the image intensity and fringe contrast changed to some extent. It was confirmed by the present computation that such almost-settled fringe patterns continue to be observed up to [\Delta {\Theta _{\rm inc}}\simeq 60''], which may be referred to as the condition in Lang topography.

[Figure 2]
Figure 2
Curves of characteristic quantities of the moiré images in Fig. 1[link], showing local variations of the mean image intensity [[{R_{og'}}(x) + {R_{gg'}}(x )]], fringe contrast Vg(x ), amplitude intensities RAAg(x ) and RBBg(x ), and PL phase [{\Phi _{{\rm{P, }}g}}(x )]. Graphs (a), (b) and (c) are related to the moiré images (a), (b) and (c), respectively, in Fig. 1[link]. For further information see text.

Fig. 2[link] shows the curves of the mean image intensity [[R_{o,g'}(x) + R_{g,g'}(x )]] unrelated to the moiré interference, fringe contrast Vg(x ), amplitude intensities RAAg(x ) and RBBg(x ) in the interference terms [in equation (1)[link]], and of the PL phase [\Phi _{{\rm{P}},g}(x )] for the moiré images in Figs. 1[link](a)–1[link](c). These characteristic quantities were computed separately from the computation of moiré images. (From here on, the position variable in these functions is denoted by x instead of [{\bf{r}}_{b'}] as in Section 2[link].) The results of the computation of the moiré images and characteristic curves agree well with each other. We can understand the characteristics of the moiré images clearly and in detail from these characteristic curves. Firstly, it can be pointed out that the phase curve [\Phi _{{\rm{P}},g}(x )] in Fig. 2[link](a) is very similar to the phase curve [[\phi _{{\rm{p}},g}({\bf{r}},u)]] of the plane-wave image shown in Fig. I-9(b); it is seen that the wave forms of oscillation in these phase curves are both slightly asymmetric and lean to the left side. This similarity between the two phase curves confirms the correctness of the similarity mentioned above for the integrated and plane-wave moiré images [Fig. 1[link](a) and Fig. I-8(a)]. The fringe-line equation obtained by setting the argument of the cosine function in equation (11)[link] to be equal to [2\pi N] is given in this case ([\Delta d/ d = 0]) by

[(2\pi \Delta \rho/ d)y = \Phi _{{\rm{P}},g}(x ) + 2\pi N. \eqno(18)]

As seen from this equation, the shape of the fringe lines is governed by the phase curve [{\Phi _{{\rm{P}},g}}(x )].

It can be seen in Figs. 2[link](a)–2[link](c) that the oscillation of moiré-fringe lines is damped rapidly with an increase in the angular width [\Delta {\Theta _{\rm inc}}], and the shape of the fringe lines approaches that of a horizontal straight line. This damping is caused by the convergence of RBBg(x ) to the zero value, which is set as the numerator in the expression of [{\Phi _{{\rm{P}},g}}(x )] [equation (13)[link]], while the denominator RAAg(x ) remains at a finite value. Although the curve of [{\Phi _{{\rm{P}},g}}(x )] in the lower graph in Fig. 2[link](c) is drawn on an enlarged scale in the vertical direction, to show that an infinitesimal oscillation still exists, it is practically a horizontal straight line in the scale of Fig. 2[link](a) or 2[link](b). It must be added that the property [| {RBB_g(x )} | \to 0] with increase in [\Delta {\Theta _{\rm inc}}] is true only when [{t_{{\rm{gap}}}} \simeq 0]; in the case of [{t_{{\rm{gap}}}} \ne 0], RBBg(x ) does not show such convergence. As seen from equations (14a)[link], (14b)[link], when [t_{\rm gap} \simeq 0],

[RAA_g(x) = (1/\Delta \Theta _{\rm inc})\textstyle\int\limits_{\Delta \theta _i}^{\Delta \theta _f} A_g(x,u)\,{\rm d}\Delta \theta , \eqno(19a)]

[RBB_g(x) = (1/\Delta \Theta _{\rm inc})\textstyle\int\limits_{\Delta \theta _i}^{\Delta \theta _f} B_g ({x,u} )\,{\rm d}\Delta \theta .\eqno(19b)]

In addition to the fringe lines, the fringe contrast Vg(x ) also shows an oscillatory variation, and its minimum value becomes zero in some cases. Such oscillatory variations in Vg(x ) can generally be understood to be governed by the oscillations in RAAg(x ) and RBBg(x ). As seen from equation (15)[link], the variation in Vg(x ) is approximately estimated by

[{V_g}(x )\infty \left \{[RAA_g(x)]^2 + [RBB_g(x)]^2\right\}^{1/2}. ]

A shallow drop in the curves of the mean image intensity [[R_{o,g'}(x ) + R_{g,g'}(x )]] around [x \simeq 9], seen in the upper graphs in Figs. 2[link](b), 2[link](c), is considered to show a secondary extinction effect in the crystal diffraction. In the images [Figs. 1[link](b), 1[link](c)] this reduction in intensity appears as a black, diffuse band image. The rear crystal becomes [{\omega}_2(x ) = 0] in this position [x \simeq 9], and the front and rear crystals become locally parallel to each other. Similar reductions in intensity are also observed with the crystal thicknesses of t1 = t2 = 1.5 and [{t_1} \simeq {t_2} \simeq 2.0] (see Figs. 5 and 8, respectively), but disappear with [{t_1},{t_2} \ge 2.5] (see Fig. 9). When both the front and rear crystals were unbent, i.e. s1 = s2 = 0, and tgap = 0, but [\Delta d/d = 0] and [\Delta \rho = d/0.045] as in the other cases described above, the moiré images were properly of horizontal straight fringes (hereafter HS fringes) with any value of [\Delta {\Theta _{\rm inc}}], since all the quantities concerned become independent of the x coordinate. Then, while the mean image intensity [[R_{o,g'}(x ) + R_{g,g'}(x )]] is considerably decreased compared with the images in Fig. 1[link], fringe contrast Vg(x ) is greatly increased (see Figs. 11 and 12 shown later).

3.1.2. Images with tgap = 0.05

When the supposed bicrystal specimen comes to have a gap of finite width, some change from images with tgap = 0 is expected to appear. However, up to about tgap = 0.01, the obtained moiré images were almost unchanged from the ones with tgap = 0, by comparison of the images (not the characteristic curves). Even when the gap width became tgap = 0.02, there was almost no change from the image with tgap = 0, although a small difference came to be found with a large value of [\Delta {\Theta _{\rm inc}}]. Though not clearly recognized in the images, the oscillation in the curve of [{\Phi _{{\rm{P}},g}}(x )] appeared to increase compared with the level when tgap = 0, and the fringes sloped slightly upwards to the right. The previous moiré images by Brádler & Lang (1968[Brádler, J. & Lang, A. R. (1968). Acta Cryst. A24, 246-247.]) and by Lang (1968[Lang, A. R. (1968). Nature, 220, 652-657.]) are considered to be taken under such a condition, but with no curvature strain in the specimens.

When the gap width was increased further, the difference from images with tgap = 0 gradually became clear. Moiré images with tgap = 0.05 and graphs of the characteristic curves concerned are shown in Figs. 3[link] and 4[link], respectively. Figs. 3[link](a), 3[link](b) are to be compared with Figs. 1[link](b), 1[link](c), respectively. When the angular width [\Delta {\Theta _{\rm inc}}] remains small ([\Delta {\Theta _{\rm inc}}\le 0.2'']), images with tgap = 0.05 are nearly the same as those with tgap = 0. However, when the angular width was increased further, the difference in images began to be noticeable. Though not clear in the images, the graphs in Fig. 4[link] show that the fringe contrast Vg(x ) is lowered considerably with tgap = 0.05, compared with the case with tgap = 0. When the angular width was increased up to [\Delta {\Theta _{\rm inc}}\simeq 4''], the fringe pattern had almost settled to the one shown in Fig. 3[link](c), and remained almost unchanged with further increase in [\Delta {\Theta _{\rm inc}}].

[Figure 3]
Figure 3
Moiré images simulation-computed under similar conditions to those for Fig. 1[link], except for the gap width of tgap = 0.05. t1 = t2 = 0.8. (a) Moiré image with [\Delta\Theta_{\rm inc}] = 0.50′′, (b) moiré image with [\Delta\Theta_{\rm inc}] = 10.0′′. For further information see text.
[Figure 4]
Figure 4
Characteristic curves of the moiré images in Fig. 3[link]. Graphs (a) and (b) are related to the moiré images (a) and (b), respectively, in Fig. 3[link].

3.2. Moiré images of bicrystals with t1 = t2 = 0.8 (μt1,2 = 1.170) and tgap = 0.24

3.2.1. General survey of resulting images

Figs. 5[link](a)–5[link](d) show some representative moiré images when the crystal thicknesses were t1 = t2 = 0.8, but the gap width was set to be tgap = 0.24. Here, computation was performed with the mid deviation angle [\Delta\theta_{\rm m}=-0.12''], which was changed from [\Delta\theta_{\rm m}=0.32''] for the images in Figs. 1[link] and 3[link]. In Paper I, the plane-wave image in Fig. I-10(a) was computed under this condition, namely [\Delta\theta_{\rm m}=-0.12'']. When the angular width [\Delta {\Theta _{\rm inc}}] was increased gradually, the integrated images did not show any obvious difference from the plane-wave images up to about [\Delta {\Theta _{\rm inc}}\simeq 0.08'']; however, sharp differences became visible from about [\Delta {\Theta _{\rm inc}}\simeq 0.12'']. Fig. 5[link](a) shows an example of such an integrated image with a relatively small [\Delta {\Theta _{\rm inc}}]. For the increased gap width of tgap = 0.24, the fringe patterns became considerably modified from the simple HS fringes. Representative examples of such moiré images when tgap = 0.24 and [\Delta {\Theta _{\rm inc}}] is of middle magnitudes are shown in Figs. 5[link](b) and 5[link](c); they will be explained in detail later. In the case of tgap = 0.24 also, a large change in the fringe pattern was not observed with angular width larger than [\Delta {\Theta _{\rm inc}}\simeq4''], and the fringe pattern had almost settled to the one shown in Fig. 5[link](d). This image is roughly similar to the one in Fig. 3[link](b), but the slope of the fringes becomes larger than in Fig. 3[link](b) and vertical streaks or bands appear more strongly. The slope angle of the fringes was 3.5° in the image shown in Fig. 5[link](d), which was measured directly on the image, while it was 1.3° in the image shown in Fig. 3[link](b). While the slope angle appears to be related to the magnitude of tgap, it is also thought to be related to the curvature strain in the crystal. To confirm this thought, computations of moiré images were made to confirm that the slope angle of the fringes is zero with s2 = 0, and increases with the value of s2.

[Figure 5]
Figure 5
Moiré images computed under similar conditions to those in Fig. 1[link], except for the gap width of tgap = 0.24 and the mid deviation angle of [\Delta\theta_{\rm m}] = −0.12′′. t1 = t2 = 0.8. The incident-wave angular widths for the computed moiré images were assumed to be: for the image (a) [\Delta\Theta_{\rm inc}] = 0.40′′, (b) [\Delta\Theta_{\rm inc}] = 1.08′′, (c) [\Delta\Theta_{\rm inc}] = 1.20′′ and (d) [\Delta\Theta_{\rm inc}] = 10.0′′. For further information see text.
3.2.2. Effect of gap phase on the formation of fringe pattern

In the following we discuss the images in Figs. 5[link](b) and 5[link](c) in detail. In the moiré image in Fig. 5[link](b) fringes slope by about 8° from the x axis, in the direct measurement on the image. Such a fringe pattern appears to be extraordinary for that of rotation-moiré fringes with [\Delta d/d = 0]. Although it is difficult to give a satisfactory explanation for this strange fringe pattern at present, the following is a partial explanation based on an incomplete study. As mentioned in Section 1[link], when the deviation angle [\Delta \theta] changes, the position of the fringe lines also changes due to the change in the gap phase [{\phi _{{\rm{gap}}}}] [see equations (8)[link] and (9b)[link]]. For a change of [\Delta \theta] in the deviation angle, the gap phase changes by

[{\phi _{{\rm{gap}}}} = {t_{{\rm{gap}}}}K\Delta \theta \sin 2{\theta _{\rm{B}}} = 38.1\Delta \theta, ]

with tgap = 0.24 (here, [\Delta \theta] is given in arcsec). The fringe positions move upwards or downwards due to this phase change. Therefore, the effective width of the Nth fringe is expanded by such fringe-position movements, and it becomes [2\pi] due to the [\Delta \theta] change in the overall width of [\delta(\Delta\theta)=0.165''], and [4\pi] for [\delta(\Delta\theta)=0.33''] with its expansion forefronts reaching the central positions of the neighbouring ([N \pm 1])-th fringes. With the incident-wave angular width of [\Delta {\Theta _{\rm inc}}=1.08''] for the image under consideration, the range of one fringe extends to the distance of 3.3 times the fringe spacing to each side of its central position. Each fringe is thus expanded, and overlaps with neighbouring fringes. The patterns of obliquely extending fringes (hereafter OE fringes or OE fringe pattern) and of broadly horizontal fringes (BH fringes or BH fringe pattern), which are shown, respectively, in Figs. 5[link](b) and 5[link](c), are considered to be produced in connection with such a complicated fringe arrangement.

While simulating the formation process of the OE and BH fringe patterns, it was seen that the resulting fringe-line configuration in the integrated image is formed by connecting the regions where plane-wave fringe lines are most densely overlapped. Then, due to a slight difference in the overlap of plane-wave fringe lines and/or a slight change in the shape of the fringe lines, the resulting fringe pattern seems to be divided into OE and BH fringe patterns. While some fringes (in integrated images) continue to lie in their initial y positions, though making an oscillation with the x position, the other fringes switch their initial vertical positions so that they are connected to a neighbouring fringe at a site that is one spacing higher or lower, in the way of the continuation along the x axis. The fringes repeat such switching of the vertical position, thereby extending obliquely to a higher or lower position. Such a difference in the continuation behaviour of fringes seems to cause the difference of the OE and BH fringe patterns. However, further investigation of these fringe patterns is not an easy task, because it requires precise simulation of complicated fringe arrangements.

3.2.3. OE and BH fringe patterns viewed in the related ΦP,g(x) phase curves

The formation of the OE and BH fringe patterns discussed above can be understood to a fair extent from the observation of [{\Phi _{{\rm{P}},g}}(x )] phase curves shown in Figs. 6[link](a), 6[link](b). In each of these figures, the `as-output [{\Phi _{{\rm{P}},g}}(x )]' phase curve is additionally given in the middle-column graph to provide a good understanding of the `corrected [{\Phi _{{\rm{P}},g}}(x )]' phase curve shown in the lowest graph. The as-output [{\Phi _{{\rm{P}},g}}(x)] curves are phase curves computed using equation (13)[link]. They have discontinuous jumps of [\pi] or [ - \pi], wherever RAAg(x) in the denominator in equation (13)[link] crosses the zero line. These discontinuities were all corrected so that the as-output [{\Phi _{{\rm{P}},g}}(x )] curves are changed to continuous curves of the corrected [{\Phi _{{\rm{P}},g}}(x )] phase, as described in Fig. I-7. The corrected [{\Phi _{{\rm{P}},g}}(x )] curve is the true phase curve, and it gives the shape of the fringe lines in the moiré image, as shown in equation (18)[link]. The moiré images are computed using equation (1)[link], so that they conform with the corrected [{\Phi _{{\rm{P}},g}}(x )] phase curve from the beginning.

[Figure 6]
Figure 6
Characteristic curves of the moiré images in Fig. 5[link]. Graph (a) is related to the moiré image (b) in Fig. 5[link], and graph (b) is related to the moiré image (c) in Fig. 5[link]. For further information see text.

It is seen in Fig. 6[link](a) for the OE fringes that the discontinuous phase change is [ - \pi] (seen from left to right) at all phase jumps in the as-output [{\Phi _{{\rm{P}},g}}(x )] curve, and the curve position is raised by [ + \pi] at every correction site. In more detail, there occur two types of phase jumps, A and B. While in the type-A jump the curves on both sides of the jumping site are connected to each other smoothly after the correction, in the type-B jump an upward bend remains in the corrected curve, to enhance the rise of the curve position. Such a rise in the vertical position of the phase curve is presumed to agree with the fringe-position rise of OE fringes described qualitatively in Section 3.2.2[link]. In contrast to the case shown in Fig. 6[link](a), discontinuous phase changes of [ + \pi] and [ - \pi] occur alternately in the as-output [{\Phi _{{\rm{P}},g}}(x )] curve shown in Fig. 6[link](b) for the BH fringes. Although an upward slope of the phase curve is also seen in this case [see the corrected [{\Phi _{{\rm{P}},g}}(x )] curve in the lowest graph], gradual rises in the curve position are cancelled by sharp downward bends at the type-B correction sites, to keep the phase curve broadly horizontal. If we want to inquire further into the behaviours of [{\Phi _{{\rm{P}},g}}(x)] phase curves in these two cases, the study of the curves of RAAg(x ) and RBBg(x ) is necessary. However, it is omitted. So far, we have described the OE and BH fringe patterns of moiré images as representatives in the middle-magnitude domain of [\Delta {\Theta _{\rm inc}}]. Roughly speaking, these two types of fringe patterns appear alternately with an increase in [\Delta {\Theta _{\rm inc}}] in the range from about 0.4′′ to about 4′′. The type of fringe patterns produced changes with [\Delta {\Theta _{\rm inc}}]. The frequency of occurrence of the two types seems to depend on the crystal thicknesses t1 and t2. Images where the two types of fringe patterns co-appear in such a way that one type of fringe pattern is connected smoothly to another type also appeared fairly frequently.

3.2.4. Low-contrast bands in the BH fringe pattern

Further comments are made below on the BH fringe pattern in Figs. 5[link](c) and 6[link](b). At the positions of the type-B phase jumps in Fig. 6[link](b), [x \cong 2.2,4.7,8.3], the fringe contrast Vg(x) locally falls to produce a streak or band pattern, accompanying a sharp bend in the fringe lines. Although similar low-contrast bands can be seen in Figs. 5[link](b) and 6[link](a) for the OE fringe pattern, the fringe bends are not so sharp in that case. As seen in the graph in Fig. 6[link](b), sharp fringe bends occur where RAAg(x) = 0 and RBBg(x ) approaches its maximum. Then, [| {RB{B_g}(x )} | \simeq 0], and the fringe contrast Vg(x ) is minimized according to equation (15)[link]. The graphs of the as-output and corrected [{\Phi _{{\rm{P}},g}}(x)] curves are similar to those of the phase curves in Fig. I-9(a) for the plane-wave image with which `abrupt fringe jump' is explained, although the wave form in the phase-curve oscillation leans to the opposite side to the lean of the wave form in the present phase curve in Fig. 6[link](b). The curves of RAAg(x ), RBBg(x ) and [{\Phi _{{\rm{P}},g}}(x )] in Fig. 6[link](b) correspond to the curves of [{A_g}({{\bf{r}},u})], [{B_g}({{\bf{r}},u} )] and [{\phi _{{\rm{P}},g}}({\bf{r}} )] in Figs. I-9(a), respectively. On the basis of this similarity, it is evident that the `low-contrast bands' in this paper are of the same origin as the `vertical bands of abrupt fringe jump' for the plane-wave image in Paper I.

The band patterns mentioned above are produced as a result of oscillation with intensities RAAg(x ) and RBBg(x), which were originally caused by the PL oscillation in the crystals. The contour or fringe patterns caused by the PL oscillation in the diffracted image of a curved crystal are called bend-extinction or equal-inclination fringes (Hirsch et al., 1965[Hirsch, P. B., Howie, A., Nicholson, R. B., Pashley, D. W. & Whelan, M. J. (1965). Electron Microscopy of Thin Crystals. London: Butterworths.]; Sugii et al., 1971[Sugii, K., Matsushita, T. & Kohra, K. (1971). J. Phys. Soc. Jpn, 31, 1295-1296.]). It is certain that the band pattern in the present moiré images is of the same kind as this PL oscillation-connected fringe pattern. From the computations made so far, it is confirmed that the number of bands produced is increased or decreased depending on the curvature strength [| {{s_1}} |] or [| {{s_2}} |]. However, in the study for the case of tgap = 0, which is described in Section 1[link], it is observed that band patterns disappear in moiré images with [\Delta\Theta_{\rm inc}\ge 1''], certainly due to smoothing of the intensity oscillation. Nevertheless, the band pattern in the moiré image in Fig. 5[link](c) is clearly observed even with an angular width of [\Delta\Theta_{\rm inc}= 1.2'']. This disagreement is considered to indicate that the band pattern in the integrated moiré images with [{t_{{\rm{gap}}}} \ne 0] is produced not only by the PL oscillation, but also by a collaboration of the PL oscillation and gap phase effect. In this respect the band pattern under discussion is not entirely similar to the intrinsic equal-inclination fringes. It is tentatively referred to as the low-contrast band (LC band).

3.2.5. Other remarks

One more comment that should be made on Figs. 5[link](c) and 6[link](b) is regarding the direction of lean of the oscillation wave form in the fringe lines or the [{\Phi _{{\rm{P}},g}}(x)] phase curve. In the phase curve in the plane-wave moiré image, the oscillation wave form always leans to the left side, as seen in Figs. I-8(a), (b) and Figs. I-9(a), (b). In the integrated images the oscillation wave form in question also leans to the left side as long as the angular width [\Delta {\Theta _{\rm inc}}] remains very small. An example of such a phase-curve oscillation wave form is seen in Fig. 2[link](a), in the case of tgap = 0. This characteristic of the leaning direction also holds in the case of tgap [\ne] 0. However, as the value of [\Delta {\Theta _{\rm inc}}] increases, the oscillation wave form leans to the right side, as seen in Fig. 4[link](a) and Fig. 6[link](b), via a neutral symmetric form in a narrow range of [\Delta {\Theta _{\rm inc}}]. As studied by the computation of moiré images and characteristic phase curves, the oscillation wave form in the integrated images when [{s_2} \,\gt\, 0] seems to generally lean to the right side except when [\Delta {\Theta _{\rm inc}}] is very small. The reason for this characteristic is yet unknown.

With reference to the images in Figs. 5[link](b), 5[link](c), the moiré images computed by changing the curvature value from s2 = 0.045′′ mm−1 to s2 = −0.045′′ mm−1 are shown in Figs. 7[link](a), 7[link](b); the other numerical conditions for this computation were kept the same as in the case of Figs. 5[link](b), 5[link](c), though the image in Fig. 7[link](a) was obtained with a slightly smaller value of [\Delta {\Theta _{\rm inc}}] than that for Fig. 5[link](b), in order to obtain the best image (see figure captions of Figs. 5[link] and 7[link]). In the OE fringes in Fig. 7[link](a) the fringe slope is observed to occur in the opposite direction to that in Fig. 5[link](b) for s2 = 0.045′′ mm−1. The slope angle was −7.3° from the x axis in the direct measurement on the image. In the BH fringes in Fig. 7[link](b) the wave form of the fringe-line oscillation is observed to lean to the left side, unlike that in the case of Fig. 5[link](c). From the study in some detail by computation, the lean of the oscillation wave form was always to this direction in the case of [{s_2} \,\lt\, 0], except when [\Delta {\Theta _{\rm inc}}] was very small.

[Figure 7]
Figure 7
Moiré images computed under similar conditions to those of the two images in Figs. 5[link](b) and 5[link](c), but with the curvature deformation of the rear crystal reversed in sign, i.e. with s2 = −0.045′′ mm−1. The image (a) computed with [\Delta\Theta_{\rm inc}] = 1.04′′ in this figure is compared with the image in Fig. 5[link](b), and the image (b) with [\Delta\Theta_{\rm inc}] = 1.20′′ is compared with the image in Fig. 5[link](c).

3.3. Moiré images of bicrystals of t1 = t2 = 1.5 (μt1,2 = 2.193) and t1 = t2 = 1.6 (μt1,2 = 2.339)

Examples of moiré images with t1 = t2 = 1.5 and t1 = t2 = 1.6 (tgap = 0.24) are shown in Figs. 8[link](a)–8[link](d). Similar to the case of t1 = t2 = 0.8 and tgap = 0.24 in Section 3.2[link], the appearance of images in this case when [\Delta {\Theta _{\rm inc}}] is small ([\Delta {\Theta _{\rm inc}}] [\le] 0.12′′) was similar to that of the plane-wave image in Fig. I-10(a) where fringe lines show a gentle oscillation. However, the oscillation amplitude of the fringe lines was much smaller in this case than those in the plane-wave image and in the integrated image of t1 = t2 = 0.8, and therefore the fringes practically looked like HS fringes. When the angular width [\Delta {\Theta _{\rm inc}}] was increased to [\Delta {\Theta _{\rm inc}}] [\ge] 0.16′′, fringe lines came to show a clearly discernible oscillation, with LC bands appearing in the image. As described in Section 3.2[link], OE and BH fringe patterns appeared roughly alternately with the increase in [\Delta {\Theta _{\rm inc}}], though in this case (t1 = t2 = 1.5) the BH fringe pattern seemed to occur more predominantly than the OE fringe pattern. Fig. 8[link](a) shows a moiré image of the BH fringe pattern when [\Delta {\Theta _{\rm inc}}] = 0.36′′, as an example in this domain of [\Delta {\Theta _{\rm inc}}]. Compared with the image with t1 = t2 = 0.8 in Fig. 5[link](c), the step-down (or step-up) of fringe lines at LC band sites (x ≃ 3.0, 8.6) becomes smaller, owing to the increase in the crystal thicknesses. Of the images presented in this paper, this image was computed under the nearest condition to the previous experimental moiré images [see Fig. I-1 in Paper I; further see Yoshimura (1996[Yoshimura, J. (1996). J. Appl. Phys. 80, 2138-2141.], 1997[Yoshimura, J. (1997). J. Appl. Phys. 82, 4697.])]. Two LC bands are observed in this image, while three LC bands appeared in the experimental image with nearly the same width in the x direction. This disagreement is considered to be due to the fact that the employed values of t1, t2 and s2 do not exactly agree with the values of the experimental images. When a further moiré image was computed by assuming that t1 = t2 = 1.6 and s2 = 0.05′′ mm−1 (s1 = 0), the number of LC bands in the resulting image agreed approximately with that in the experimental image.

[Figure 8]
Figure 8
Moiré images computed for bicrystal models of t1 = t2 = 1.5 and t1 = t2 = 1.6. The thicknesses of the component crystals and the incident-wave angular widths for the computed images are as follows: (a) t1 = t2 = 1.5, [\Delta\Theta_{\rm inc}] = 0.36′′; (b) t1 = t2 = 1.5, [\Delta\Theta_{\rm inc}] = 0.90′′; (c) t1 = t2 = 1.6, [\Delta\Theta_{\rm inc}] = 0.90′′; (d) t1 = t2 = 1.6, [\Delta\Theta_{\rm inc}] = 1.50′′. Other numerical conditions for the computation are the same as for Fig. 5[link].

When [\Delta {\Theta _{\rm inc}}] was increased to [\Delta {\Theta _{\rm inc}}] ≃ 0.9′′, the resulting moiré image was the one shown in Fig. 8[link](b). The moiré image again became very close to the HS fringe pattern, with a fringe-line oscillation of a very small amplitude and an almost disappearing LC band pattern. The fringe pattern did not change significantly from the one in Fig. 8[link](b) with further increase in [\Delta {\Theta _{\rm inc}}], though some exceptional images of low contrast were produced at exceptional values of [\Delta {\Theta _{\rm inc}}] such as [\Delta {\Theta _{\rm inc}}] = 1.0′′, 1.5′′ etc. With the gap width of tgap = 0.24 and the crystal thicknesses of t1 = t2 = 1.5, the fringe contrast decreased rapidly with increase in [\Delta {\Theta _{\rm inc}}], and almost disappeared with [\Delta {\Theta _{\rm inc}}] ≃ 1.5′′ (see fringe contrast diagram in Fig. 12).

As has been shown so far, the fringe contrast Vg(x ) varies in an oscillating way with the x coordinate in the image, and also varies in an oscillating way with the angular width [\Delta {\Theta _{\rm inc}}]; the fringe pattern also varies with [\Delta {\Theta _{\rm inc}}] in an oscillating way when [\Delta {\Theta _{\rm inc}}] is in the middle-magnitude domain. In addition to such dependences on x and [\Delta {\Theta _{\rm inc}}], the fringe contrast and the fringe pattern seem to vary in an oscillating way with the crystal thicknesses t1 and t2, as seen from the results of this trial simulation research. The variations in RAAg(x ), RBBg(x ) and [{\Phi _{{\rm{P}},g}}(x )] with t1 = t2 = 1.5 seem to be more gentle than with t1 = t2 = 1.4 and t1 = t2 = 1.6. Figs. 8[link](c) and 8[link](d) show moiré images with t1 = t2 = 1.6 as an example of the case where a clear variation of fringe pattern continues up to a larger value of [\Delta {\Theta _{\rm inc}}] than in the case of t1 = t2 = 1.5. In Fig. 8[link](c), which is compared with Fig. 8[link](b) with t1 = t2 = 1.5, an OE fringe pattern is produced even with [\Delta {\Theta _{\rm inc}}] ≃ 0.9′′. The fringe pattern continued to vary with further increase in [\Delta {\Theta _{\rm inc}}], and almost settled to the pattern shown in Fig. 8[link](d) with [\Delta {\Theta _{\rm inc}}] [\ge] 1.5′′. It may be noted that this fringe pattern is similar to that in Fig. 5[link](d) with t1 = t2 = 0.8, apart from the differences in the position where LC bands occur and in the fringe contrast. The slope angle of the fringes is about 3.5° in this case also.

3.4. Moiré images of bicrystals of t1 = t2 ≥ 2.5 (μt1,2 ≥ 3.655)

When the crystal thicknesses were increased from t1 = t2 = 1.5 and t1 = t2 = 1.6, the fringe patterns that were produced and their variations with [\Delta {\Theta _{\rm inc}}] were basically analogous to the results for t1 = t2 = 0.8, t1 = t2 = 1.5 and t1 = t2 = 1.6, until the thickness was increased to t1 = t2 ≃ 2.5. Figs. 9[link] and 10[link] show examples of moiré images with t1 = t2 [\ge] 2.5 and graphs of the characteristic curves concerned. Figs. 9[link](a) and 10(a) show, respectively, a moiré image and the characteristic curves with t1 = t2 = 2.0, to be compared with those with t1 = t2 [\ge] 2.5. The LC bands and the fringe-line bends are still observed at the positions x ≃ 1.8, 5.4, 9.0 (indicated by arrows), though the fall in fringe contrast is considerably moderated. When the crystal thicknesses were increased further, a similar aspect of the produced fringe patterns continued up to t1 = t2 = 2.4. However, with t1 = t2 ≃ 2.5, the OE fringe patterns and BH fringe patterns with LC bands disappeared almost completely, and all the fringe patterns became of HS fringes as shown in Fig. 9[link](b), except for some exceptional cases. From the graph of the characteristic curves in Fig. 10[link](b), it is seen that this transition to HS fringes results from the convergence of RAAg(x ) and RBBg(x ) to their respective finite values, unlike the case of Figs. 1[link](c) and 2[link](c) with t1 = t2 = 0.8 and tgap = 0.

[Figure 9]
Figure 9
Moiré images computed for bicrystal models of t1 = t2 [\ge] 2. The thicknesses of the component crystals and the incident-wave angular width for the computed images are as follows: (a) t1 = t2 = 2.0, [\Delta\Theta_{\rm inc}] = 1.08′′; (b) t1 = t2 = 2.5, [\Delta\Theta_{\rm inc}] = 0.90′′; (c) t1 = t2 = 3.0, [\Delta\Theta_{\rm inc}] = 0.90′′; (d) t1 = t2 = 3.0, [\Delta\Theta_{\rm inc}] = 1.0′′. Other numerical conditions for the computation are the same as for Fig. 5[link].
[Figure 10]
Figure 10
Characteristic curves of the moiré images in Fig. 9[link]. Graphs (a), (b) and (c) are related to the moiré images (a), (b) and (c), respectively, in Fig. 9[link]. For further information see text.

When the crystal thicknesses became t1 = t2 = 2.6, OE and BH fringe patterns appeared again, though with weak LC bands, for the reverse tendency in the dependence on the crystal thicknesses as shown in Figs. 8[link](c) and 8[link](d). However, when the crystal thicknesses were further increased, the fringe patterns produced again became of HS fringes as seen when t1 = t2 = 2.5. For the crystal thicknesses approximately t1 = t2 [\ge] 3, the fringe patterns produced were always nearly of HS fringe type, in spite of the increase in t1 and t2. This result indicates that the minimum specimen thickness for obtaining a good moiré image without disturbance from PL oscillation and crystal strain is t1 = t2 ≃ 3 ([\mu {t_1}] = [\mu {t_2}] ≃ 4.4). Figs. 9[link](c) and 10[link](c) show, respectively, a moiré image with t1 = t2 ≃ 3 and a graph of the characteristic curves concerned. The values of the angular width [\Delta {\Theta _{\rm inc}}], with which the produced fringe patterns were not of HS fringes, were [\Delta\Theta_{\rm inc}] = 0.5′′, 1.0′′, 1.5′′ etc. with any value of t1 or t2. With these values of [\Delta {\Theta _{\rm inc}}], the image intensity and fringe contrast momentarily fell markedly, and OE or BH fringe patterns appeared in a disordered manner in a low-intensity image. Fig. 9[link](d) shows an example of such a specific moiré image with t1 = t2 = 3.0 and [\Delta\Theta_{\rm inc}] = 1.0′′. The reason for this specificity is as yet unknown.

3.5. Diffraction intensity curves for the presented moiré images and wide-view diagram of the fringe contrast variation

3.5.1. Curves of mean image intensity versus Δθm

Fig. 11[link] gives diffraction intensity curves that show the dependence of the mean image intensity of moiré images on the mid deviation angle [\Delta {\theta _{\rm m}}], for the moiré images presented in this paper so far. The mean image intensity [[R_{o,g'}(x) + R_{g,g'}(x )]] averaged over the represented range of moiré images [0.3 \le x \le 11.5] is plotted. The value of intensity is represented as the ratio of the integrated image intensity to the total incident intensity [{I_{\rm inc}} \times \Delta {\Theta _{\rm inc}}] [see equations (12a)[link], (12b)[link]]. These intensity curves would approximate well real rocking curves measured experimentally. The curves (B)–(H) are related to the presented moiré images as follows: (B) Figs. 1[link](b) and 3[link](a), (C) Fig. 5[link](c), (E) Fig. 8[link](a), (F) Fig. 8[link](b) and (H) Fig. 9[link](b). These intensity curves are the same for all values of tgap, including the case of tgap = 0. This is because the factor tgap is not involved in the functions [{I_{o,g'}}({x,u} )] and [{I_{g,g'}}({x,u} )] which are the integrands in [{R_{o,g'}}(x )] and [{R_{g,g'}}(x)] [see equations (I-22a,b)]. Furthermore, the intensity curves do not depend on the value of [\Delta \rho], which gives the spacing of moiré fringes (in the case of rotation moiré), and are the same for all values of [\Delta \rho]. The reason for this is also that the factor [\Delta \rho] is not involved in [{I_{o,g'}}({x,u} )] and [{I_{g,g'}}({x,u} )].

[Figure 11]
Figure 11
Plots of the mean intensity of moiré images [[{R_{og'}}(x ) + {R_{gg'}}(x )]] versus the mid deviation angle [\Delta {\theta _{\rm m}}]. Numerical conditions for the computation are as follows: curve (A) t1 = t2 = 0.8, [\Delta\Theta_{\rm inc}] = 0.5′′, s2 = 0; (B) t1 = t2 = 0.8, [\Delta\Theta_{\rm inc}] = 0.5′′; (C) t1 = t2 = 0.8, [\Delta\Theta_{\rm inc}] = 1.2′′; (D) t1 = t2 = 0.8, [\Delta\Theta_{\rm inc}] = 2.0′′; (E) t1 = t2 = 1.5, [\Delta\Theta_{\rm inc}] = 0.36′′; (F) t1 = t2 = 1.5, [\Delta\Theta_{\rm inc}] = 0.9′′; (G) t1 = t2 = 1.5, [\Delta\Theta_{\rm inc}] = 2.0′′; (H) t1 = t2 = 2.5, [\Delta\Theta_{\rm inc}] = 0.9′′; (I) t1 = t2 = 2.5, [\Delta\Theta_{\rm inc}] = 2.0′′; (J) t1 = t2 = 2.5, [\Delta\Theta_{\rm inc}] = 4.0′′. The curvature strength of the rear crystal was s2 = 0.045′′ mm−1 for all the curves from (B) to (J). The half widths (FWHM) of these intensity curves are: 2.18′′, 2.15′′, 2.55′′ and 2.79′′ in curves (A), (B), (C) and (D), respectively; 1.42′′, 1.57′′, 2.27′′ in curves (E), (F) and (G), respectively; 1.06′′, 2.09′′, 4.01′′ in curves (H), (I) and (J), respectively.

The curve (A) is the intensity curve when the curvature deformation is not induced, i.e. s1 = s2 = 0, with t1 = t2 = 0.8. Unlike the other curves (B)–(J), this curve has a symmetric shape with respect to the coordinate origin [\Delta {\theta _{\rm m}}] = 0. This characteristic is understood as follows: though the integrand functions [{I_{o,g'}}({x,u} )] and [{I_{g,g'}}({x,u} )] in [{R_{o,g'}}(x )] and [{R_{g,g'}}(x)] are both asymmetric functions, the function [[I_{o,g'}({x,u} ) + I_{g,g'}({x,u} )]] becomes symmetric with respect to u, i.e. [\Delta \theta] (u = [K\Delta \theta \sin 2{\theta _{\rm{B}}}]), since the relationship [{I_{o,g'}}({x,u} )] = [{I_{g,g'}}(x, - u)] holds when [(\Delta{\bf g}_1\cdot \hat{\bf K}_g)] = [(\Delta{\bf g}\cdot\hat{\bf K}_g)] = [(\Delta{\bf g}\cdot\hat{\bf K}_o)] = 0 [i.e. [({\Delta d} / d)_1] = [({\Delta d} / d)] = 0, [\Delta {\omega _1}] = [\Delta \omega] = 0 (s1 = s = 0)] and [{\gamma _o}] = [{\gamma _g}], t1 = t2. When these conditions do not hold, even in part, the function [[I_{o,g'}({x,u} ) + I_{g,g'}({x,u} )]] becomes asymmetric. Based on this understanding, it is considered to be due to the failure of the condition s = 0 that the curves (B), (C) etc. are asymmetric as a whole, with their peak positions slightly shifted to the low-angle side. In fact, when the computation was performed on the assumption that [s \,\lt\, 0], or when the x position corresponding to the centre of curvature was moved to the left side with [s \,\gt\, 0], the peak position of the intensity curves was displaced to the high-angle side.

As to the half width (FWHM) of the intensity curves, that of the curve (A) was [\Delta\Theta_{\rm hw}] = 2.18′′, while the half width of the intrinsic diffraction curve with no absorption is [\Delta\Theta_{\rm hw}] = 2.29′′. The half widths of the other curves were as described in the figure caption. The image intensity decreases fairly rapidly with increase in the crystal thicknesses, due to the increased absorption of X-rays. When the angular width [\Delta {\Theta _{\rm inc}}] is increased, the peak of the intensity curves is lowered under the present definition of diffracted intensity, while their half widths increase. When [\Delta {\Theta _{\rm inc}}] becomes larger than circa 2′′, the increase in the half width with [\Delta {\Theta _{\rm inc}}] proceeds relatively rapidly. The curve (J) (t1 = t2 = 2.5) is an example of a curve of a large half width ([\Delta\Theta_{\rm hw}] = 4.01′′) with [\Delta {\Theta _{\rm inc}}] = 4′′. When t1 and t2 took values other than those in this case, intensity curves with [\Delta {\Theta _{\rm inc}}] = 4′′ were also of a similar wide flat-top shape but with a slight difference in the top height; this was also the case when s1 = s2 = 0. Such shapes with a wide flat top are rarely seen for X-ray rocking curves. The intensity curve with [\Delta {\Theta _{\rm inc}}] = 10′′, which is related to the moiré images in Figs. 1[link](c), 3[link](b) and 5[link](d), would be like a low-height horizontal line if it is plotted in the same graph of Fig. 11[link].

3.5.2. Curves of mean fringe contrast versus ΔΘinc

Fig. 12[link] shows the characteristic curves of the fringe contrast variation with the angular width [\Delta {\Theta _{\rm inc}}]. Similar to the image intensity curves in Fig. 11[link], fringe contrasts Vg(x ) averaged over the represented range of moiré images [0.3 \le x \le 11.5] are plotted. The curve (A) in the inset shows the variation of the fringe contrast related to the moiré image when s1 = s2 = 0 with t1 = t2 = 0.8 and tgap = 0. The other curves (B), (D) etc. are related to the moiré images presented so far as follows: (B) Fig. 1[link], (D) Fig. 3[link], (F) Fig. 5[link], (G) Figs. 8[link](a), 8[link](b) and (H) Fig. 9[link](b) [s2 = 0.045′′ mm−1 in all the curves of (B) to (J)]. These fringe contrast curves change with the curvature values s1 and s2, but do not change with the value of [\Delta \rho]. The reason is that [\Delta \rho] is not involved in the expression of Vg(x ) in equation (15)[link]. While it has been observed in Figs. 2[link], 4[link], 6[link] etc. that the fringe contrast oscillates with the x coordinate in the image, the oscillation with the [\Delta {\Theta _{\rm inc}}] value in the mean fringe contrast Vg is seen understandably in this graph. The Vg oscillation occurs in the domain of relatively small [\Delta {\Theta _{\rm inc}}], and nearly converges in the domain of [\Delta {\Theta _{\rm inc}}\ge 2''] (though a very weak oscillation continues further). The curves of Vg variation gradually descend with [\Delta {\Theta _{\rm inc}}], and almost settle to the given values. Like the Vg(x ) oscillation with the x coordinate, the oscillation with [\Delta {\Theta _{\rm inc}}] in the mean fringe contrast Vg would also be related to the oscillations in RAAg(x ) and RBBg(x ), which arise from the combined effect of the PL oscillation and gap phase difference.

[Figure 12]
Figure 12
Plots of the mean fringe contrast of computed moiré images versus the incident-wave angular width [\Delta {\Theta _{\rm inc}}]. Numerical conditions for the computation are as follows: curve (A) t1 = t2 = 0.8, tgap = 0.0, s2 = 0; (B) t1 = t2 = 0.8, tgap = 0.0; (C) t1 = t2 = 0.8, tgap = 0.02; (D) t1 = t2 = 0.8, tgap = 0.05; (E) t1 = t2 = 0.8, tgap = 0.2; (F) t1 = t2 = 0.8, tgap = 0.24; (G) t1 = t2 = 1.5, tgap = 0.24; (H) t1 = t2 = 2.5, tgap = 0.24; (I) t1 = t2 = 1.5, tgap = 0.0; (J) t1 = t2 = 2.5, tgap = 0.0. The curvature strength of the rear crystal was s2 = 0.045′′ mm−1 for all the curves from (B) to (J).

With the increase in the crystal thicknesses t1, t2 and the gap width tgap, the curves of Vg variation descend rapidly with increase in [\Delta {\Theta _{\rm inc}}] [curves (D)–(H)]. In the domain of small [\Delta {\Theta _{\rm inc}}], the Vg curves oscillate with a large amplitude for a small change in [\Delta {\Theta _{\rm inc}}]. In the case of t1 = t2 = 2.5 [curve (H)], changes in the fringe contrast corresponding to the minima and maxima of the Vg curve were recognized clearly among the moiré images, although it was not so clear in the case of t1 = t2 = 1.5 [curve (G)]. It must be added that no other co-occurring change was found in these moiré images, in spite of the large, rapid change in the fringe contrast. By adjusting the magnitude of [\Delta {\Theta _{\rm inc}}], it would be possible to adjust the fringe contrast in the obtained moiré images to a fair extent. This would be useful knowledge when planning future moiré experiments.

To expand on the case of low fringe contrast, the moiré fringe pattern can be fairly well observed with the contrast [{V_g} \ge 0.01], so far as this computer-simulation study (256 graduations) is concerned, aside from the answer given by the experiment. However, with [{V_g} \le 0.005] [at [\Delta {\Theta _{\rm inc}}] ≃ 1.0′′, 1.5′′ on the curves (G) and (H)] the fringe pattern became almost invisible even in this simulation study. In the following we assume that Vg = 0.01 is the visibility limit of the fringe pattern. According to this graph in Fig. 12[link], the fringe contrast becomes practically zero with [\Delta {\Theta _{\rm inc}} \ge 2''] in the integrated images with t1 = t2 = 1.5 and t1 = t2 = 2.5 (tgap = 0.24). This gives an accurate theoretical explanation for the known experimental fact that the moiré image from a bicrystal specimen (tgap ≃ 0.25) which appears for the quasi-plane-wave incidence ([\Delta {\Theta _{\rm inc}} \ll 2'']) becomes invisible in Lang traverse topography ([\Delta {\Theta _{\rm inc}}\ge 60'']) (though visible in section topography). Furthermore, the graph in Fig. 12[link] suggests that the moiré image could be observed with some fringe contrast even by Lang topography, if the sample crystal was somewhat thinner than that mentioned. In the results of the computation, the fringe contrast was Vg = 0.015 with [\Delta {\Theta _{\rm inc}}] = 60′′, when t1 = t2 = 0.8 and tgap = 0.24 [curve (F)]. Additionally, when the gap width was increased with the crystal thicknesses kept at t1 = t2 = 1.5, the fringe contrasts obtained by the computation were Vg = 0.042 for tgap = 0.5, [\Delta {\Theta _{\rm inc}}] = 0.5′′, Vg = 0.026 for tgap = 0.5, [\Delta {\Theta _{\rm inc}}] = 1.0′′, and Vg = 0.021 for tgap = 1.0, [\Delta {\Theta _{\rm inc}}] = 0.5′′. These estimates suggest that moiré fringes can be observed even with large gap widths if the angular width of the incident wave is suitably narrowed.

Curves (I) and (J) (in the inset) show the Vg variation curves with t1 = t2 = 1.5, tgap = 0 and with t1 = t2 = 2.5, tgap = 0 for reference. The curves (G) and (H), which give a low fringe contrast with tgap = 0.24, change to such curves as above, giving a high fringe contrast with tgap = 0. In these cases, while the image intensity [[{R_{o,g'}}(x ) + {R_{g,g'}}(x )]] is decreased considerably, the fringe contrast Vg is increased greatly.

4. Concluding remarks

(i) As described above, in connection with the question of what type of moiré images are given by the X-ray moiré-fringe theory in Paper I when practical experimental conditions are applied, we have simulation-computed and observed many integrated moiré images for an assumed bicrystal specimen; this was done by changing the crystal thicknesses and incident-wave angular width over a wide range, and by setting the width of the interspacing gap to different values. It was shown that the interference pattern of intrinsic moiré images is considerably modified by the combined effect of PL oscillation and gap phase difference, related to the effect of strain in the specimen. As examples of such modified and peculiar fringe patterns in the integrated images, the OE fringe pattern and BH fringe pattern were observed. While the BH fringe pattern has been observed in previous experiments, occurrence of the OE fringe pattern in this study is quite a new finding. The experimental verification of it is hoped to be performed in the near future. Furthermore, not only the images of this newly found fringe pattern, but also all the moiré images shown by the present computations are hoped to be really observed and verified in future experiments, although there is no particular question on the correctness of the computed images..

(ii) This theory of moiré fringes was developed by assuming a gapped bicrystal as the specimen for moiré images. It is related to previous experimental studies by the author (e.g. Yoshimura, 1996[Yoshimura, J. (1996). J. Appl. Phys. 80, 2138-2141.], 1997[Yoshimura, J. (1997). J. Appl. Phys. 82, 4697.]), where gapped bicrystals were used as the specimen. However, studies of X-ray moiré fringes are made mainly using an X-ray interferometer at present, and therefore this theory of moiré fringes should be extended so as to be applicable for interferometer moiré fringes. Nevertheless, the description in this paper has been made with the intention of understanding fully the properties of bicrystal moiré fringes. As a result of such moiré-fringe study, the relationship of fringe contrast versus angular width of the incident wave was extensively studied as shown in the graph in Fig. 12[link]. As a discussion of this graph, the lower limits of the visibility of bicrystal moiré fringes were described, and it was suggested that if the angular width of the incident wave is narrowed sufficiently, moiré fringes could be observed even with a gap width of 0.5, 1.0 or more (mm). This result would be worthy of attention although its validity must be further studied by experiment. To the best of the author's knowledge, this theory of moiré fringes would be the first theory applicable to the exact and detailed theoretical understanding of moiré images observed in experiments. Though the topic of X-ray moiré fringes does not attract much interest at present, the theory could be useful for future studies of diffraction moiré fringes, and for the development of related techniques.

References

First citationAuthier, A., Milne, A. D. & Sauvage, M. (1968). Phys. Status Solidi, 26, 469–484.  CrossRef CAS Web of Science Google Scholar
First citationBrádler, J. & Lang, A. R. (1968). Acta Cryst. A24, 246–247.  CrossRef IUCr Journals Web of Science Google Scholar
First citationHart, M. & Milne, A. D. (1970). Acta Cryst. A26, 223–229.  CrossRef CAS IUCr Journals Web of Science Google Scholar
First citationHirsch, P. B., Howie, A., Nicholson, R. B., Pashley, D. W. & Whelan, M. J. (1965). Electron Microscopy of Thin Crystals. London: Butterworths.  Google Scholar
First citationLang, A. R. (1968). Nature, 220, 652–657.  CrossRef CAS Google Scholar
First citationPinsker, Z. G. (1978). Dynamical Scattering of X-rays in Crystals. Berlin: Springer-Verlag.  Google Scholar
First citationSugii, K., Matsushita, T. & Kohra, K. (1971). J. Phys. Soc. Jpn, 31, 1295–1296.  CAS Google Scholar
First citationYoshimura, J. (1991). Phys. Status Solidi A, 125, 429–440.  CrossRef Web of Science Google Scholar
First citationYoshimura, J. (1996). J. Appl. Phys. 80, 2138–2141.  CrossRef CAS Google Scholar
First citationYoshimura, J. (1997). J. Appl. Phys. 82, 4697.  CrossRef Google Scholar
First citationYoshimura, J. (2015). Acta Cryst. A71, 368–381.  CrossRef IUCr Journals Google Scholar

This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.

Journal logoFOUNDATIONS
ADVANCES
ISSN: 2053-2733
Follow Acta Cryst. A
Sign up for e-alerts
Follow Acta Cryst. on Twitter
Follow us on facebook
Sign up for RSS feeds