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

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775

Dynamic crystal rotation resolved by high-speed synchrotron X-ray Laue diffraction

CROSSMARK_Color_square_no_text.svg

aThe Peac Institute of Multiscale Sciences, Chengdu, Sichuan 610031, People's Republic of China, bKey Laboratory of Advanced Technologies of Materials, Ministry of Education, Southwest Jiaotong University, Chengdu, Sichuan 610031, People's Republic of China, and cAdvanced Photon Source, Argonne National Laboratory, Argonne, IL 60439, USA
*Correspondence e-mail: sluo@pims.ac.cn

Edited by A. F. Craievich, University of São Paulo, Brazil (Received 24 December 2015; accepted 4 February 2016; online 30 March 2016)

Dynamic compression experiments are performed on single-crystal Si under split Hopkinson pressure bar loading, together with simultaneous high-speed (250–350 ns resolution) synchrotron X-ray Laue diffraction and phase-contrast imaging. A methodology is presented which determines crystal rotation parameters, i.e. instantaneous rotation axes and angles, from two unindexed Laue diffraction spots. Two-dimensional translation is obtained from dynamic imaging by a single camera. High-speed motion of crystals, including translation and rotation, can be tracked in real time via simultaneous imaging and diffraction.

1. Introduction

Grain or crystal rotation is a key mechanism of accommodating elastic or plastic deformation in crystalline solids (Chen et al., 2013[Chen, P., Mao, S. C., Liu, Y., Wang, F., Zhang, Y. F., Zhang, Z. & Han, X. D. (2013). Mater. Sci. Eng. A, 580, 114-124.]). In quasi-static experiments, crystal rotation is normally resolved with electron backscattering diffraction (Zaafarani et al., 2008[Zaafarani, N., Raabe, D., Roters, F. & Zaefferer, S. (2008). Acta Mater. 56, 31-42.]) or X-ray Laue diffraction (Maaß et al., 2007[Maaß, R., Van Petegem, S., Van Swygenhoven, H., Derlet, P. M., Volkert, C. A. & Grolimund, D. (2007). Phys. Rev. Lett. 99, 145505.], 2008[Maaß, R., Van Petegem, S., Grolimund, D., Van Swygenhoven, H., Kiener, D. & Dehm, G. (2008). Appl. Phys. Lett. 92, 071905.]). Single-shot multiframe in situ real-time measurements on crystal rotation during dynamic loading have been an experimental challenge, especially for high strain rate loading in the nano- to micro-second regimes with a split Hopkinson pressure bar (SHPB) (Kolsky, 1949[Kolsky, H. (1949). Proc. Phys. Soc. B, 62, 676-700.]) or a gas gun (Fowles et al., 1970[Fowles, G. R. (1970). Rev. Sci. Instrum. 41, 984-996.]). In powder technology and fluid mechanics, rotation of particles or granules plays a significant role in the gas–solid interactions, including two-phase flows, heat transfer and coal combustion (Kajishima, 2004[Kajishima, T. (2004). Int. J. Heat Fluid Flow, 25, 721-728.]; Goldschmidt et al., 2004[Goldschmidt, M. J. V., Beetstra, R. & Kuipers, J. A. M. (2004). Powder Technol. 142, 23-47.]; Sun & Battaglia, 2006[Sun, J. & Battaglia, F. (2006). Chem. Eng. Sci. 61, 1470-1479.]). Detection of high-speed particle rotation is in urgent need. For example, Wu et al. (2008[Wu, X., Wang, Q., Luo, Z., Fang, M. & Cen, K. (2008). Powder Technol. 181, 21-30.]) used optical imaging to measure particle rotation at about 500 revolutions per second (r s−1) in a circulating fluidized bed. X-ray diffraction is a useful complement for rotation measurements given its penetration capability and higher resolution. X-ray diffraction was used to study the quasi-static packing and deformation characteristics of sands (Hall et al., 2011[Hall, S. A., Wright, J., Pirling, T., Andò, E., Hughes, D. J. & Viggiani, G. (2011). Granul. Matter, 13, 251-254.]). Other applications may involve measuring remotely translation and rotation of a fast-moving object.

Synchrotron X-ray Laue diffraction possesses unique advantages for lattice-level measurements under dynamic loading and has been under intensive development along with synchrotron X-ray imaging (Luo et al., 2012[Luo, S. N., Jensen, B. J., Hooks, D. E., Fezzaa, K., Ramos, K. J., Yeager, J. D., Kwiatkowski, K. & Shimada, T. (2012). Rev. Sci. Instrum. 83, 073903.]; Hudspeth et al., 2013[Hudspeth, M., Claus, B., Dubelman, S., Black, J., Mondal, A., Parab, N., Funnell, C., Hai, F., Qi, M. L., Fezzaa, K., Luo, S. N. & Chen, W. (2013). Rev. Sci. Instrum. 84, 025102.]; Lambert et al., 2014[Lambert, P. K., Hustedt, C. J., Vecchio, K. S., Huskins, E. L., Casem, D. T., Gruner, S. M., Tate, M. W., Philipp, H. T., Woll, A. R., Purohit, P., Weiss, J. T., Kannan, V., Ramesh, K. T., Kenesei, P., Okasinski, J. S., Almer, J., Zhao, M., Ananiadis, A. G. & Hufnagel, T. C. (2014). Rev. Sci. Instrum. 85, 093901.]; Turneaure et al., 2009[Turneaure, S. J., Gupta, Y., Zimmerman, K., Perkins, K., Yoo, C. & Shen, G. (2009). J. Appl. Phys. 105, 053520.]; Wei et al., 2004[Wei, Q., Jiao, T., Ramesh, K. & Ma, E. (2004). Scr. Mater. 50, 359-364.]) for studying high-strain-rate deformation, fracture and phase change of a variety of materials. In particular, a real-time in situ simultaneous X-ray imaging and diffraction technique has been demonstrated recently at the Advanced Photon Source (Fan et al., 2014[Fan, D., Lu, L., Li, B., Qi, M. L., Zhao, F., Sun, T., Fezzaa, K., Chen, W. & Luo, S. N. (2014). Rev. Sci. Instrum. 85, 113902.]; Hudspeth et al., 2015[Hudspeth, M., Sun, T., Parab, N., Guo, Z., Fezzaa, K., Luo, S. & Chen, W. (2015). J. Synchrotron Rad. 22, 49-58.]). In this work, we choose single-crystal Si, brittle and of high modulus, for dynamic loading with a SHPB, and perform simultaneous X-ray imaging and Laue diffraction, in order to establish an illustrative case for quantitative rotation and translation measurements. A theoretical procedure is introduced to resolve three-dimensional high-speed rotation of a single crystal with two unindexed diffraction spots. Three-dimensional crystal rotation at approximately 300 r s−1 is resolved. High-speed motion of crystals, including two-dimensional translation and three-dimensional rotation, can be tracked in real time from simultaneous imaging and diffraction.

2. Methodology

The experimental methodology for simultaneous synchrotron X-ray imaging and diffraction has been established at beamline 32-ID-B of the Advanced Photon Source (Fan et al., 2014[Fan, D., Lu, L., Li, B., Qi, M. L., Zhao, F., Sun, T., Fezzaa, K., Chen, W. & Luo, S. N. (2014). Rev. Sci. Instrum. 85, 113902.]; Hudspeth et al., 2015[Hudspeth, M., Sun, T., Parab, N., Guo, Z., Fezzaa, K., Luo, S. & Chen, W. (2015). J. Synchrotron Rad. 22, 49-58.]). The schematic setup for simultaneous diffraction and phase contrast imaging is shown in Fig. 1[link], along with the coordinate system. We perform high-strain-rate loading on a brittle Si single-crystal with SHPB. Prior to loading, the crystallographic directions [0[\bar{1}]1], [2[\bar{1}][\bar{1}]] and [111] are oriented approximately along the x, y and z directions, respectively. Dimensions of the sample perpendicular to the X-ray direction ([O^{\,\prime}B]) are 3 mm × 3 mm and the thickness along the X-ray direction is 1 mm. After the gas gun of the SHPB is fired, a compression pulse is imposed on the sample. A tensile pulse ensues upon reflection between the bar–sample interface, and the incident bar is then separated from the sample. The sample retains its overall integrity and is subjected to rigid-body motion, including translation and rotation, when a stress is below the fracture strength of Si (∼1 GPa).

[Figure 1]
Figure 1
Experimental setup for dynamic X-ray diffraction with simultaneous imaging and the coordinate system (Oxyz). Images form on the [O^{{\,\prime\prime\prime}}D] plane along the [O^{\,\prime}O] direction. [O^{{\,\prime\prime}}A] is the diffraction detector plane. [OA^{{\prime}}] is the virtual detector for diffraction. SHPB: split Hopkinson pressure bar (loading device); [O^{\,\prime}]: sample position; [BO^{\,\prime}]: incident X-rays; [O^{\,\prime}A^{{\prime}}]: diffracted X-rays; [A^{{\prime}}]: a diffraction spot on the virtual detector (the xy plane). Also see text.

We perform simultaneous X-ray imaging and diffraction measurements on the sample. The probe X-rays are from an APS `undulator A' light source and the undulator gap is 25 mm. The spectral flux–photon energy curves of the undulator source are presented elsewhere (Luo et al., 2012[Luo, S. N., Jensen, B. J., Hooks, D. E., Fezzaa, K., Ramos, K. J., Yeager, J. D., Kwiatkowski, K. & Shimada, T. (2012). Rev. Sci. Instrum. 83, 073903.]). The probe X-rays ([O^{\,\prime}O]) illuminate the sample perpendicular to the loading direction. The transmitted X-rays are collected by an imaging scintillator, Ce-doped Lu3Al5O12 (LuAG:Ce), and recorded with a Photron FastCam SA-Z ([O^{{\,\prime\prime\prime}}D]). The exposure time is 0.35 µs and the frame interval is 10 µs. At the same time, the scattered light is collected by a diffraction scintillator, Ce-doped Lu2–2xY2xSiO5 (LYSO:Ce), and is recorded with another SA-Z after passing a microchannel plate (MCP, Quantum Leap E, Stanford Computer Optics, Inc.), synchronized with the imaging camera. The exposure time is 0.25 µs and the frame interval is 10 µs. The number of activated pixels is 384×408. The diffraction camera is mounted on a rotation stage for continuous rotation in the circumferential direction. A Huber 410 goniometer is used to control the rotation stage with a resolution of 0.001°. The scintillator, MCP and SA-Z together are equivalent to an area detector ([O^{{\,\prime\prime}}A]). The pixel size is calibrated to be 60 µm × 60 µm for the current geometry. The sample-to-diffraction detector distance (d) is 230 mm and the angle [\omega = \angle OO^{\,\prime}O^{{\,\prime\prime}}] between the incident beam and the normal of the diffraction detector plane is 24.5°.

Each Laue diffraction spot corresponds to lattice planes with certain spatial orientation. Lattice deformation and crystal rotation can both contribute to its shift. Laue diffraction pattern analysis has been widely applied to obtain strain tensor (Li et al., 2015[Li, Y., Qian, D., Xue, J., Wan, J., Zhang, A., Tamura, N., Song, Z. & Chen, K. (2015). Appl. Phys. Lett. 107, 181902.]; Liu et al., 2014[Liu, X., Thadesar, P. A., Taylor, C. L., Oh, H., Kunz, M., Tamura, N., Bakir, M. S. & Sitaraman, S. K. (2014). Appl. Phys. Lett. 105, 112109.]; Wang et al., 2011[Wang, L., Barabash, R. I., Yang, Y., Bieler, T. R., Crimp, M. A., Eisenlohr, P., Liu, W. & Ice, G. E. (2011). Metall. Mater. Trans. A, 42, 626-635.]; Barabash et al., 2010[Barabash, O. M., Santella, M., Barabash, R. I., Ice, G. E. & Tischler, J. (2010). JOM, 62, 29-34.]) and grain rotation (Maaß et al., 2007[Maaß, R., Van Petegem, S., Van Swygenhoven, H., Derlet, P. M., Volkert, C. A. & Grolimund, D. (2007). Phys. Rev. Lett. 99, 145505.], 2008[Maaß, R., Van Petegem, S., Grolimund, D., Van Swygenhoven, H., Kiener, D. & Dehm, G. (2008). Appl. Phys. Lett. 92, 071905.]); there are many dedicated software packages available for this analysis (Tamura, 2014[Tamura, N. (2014). Strain and Dislocation Gradients from Diffraction: Spatially Resolved Local Structure and Defects, edited by R. Barabash and G. Ice. Imperial College Press.]; Micha & Robach, 2014[Micha, J. S. & Robach, O. (2014). Lauetools, https://sourceforge.net/projects/lauetools/.]; Huang, 2010[Huang, X. R. (2010). J. Appl. Cryst. 43, 926-928.]). A method based on digital image correlation using Laue spots as speckles has been developed recently for strain analysis (Borbély, 2015[Borbély, A. (2015). J. Appl. Cryst. 48, 1614-1616.]). However, the effective area of detectors available for dynamic events is very limited and the sample–detector distance has to be sufficiently large in order to accommodate the imaging camera as well as protective components for impact loading. As a result, the number of recorded diffraction spots is small (e.g. 1 or 2). Meanwhile, a conventional Hopkinson bar setup is not intended for high-accuracy alignment as required for Laue diffraction indexing. So indexing the spots and orienting a crystal are currently impractical, and the aforementioned Laue analysis tools are not applicable. In our experiments, the sample moves stress-free, so the remaining lattice strain after low-speed impact is small and its contribution to diffraction spot movement is negligible compared with rotation. Therefore, we establish a quantitative relationship between three-dimensional crystal rotation and two-dimensional movement of two or more diffraction spots, which is elaborated in Appendix A[link]. In this method, neither crystal orientation nor indices of diffraction spots are needed to obtain rotation parameters (rotation axes and rotation angles).

3. Results and discussion

The dynamic loading experiment was conducted at a projectile velocity of 20 m s−1. Results from simultaneous imaging and diffraction are presented in Figs. 2[link] and 3[link].

[Figure 2]
Figure 2
Selected phase-contrast images of the sample obtained simultaneously with diffraction shown in Fig. 3[link]. The black stripes on the left refer to the recoiling steel bar after impact. fN denotes the Nth frame. The frame interval is 10 µs. Loading direction is from left to right.
[Figure 3]
Figure 3
Time-resolved synchrotron X-ray diffraction. Columns 1 and 3 are two-dimensional diffraction patterns, and columns 2 and 4 are the corresponding [2\theta]γ plots (filled circles). Empty circles denote the diffraction spots predicted from the rotation parameters. V: vertical; H: horizontal.

During recoil, the incident bar separates from the Si sample at 1 m s−1. The sample remains essentially intact except that minor fragments form from stress concentrations during impact and flies in free rigid-body motion along the loading direction (Fig. 2[link]). The two-dimensional phase-contrast images supply information on planar translation and rotation of the sample, but not on out-of-plane motion. The translation velocities are approximately 5 m s−1 and 1 m s−1 in the x and y directions, respectively, and the in-plane rotation speed is small.

Eight representative frames of diffraction patterns containing two diffraction spots are chosen for analysis (Fig. 3[link], columns 1 and 3). Two diffraction spots shift simultaneously from the lower right to upper left of the screen, without change in spot shape. This shift is attributed to sample or crystal rotation, rather than lattice compression or tension, because the sample has already been subjected to unloading and detached from the incident bar during the time window of interest here. The minute residual lattice deformation (Noyan & Cohen, 2013[Noyan, I. C. & Cohen, J. B. (2013). Residual Stress: Measurement by Diffraction and Interpretation. Berlin: Springer-Verlag.]) cannot induce such a pronounced shift. Note that there is only one diffraction spot recorded prior to frame 1 (Fig. 3[link]). As the sample undergoes rotation during impact (before frame 1), indexing the spots before and after frame 1 is not possible even though its initial orientation is known. However, a quantitative analysis of crystal rotation can be conducted from two Laue spots with the scheme detailed in Appendix A[link].

As each harmonic of the undulator radiation has a finite bandwidth (∼2% for the APS undulator A light source) (Fan et al., 2014[Fan, D., Lu, L., Li, B., Qi, M. L., Zhao, F., Sun, T., Fezzaa, K., Chen, W. & Luo, S. N. (2014). Rev. Sci. Instrum. 85, 113902.]), some diffraction spots do not vanish from the field of view during rotation, within the time window of observation. As discussed above, relative translation and rotation of the sample between adjacent frames on the Oxy plane is tiny due to short frame intervals (Fig. 3[link]), so the small-angle assumption in equation (8)[link] of Appendix A[link] is reasonable.

Each pixel on the diffraction detector is mapped onto the [2\theta]γ plane. Given the finite size of a diffraction spot, we use Gaussian fitting to obtain the coordinates of its center ([2\theta], γ). The diffraction spots are plotted on the [2\theta]γ plane in columns 2 and 4 of Fig. 3[link]. Then, the relative rotation axes (expressed in three projection quantities) and the rotation angles between adjacent frames are calculated (Table 1[link]). The rotation directions during the free motion of the sample appear to be random. The rotation angles are extremely small between two neighboring frames (0.05–0.13°), consistent with the small-angle assumption. The sample mainly undergoes translation along the loading (x) direction with concurrent, small and irregular rotation. The high angular resolution of dynamic Laue diffraction is advantageous for such applications as determining grain rotation in a polycrystalline solid.

Table 1
Rotation parameters at the instants noted in Fig. 3[link]

  Rotation axis  
Frame x y z Rotation angle (°)
1–2 −0.5551 0.4298 0.7121 0.1149
2–3 −0.7023 0.7039 −0.1065 0.0744
3–4 −0.7398 0.5853 0.3319 0.0988
4–5 −0.2884 0.8321 0.4737 0.0626
5–6 −0.4736 0.7842 0.4010 0.0689
6–7 0.0135 0.8960 0.4439 0.0523
7–8 −0.8573 0.3840 0.3429 0.1262

As a self-consistency check, we conduct forward calculations using the rotation parameters in Table 1[link] to compute the diffraction spot positions on the diffraction detector. Given the rotation axes and angles as well as the initial positions of diffraction spots, their positions at any instant can be obtained as elaborated in Appendix A[link]. The simulated diffraction spot positions are listed in Table 2[link] and plotted as empty circles in columns 1 and 3 of Fig. 3[link]. The comparison shows a complete coincidence between the simulations and the measurement.

Table 2
Measured and simulated positions of the two diffraction spots on the detector

V and H denote vertical and horizontal positions in pixels, respectively. Subscripts 1 and 2 denote Laue diffraction spots 1 and 2, respectively.

  Experimental Simulated
Frame V1 H1 V2 H2 V1 H1 V2 H2
1 86 40 341 353 86 40 341 353
2 81 39 335 348 80 38 336 349
3 73 40 330 348 73 39 331 349
4 66 40 324 346 65 39 325 347
5 59 39 319 345 59 39 319 345
6 53 39 313 343 52 39 313 344
7 46 39 308 343 46 39 308 342
8 39 39 303 340 39 39 304 340

The two-dimensional translation and full rotation parameters as a function of time (Table 1[link]) are illustrated in Fig. 4[link], using a schematic object. Therefore, rigid-body motion, including translation and three-dimensional rotation, of the sample can be strictly mapped through simultaneous two-dimensional X-ray imaging and diffraction. Translation mapping can be extended into three dimensions if an extra camera is used at a different view angle. Simultaneous high-speed tracking of translation and rotation in three dimensions offers great potential in a multitude of applications.

[Figure 4]
Figure 4
Illustrations of the evolution of sample position and orientation, corresponding to Figs. 2[link] and 3[link]. The rotation angles are magnified by ten for clarity. Red arrows: rotation axis; blue elliptic arrows: rotation direction.

We present the simplest case in the above discussion. For strong shocks, both lattice strain and grain rotation contribute to the movement of a diffraction spot. Strain can be evaluated from conventional methods such as known equation of state, and it is then decoupled from rotation to a certain extent. However, more diffraction spots are certainly desirable.

4. Conclusion

We have performed high-speed simultaneous synchrotron X-ray Laue diffraction and phase-contrast imaging measurements on single-crystal Si under SHPB loading, and developed a methodology for quantifying rotation parameters, i.e. rotation axis and rotation angles. The exposure time is 250–350 ns and the frame interval is 10 µs. High-speed crystal rotation at about 300 r s−1 is resolved. High-speed motion of crystals, including two-dimensional translation and three-dimensional rotation, can be tracked in real time from simultaneous imaging and diffraction.

APPENDIX A

Determining rotation parameters from two Laue diffraction spots

The geometry of diffraction is defined in Figs. 5[link](a) and 5(b). Incident X-rays ([O^{\,\prime}B]) are scattered by lattice planes ([O^{\,\prime}C]) and the scattered X-rays ([O^{\,\prime}A^{{\prime}}]) yield a diffraction spot ([A^{{\prime}}]) on a diffraction detector ([OA^{{\prime}}]) located at an arbitrary position. The origin, O, is the transmission spot on the detector plane of the incident X-ray beam propagating along the −z direction. For convenience, we introduce the virtual detector plane containing O, i.e. the Oxy plane, perpendicular to the incident X-ray beam. A diffraction spot ([A^{{\prime}}]) on the detector is mapped onto the virtual detector plane or vice versa. Discussions refer to the virtual detector plane ([OA^{{\prime}}]) unless stated otherwise.

[Figure 5]
Figure 5
(a) Diffraction geometry. (b) Definitions of lattice plane normal n, and angles of γ and η. Symbols are the same as those in Fig. 1.[link]

The diffraction angle [\theta = \angle\,BO^{\,\prime}C] and [2\theta = \angle\,A^{{\prime}}O^{\,\prime}O]. The azimuthal angle of a diffraction spot, γ, is the angle formed between [A^{{\prime}}O] and the x axis. The distance between the sample and the origin is [d^{\,{\prime}} = OO^{\,\prime}]. In terms of [d^{\,{\prime}}] and the mapped coordinates (x,y) of a diffraction spot:

[2\theta = \arctan\Big[{\big({x^{2}+y^{2}\big)^{1/2}}/ {d^{\,\prime}}}\Big], \eqno(1)]

[\gamma = \arctan({y}/{x}).\eqno(2)]

The unit vector n is the normal of certain parallel lattice planes. The angle between n and [OA^{{\prime}}] is defined as η, and [\eta = \theta]. Therefore,

[n_{x} = \cos\eta\cos\gamma = \cos\theta\cos\gamma,\eqno(3)]

[n_{y} = \cos\eta\sin\gamma = \cos\theta\sin\gamma,\eqno(4)]

[n_{z} = \sin\eta = \sin\theta.\eqno(5)]

n is expressed in terms of γ and θ obtained from two-dimensional diffraction patterns, so the instantaneous orientation of the lattice planes is determined.

The rotation axis and rotation angle are chosen to describe sample rotation, and are denoted by a unit vector [{\bf u} = (u_{x},u_{y},u_{z})^{{\rm T}}] and α, respectively. Here, superscript T denotes the transpose. The corresponding rotation matrix [{\bb R}] can be expressed in terms of u and α as

[{\left[\matrix{\matrix{\cos\alpha \,\,+\cr u_{x}^{2}(1-\cos\alpha_{\vphantom{\Big|}})} &\matrix {u_{x}u_{y}(1-\cos\alpha)\cr \,-\, u_{z}\sin\alpha_{\vphantom{\Big|}}} &\matrix{u_{x}u_{z}(1-\cos\alpha)\cr \,+ \, u_{y}\sin\alpha_{\vphantom{\Big|}}} \cr \cr \matrix{u_{y}u_{x}(1-\cos\alpha)\cr \,+\, u_{z}\sin\alpha_{\vphantom{\Big|}}} &\matrix{\cos\alpha \,\, +\cr u_{y}^{2}(1-\cos\alpha_{\vphantom{\Big|}})} &\matrix{u_{y}u_{z}(1-\cos\alpha)\cr \,-\, u_{x}\sin\alpha_{\vphantom{\Big|}}} \cr \cr \matrix{u_{z}u_{x}(1-\cos\alpha)\cr \,-\, u_{y}\sin\alpha} &\matrix {u_{z}u_{y}(1-\cos\alpha)\cr \,+\, u_{x}\sin\alpha } &\matrix{\cos\alpha \,\, +\cr u_{z}^{2}(1-\cos\alpha)}}\right]} .\eqno(6)]

After sample rotation, instantaneous normal vector [{\bf n}_{{\rm 1}} = (n_{{x_1}},n_{{y_1}},n_{{z_1}})^{{\rm T}}] at time t1 changes to [{\bf n}_{{\rm 2}}] = (nx2,ny2,nz2)T at t2. The rotation matrix, [{\bb R}], relates [{\bf n}_{1}] and [{\bf n}_{2}] via

[{\bf n}_{{\rm 2}} = {\bb R}{\bf n}_{1}. \eqno(7)]

The objective is to solve [{\bf u}] and α from equation (7)[link], which fully characterizes the rigid-body rotation of the sample. This equation can be solved numerically. For small rotation angles, it can be simplified with [\cos\alpha\approx 1] and [\sin\alpha\approx\alpha]. Equation (7)[link] is then reduced to

[\pmatrix{n_{{x_{2}}} \cr n_{{y_{2}}} \cr n_{{z_{2}}}} = \pmatrix{1 &-u_{z}\alpha&u_{y}\alpha \cr u_{z}\alpha &1 &-u_{x}\alpha\cr -u_{y}\alpha&u_{x}\alpha&1 } \pmatrix{n_{{x_{1}}} \cr n_{{y_{1}}} \cr n_{{z_{1}}}}. \eqno(8)]

Subtracting [{\bf n}_{{\rm 1}}] on both sides of equation (8)[link] yields

[\pmatrix{\Delta n_{x}\cr \Delta n_{y}\cr \Delta n_{z} } = \alpha \pmatrix{ 0&-u_{z}&u_{y}\cr u_{z}&0&-u_{x}\cr -u_{y}&u_{x}&0} \pmatrix{n_{{x_{1}}}\cr n_{{y_{1}}}\cr n_{{z_{1}}} }, \eqno(9)]

where [\Delta n_{{x}} = n_{{x_{2}}}-n_{{x_{1}}}], [\Delta n_{y} = n_{{y_{2}}}-n_{{y_{1}}}] and [\Delta n_{z} = n_{{z_{2}}}-n_{{z_{1}}}].

For convenience, the Euler vector E is defined as

[{\bf E} = (E_{x},E_{y},E_{z})^{{\rm T}}\,\, \doteq \,\, \alpha(u_{x},u_{y},u_{z})^{{\rm T}}. \eqno(10)]

It follows that

[\eqalign{\Delta n_{x}& = -n_{{y_{1}}}E_{z}+n_{{z_{1}}}E_{y},\cr \Delta n_{y}& = -n_{{z_{1}}}E_{x}+n_{{x_{1}}}E_{z},\cr \Delta n_{z}& = -n_{{x_{1}}}E_{y}+n_{{y_{1}}}E_{x}.} \eqno(11)]

There are three unknowns, Ex, Ey, Ez. But these three equations are linearly dependent or contradictory. Thus, in order to determine [{\bf E}], another diffraction spot is needed. Given two diffraction spots from different sets of lattice planes, this problem becomes overdetermined (six equations versus three unknowns) and the rotation parameters can be fully solved as

[\alpha = \Big({E_{x}^{\,2}+E_{y}^{\,2}+E_{z}^{\,2}}\Big)^{1/2},\eqno(12)]

[{\bf u} = {{\bf E} /\alpha}.\eqno(13)]

Note that even more diffraction spots can be used and the uncertainties are reduced.

For each diffraction spot, a set of [(\theta,\gamma)] and thus a set of (nx,ny,nz) are obtained. Two independent spots yield two sets of (nx,ny,nz). The rotation parameters, α and u, are solved from four diffraction spots for frames 1 and 2. Repeating such a process between two adjacent frames yields the rotation history of the sample. A MATLAB code is developed to perform such analysis in a batch mode.

To validate the above derivations, we perform forward analysis, i.e. calculate the diffraction spots on the detector from rotation parameters α and u obtained from the diffraction spots of frames 1 and 2 as described above. [{\bf n}_{1}] is calculated from the initial positions of a diffraction spot at frame 1. After rotation, [{\bf n}_{1}] changes to [{\bf n}_{2}], which is calculated with the Rodrigues' rotation formula (Rodrigues, 1840[Rodrigues, O. (1840). J. Math. Pures Appl. 5, 380-440.])

[{\bf n}_{2} = {\bf n}_{1}\cos\alpha+({ \bf u}\times {\bf n}_{1})\sin\alpha+{\bf u}({\bf u}\cdot{\bf n}_{1})(1-\cos\alpha). \eqno(14)]

Then,

[\gamma_{2} = \arctan\Big({ n_{y_{2}}/ {n_{x_{2}}} }\Big), \eqno(15)]

[\eta _{2} = \arcsin { {n_{{z_{2}}}}\Big({n_{{x_{2}}}^{2}+n_{{y_{2}}}^{2}}}\Big)^{-1/2}.\eqno(16)]

On the virtual detector plane,

[x_{2} = d^{\,{\prime}}\tan 2\theta _{2}\cos\gamma _{2} = d^{\,{\prime}}\tan 2\eta _{2}\cos\gamma _{2}, \eqno(17)]

[y_{2} = d^{\,{\prime}}\tan 2\theta _{2}\sin\gamma _{2} = d^{\,{\prime}}\tan 2\eta _{2}\sin\gamma _{2}.\eqno(18)]

The simulated diffraction spot (x2,y2) is then mapped onto the detector and should coincide with the directly measured spot at frame 2.

Acknowledgements

The MATLAB program developed in this work has benefited from the HisPoD diffraction simulation code. We are grateful to J. Wang at APS, B. X. Bie, D. Fan and L. Lu at PIMS for their help with the experiments. This work was partially supported by the 973 project (2014CB845904) and NSF (11472253) of China. Use of the Advanced Photon Source, an Office of Science User Facility operated for the US Department of Energy (DOE) Office of Science by Argonne National Laboratory, was supported by the US DOE under contract No. DE-AC02-06CH11357.

References

First citationBarabash, O. M., Santella, M., Barabash, R. I., Ice, G. E. & Tischler, J. (2010). JOM, 62, 29–34.  Web of Science CrossRef CAS Google Scholar
First citationBorbély, A. (2015). J. Appl. Cryst. 48, 1614–1616.  Google Scholar
First citationChen, P., Mao, S. C., Liu, Y., Wang, F., Zhang, Y. F., Zhang, Z. & Han, X. D. (2013). Mater. Sci. Eng. A, 580, 114–124.  Web of Science CrossRef CAS Google Scholar
First citationFan, D., Lu, L., Li, B., Qi, M. L., Zhao, F., Sun, T., Fezzaa, K., Chen, W. & Luo, S. N. (2014). Rev. Sci. Instrum. 85, 113902.  Web of Science CrossRef PubMed Google Scholar
First citationFowles, G. R. (1970). Rev. Sci. Instrum. 41, 984–996.  CrossRef Web of Science Google Scholar
First citationGoldschmidt, M. J. V., Beetstra, R. & Kuipers, J. A. M. (2004). Powder Technol. 142, 23–47.  Web of Science CrossRef CAS Google Scholar
First citationHall, S. A., Wright, J., Pirling, T., Andò, E., Hughes, D. J. & Viggiani, G. (2011). Granul. Matter, 13, 251–254.  Web of Science CrossRef Google Scholar
First citationHuang, X. R. (2010). J. Appl. Cryst. 43, 926–928.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationHudspeth, M., Claus, B., Dubelman, S., Black, J., Mondal, A., Parab, N., Funnell, C., Hai, F., Qi, M. L., Fezzaa, K., Luo, S. N. & Chen, W. (2013). Rev. Sci. Instrum. 84, 025102.  Web of Science CrossRef PubMed Google Scholar
First citationHudspeth, M., Sun, T., Parab, N., Guo, Z., Fezzaa, K., Luo, S. & Chen, W. (2015). J. Synchrotron Rad. 22, 49–58.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationKajishima, T. (2004). Int. J. Heat Fluid Flow, 25, 721–728.  Web of Science CrossRef Google Scholar
First citationKolsky, H. (1949). Proc. Phys. Soc. B, 62, 676–700.  CrossRef Google Scholar
First citationLambert, P. K., Hustedt, C. J., Vecchio, K. S., Huskins, E. L., Casem, D. T., Gruner, S. M., Tate, M. W., Philipp, H. T., Woll, A. R., Purohit, P., Weiss, J. T., Kannan, V., Ramesh, K. T., Kenesei, P., Okasinski, J. S., Almer, J., Zhao, M., Ananiadis, A. G. & Hufnagel, T. C. (2014). Rev. Sci. Instrum. 85, 093901.  Web of Science CrossRef PubMed Google Scholar
First citationLi, Y., Qian, D., Xue, J., Wan, J., Zhang, A., Tamura, N., Song, Z. & Chen, K. (2015). Appl. Phys. Lett. 107, 181902.  Web of Science CrossRef Google Scholar
First citationLiu, X., Thadesar, P. A., Taylor, C. L., Oh, H., Kunz, M., Tamura, N., Bakir, M. S. & Sitaraman, S. K. (2014). Appl. Phys. Lett. 105, 112109.  Web of Science CrossRef Google Scholar
First citationLuo, S. N., Jensen, B. J., Hooks, D. E., Fezzaa, K., Ramos, K. J., Yeager, J. D., Kwiatkowski, K. & Shimada, T. (2012). Rev. Sci. Instrum. 83, 073903.  Web of Science CrossRef PubMed Google Scholar
First citationMaaß, R., Van Petegem, S., Grolimund, D., Van Swygenhoven, H., Kiener, D. & Dehm, G. (2008). Appl. Phys. Lett. 92, 071905.  Google Scholar
First citationMaaß, R., Van Petegem, S., Van Swygenhoven, H., Derlet, P. M., Volkert, C. A. & Grolimund, D. (2007). Phys. Rev. Lett. 99, 145505.  Web of Science PubMed Google Scholar
First citationMicha, J. S. & Robach, O. (2014). Lauetools, https://sourceforge.net/projects/lauetools/Google Scholar
First citationNoyan, I. C. & Cohen, J. B. (2013). Residual Stress: Measurement by Diffraction and Interpretation. Berlin: Springer-Verlag.  Google Scholar
First citationRodrigues, O. (1840). J. Math. Pures Appl. 5, 380–440.  Google Scholar
First citationSun, J. & Battaglia, F. (2006). Chem. Eng. Sci. 61, 1470–1479.  Web of Science CrossRef CAS Google Scholar
First citationTamura, N. (2014). Strain and Dislocation Gradients from Diffraction: Spatially Resolved Local Structure and Defects, edited by R. Barabash and G. Ice. Imperial College Press.  Google Scholar
First citationTurneaure, S. J., Gupta, Y., Zimmerman, K., Perkins, K., Yoo, C. & Shen, G. (2009). J. Appl. Phys. 105, 053520.  Web of Science CrossRef Google Scholar
First citationWang, L., Barabash, R. I., Yang, Y., Bieler, T. R., Crimp, M. A., Eisenlohr, P., Liu, W. & Ice, G. E. (2011). Metall. Mater. Trans. A, 42, 626–635.  Web of Science CrossRef CAS Google Scholar
First citationWei, Q., Jiao, T., Ramesh, K. & Ma, E. (2004). Scr. Mater. 50, 359–364.  Web of Science CrossRef CAS Google Scholar
First citationWu, X., Wang, Q., Luo, Z., Fang, M. & Cen, K. (2008). Powder Technol. 181, 21–30.  Web of Science CrossRef CAS Google Scholar
First citationZaafarani, N., Raabe, D., Roters, F. & Zaefferer, S. (2008). Acta Mater. 56, 31–42.  Web of Science CrossRef CAS 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 logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775
Follow J. Synchrotron Rad.
Sign up for e-alerts
Follow J. Synchrotron Rad. on Twitter
Follow us on facebook
Sign up for RSS feeds