research papers
Dynamic crystal rotation resolved by high-speed synchrotron X-ray Laue diffraction
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
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.
Keywords: synchrotron X-ray Laue diffraction; crystal rotation.
1. Introduction
Grain or crystal rotation is a key mechanism of accommodating elastic or plastic deformation in crystalline solids (Chen et al., 2013). In quasi-static experiments, crystal rotation is normally resolved with electron backscattering diffraction (Zaafarani et al., 2008) or X-ray Laue diffraction (Maaß et al., 2007, 2008). 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) or a gas gun (Fowles et al., 1970). 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; Goldschmidt et al., 2004; Sun & Battaglia, 2006). Detection of high-speed particle rotation is in urgent need. For example, Wu et al. (2008) used optical imaging to measure particle rotation at about 500 revolutions per second (r s−1) in a circulating 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). 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; Hudspeth et al., 2013; Lambert et al., 2014; Turneaure et al., 2009; Wei et al., 2004) 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; Hudspeth et al., 2015). 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; Hudspeth et al., 2015). The schematic setup for simultaneous diffraction and phase contrast imaging is shown in Fig. 1, 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 [01], [2] and [111] are oriented approximately along the x, y and z directions, respectively. Dimensions of the sample perpendicular to the X-ray direction () 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).
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). The probe X-rays () 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 (). 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 (). 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 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; Liu et al., 2014; Wang et al., 2011; Barabash et al., 2010) and grain rotation (Maaß et al., 2007, 2008); there are many dedicated software packages available for this analysis (Tamura, 2014; Micha & Robach, 2014; Huang, 2010). A method based on digital image correlation using Laue spots as speckles has been developed recently for strain analysis (Borbély, 2015). 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. 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 and 3.
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). 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, 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) cannot induce such a pronounced shift. Note that there is only one diffraction spot recorded prior to frame 1 (Fig. 3). 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.
As each harmonic of the undulator radiation has a finite bandwidth (∼2% for the APS undulator A light source) (Fan et al., 2014), 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), so the small-angle assumption in equation (8) of Appendix A is reasonable.
Each pixel on the diffraction detector is mapped onto the –γ plane. Given the finite size of a diffraction spot, we use Gaussian fitting to obtain the coordinates of its center (, γ). The diffraction spots are plotted on the –γ plane in columns 2 and 4 of Fig. 3. Then, the relative rotation axes (expressed in three projection quantities) and the rotation angles between adjacent frames are calculated (Table 1). 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.
|
As a self-consistency check, we conduct forward calculations using the rotation parameters in Table 1 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. The simulated diffraction spot positions are listed in Table 2 and plotted as empty circles in columns 1 and 3 of Fig. 3. The comparison shows a complete coincidence between the simulations and the measurement.
|
The two-dimensional translation and full rotation parameters as a function of time (Table 1) are illustrated in Fig. 4, 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.
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(a) and 5(b). Incident X-rays () are scattered by lattice planes () and the scattered X-rays () yield a diffraction spot () on a diffraction detector () 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 () on the detector is mapped onto the virtual detector plane or vice versa. Discussions refer to the virtual detector plane () unless stated otherwise.
The diffraction angle and . The azimuthal angle of a diffraction spot, γ, is the angle formed between and the x axis. The distance between the sample and the origin is . In terms of and the mapped coordinates (x,y) of a diffraction spot:
The unit vector n is the normal of certain parallel lattice planes. The angle between n and is defined as η, and . Therefore,
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 and α, respectively. Here, superscript T denotes the transpose. The corresponding rotation matrix can be expressed in terms of u and α as
After sample rotation, instantaneous normal vector at time t1 changes to = (nx2,ny2,nz2)T at t2. The rotation matrix, , relates and via
The objective is to solve and α from equation (7), 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 and . Equation (7) is then reduced to
Subtracting on both sides of equation (8) yields
where , and .
For convenience, the Euler vector E is defined as
It follows that
There are three unknowns, Ex, Ey, Ez. But these three equations are linearly dependent or contradictory. Thus, in order to determine , 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
Note that even more diffraction spots can be used and the uncertainties are reduced.
For each diffraction spot, a set of 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. is calculated from the initial positions of a diffraction spot at frame 1. After rotation, changes to , which is calculated with the Rodrigues' rotation formula (Rodrigues, 1840)
Then,
On the virtual detector plane,
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
Barabash, O. M., Santella, M., Barabash, R. I., Ice, G. E. & Tischler, J. (2010). JOM, 62, 29–34. Web of Science CrossRef CAS Google Scholar
Borbély, A. (2015). J. Appl. Cryst. 48, 1614–1616. Google Scholar
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. Web of Science CrossRef CAS Google Scholar
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. Web of Science CrossRef PubMed Google Scholar
Fowles, G. R. (1970). Rev. Sci. Instrum. 41, 984–996. CrossRef Web of Science Google Scholar
Goldschmidt, M. J. V., Beetstra, R. & Kuipers, J. A. M. (2004). Powder Technol. 142, 23–47. Web of Science CrossRef CAS Google Scholar
Hall, 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
Huang, X. R. (2010). J. Appl. Cryst. 43, 926–928. Web of Science CrossRef CAS IUCr Journals Google Scholar
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. Web of Science CrossRef PubMed Google Scholar
Hudspeth, 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
Kajishima, T. (2004). Int. J. Heat Fluid Flow, 25, 721–728. Web of Science CrossRef Google Scholar
Kolsky, H. (1949). Proc. Phys. Soc. B, 62, 676–700. CrossRef Google Scholar
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. Web of Science CrossRef PubMed Google Scholar
Li, 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
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. Web of Science CrossRef Google Scholar
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. Web of Science CrossRef PubMed Google Scholar
Maaß, R., Van Petegem, S., Grolimund, D., Van Swygenhoven, H., Kiener, D. & Dehm, G. (2008). Appl. Phys. Lett. 92, 071905. Google Scholar
Maaß, 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
Micha, J. S. & Robach, O. (2014). Lauetools, https://sourceforge.net/projects/lauetools/. Google Scholar
Noyan, I. C. & Cohen, J. B. (2013). Residual Stress: Measurement by Diffraction and Interpretation. Berlin: Springer-Verlag. Google Scholar
Rodrigues, O. (1840). J. Math. Pures Appl. 5, 380–440. Google Scholar
Sun, J. & Battaglia, F. (2006). Chem. Eng. Sci. 61, 1470–1479. Web of Science CrossRef CAS Google Scholar
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. Google Scholar
Turneaure, S. J., Gupta, Y., Zimmerman, K., Perkins, K., Yoo, C. & Shen, G. (2009). J. Appl. Phys. 105, 053520. Web of Science CrossRef Google Scholar
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. Web of Science CrossRef CAS Google Scholar
Wei, Q., Jiao, T., Ramesh, K. & Ma, E. (2004). Scr. Mater. 50, 359–364. Web of Science CrossRef CAS Google Scholar
Wu, X., Wang, Q., Luo, Z., Fang, M. & Cen, K. (2008). Powder Technol. 181, 21–30. Web of Science CrossRef CAS Google Scholar
Zaafarani, 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.