# The effect of crystal defects on 3D high resolution diffraction peaks: a FFT-based method

Komlavi Eloh, Alain Jacques, Gabor Ribarik, Stéphane Berbenni

► To cite this version:

Komlavi Eloh, Alain Jacques, Gabor Ribarik, Stéphane Berbenni. The effect of crystal defects on 3D high resolution diffraction peaks: a FFT-based method. Materials, 2018, 11 (9), pp.1669. 10.3390/ma11091669. [https://doi.org/10.3390/ma11091669](https://doi.org/10.3390/ma11091669)

HAL Id: hal-02360465

https://hal.science/hal-02360465

Submitted on 12 Nov 2019

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# The effect of crystal defects on 3D high resolution diffraction peaks: a FFT-based method

K. Eloh<fnref n="1" />,<fnref n="2" />,<fnref n="3" />, A. Jacques<fnref n="1" />,<fnref n="3" />,*, G. Ribarik<fnref n="3" />,<fnref n="4" />, S. Berbenni<fnref n="2" />,<fnref n="3" />

**Abstract:** Forward modeling of diffraction peaks is a potential way to compare the results of theoretical mechanical simulations and experimental X-Ray Diffraction data recorded during in situ experiments. As the input data are the strain or displacement field within a representative volume of the material containing dislocations, a computer-aided efficient and accurate method to generate these fields is necessary. With this aim, a current and promising numerical method is based on the use of the Fast Fourier Transform (FFT) method. However, classic FFT-based methods present some numerical artifacts due to the Gibbs phenomenon or "aliasing" and to "voxelization" effects. Here, we propose several improvements: first, a consistent discrete Green operator to remove "aliasing" effects and second, a method to minimize the voxelization artifacts generated by dislocation loops inclined with respect to the computational grid. Then we then show the effect of these improvements on theoretical diffraction peaks.

**Keywords:** Dislocations; diffraction; FFT-based method; Discrete Green operator; voxelization artifacts; sub-voxel method; simulated diffraction peaks; scattered Intensity

## 1. Introduction

X-Ray diffraction is one of the most powerful non-destructive tools to investigate materials, as their wavelength is commensurate with the distance between atoms within a crystal [1–9]. Successive improvements of both the X-Ray sources (from X-Ray tubes to third generation synchrotrons) and detectors (from photographic plates and gas counters to fast two-dimensional arrays) have led to a tremendous increase in the quantity of data recorded per unit time, allowing real time *in situ* or *in operando* measurements [10,11]. It is now possible to determine the 3D grain microstructure of a bulk material with a submicron resolution (using Topo-Tomography), to follow the evolution of the elastic strain state of the grains of a polycrystal during mechanical tests (3D-XRD, far field diffractometry), or to measure the distribution of strains within a few grains in real time (2D diffractometry) [12,13]. Such experiments result in terabytes of data recorded within a few days, which need to be efficiently analyzed. In fact, only a low fraction of those data is actually treated because scientists lack both time and software for further analysis [14].

The classical techniques used to analyze the 1D or 2D diffraction patterns recorded during tests performed on polycrystalline specimens such as the Rietveld method, the square sines method to measure internal stresses, or CMWP fitting for dislocations content often rely on simplified and mathematically tractable models of a microstructure. Calculations which may involve simplifying hypothesis lead to a general formula which can be used to fit one or several parameters of the

<footnote n="1">¹Université de Lorraine, CNRS, IJL, F-54000 Nancy, France; komlavi-senyo.eloh@univ-lorraine.fr</footnote>

<footnote n="2">²Université de Lorraine, CNRS, Arts et Métiers Paris Tech, LEM3, F-57000 Metz, France; stephane.berbenni@univ-lorraine.fr</footnote>

<footnote n="3">³Laboratory of Excellence on Design of Alloy Metals for low-mAss Structures (DAMAS), Université de Lorraine, France.</footnote>

<footnote n="4">⁴Department of Materials Physics, Eötvös University, Budapest, Hungary.</footnote>

<footnote>* Correspondence: alain.jacques@univ-lorraine.fr; Tel.: +03 72 74 26 79</footnote>

<footnote>Received: date; Accepted: date; Published: date</footnote>
45

46

microstructure (dislocation densities and type, internal stress tensor...) to the diffraction pattern (peak profiles, variation of the 2θ_B angle with orientation...)

47

48

49

50

51

52

53

54

55

56

57

58

During the last ten years, several authors proposed the opposite approach: forward modeling [15–22]. This requires the design of a microstructure and the simulation of its behavior (often under process or thermo-mechanical solicitation), the computation of the elastic strain field or the displacement field. The last step is the generation through a ‘virtual diffractometer’ of a theoretical diffraction pattern (different G vectors and different orientations of the lattice planes), which can be compared with the experimental one. Depending on the size of the simulated representative volume of matter and the experimental conditions such as the X-Ray beam coherence, different assumptions can be done such as a coherent beam (where the amplitudes scattered by different points add) or an incoherent beam (scattered intensities add), or for a partially coherent beam where a full calculation may be necessary. Such modeling can be quite successful and can be used to validate the different steps involved, mainly the microstructure and the constitutive law used to simulate the material’s behavior.

59

60

61

62

63

64

However, as diffraction peaks contain information on different scales of a specimen: from average quantities such as Type I (average) stresses related to the peaks' positions, Type II (at grain level) stresses related to its width, and Type III stresses (near the core of defects such as dislocations) related to the peaks' tails, a realistic simulation of a diffraction peak requires a description of a material's Representative Volume Element with a very fine mesh, i.e. a huge amount of CPU time with classical methods used for simulations such as the Finite Element Method.

65

66

67

68

69

70

71

72

Numerical approaches based on the FFT for calculating the stress and strain fields within a composite material received a surge of interest since the pioneering work of Moulinec and Suquet [23,24]. They were first developed to compute effective properties and mechanical field of linear elastic composites[23–26] and were extended to heterogeneous materials with eigenstrains (dislocations, thermal strains...) [14,27–30]. They are also used for conductivity problems [31], non-linear materials [25,27], viscoplastic or elastoviscoplastic polycrystals [32–36]. Today, FFT based approaches represent alternative to finite element method because they are rather attractive in terms of computation time [32].

73

74

75

76

77

78

79

However initial tests indicate that the displacement field computed (essential for diffraction pattern generation) with FFT algorithms presents some numerical artifacts. These numerical artifacts are due to Gibbs phenomenon or “aliasing” and to voxelization. The accuracy of the calculated strain or displacement field is strongly influenced by these shortcomings and the simulated peaks may provide wrong information on mechanical behavior or material characteristics. Therefore, it is important to control these artifacts to simulate correct diffraction pattern in the case of a microstructure containing different phases, grains, and crystal defects.

80

81

82

83

84

85

86

87

88

89

The aim of this paper is to improve the accuracy of the displacement field for diffraction peak generation. This improvement is based on the introduction of a consistent discrete periodized Green operator associated with the displacement field in order to take explicitly into account the discreteness of the discrete Fourier Transform method [37]. The improvement of the voxelization in FFT-method is performed through a sub-voxelization method will be described for inclined dislocation loops. These improvements are reported and discussed. In the section 2, the FFT-based method to compute the displacement field in a periodic medium is described. In section 3, the treatment of voxelization problems in FFT-based approaches by a sub-voxelization method is detailed in the case of slip plane not conforming to FFT grid. In the section 4, simulation of diffraction peaks is reported and discussed.

90

## 2. Displacement field

91

### 2.1. FFT-based algorithm and mechanical fields

92

93

94

Let us consider a homogeneous elastic medium with eigenstrain assuming a periodic unit cell discretized in $N \times N \times N$ voxels and subjected to an uniform overall strain tensor denoted **E**. Here, this overall strain is the spatial average of the strain field in the unit cell (with external loading and a
given eigenstrain field). The unit cell contains dislocations (line defects) which are modeled with an eigenstrain tensor. These dislocations create a displacement field during their motion along various slip systems and thus generates strain/stress field [38].

The displacement vector is denoted **u** and in the forthcoming equations, **x** denotes any position vector within the unit cell. All vectorial and tensorial fields will be written using bold characters.

Starting from the equation for mechanical equilibrium, $div \sigma(\mathbf{x}) = 0$, and using elastic equations, the displacement field is given at every positions by [39]:

$$ \mathbf{u}(\mathbf{x}) = (\mathbf{B} * \mathbf{c}^0: \boldsymbol{\varepsilon}^*)(\mathbf{x}) \quad (1) $$

Where the symbol * denotes the spatial convolution product, $\mathbf{c}^0$ is the homogeneous linear elastic stiffness, $\boldsymbol{\varepsilon}^*$ is the eigenstrain field and $\mathbf{B}$ is a third order Green operator defined in Fourier space as:

$$ \hat{\mathbf{B}}_{ijk}(\xi) = \frac{i}{2} (\hat{G}_{ij} \xi_k + \hat{G}_{ik} \xi_j) \quad (2) $$

in which $\hat{\mathbf{B}}$ is the Fourier transform of $\mathbf{B}$ and $\hat{\mathbf{G}}$ is the Fourier transform of the elastic Green tensor [39]. Therefore, using the Fourier transform of spatial convolution product, Equation (1) can be written in Fourier space as:

$$ \hat{\mathbf{u}}(\xi) = \hat{\mathbf{B}}(\xi):\mathbf{c}^0: \boldsymbol{\varepsilon}^*(\xi) \quad (3) $$

Several numerical results showed that the use of the third order operator $\hat{\mathbf{B}}$ derived from the classic Green $\hat{\mathbf{G}}$ leads to spurious oscillations on the computed displacement field near materials discontinuities. The Discrete Fourier transform (DFT) used in this algorithm indeed transforms a periodic function in real space into a periodic function in reciprocal space. However, the operator $\hat{\mathbf{B}}$ commonly used is the continuous analytic operator truncated to the size of the unit cell of the reciprocal space: it is not periodic function. To fix this problem, we a periodized consistent discrete Green operator using the DFT. The Fourier transform of this new operator $\hat{\mathbf{B}}'$ is written as function of $\hat{\mathbf{B}}$ (the mathematical details about its derivation are given in [37]) and reads:

$$ \hat{\mathbf{B}}'(\xi_{ijk}) = A_{ijk} \sum_{m,n,p=-\infty}^{+\infty} \frac{(-1)^{m+n+p}}{(mN+i)} \frac{1}{(nN+j)} \frac{1}{(pN+k)} \hat{\mathbf{B}}(\xi_{mN+i,nN+j,pN+k}) $$

$$ \text{With} \quad A_{ijk} = \left(\frac{N}{\pi}\right)^3 \sin\left(\frac{i\pi}{N}\right) \sin\left(\frac{j\pi}{N}\right) \sin\left(\frac{k\pi}{N}\right) \quad (4) $$

Discrete frequency appearing in this equation is given when N is even by (T is the unit cell period):

$$ \xi = \left(-\frac{N}{2} + 1\right) \frac{1}{T}, \left(-\frac{N}{2} + 2\right) \frac{1}{T}, \dots, \frac{1}{T}, 0, \frac{1}{T}, \dots, \left(\frac{N}{2} - 1\right) \frac{1}{T}, \left(\frac{N}{2}\right) \frac{1}{T} $$

Here the sum on the $\hat{\mathbf{B}}$ operator is extended to the whole reciprocal space (in practice for m, n, p up to a few tens) and folded up onto the unit cell of the DFT with suitable coefficients. The inverse transform of $\hat{\mathbf{u}}(\xi)$ gives the displacement field at the center of each voxel.

We can also compute the displacement field at each voxel's corner with a shifted operator using the shift theorem:

$$ \hat{\mathbf{B}}''(\xi_{ijk}) = A_{ijk} e^{\frac{p\pi i+j+k}{N}} \sum_{m,n,p=-\infty}^{+\infty} \frac{1}{(mN+i)(nN+j)(pN+k)} \hat{\mathbf{B}}(\xi_{mN+i,nN+j,pN+k}) \quad (5) $$

## 2.2. Numerical examples

Let us consider a homogeneous material with isotropic elastic constant: Young's modulus $E = 333.4$ GPa and the Poisson ratio $\vartheta = 0.26$. This approximately corresponds to the room temperature elastic constants of single crystalline Ni-based Superalloys. The unit cell (Figure 1a) is discretized in 128 × 128 × 128 voxels and contains a square-shaped inclusion discretized in
32 × 32 × 1 voxels corresponding to an Eshelby-like square prismatic loop perpendicular to the z-axis. In order to generate a shift of the upper surface of the inclusion relative to its lower surface by a Burgers vector $b(0, 0, b_3)$ the voxels within the inclusion are submitted to an eigenstrain: $\epsilon_{ij}^* = 0$ except $\epsilon_{33}^* = 1$. We thus have $b_3 = t \times \epsilon_{33}^*$ where t the thickness of the inclusion in the z-direction (i.e. the voxel size). This displacement field computed with the FFT algorithm using the different periodized discrete Green operators is represented along z-axis in figure 1.

<figure><img src="image_3.png" /><figcaption>**Figure 1:** (a) Simulation of a squared dislocation in plane (001) by a platelet with eigenstrain; (b) Component $u_3$ of the displacement field (normalized by $b_3$) along the z axis (arrow) computed with the Green operator **B** and showing oscillations; (c) Same component correct $u_3$ computed with **B'**. The displacement at voxel (64,64,64) is zero in the center of the inclusion (c) and $b_3/2$ on its surface (c). (d) Same component correct $u_3$ computed with **B''**.</figcaption></figure>

When computed along a line crossing a dislocation loop, the displacement field exhibits a discontinuity with a jump equal to Burgers vector **b**. This is indeed observed in figure 1. However, the displacement field computed with the usual Green operator **B** (figure 1b) also shows oscillations, while which are not observed with the periodized operators **B'** and **B''**. An artificial damping of the oscillations in figure 1b (such as a low pass filtering) might smooth these oscillations, but it would also smooth the discontinuity.

## 2.3. Voxelization effect on the displacement field.

While the displacement field computed for dislocation loops having their planes parallel to the faces of the simulated volume, artifacts appear for inclined loops, as shown in figure 2 for a dislocation loop with a [011] Burgers vector lying in a (111) slip plane of a fcc crystal. The eigenstrain tensor is constrained in the region occupied by the dislocation loop (transformed voxels) and is given by:
$$ \varepsilon_{ij}^* = \frac{A_s}{2V} (n_i b_j + n_j b_i) \qquad (6) $$

where $A_s$ is the area on which planes with normal $n(n_1, n_2, n_3)$ has slipped by a relative amount $b(b_1, b_2, b_3)$ and $V$ is the volume occupied by the loop [40,41]. As before, the dislocation loop is 32 voxels wide in the x and y directions, and 1 voxel thick but now with a z position such that $x + y + z = constant$. The displacement has been computed at the center of voxels with the periodized operator $B'$ along z (figure 2b). As in figure 1c, the displacement in the center of a voxel belonging to the loop plane (black dot in the reddish transformed voxel in figure 2b) is zero. The displacement in the first neighboring voxels (red dot in figure 2b) are shifted relative to the expected position, so that the displacement difference between these voxels is significantly lower than $b$, see figure 2c. It can be checked in figure 2b that each of these voxels shares three faces with a transformed voxel. A more detailed analysis shows that the second neighbors (which share three edges with transformed voxels) are also slightly shifted in the opposite direction. The result is shown in figure 2d with: a strong localized oscillation of the phase (taken here as the displacement modulo $b$).

Although the amplitude of this shift is small (less than 10% of the Burgers vector) it has unwanted consequences on the diffraction peak simulation:

* The dislocation loops are surrounded by four impaired layers of voxels: As the scattered X-Ray amplitude is proportional to the Fourier transform of $G \cdot u$ (see equation (8) in section 4) we can expect a phantom streak in the intensity in a direction perpendicular to the loop plane.

* The displacement field near the edges of the loop (near the dislocation line) will be quite different from its expected value, and the strain field will not vary with the distance $r$ to the dislocation line as $1/r$. This will strongly affect the tails of the diffraction peaks

<figure><img src="image_4.png" /><figcaption>**Figure 2.** (a) Modeling of a squared dislocation loop in a (111) plane as a layer of voxels with eigenstrain; (b) position of the computed points relative to the transformed voxels with eigenstrains; (c) Plot of the displacement field $u_3$ (normalized by $b_3$) along the z direction for dislocation loop</figcaption></figure>
186 illustrated on figure 2. (d) Local oscillation of the phase due to the voxelization of the dislocation loop
187 (the representation is made for 32 voxels centered in the unit cell along z direction). The red line is
188 approximately equal to the phase expected for this displacement field.

189 **3. Sub-voxelization method to correct voxelization effects.**

190 **3.1. Sub-voxelization method**

191 The (conceptually) simplest way to remove this voxelization artifact would be to work on a
192 multiple grid (to multiply the number of voxels along each direction by 2, 4, or more), then to
193 downsample the displacement field data. In that case, FFT algorithms would lose much of their
194 interest due to these more demanding computational efforts. We show below that this can be done
195 in a more economical way by applying a patch to the FFT-computed displacement field. The basic
196 method is to compute on the same grid the difference vector:

$$
197 \Delta_i(\mathbf{x}) = u_i^{sub}(\mathbf{x}) - u_i^{hom}(\mathbf{x}) \qquad (7)
$$

198 where $u_i^{sub}(\mathbf{x})$ is the displacement vector calculated for voxels where this eigenstrain is
199 concentrated on a single plane of sub voxels (figure 3b) and $u_i^{hom}$ the displacement field in direction
200 i of voxels with a uniform eigenstrain (figure 3a). For the sake of clarity, we use 2D diagrams in Figure
201 3, but here the technique is applied to real 3D problems.

<figure><img src="image_5.png" /><figcaption>Figure 3. 2D representation of a dislocation loop in a tilted plane on a (8 × 8) FFT grid. (a) With a homogeneous eigenstrain in the voxels occupied by the dislocation loop. (b) With each voxel subdivided into 4\*4 sub-voxels, only 4 of which are eigenstrained.</figcaption></figure>

205 In order to compute the displacement due to sub-voxels, we use a $N \times N$ ($N \times N \times N$) grid for
206 2D (resp. 3D) problems where each voxel can be subdivided into $n \times n$ ($n \times n \times n$) sub voxels. Only
207 $n$ ($n \times n$) sub voxels are submitted to an eigenstrain field. At a point $\mathbf{A}$ of the grid (black dots, figure
208 4a), we need to compute the sum of the displacements $\mathbf{u}_i^j$ due to the $n$ ($n \times n$) sub voxels $j$ (center
209 $\mathbf{B}_j$) within a voxel centered at point $\mathbf{O}$. This sum is equivalent to the sum of the displacements due
210 to a strained sub voxel at point $\mathbf{O}$ on the grid points $\mathbf{A}_j$ such as $\mathbf{O}\mathbf{A}_j = \mathbf{B}_j\mathbf{A}$ (figure 4b). It is also
211 equivalent to the sum of the displacements $\mathbf{u}'_i^j$ due to a full voxel at point $\mathbf{O}$ on the initial grid on
212 points $\mathbf{A}'_j$ such as $\mathbf{O}\mathbf{A}'_j = n\mathbf{B}_j\mathbf{A}$ (figure 4c). The only difference between these last two sums is due to
213 the long-range strain field, and approximately results in a linear drift of the displacement. As the end
214 of the vectors $\mathbf{O}\mathbf{A}'_j$ does not lie on the grid points (voxel centers) but on the corners of the voxels, the
215 $\mathbf{u}'_i^j$ displacements must be calculated with the shifted operator $\mathbf{B}''$ (equation (5)). A last point is the
216 scaling of the $\mathbf{u}_i^j$ and $\mathbf{u}'_i^j$ sums during the operations of figure 4. To keep the one Burgers vector
217 jump between both sides of the sub voxels plane in Figure 4a, the eigenstrain in the sub voxels must
218 be multiplied by $n$. The backwards change of scale requires a division by $n$: there is no scaling factor
219 between $u_i^{hom}$ and $u_i^{sub} = \sum \mathbf{u}'_i^j$.
<figure><img src="image_6.png" /><figcaption>**Figure 4.** (a). 2D representation of the computational grid. Black dots correspond to the voxels center. A voxel with center O is discretized in 4 × 4 in 2D (4 × 4 × 4 in 3D) sub-voxels. Some sub-voxels contain an eigenstrain (red sub-voxels). We want to compute the displacement field at voxel centered at point A, due to these deformed sub-voxels centered at $B_j$. (b) Effect of one deformed sub-voxel centered at O on a row of sub-voxels centered at $B_j$ such as $OA_j = B_jA$. The sum of these effects is equal to the previous displacement field. (c) Effect of deformed voxels centered at O on a row of voxels (computed at corners $A'_j$ using Green operator $B''$) such as $OA'_j = nB_jA$. This sum is equal to previous wanted sum.</figcaption></figure>

We need to compute $\Delta_{ij\_pl}(\mathbf{x})$ the difference in displacement in direction i due to a voxel which belongs to the plane "pl" (for fcc "pl" is equal to (111), (\bar{1}11), (1\bar{1}1), (11\bar{1}) ) of a dislocation loop with a Burgers vector j at a position x relative to the transformed voxel. In practice, in a material with cubic symmetry, it is sufficient to compute $\Delta_{13\_\bar{(1}1\bar{1)}}(\mathbf{x})$ and $\Delta_{3\bar{3}\_\bar{(1}1\bar{1)}}(\mathbf{x})$, and to use the symmetries of the cube (fourfold [001] axis, threefold [1\bar{1}1] axis, and (1\bar{1}0) symmetry plane) (and suitable exchanges of the components of $\mathbf{x}$) to obtain the required components. As it can be seen in figure 2b, $\Delta_{ij\_pl}(\mathbf{x})$ is non-zero only for the neighbors of the transformed voxel, except the drift due to the long range strain alluded above. The final recipe to compute $\Delta_{ij\_pl}(\mathbf{x})$ and use the patch becomes:

* Compute the field $c^0: \varepsilon^*$ defined in equation (1) for an isolated voxel with the eigenstrain associated to a dislocation loop (equation (6)) with a Burgers vector [001] in a (111) plane (see figure 2).

* Compute the displacement field in directions x($u_1^{hom}$) and z ($u_3^{hom}$) at the voxels' center around the transformed voxel by convolution with the discrete periodized operator $B'$ (equation (4))

* Compute the displacement field in directions x and z at the voxels' corners around the transformed voxel by convolution with the shifted operator $B''$ (equation (5))

* Calculate the $u_1^{sub} = \sum u_1^j$ and $u_3^{sub} = \sum u_3^j$ sums ($n \times n$ terms for each sum) as in figure 4c, then the raw $\Delta_{13(111)}(\mathbf{x})$ and $\Delta_{3\bar{3}(111)}(\mathbf{x})$ for ($x_1, x_2, x_3$) going from -3 to 3 times the voxel size t.

* Use the farthest voxels to correct the drift of the components so that all terms for large $\mathbf{x}$ are zero, and keep only the terms for the first three neighbors non zero.

The patch can then be applied on the raw (FFT-based) displacement field by adding the convolution of all transformed voxels of the different slip systems by the relevant $\Delta_{ij\_pl}(\mathbf{x})$.
## 3.2. Results

For tests, we used the same 128 × 128 × 128 grid as above, and the transformed voxel was divided into 8*8*8 sub voxels (using a reference medium with the same elastic constants as before). Only the final values in units of b (after drift correction) of $\Delta_{13_{(111)}}(\mathbf{x})$ and $\Delta_{33_{(111)}}(\mathbf{x})$ are used and others components are obtained by symmetries (Appendix A).

<figure><img src="image_7.png" /><figcaption>Figure 5. (a) Plot of the displacement field $u_3$ (normalized by $b_3$) along the z direction for dislocation loop illustrated on figure 2. Artifacts are removed by sub-voxel method described above. (b) the phase (i.e. the displacement modulo a Burgers vector). With this correction, the phase is almost continuous.</figcaption></figure>

The patch was used on the same configuration as in figure 2. Figure 5a shows the resulting displacement field and figure 5b the phase (i.e. the displacement modulo a Burgers vector) in Burgers vector units. As it can be observed from figure 5a, the artifacts of the displacement field have disappeared. In addition, the resulting phase varies smoothly even during the crossing of the dislocation loop, see figure 5b.

## 4. Application on diffraction peak simulation

In this section, we show simulated diffraction peaks in order to point the effects of voxelization artifacts and of the patch on numerical results. Under kinematical conditions and assuming a coherent beam, the amplitude of a diffracted wave at a position $\mathbf{q}$ in the vicinity of a reciprocal $\mathbf{G}$ lattice vector is [14,18,42–44]:

$$A(\mathbf{q}) = \mathrm{FT}[A_0(\mathbf{r}) \times F(\mathbf{G}, \mathbf{x}) \times \exp(-2i\pi \mathbf{G} \cdot \mathbf{u}(\mathbf{x}))] \quad (8)$$

where $\mathbf{x}$ is the position of the scattering atom, $A_0(\mathbf{x})$ is the amplitude of the incidence wave, $F(\mathbf{G}, \mathbf{x})$ is the local structure factor and $\mathbf{u}(\mathbf{x})$ the displacement field. The scattered intensity is $I(\mathbf{q})=|A(\mathbf{q})|^2$. For a face-centered cubic crystal, this intensity is non zero when $\mathbf{G}$ (h, k, l) is such as h, k, and l have the same parity. Here two diffraction vectors $\mathbf{G}$ (200) and $\mathbf{G}$ (002) are used. They respectively correspond to $\mathbf{G} \cdot \mathbf{b} = 0$ and $\mathbf{G} \cdot \mathbf{b} = 1$. The 3D diffracted intensity has been calculated using the FFT instead of the continuous Fourier Transform, then summed in the planes perpendicular to the $\mathbf{G}$ vector to obtain a linear plot along $\mathbf{G}$ equivalent to a $I(2\theta)$ plot. In figures 6a ($\mathbf{G}$ (200)) and 6c ($\mathbf{G}$ (002)), we show the diffracted intensity (logarithmic scale) as a function of the pixel position $i$, and in figures 6b and 6d a logarithmic/logarithmic plot of the intensity vs. $|i - i_0|$ where $i_0$ is the center of the peak. In order to only study the effect of the displacement fields, we set $A_0(\mathbf{x}) = 1$ and $F(\mathbf{G}, \mathbf{x}) = 1$ for these simulations.
<figure><img src="image_8.png" /><figcaption>Figure 6. Simulated diffracted intensity as a function of the pixel position (logarithmic scale). 3D configuration is represented in 1D plot by making sum in each plane along an x-axis. Different way for computing the displacement fields are studied for a dislocation loop with a [011] Burgers vector lying in a (111) slip plane. (a) Diffracted vector studied is **G** (200) corresponding to **G**·**b** = 0. (b) log/log representation of the intensity vs. |**i** − **i**₀|. (c) and (d) Same as (a) et (b) but the studied diffracted vector is **G** (002).</figcaption></figure>

The peak shape near the top of the peaks is the same for both computing methods. It is perfectly symmetric in the **G**·**b** = 0 case and exhibits a bump on the right side for **G**·**b** = 1. Their long-range behavior is however quite different. When the displacement field has been calculated with the usual truncated operator (black line), a phantom peak is observed at large |**i** − **i**₀| (at large **q**), which is due to the short period oscillations near the displacement field discontinuity (Figure 1a). The behavior of the peak calculated with the modified Green operator (blue curve) is only slightly better: the intensity at large **q** is underestimated in one case and overestimated in the other. When the intensity has been calculated with the sub voxel patch (red curve) the long range intensity follows the expected I₀|**i** − **i**₀|⁻³ law [45,46]: the peak tails are indeed related to the highly distorted zones near the dislocations' cores. However, the red curve saturates at very large **q**. We suppose this is due to the use of the Fast Fourier Transform instead of the continuous Fourier Transform in the calculation of the scattered amplitude (equation (8)): The plot of figure 6 represents only one period in Fourier space, and is repeated over and over on all Fourier space. We can now calculate the intensity of the tails of these repetitions:

$$I_{neib.} = \sum I_0 |i - i_0 - 128m|^{-3} \qquad (9)$$

where *m* varies from -5 to 5 (zero excluded). We obtain the pink curve at the bottom of figures 5a and 5c. If we now plot the difference between the red intensity curve and this pink background line, we obtain the green curve. On the log./log. plots, figure 6b and 6d it can be checked that this curve follows the I₀|**i** − **i**₀|⁻³ law to the end. Thus, the residual error in the intensity computed by FFT results of the FFT itself, and not from a residual error on the sub voxel corrected displacement field. If the
310 number of voxels is increased to 512³ or 1024³ while keeping the physical size of the Representative
311 Volume constant this residual error should fall down to undetectable levels.

312 **5. Conclusions**

313 In this paper, we showed that although the use of a periodized Green operator in the FFT-based
314 method improves the final displacement field solution in a Representative Volume containing
315 discontinuities (dislocation loops), artifacts due to the voxelization of the dislocation loop planes are
316 still present with respect to analytical solutions. These artifacts have unwanted consequences on the
317 tails of diffraction peaks simulated by using this displacement field as input data.

318 We have introduced a patch which corrects these artifacts by simulating the displacement field
319 which would be obtained with a much finer voxelization without need to do the computations on a
320 finer grid. A simple construction method for this patch has been given and the patch can be used in
321 a single post-processing step to modify the initial FFT-based displacement field.

322 The modified displacement field has been used to simulate one-dimensional diffraction peaks.
323 The procedure strongly improves the shape of the peaks' tails, i.e. it gives a good description of the
324 displacement field near the dislocation lines.

325 Acknowledgments: This work was supported by the French State through the program "Investment in the
326 future" operated by the National Research Agency (ANR) and referenced by ANR-11-LABX-0008-01
327 (LabEx DAMAS). Gabor Ribarik greatly acknowledge the support of the Janos Bolyai Research Fellowship
328 of the Hungarian Academy of Sciences.

329 Author Contributions: For the present paper, Komlavi S. Eloh performed the simulation work and wrote
330 the first draft of the paper; Alain Jacques analyzed numerical results and helped the simulation part; Gabor
331 Ribarik helped in constructing the computer code and helped in analyzing results; Stéphane Berbenni
332 helped in developing the first theoretical part of the paper. The paper was written by Komlavi S. Eloh,
333 Alain Jacques and Stéphane Berbenni.

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

335 **Appendix A**

336 $\Delta_{33\_\langle 111 \rangle}(\mathbf{x}_1, \mathbf{x}_2, \mathbf{x}_3)$ and $\Delta_{13\_\langle 111 \rangle}(\mathbf{x}_1, \mathbf{x}_2, \mathbf{x}_3)$ have these symmetries due to the permutation
337 properties of the plan (111)

$$
338 \Delta_{33\_\langle 111 \rangle}(x_1, x_2, x_3) = \Delta_{33\_\langle 111 \rangle}(x_2, x_1, x_3)
$$

$$
339 \Delta_{33\_\langle 111 \rangle}(x_1, x_2, x_3) = -\Delta_{33\_\langle 111 \rangle}(-x_1, -x_2, -x_3)
$$

$$
340 \Delta_{13\_\langle 111 \rangle}(x_1, x_2, x_3) = -\Delta_{13\_\langle 111 \rangle}(-x_1, -x_2, -x_3)
$$

341 The others values of $\Delta_{ij\_\langle 111 \rangle}(\mathbf{x})$ are given as function of $\Delta_{33\_\langle 111 \rangle}(x_1, x_2, x_3)$ and $\Delta_{13\_\langle 111 \rangle}(x_1, x_2, x_3)$:

$$
342 \Delta_{11\_\langle 111 \rangle}(x_1, x_2, x_3) = \Delta_{33\_\langle 111 \rangle}(x_3, x_2, x_1) \qquad \Delta_{12\_\langle 111 \rangle}(x_1, x_2, x_3) = \Delta_{13\_\langle 111 \rangle}(x_1, x_3, x_3)
$$

$$
343 \Delta_{22\_\langle 111 \rangle}(x_1, x_2, x_3) = \Delta_{33\_\langle 111 \rangle}(x_1, x_3, x_2) \qquad \Delta_{31\_\langle 111 \rangle}(x_1, x_2, x_3) = \Delta_{13\_\langle 111 \rangle}(x_3, x_2, x_1)
$$

$$
344 \Delta_{32\_\langle 111 \rangle}(x_1, x_2, x_3) = \Delta_{13\_\langle 111 \rangle}(x_3, x_1, x_2) \qquad \Delta_{21\_\langle 111 \rangle}(x_1, x_2, x_3) = \Delta_{13\_\langle 111 \rangle}(x_1, x_3, x_2)
$$

$$
345 \Delta_{23\_\langle 111 \rangle}(x_1, x_2, x_3) = \Delta_{13\_\langle 111 \rangle}(x_2, x_1, x_3)
$$

346 The value of $\Delta_{33}(x_1, x_2, x_3)$ and $\Delta_{13}(x_1, x_2, x_3)$ for the remaining plane $(\bar{1}11), (1\bar{1}1), (11\bar{1})$ are
347 obtained using these symmetries:

$$
348 \Delta_{33\_\langle \bar{1}11 \rangle}(x_1, x_2, x_3) = \Delta_{33\_\langle 111 \rangle}(x_2, -x_1, x_3) \qquad \Delta_{13\_\langle \bar{1}11 \rangle}(x_1, x_2, x_3) = \Delta_{13\_\langle 111 \rangle}(-x_1, x_2, x_3)
$$

$$
349 \Delta_{33\_\langle \bar{1}11 \rangle}(x_1, x_2, x_3) = \Delta_{33\_\langle 111 \rangle}(-x_2, x_1, x_3) \qquad \Delta_{13\_\langle \bar{1}11 \rangle}(x_1, x_2, x_3) = \Delta_{13\_\langle 111 \rangle}(x_1, -x_2, x_3)
$$

$$
350 \Delta_{33\_\langle 11\bar{1} \rangle}(x_1, x_2, x_3) = -\Delta_{33\_\langle 111 \rangle}(x_2, x_1, -x_3) \qquad \Delta_{13\_\langle 11\bar{1} \rangle}(x_1, x_2, x_3) = -\Delta_{13\_\langle 111 \rangle}(x_1, x_2, -x_3)
$$

351 **References**

352 1. Graverend, J.-B. L.; Dirand, L.; Jacques, A.; Cormier, J.; Ferry, O.; Schenk, T.; Gallerneau, F.; Kruch, S.;
353 Mendez, J. In Situ Measurement of the γ/γ' Lattice Mismatch Evolution of a Nickel-Based Single-Crystal
354 Superalloy During Non-isothermal Very High-Temperature Creep Experiments. *Metall and Mat Trans A* **2012**, *43*, 3946–3951, doi:10.1007/s11661-012-1343-x.

355 2. Robinson, I.; Harder, R. Coherent X-ray diffraction imaging of strain at the nanoscale. *Nat Mater* **2009**, *8*, 291–298, doi:10.1038/nmat2400.

358 3. Pfeifer, M. A.; Williams, G. J.; Vartanyants, I. A.; Harder, R.; Robinson, I. K. Three-dimensional mapping of a deformation field inside a nanocrystal. *Nature* **2006**, *442*, 63.

360 4. Ungár, T. Strain Broadening Caused by Dislocations Available online: https://www.scientific.net/MSF.278-281.151 (accessed on Jun 16, 2018).

362 5. Ungár, T.; Gubicza, J.; Ribárik, G.; Borbély, A. Crystallite size distribution and dislocation structure determined by diffraction profile analysis: principles and practical application to cubic and hexagonal crystals. *Journal of applied crystallography* **2001**, *34*, 298–310.

365 6. Ribárik, G.; Ungár, T.; Gubicza, J. MWP-fit: a program for multiple whole-profile fitting of diffraction peak profiles by ab initio theoretical functions. *Journal of Applied Crystallography* **2001**, *34*, 669–676.

367 7. Ribárik, G.; Gubicza, J.; Ungár, T. Correlation between strength and microstructure of ball-milled Al–Mg alloys determined by X-ray diffraction. *Materials science and engineering: A* **2004**, *387*, 343–347.

368 8. Balogh, L.; Ribárik, G.; Ungár, T. Stacking faults and twin boundaries in fcc crystals determined by x-ray diffraction profile analysis. *Journal of applied physics* **2006**, *100*, 023512.

371 9. Groma, I. X-ray line broadening due to an inhomogeneous dislocation distribution. *Physical Review B* **1998**, *57*, 7535.

373 10. Tréhorel, R.; Ribarik, G.; Schenk, T.; Jacques, A. Real time study of transients during high temperature creep of a Ni-based superlloy by far field high energy synchrotron X-rays diffraction. *Journal of applied crystallography (under review)*.

376 11. Tréhorel, R. Comportement mécanique haute température du superalliage monocristallin AM1: Etude in-situ par une nouvelle technique de diffraction en rayonnement synchrotron, Université de Lorraine: Nancy, France, 2018.

379 12. Bernier, J. V.; Park, J.-S.; Pilchak, A. L.; Glavicic, M. G.; Miller, M. P. Measuring Stress Distributions in Ti-6Al-4V Using Synchrotron X-Ray Diffraction. *Metallurgical and Materials Transactions A* **2008**, *39*, 3120–3133, doi:10.1007/s11661-008-9639-6.

382 13. Miller, M. P.; Bernier, J. V.; Park, J.-S.; Kazimirov, A. Experimental measurement of lattice strain pole figures using synchrotron x rays. *Review of Scientific Instruments* **2005**, *76*, 113903, doi:10.1063/1.2130668.

384 14. Jacques, A. From Modeling of Plasticity in Single-Crystal Superalloys to High-Resolution X-rays Three-Crystal Diffractometer Peaks Simulation. *Metallurgical and Materials Transactions A* **2016**, *47*, 5783–5797, doi:10.1007/s11661-016-3793-z.

387 15. Weisbrook, C. M.; Gopalaratnam, V. S.; Krawitz, A. D. Use of finite element modeling to interpret diffraction peak broadening from elastic strain distributions. *Materials Science and Engineering: A* **1995**, *201*, 134–142, doi:10.1016/0921-5093(95)09757-0.

390 16. Miller, M. P.; Dawson, P. R. Understanding local deformation in metallic polycrystals using high energy X-rays and finite elements. *Current Opinion in Solid State & Materials Science* **2014**, *5*, 286–299, doi:10.1016/j.cossms.2014.09.001.

393 17. Demir, E.; Park, J.-S.; Miller, M. P.; Dawson, P. R. A computational framework for evaluating residual stress distributions from diffraction-based lattice strain data. *Computer Methods in Applied Mechanics and Engineering* **2013**, *265*, 120–135, doi:10.1016/j.cma.2013.06.002.

396 18. Vaxelaire, N.; Proudhon, H.; Labat, S.; Kirchlechner, C.; Keckes, J.; Jacques, V.; Ravy, S.; Forest, S.; Thomas, O. Methodology for studying strain inhomogeneities in polycrystalline thin films during in situ thermal loading using coherent x-ray diffraction. *New J. Phys.* **2010**, *12*, 035018, doi:10.1088/1367-2630/12/3/035018.

399 19. Song, X.; Xie, M.; Hofmann, F.; Illston, T.; Connolley, T.; Reinhard, C.; Atwood, R. C.; Connor, L.; Drakopoulos, M.; Frampton, L. Residual stresses and microstructure in powder bed direct laser deposition (PB DLD) samples. *International Journal of Material Forming* **2015**, *8*, 245–254.

402 20. Hofmann, F.; Song, X.; Jun, T.-S.; Abbey, B.; Peel, M.; Daniels, J.; Honkimäki, V.; Korsunsky, A. M. High energy transmission micro-beam Laue synchrotron X-ray diffraction. *Materials Letters* **2010**, *64*, 1302–1305.

404 21. Hofmann, F.; Abbey, B.; Liu, W.; Xu, R.; Usher, B. F.; Balaur, E.; Liu, Y. X-ray micro-beam characterization of lattice rotations and distortions due to an individual dislocation. *Nature communications* **2013**, *4*, 2774.
406

407

408

22. Suter, R. M.; Hennessy, D.; Xiao, C.; Lienert, U. Forward modeling method for microstructure reconstruction using x-ray diffraction microscopy: Single-crystal verification. *Review of Scientific Instruments* **2006**, 77, 123905, doi:10.1063/1.2400017.

409

410

411

23. Moulinec, H.; Suquet, P. A numerical method for computing the overall response of nonlinear composites with complex microstructure. *Computer Methods in Applied Mechanics and Engineering* **1998**, *157*, 69–94, doi:10.1016/S0045-7825(97)00218-1.

412

413

24. Moulinec, H.; Suquet, P. Fast numerical method for computing the linear and nonlinear properties of composites. *Comptes Rendus de l'Académie des Sciences. Série II* **1994**, 318.

414

415

25. Michel, J. C.; Moulinec, H.; Suquet, P. A computational scheme for linear and non-linear composites with arbitrary phase contrast. *International Journal for Numerical Methods in Engineering* **2001**, *52*, 139–160.

416

417

26. Müller, W. Mathematical vs. Experimental Stress Analysis of Inhomogeneities in Solids. *Journal de Physique IV Colloque* **1996**, *06*, C1-139-C1-148, doi:10.1051/jp4:1996114.

418

419

27. Vinogradov, V.; Milton, G. W. An accelerated FFT algorithm for thermoelastic and non-linear composites. *International Journal for Numerical Methods in Engineering* **2008**, *76*, 1678–1695, doi:10.1002/nme.2375.

420

421

422

28. Anglin, B. S.; Lebensohn, R. A.; Rollett, A. D. Validation of a numerical method based on Fast Fourier Transforms for heterogeneous thermoelastic materials by comparison with analytical solutions. *Computational Materials Science* **2014**, 209–217.

423

424

29. Graham, J. T.; Rollett, A. D.; LeSar, R. Fast Fourier transform discrete dislocation dynamics. *Modelling Simul. Mater. Sci. Eng.* **2016**, *24*, 085005, doi:10.1088/0965-0393/24/8/085005.

425

426

427

30. Berbenni, S.; Taupin, V.; Djaka, K. S.; Fressengeas, C. A numerical spectral approach for solving elasto-static field dislocation and g-disclination mechanics. *International Journal of Solids and Structures* **2014**, *51*, 4157–4175, doi:10.1016/j.ijsolstr.2014.08.009.

428

429

31. Eyre, D. J.; Milton, G. W. A fast numerical scheme for computing the response of composites using grid refinement. *Eur. Phys. J. AP* **1999**, *6*, 41–47, doi:10.1051/epjap:1999150.

430

431

432

32. Prakash, A.; Lebensohn, R. A. Simulation of micromechanical behavior of polycrystals: finite elements versus fast Fourier transforms. *Modelling Simul. Mater. Sci. Eng.* **2009**, *17*, 064010, doi:10.1088/0965-0393/17/6/064010.

433

434

33. A. Lebensohn, R. N-site modeling of a 3D viscoplastic polycrystal using Fast Fourier Transform. *Acta Materialia* **2001**, *49*, 2723–2737, doi:10.1016/S1359-6454(01)00172-0.

435

436

34. Lebensohn, R. A.; Rollett, A. D.; Suquet, P. Fast fourier transform-based modeling for the determination of micromechanical fields in polycrystals. *JOM* **2011**, *63*, 13–18, doi:10.1007/s11837-011-0037-y.

437

438

439

35. Lebensohn, R. A.; Kanjarla, A. K.; Eisenlohr, P. An elasto-viscoplastic formulation based on fast Fourier transforms for the prediction of micromechanical fields in polycrystalline materials. *International Journal of Plasticity* **2012**, 32–33, 59–69, doi:10.1016/j.ijplas.2011.12.005.

440

441

442

36. Suquet, P.; Moulinec, H.; Castelnau, O.; Montagnat, M.; Lahellec, N.; Grennerat, F.; Duval, P.; Brenner, R. Multi-scale modeling of the mechanical behavior of polycrystalline ice under transient creep. *Procedia IUTAM* **2012**, *3*, 76–90, doi:10.1016/j.piutam.2012.03.006.

443

444

445

37. Eloh, K. S.; Jacques, A.; Berbenni, S. Development of a new consistent discrete Green operator for FFT-based methods to solve heterogeneous problems with eigenstrain. *International Journal of Plasticity*(submitted) **2018**.

446

447

448

38. Hirth, J. P.; Lothe, J. *Theory of Dislocations*; Krieger Publishing Company, 1982; ISBN 978-0-89464-617-1.

449

450

39. Mura, T. *Micromechanics of Defects in Solids*; Mechanics of Elastic and Inelastic Solids; 2nd ed.; Springer Netherlands, 1987; ISBN 978-90-247-3256-2.

451

452

453

40. Li, Q.; Anderson, P. M. A Compact Solution for the Stress Field from a Cuboidal Region with a Uniform Transformation Strain. *Journal of Elasticity* **2001**, *64*, 237–245, doi:10.1023/A:1015203721914.

454

455

41. Anderson, P. M. Crystal-based plasticity. *Fundamentals of Metals Forming* 1982, 277–279.

456

457

42. Takagi, S. A Dynamical Theory of Diffraction for a Distorted Crystal. *Journal of the Physical Society of Japan* **1969**, *26*, 1239–1253, doi:10.1143/JPSJ.26.1239.

458

459

43. Vartanyants, I. A.; Yefanov, O. M. Coherent X-ray Diffraction Imaging of Nanostructures. *arXiv:1304.5335 [cond-mat]* **2013**.

44. Takagi, S. Dynamical theory of diffraction applicable to crystals with any kind of small distortion. *Acta Crystallographica* **15**, 1311–1312, doi:10.1107/S0365110X62003473.

45. Ungár, T. Microstructural parameters from X-ray diffraction peak broadening. *Scripta Materialia* **2004**, *51*, 777–781, doi:10.1016/j.scriptamat.2004.05.007.
460
461

46. Krivoglaz, M. A. Theory of X-Ray and Thermal Neutron Scattering by Real Crystals; Springer US, 1969; ISBN 978-1-4899-5584-5.

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

© 2018 by the authors. Submitted for possible open access publication under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).