The Convolutional Multiple Whole Proﬁle (CMWP) Fitting Method, a Global Optimization Procedure for Microstructure Determination

: The analysis of line broadening in X-ray and neutron di ﬀ raction patterns using proﬁle functions constructed on the basis of well-established physical principles and TEM observations of lattice defects has proven to be a powerful tool for characterizing microstructures in crystalline materials. These principles are applied in the convolutional multiple-whole-proﬁle (CMWP) procedure to determine dislocation densities, crystallite size, stacking fault and twin boundary densities, and intergranular strains. The di ﬀ erent lattice defect contributions to line broadening are separated by considering the hkl dependence of strain anisotropy, planar defect broadening and peak shifts, and the defect dependent proﬁle shapes. The Levenberg–Marquardt (LM) peak ﬁtting procedure can be used successfully to determine crystal defect types and densities as long as the di ﬀ raction patterns are relatively simple. However, in more complicated cases like hexagonal materials or multiple-phase patterns, using the LM procedure alone may cause uncertainties. Here, we extended the CMWP procedure by including a Monte Carlo statistical method where the LM and a Monte Carlo algorithm were combined in an alternating manner. The updated CMWP procedure eliminated uncertainties and provided global optimized parameters of the microstructure in good correlation with electron microscopy methods. in a hydrated blister on a Zircaloy-4 specimen. The two-dimensional lateral and depth proﬁles of the dislocation densities, the stacking fault densities, and the subgrain size were determined by evaluating position-sensitive high-resolution synchrotron X-ray di ﬀ raction patterns. The results were in good correlation with the mechanical and electron microscopy analysis done on the same specimen.

The Levenberg-Marquardt optimization procedure has been used successfully in innumerable cases, cf.[5][6][7]19,[25][26][27][28][29][30][31][32][33][34][35] to determine crystal defect types and densities.In the present work, we developed a more robust procedure where the Levenberg-Marquardt and the Monte Carlo methods are applied consecutively to provide the global optimum values of the physical parameters characterizing microstructures.In a previous short letter-type report [15], the basics of the method were outlined briefly.Here, we provide a detailed and more elaborate description of the procedure and show its power by analyzing the depth dependence of the dislocation density, the planar defect density, and the grain size in a Zr matrix and in ZrH precipitates in a hydrated Zr specimen.

Fundamental Principles of the Convolutional Multiple-Whole-Profile (CMWP) Optimization Procedure
The convolutional multiple-whole-profile (CMWP) optimization method is based on physically modeled profile functions of different microstructure elements [12][13][14][15].The two fundamental microstructure elements are size and strain [1].In diffraction patterns they combine as convolution, and in a particular hkl diffraction peak they are where I S hkl (s) and I D hkl (s) are the size and strain profiles.The variable s is where g hkl is the fundamental reciprocal lattice vector of the hkl reflection and K is an arbitrary reciprocal space vector.Diffraction peaks are three-dimensional in reciprocal space [35].In powder diffraction experiments, the intensity distributions are integrated along the surface perpendicular to the diffraction vector, g hkl , and Equations (1) and (2) reduce to one dimension in the radial direction.
I hkl (s) = I S hkl (s) * I D hkl (s) where s = 2(sinθ − sinθ B )/λ, θ and θ B are the diffraction angle and the exact Bragg angle of the hkl peak, λ is the wavelength of radiation and * indicates convolution.The equivalent of Equation (1) in Fourier space is where L is the Fourier variable.If planar defects, intergranular strains, and instrumental effects become substantial, the above equation is extended: where A PD hkl (L), A IGS hkl (L), and A Inst hkl (L) are the Fourier transforms of the profile functions of planar defects, intergranular strains, and instrumental effects.The Fourier transforms A S hkl (L), A D hkl (L), A PD hkl (L), and A IGS hkl (L) are modeled using physical properties of these lattice defects [15,25,[30][31][32][33][34], whereas A Inst hkl (L) is determined by using measured patterns of defect-free standard materials.The latter can be determined by measuring the patterns of standard specimens, e.g., of Si, CeO 2 , diamond, or LaB 6 standards.The diffraction pattern is calculated from the inverse Fourier transforms of A hkl (L): Crystals 2020, 10, 623 where θ hkl 0 and I hkl 0 are the exact Bragg position and the intensity of the hkl reflection.The least-squares optimization is made by minimizing the weighted sum of squared residuals (WSSR): where BG(2θ i ) is the background intensity at 2θ i and w i are weights applied to the i th measured data.
The size profile is calculated by taking into account the shape and size distribution of coherently scattering domains [17].Assuming log-normal size distribution, the Fourier transform of size profiles can be written as [12,19] A S (L) = m 3 exp(4.5σ 2 ) where L is the Fourier variable, m and σ are the median and variance of the log-normal size distribution density function, and erfc is the complementary error function.Equation ( 8) is valid for spherical crystallites.The CMWP procedure yields the option to evaluate elliptical flat of long crystallites using the ellipticity parameter e which is the ratio of the long and short axis of rotational ellipsoids.With m and σ, the arithmetic-, area-, and volume-weighted mean crystallite diameters can be calculated [36].
where k = 0.5, 2.5, and 3.5 for arithmetic-, area-, and volume-weighted means, respectively, and j stands for these averages.
The Fourier coefficients of the strain profile are [1] where ε 2 g,L is the mean square strain (mss).Krivoglaz showed that strain broadening can only be caused by one-dimensional linear defects [20,21], and calculated ε 2 g,L for dislocations at small L values: where ρ, C, and b are the density, the contrast, and the Burgers vector of dislocations, and D is the size of the crystal.This expression is logarithmically singular with increasing crystal size and is only valid for small L values.Wilkens [2,3] corrected the expression by replacing the crystal size with the effective outer cut-off radius of dislocations, R e , and calculated ε 2 g,L in the entire L range: where η=L/R e .The f (η) function is logarithmic at small L values and hyperbolic at large L values.The explicit form of f (η) is given in Equations (A6) to (A8) in Reference [2] and is shown in Figure 1.
The figure shows f (η) as a solid blue line, and the logarithmic and hyperbolic components are shown as dash-dot-dot and dashed lines, respectively.The profile functions of faulted and twinned crystals consist of several sub-profiles [1].Typical sub-profiles of the 311 reflection of copper containing 4% intrinsic stacking faults are shown in Figure 2 [30].The number, positions and breadths of sub-profiles depend on the hkl indices of the fundamental Bragg peaks [1,[27][28][29][30][31][32].In the CMWP procedure, the shifts and breadths of the sub-profiles are parameterized as a function of the density of specific planar faults or twin boundaries [30,31].The fractions of the sub-profiles are based on the multiplicities within the hkl reflection.The shifts and breadths of the sub-profiles are given by fifth-order polynomials of the planar fault density.The coefficients of the polynomials, along with the fractions of sub-profiles, are listed in the materials' specific parameter tables.CMWP uses these tables to evaluate the density of planar defects.The parameter tables are free to edit by users.Based on theoretical considerations, it was shown that the sub-profiles are the sum of a symmetrical and an anti-symmetrical Lorentzian function [30,31].The theoretical derivations were verified by numerical simulations using the DIFFaX software [32].The anti-symmetrical component of the profile function is produced by interference between two overlapping sub-reflections in reciprocal space.One of these sub-reflections corresponds to the parent, whereas the other corresponds to the twin crystal.The symmetrical and the anti-symmetrical parts of the sub-reflections are correlated with each other and can be characterized by a breadth and an anti-symmetry parameter [27].These two values are also parameterized as a function of stacking fault and twin density, and are included in the parameter tables freely available through the web [30,31].
stands for these averages.
The Fourier coefficients of the strain profile are [1] where 〈 , 2 〉 is the mean square strain (mss).Krivoglaz showed that strain broadening can only be caused by one-dimensional linear defects [20,21], and calculated 〈 , 2 〉 for dislocations at small L values: where , C, and b are the density, the contrast, and the Burgers vector of dislocations, and D is the size of the crystal.This expression is logarithmically singular with increasing crystal size and is only valid for small L values.Wilkens [2,3] corrected the expression by replacing the crystal size with the effective outer cut-off radius of dislocations, Re, and calculated 〈 , 2 〉 in the entire L range: where =L/Re.The f() function is logarithmic at small L values and hyperbolic at large L values.The explicit form of f() is given in Equations (A6) to (A8) in Reference [2] and is shown in Figure 1.The figure shows f() as a solid blue line, and the logarithmic and hyperbolic components are shown as dash-dot-dot and dashed lines, respectively.(blue curve).The dot-dot-dash line is the logarithmic part, as in the Krivoglaz approximation [21].The dash curve is the hyperbolic part.(Copyright by courtesy of Ribárik [12].) The profile functions of faulted and twinned crystals consist of several sub-profiles [1].Typical sub-profiles of the 311 reflection of copper containing 4% intrinsic stacking faults are shown in Figure 2 [30].The number, positions and breadths of sub-profiles depend on the hkl indices of the fundamental Bragg peaks [1,[27][28][29][30][31][32].In the CMWP procedure, the shifts and breadths of the subprofiles are parameterized as a function of the density of specific planar faults or twin boundaries [30,31].The fractions of the sub-profiles are based on the multiplicities within the hkl reflection.The shifts and breadths of the sub-profiles are given by fifth-order polynomials of the planar fault density.The coefficients of the polynomials, along with the fractions of sub-profiles, are listed in the materials' specific parameter tables.CMWP uses these tables to evaluate the density of planar defects.The parameter tables are free to edit by users.Based on theoretical considerations, it was shown that the sub-profiles are the sum of a symmetrical and an anti-symmetrical Lorentzian function [30,31].The theoretical derivations were verified by numerical simulations using the DIFFaX software [32].The anti-symmetrical component of the profile function is produced by interference between two overlapping sub-reflections in reciprocal space.One of these sub-reflections corresponds to the parent, whereas the other corresponds to the twin crystal.The symmetrical and the anti-symmetrical parts of the sub-reflections are correlated with each other and can be characterized by a breadth and an anti-symmetry parameter [27].These two values are also parameterized as a function of stacking fault and twin density, and are included in the parameter tables freely available through the web [30,31].

Physical and Secondary Parameters in the CMWP Procedure
The aim of the CMWP procedure is the qualitative and quantitative characterization of microstructures in crystalline materials.The crystallographic structure of the different components in the materials is assumed to be known.In this sense, CMWP is different, but complementary to Rietveld methods.The parameters used in the optimization procedure are divided into two groups of (i) physical and (ii) secondary ones.Physical parameters include the size distribution parameters, m and , the dislocation parameters,  and M, the planar defect densities,  or  for stacking faults of twin boundaries, and the elastic intergranular strain,  IGS .The peak positions and peak intensities are also fitted during the optimization procedure; however, these parameters are only used to improve the match between the measured and calculated diffraction patterns.

Physical and Secondary Parameters in the CMWP Procedure
The aim of the CMWP procedure is the qualitative and quantitative characterization of microstructures in crystalline materials.The crystallographic structure of the different components in the materials is assumed to be known.In this sense, CMWP is different, but complementary to Rietveld methods.The parameters used in the optimization procedure are divided into two groups of (i) physical and (ii) secondary ones.Physical parameters include the size distribution parameters, m and σ, the dislocation parameters, ρ and M, the planar defect densities, α or β for stacking faults of twin boundaries, and the elastic intergranular strain, ε IGS .The peak positions and peak intensities are also fitted during the optimization procedure; however, these parameters are only used to improve the match between the measured and calculated diffraction patterns.

The Density and the Arrangement Parameter of Dislocations
Krivoglaz and coworkers [20,21] showed that the Fourier transform of the strain profile is a logarithmic function of K as long as D/L is larger than unity, where D is a number of length dimensions.In References [20,21] D was suggested to be the grain size.The problem with this definition of D is the same as in the elastic stored energy of dislocations, which becomes logarithmically singular if the upper integration limit is set to be the crystal size.In the case of elastic stored energy, this problem is solved by setting the upper integration limit by the average distance of dislocations [37,38].Wilkens realized the correlation between the elastic stored energy of dislocations and strain broadening of diffraction peaks, and replaced D in the Krivoglaz formula with the effective outer cut-off radius of dislocations, R e [2,3,39].Wilkens introduced the concept of "restrictedly random" dislocation distributions and calculated the strain function, f (η), in the entire L range [3].The concept assumes that the equal number of plus and minus straight parallel screw dislocations are randomly distributed within cylinders of R e radii, where the cylinders cover the whole crystal.Since the value of R e depends on the actual dislocation density, Wilkens introduced the dimensionless parameter M = R e √ ρ to describe the dislocation arrangements [3].The value of M is large when the dislocations are uncorrelated and arranged randomly and the related strain fields are of long-range character.In contrast, M is small when opposite-sign dislocations are in strong correlation close to each other and the related strain fields, due to screening, are of short-range character.An intuitive schematic illustration is shown in in Figure 2, where 12 dislocations are randomly distributed in Figure 2a and in a strong dipole configuration in Figure 2b.In the case of random distribution, the strain fields reach out much further than the average dislocation distance, d disloc , (the red circle is R e and the blue double arrow the average dislocation distance), whereas in the case of strong dipole configuration, the strain fields are strongly screened and R e becomes shorter than d disloc .The M value is in direct correlation with the average dislocation distance: Random or strong correlated dislocation arrangements can be called weak or strong dipole characters.In the first case, M is larger than unity, M 1, whereas in the second case, it is close to unity or even smaller: M ≤ 1. Due to the reciprocity of length scales in crystal and reciprocal space, Crystals 2020, 10, x FOR PEER REVIEW 5 of 18 Krivoglaz and coworkers [20,21] showed that the Fourier transform of the strain profile is a logarithmic function of K as long as D/L is larger than unity, where D is a number of length dimensions.In References [20,21] D was suggested to be the grain size.The problem with this definition of D is the same as in the elastic stored energy of dislocations, which becomes logarithmically singular if the upper integration limit is set to be the crystal size.In the case of elastic stored energy, this problem is solved by setting the upper integration limit by the average distance of dislocations [37,38].Wilkens realized the correlation between the elastic stored energy of dislocations and strain broadening of diffraction peaks, and replaced D in the Krivoglaz formula with the effective outer cut-off radius of dislocations, Re [2,3,39].Wilkens introduced the concept of "restrictedly random" dislocation distributions and calculated the strain function, f(), in the entire L range [3].The concept assumes that the equal number of plus and minus straight parallel screw dislocations are randomly distributed within cylinders of Re radii, where the cylinders cover the whole crystal.Since the value of Re depends on the actual dislocation density, Wilkens introduced the dimensionless parameter M = Re√ to describe the dislocation arrangements [3].The value of M is large when the dislocations are uncorrelated and arranged randomly and the related strain fields are of long-range character.In contrast, M is small when opposite-sign dislocations are in strong correlation close to each other and the related strain fields, due to screening, are of short-range character.An intuitive schematic illustration is shown in in Figure 2, where 12 dislocations are randomly distributed in Figure 2a and in a strong dipole configuration in Figure 2b.In the case of random distribution, the strain fields reach out much further than the average dislocation distance, ddisloc, (the red circle is Re and the blue double arrow the average dislocation distance), whereas in the case of strong dipole configuration, the strain fields are strongly screened and Re becomes shorter than ddisloc.The M value is in direct correlation with the average dislocation distance: Random or strong correlated dislocation arrangements can be called weak or strong dipole characters.In the first case, M is larger than unity, M ≫ 1, whereas in the second case, it is close to unity or even smaller: M  1. Due to the reciprocity of length scales in crystal and reciprocal space, Intensity distributions of the {110} reflections of a V-5Cr-5Ti alloy deformed by highpressure torsion at 4 GPa pressure to the shear strain of  = 37 (red curve) and after heating to 800 C (blue curve) with logarithmic intensity scale normalized to the maximum intensities and the FWHM.More details about material processing and the heat treatments can be found in Reference [40].
When the strain fields are of long-or short-range character, the tail regions of diffraction peaks decay faster or slower.Figure 3 shows the intensity distributions of the {110} reflections of a V-5Cr-5Ti [40] Figure 3. Intensity distributions of the {110} reflections of a V-5Cr-5Ti alloy deformed by high-pressure torsion at 4 GPa pressure to the shear strain of γ = 37 (red curve) and after heating to 800 • C (blue curve) with logarithmic intensity scale normalized to the maximum intensities and the FWHM.More details about material processing and the heat treatments can be found in Reference [40].When the strain fields are of long-or short-range character, the tail regions of diffraction peaks decay faster or slower.Figure 3 shows the intensity distributions of the {110} reflections of a V-5Cr-5Ti [40] alloy deformed by high-pressure torsion at 4 GPa pressure to the shear strain of γ = 37 (blue curve) and after.
Heating to 800 • C (red curve) with a logarithmic intensity scale normalized to the maximum and the FWHM.The dislocation density did not decrease after heat treatment at 800 • C for 30 min, but the dislocations rearranged into narrow dipole configurations [40].This rearrangement of dislocations is a kind of recovery process before recrystallization, when the dislocations annihilate and their density drops to low values.In the as-deformed state (blue curve), the tail of the peak decays much faster than in the heat-treated (red curve) recovery state.In the recovery state, the plus-minus dislocations form close pairs of dipoles for which the efficiently screened strain fields are of short-range character.The change of the M parameter from 5.8 to 1.6 is the quantitative measure of rearrangement of the dislocation distribution during the recovery process.When the strain fields are of long-range character, i.e.M 1, the profile tails will be short and the profiles will have a bell-shaped character.Meanwhile, when the strain fields are of short-range character, the profile tails will be long and the profiles will be more of Lorentzian character.This is one of the reasons why attributing Gaussian or Lorentzian components in a pseudo-Voight profile function to strain or size broadening is physically incorrect.

The Contrast Factor of Dislocations
The mean square strain is proportional to the contrast, C, of dislocations.In X-ray or neutron diffraction, the contrast has the same physical meaning and origin as in TEM.It depends on the relative orientation between the Burgers and line vectors, b and l, the particular diffraction vector, g, and the elastic constants of the crystal, c ijkl : C = C(b,l,g,c ijkl ) [2,3,19].In a polycrystal, or if all possible Burgers vectors are randomly populated, C can be averaged either over the permutations of the hkl indices or all possible Burgers vectors [41,42].In the case of cubic and hexagonal crystals, the average contrast factors, C, can be written as [5,43,44] where C h00 and C hk.0 are the average contrast factors for of the h00 and hk.0 type reflections, 2 , and x = (2/3)(l/ga) 2 (a is the basal lattice constant of the hexagonal crystal).
In the CMWP code, Equation ( 15) is included in a more direct manner as Equation (15a) [12]: where . The mss in Equation ( 12) consists of ρ, C, and b in the product.Therefore, only ρ and the hkl dependent q or q 1 and q 2 can be determined independently via the CMWP procedure.The scaling factors C h00 and C hk.0 can only be refined after the q or q 1 and q 2 values have been determined.The same is true for the Burgers vector values of <a>, <c>, and <c+a> slip modes in hcp materials [15,45].The effect of changing dislocation density on the elastic constants and the C h00 and C hk.0 scaling constants was investigated in References [46,47].In Reference [47], it was shown that the effects could change the absolute values of dislocation densities up to 40%.

Determination of Slip Modes in hcp Crystals
Once the a 1 and a 2 parameters are determined, the true dislocation densities, ρ true , and the <a>, <c>, and <c + a> slip mode fractions in hcp materials can be determined [15,45].A simpler and more straightforward method than in References [15,45] is briefly described here.The mss can be written as the sum of the partial mss values, ε 2 i,g,L , related to the different slip modes: where i stands for the different slip modes.We can assume that the strain function f (η) is global for all hk.l peaks: where ρ* and C e f f are the dislocation density and the contrast factor values given by the CMWP procedure.Here, , where C* hk.0 is the input value of the scaling factor and q 1 * and q 2 * are the contrast factor parameters given by the CMWP procedure.The fractions, ξ a , ξ c and ξ c+a of the partial mss values are where i ξ i = 1.The partial dislocation densities will be Since Equations ( 5) must hold for all hk.l values, ξ i can be obtained using the least-squares method.
The CMWP procedure allows the evaluation of diffraction patterns consisting of more than one phase.In such cases, the number of physical parameters increases.In the simplest case of a cubic single-phase material with a spherical crystallite shape and no planar defects, the number of physical parameters will be five: ρ, M, and q for the strain profile and m and σ for the size profile.In an hcp specimen, containing two different phases, the number of physical parameters becomes twelve: ρ, M, m, σ, and q 1 , q 2 , and ρ 1 , M 1 , m 1 , σ 1 , and q 1 1 , q 1 2 , where the upper index '1' is for the second phase.(Note that the major phase is not labeled.)Since line profile analysis evaluates only the shape and broadening of diffraction peaks, the peak positions and maxima are treated in CMWP as secondary parameters.The peak shifts caused by faulting or twinning are taken into account by the shifts of the sub-profiles included in the parameter files for planar defects [30,31].

Algorithms Used for Solving Equation (7)
Equation ( 7) requires a non-linear problem with correlations between the different parameters to be solved.The possibilities are (a) nonlinear least-squares algorithms [48], (b) direct search methods [49], and (c) statistical methods [50].The most frequently used nonlinear least-squares algorithms are the Newton method, the steepest descent or conjugate gradient method, and the LM [51,52] procedure.The nonlinear least-squares algorithms find the local minimum of the WSSR in the parameter space, which does not necessarily coincide with the global minimum of the WSSR.Statistical methods [37,39] can be used to find the global minimum of the WSSR in the parameter space.In order to find the global minimum in the CMWP procedure, we applied a special MC statistical algorithm and the LM nonlinear least-squares algorithm in a multiple successive process.The MC method only optimizes the physical parameters.Based on physical considerations, each parameter is restricted between a minimum and maximum value which cannot be bypassed; however, these can be edited by the user.In consecutive iterations, the new parameter values are searched in the proximity of the previous ones: Crystals 2020, 10, 623 where i is the parameter index and n is the current iteration number.∆ i n is defined as The value of ∆ i n is the same during the first n 0 iterations.During further iterations, ∆ i n decreases exponentially to ∆ i min .Once ∆ i n reaches ∆ i min , when n reaches n 1 , the ∆ i n values remain constant.For example, if ∆ i 0 / i min = 50, the value of ∆ i min is reached after about 4000 iterations, as shown in Figure 4.In the special MC algorithm developed for CMWP, we applied a bias towards the prevailing parameter values in the consecutive iterations.The consecutive parameter values follow a cubic probability function vs. x i n generated by a random number generator: where The star in a i * n indicates that this consecutive value has not yet been accepted as a better value than a i n−1 .The biased cubic probability function is shown schematically in Figure 4b.
Crystals 2020, 10, x FOR PEER REVIEW 8 of 18 shown in Figure 4.In the special MC algorithm developed for CMWP, we applied a bias towards the prevailing parameter values in the consecutive iterations.The consecutive parameter values follow a cubic probability function vs.    generated by a random number generator: where .The biased cubic probability function is shown schematically in Figure 4b.Due to the biased probability function, the parameter space is sampled on a finer scale around preceding values.However, as Figure 4b shows, there is also a good chance of sampling remote regions of the parameter space.The condition for accepting the new parameters,    is  The errors of the parameter values are determined in terms of p% fractions of the WSSR [53].If a particular WSSR value is larger than the best value, WSSRbest, found until the actual iteration, but within a certain p% fraction of WSSRbest, i.e., then the parameter values corresponding to this WSSR value are stored in an array.At the end of the fitting procedure, the largest and the smallest of each parameter value in the array will be considered as the plusminus absolute errors of the respective parameters with the given p% number.The relative errors are calculated and listed accordingly.The convergence criteria of the special MC procedure are that (i) the number of steps reaches a user-defined value and (ii) that the 1 + p% statistics reaches at least Np, where Np is an editable minimum number of MC iterations in which the WSSR was not larger than the last WSSR value.The systematic testing of the combined LM + MC algorithm showed that n0 = 2000, p = 3.5% and Np = 100 gave excellent pattern fittings and good physical parameter values even if the diffraction patterns were complicated.Our experience showed that the MC procedure outlined above is ideally suited in the case of line profile analysis.In agreement with the copyright conditions of CMWP, the source code of this MC is accessible in the software package, which is free to use for academia purposes.

Organization of the Combined LM and MC Algorithms
The total number of parameters can become quite large when more than ten diffraction peaks and more than one phase comprise the diffraction pattern.In such cases, the LM procedure can get Due to the biased probability function, the parameter space is sampled on a finer scale around preceding values.However, as Figure 4b shows, there is also a good chance of sampling remote regions of the parameter space.The condition for accepting the new parameters, a i n is The errors of the parameter values are determined in terms of p% fractions of the WSSR [53].If a particular WSSR value is larger than the best value, WSSR best , found until the actual iteration, but within a certain p% fraction of WSSR best , i.e., WSSR > WSSR best AND WSSR < (1 + p%)•WSSR best , (24) then the parameter values corresponding to this WSSR value are stored in an array.At the end of the fitting procedure, the largest and the smallest of each parameter value in the array will be considered as the plus-minus absolute errors of the respective parameters with the given p% number.The relative errors are calculated and listed accordingly.The convergence criteria of the special MC procedure are that (i) the number of steps reaches a user-defined value and (ii) that the 1 + p% statistics reaches at least N p , where N p is an editable minimum number of MC iterations in which the WSSR was not larger than the last WSSR value.The systematic testing of the combined LM + MC algorithm showed that n 0 = 2000, p = 3.5% and N p = 100 gave excellent pattern fittings and good physical parameter values even if the diffraction patterns were complicated.Our experience showed that the MC procedure outlined above is ideally suited in the case of line profile analysis.In agreement with the copyright conditions of CMWP, the source code of this MC is accessible in the software package, which is free to use for academia purposes.

Organization of the Combined LM and MC Algorithms
The total number of parameters can become quite large when more than ten diffraction peaks and more than one phase comprise the diffraction pattern.In such cases, the LM procedure can get frustrated and the optimization procedure can get stuck.This problem is overcome by sectioning the LM procedure.One full cycle of the optimization procedure consists of the following sections: (1) MC optimization, (2) LM optimization of only the peak positions, (3) LM optimization of only the peak heights, (4) LM optimization of only the physical parameters, (5) LM optimization of all parameters including the physical parameters, the peak positions, and peak heights, and (6) readjustment of the BG.The cycle of six steps is repeated until convergence is reached.

Systematic Comparison of the Performance of the LM and MC Procedures
Though the LM optimization procedure is a well-established and widely used analytical method, in nonlinear problems, the parameter space can have more than a single minimum and LM can get stuck in one which is not the absolute minimum.The implementation of the MC procedure combined with LM aimed to find the global minimum in the parameter space.A systematic analysis was carried out to check the performance of both the LM procedure alone and the combined alternative with LM plus MC.The most critical parameters were the dislocation density concomitant with the dipole character number, i.e., ρ and M. One of the Zircaloy-4 X-ray diffraction patterns containing about 10% δ-hydride, discussed in detail in the next paragraph, was evaluated by varying the starting values of ρ over 5 orders of magnitude, from 10 −3 × 10 14 m −2 to 100 × 10 14 m −2 .The results are shown in Figure 5. CMWP is governed by optimizing the parameters for obtaining the smallest weighted sum of squared residuals (WSSR).The WSSR values are shown in Figure 5a vs. the starting values of ρ when LM alone (open blue triangles) and LM combined with MC (open red circles) were applied.The figure shows that LM alone could not find the minimum of WSSR, especially when the initial ρ values were larger than the optimum, which was about 0.9(± 0.1) × 10 14 m −2 .The combined application of LM + MC always found the optimum WSSR whatever the initial values of ρ were.The dislocation density and the dipole character number, ρ and M, are shown vs. the initial ρ values in Figure 5b.The figure shows that both ρ and M varied up to about two orders of magnitude when only LM was applied: 0.7 ≤ ρ ≤ 100 × 10 14 m −2 , 0.4 ≤ M ≤ 60.With the combination of LM + MC, the optima of both ρ and M were found within reasonable error margins: 0.7 ≤ ρ ≤ 1.13 × 10 14 m −2 , 10 ≤ M ≤ 60.The relatively large variation of M does not mean that it could not be determined.Whenever the value of M is of the order of 10 or larger, it only means that the dipole character of dislocations is very weak.Or, in other words, the dislocation arrangement is "random" [2,3].

Figure 5b
. The figure shows that both  and M varied up to about two orders of magnitude when only LM was applied: 0.7   10010 14 m −2 , 0.4  M  60.With the combination of LM+MC, the optima of both  and M were found within reasonable error margins: 0.7    1.13  10 14 m −2 , 10  M  60.The relatively large variation of M does not mean that it could not be determined.Whenever the value of M is of the order of 10 or larger, it only means that the dipole character of dislocations is very weak.Or, in other words, the dislocation arrangement is "random" [2,3].

Extension of CMWP for Handling Satellites or Diffuse Scattering
Guinier-Preston (GP) zones [54,55], small dislocation loops [56,57], or diffuse scattering of solute atoms [58] can produce satellites or diffuse scattering peaks around or below the fundamental Bragg

Extension of CMWP for Handling Satellites or Diffuse Scattering
Guinier-Preston (GP) zones [54,55], small dislocation loops [56,57], or diffuse scattering of solute atoms [58] can produce satellites or diffuse scattering peaks around or below the fundamental Bragg peaks.Typical satellites in the lower intensity ranges of diffraction peaks of a proton-irradiated Zircaloy-2 specimen are shown in Figure 6a.The vertical arrows point to the satellite peaks produced either by vacancy-(on the left side of the peaks) or interstitial-type dislocation loops (on the right side of the 00.2 reflection).The interstitial-type satellite on the 10.1 peak cannot be seen because it is blending in with the broader main peak.The strain profile given by Equations ( 10)-( 12) relies strongly on the tail regions of Bragg peaks.Satellites or diffuse scattering peaks, not related directly to the strain profiles, may distort the physical parameter values defining the true strain profiles and the true strain parameters.The option to handle satellites or diffuse scattering peaks is also implemented in the CMWP procedure [25].Irradiation-induced dislocation loops evolve from knocked out vacancies and interstitials in irradiated materials [56,57,[59][60][61][62]. Larson and Young [61] and Mason et al. [62] showed that irradiation-induced dislocation loops have a very wide size distribution, ranging from a few tenths to a few hundred nanometers.The strain fields of small dislocation loops overlap and strengthen each other inside the loop regions.As a consequence, within the strained regions of the small vacancyor small interstitial-type loops, the lattice constants increase or decrease relative to the matrix lattice constants, respectively [61,63].The strained regions around small loops produce satellite peaks around the fundamental Bragg reflections in proton-, neutron-, or ion-irradiated materials [25,57,62].Based on the theoretical description of satellites and diffuse scattering of point defect clusters and small dislocation loops [61,63,64], we developed a procedure to handle satellites in the CMWP procedure [25].Taking into account the small size and wide size distribution of the strained volumes of small loops, the satellite profiles, I SAT hkl (s + ∆s), were modeled as size profiles according to Equation ( 8), where ∆s is the shift of the satellite peaks relative to the main hkl reflection.In References [64-66], it was shown that although the strained volumes are coherent with the matrix, the intensities scattered by the strained volumes are incoherent; therefore, the intensities of the main and the satellite peaks are additive.The satellite intensities are usually a few percent of the intensities of the main peaks; therefore, the optimization of the main diffraction pattern and the satellite peaks can be done in separate steps.The two steps can be repeated subsequently until the procedure converges.Figure 6b,c shows zoomed parts of the measured (open circles) and CMWP-calculated (red lines) diffraction patterns of a Zircaloy-2 specimen neutron irradiated to the fluence of 13.1 × 1025 n/m2 at about 300(± 25) • C [58,67].
The satellites are shown as blue lines.The shift of satellites is s = s 0 Cg, where s0 is a global shift parameter for the entire pattern and ∆s scales for the different reflections, with the average contrast factor of <a>-type dislocation loops [25].The intensities of the satellites were adjusted to the heights of the related main peaks.Although the procedure was conceived for treating satellites produced by irradiation-induced dislocation loops, it is a general option in the CMWP procedure which can be used to handle diffuse scattering satellites of any kind around fundamental Bragg peaks.
blending in with the broader main peak.The strain profile given by Equations ( 10)-( 12) relies strongly on the tail regions of Bragg peaks.Satellites or diffuse scattering peaks, not related directly to the strain profiles, may distort the physical parameter values defining the true strain profiles and the true strain parameters.The option to handle satellites or diffuse scattering peaks is also implemented in the CMWP procedure [25].Irradiation-induced dislocation loops evolve from knocked out vacancies and interstitials in irradiated materials [56,57,[59][60][61][62]. Larson and Young [61] and Mason et al. [62] showed that irradiation-induced dislocation loops have a very wide size distribution, ranging from a few tenths to a few hundred nanometers.The strain fields of small dislocation loops overlap and strengthen each other inside the loop regions.As a consequence, within the strained regions of the small vacancy-or small interstitial-type loops, the lattice constants increase or decrease relative to the matrix lattice constants, respectively [61,63].The strained regions around small loops produce satellite peaks around the fundamental Bragg reflections in proton-, neutron-, or ion-irradiated materials [25,57,62].Based on the theoretical description of satellites and diffuse scattering of point defect clusters and small dislocation loops [61,63,64], we developed a procedure to handle satellites in the CMWP procedure [25].Taking into account the small size and wide size distribution of the strained volumes of small loops, the satellite profiles,  ℎ  ( + ), were modeled as size profiles according to Equation (8), where s is the shift of the satellite peaks relative to the main hkl reflection.
In References [64-66], it was shown that although the strained volumes are coherent with the matrix, the intensities scattered by the strained volumes are incoherent; therefore, the intensities of the main and the satellite peaks are additive.The satellite intensities are usually a few percent of the intensities of the main peaks; therefore, the optimization of the main diffraction pattern and the satellite peaks can be done in separate steps.The two steps can be repeated subsequently until the procedure converges.Figure 6b,c shows zoomed parts of the measured (open circles) and CMWP-calculated (red lines) diffraction patterns of a Zircaloy-2 specimen neutron irradiated to the fluence of 13.1  1025 n/m2 at about 300( 25) C [58,67].The satellites are shown as blue lines.The shift of satellites is  =  0 √  ̅ , where s0 is a global shift parameter for the entire pattern and s scales for the different reflections, with the average contrast factor of <a>-type dislocation loops [25].The intensities of the satellites were adjusted to the heights of the related main peaks.Although the procedure was conceived for treating satellites produced by irradiation-induced dislocation loops, it is a general The actual value of the effective outer cut-off radius of dislocations, R e , as it was shown in Equation (12), depends on both the dislocation density and the dipole character of the dislocation arrangement.
A systematic investigation of the f (η) function revealed that when M drops into the range around or below unity along with large dislocation densities in the range of 10 16 × m −2 , the effective outer cut-off radius, R e , can reach small values of a few nm, approaching the lower limit of the continuum approach [66,67].The dislocation density in neutron- [57] or proton-irradiated [68,69] Zr alloys can reach a few times 10 16 •m -2 while the M value is M ≤ 1.In these diffraction patterns, the profiles broaden mainly in the tail regions of the profiles, while the FWHM change much less.The diffraction profiles of the 10.3 Bragg reflections of Zircaloy-2 proton irradiated to 2.3 dpa at 280, 350, and 450 • C are shown in Figure 8a in Reference [68].The figure shows that the tails are short when irradiation is at a high temperature, whereas they are very long when the irradiation is at a low temperature.TEM images of the same specimens are shown in Figure 3 in Reference [68]. Figure 3a in Reference [68] shows that at low temperatures, small dislocation loops form with very high density, whereas at high-temperature irradiation (see Figure 3c in Reference [68]) only a few large dislocation loops are present.When ρ is in the range of 10 16 •m -2 and M is around unity or smaller, R e can approach the lower limit of the continuum approach.During the optimization procedure, M can become uncertain in the sense that R e can approach extremely small, and ρ extremely large values.In order to avoid this kind of uncertainty, we introduced a lower limit, R c , for R e by replacing η = L/R e with η = L/(R e + R c ).
where R c is in the range of the lower limit of the length scale of the continuum approach [68,69].A systematic investigation of the effect of R c on the physical parameters is shown in Figure 7.We carried out this analysis on the diffraction pattern of a Zircaloy-2 specimen proton-irradiated to 2.3 dpa at 280 • C, in which the dislocation density was very large and the M parameter about unity [68].The figure shows ρ and M vs. R c in the range of 0 < R c < 30 nm.With R c decreasing to about R c 8 nm, M decreased, whereas ρ increased.Below R c 8 nm, there was a range of R c between about 6 ≤ R c ≤ 8 nm in which both ρ and M were stable around ρ 9 × 10 15 m −2 and M 0.75.The introduction of R c in the range of about 6 to 8 nm stabilized the physical parameters, especially the values of ρ and M. When R e is significantly larger than R c , then its value has no effect on the physical parameters.
radius, Re, can reach small values of a few nm, approaching the lower limit of the continuum approach [66,67].The dislocation density in neutron- [57] or proton-irradiated [68,69] Zr alloys can reach a few times 10 16 m -2 while the M value is M  1.In these diffraction patterns, the profiles broaden mainly in the tail regions of the profiles, while the FWHM change much less.The diffraction profiles of the 10.3 Bragg reflections of Zircaloy-2 proton irradiated to 2.3 dpa at 280, 350, and 450 C are shown in Figure 8a in Reference [68].The figure shows that the tails are short when irradiation is at a high temperature, whereas they are very long when the irradiation is at a low temperature.TEM images of the same specimens are shown in Figure 3 in Reference [68]. Figure 3a in Reference [68] shows that at low temperatures, small dislocation loops form with very high density, whereas at high-temperature irradiation (see Figure 3c in Reference [68]) only a few large dislocation loops are present.When  is in the range of 10 16 m -2 and M is around unity or smaller, Re can approach the lower limit of the continuum approach.During the optimization procedure, M can become uncertain in the sense that Re can approach extremely small, and  extremely large values.In order to avoid this kind of uncertainty, we introduced a lower limit, Rc, for Re by replacing  = L/Re with  = L/(Re + Rc).
where Rc is in the range of the lower limit of the length scale of the continuum approach [68,69].A systematic investigation of the effect of Rc on the physical parameters is shown in Figure 7.We carried out this analysis on the diffraction pattern of a Zircaloy-2 specimen proton-irradiated to 2.3 dpa at 280 C, in which the dislocation density was very large and the M parameter about unity [68].The figure shows  and M vs. Rc in the range of 0 < Rc < 30 nm.With Rc decreasing to about Rc 

Dislocation Density and Crystallite Size in Zr Matrix and Zr Hydrides in a Hydrated Zircaloy-4 Sheet Material
The performance of the CMWP procedure with the combination of the MC and the LM optimization methods was shown by the evaluation of the dislocation structure in a hydrated Zircaloy-4 specimen.The dislocation and stacking fault densities and the subgrain size were determined within the two coexisting Zr and ZrH phases.The spatial distribution of the volume fraction of ZrH in a hydrated Zircaloy-4 specimen was investigated using high-resolution scanning synchrotron experiments in References [70,71].The diffraction patterns were re-evaluated using the CMWP procedure described in Section 2.6.

Experiment
A hydride blister was grown on a Zircaloy-4 alloy platelet sectioned from a CANDU pressure tube.The alloy composition was Zr-1.49%Sn-0.2%Fe-0.11%Cr-2.5%Nb,all in wt %. [71].For the hydration procedure, see References [70,72].The schematic configuration of the hydrated platelet within the pressure tube, along with an optical micrograph of the hydrated blister, is shown in Figure 8a,b.High-angular-resolution X-ray microdiffraction was carried out at the 1-ID beamline at the Advanced Photon Source synchrotron at the Argonne National Laboratory [71].The schematic positions of the X-ray diffraction experiments on the hydrated platelet are shown in Figure 8c.X-ray diffraction scans were taken along three linear traces, one horizontal (red line) and two vertical (center green line and right blues line), with 16 and two times 11 positions, respectively.The beam size on the specimen was 300 × 300 µm 2 , illustrated approximately by the small red dot in the center vertical line in Figure 8c.
A hydride blister was grown on a Zircaloy-4 alloy platelet sectioned from a CANDU pressure tube.The alloy composition was Zr-1.49%Sn-0.2%Fe-0.11%Cr-2.5%Nb,all in wt %. [71].For the hydration procedure, see References [70,72].The schematic configuration of the hydrated platelet within the pressure tube, along with an optical micrograph of the hydrated blister, is shown in Figure 8a,b.High-angular-resolution X-ray microdiffraction was carried out at the 1-ID beamline at the Advanced Photon Source synchrotron at the Argonne National Laboratory [71].The schematic positions of the X-ray diffraction experiments on the hydrated platelet are shown in Figure 8c.X-ray diffraction scans were taken along three linear traces, one horizontal (red line) and two vertical (center green line and right blues line), with 16 and two times 11 positions, respectively.The beam size on the specimen was 300  300 m 2 , illustrated approximately by the small red dot in the center vertical line in Figure 8c.The diffraction patterns were recorded in a solid-state detector of 200 m pixel size with 2048 × 2048 pixels at a distance of 1.948 m from the sample.In order to obtain more reflections, the detector was shifted so that the beam hit the lower left corner of the detector, as shown in Figure 2a in Reference [71].The Debye-Scherrer rings of the reflections were integrated along azimuth sections of arcs in  2.5 degrees.The diffraction patterns were recorded in a solid-state detector of 200 µm pixel size with 2048 × 2048 pixels at a distance of 1.948 m from the sample.In order to obtain more reflections, the detector was shifted so that the beam hit the lower left corner of the detector, as shown in Figure 2a in Reference [71].The Debye-Scherrer rings of the reflections were integrated along azimuth sections of arcs in ±2.5 degrees.

Results
The diffraction patterns were evaluated for the dislocation density, the subgrain size, the dislocation arrangement parameter, and the stacking fault density, ρ, <x> area , M, and β in the fcc δ−ZrH and in the Zr matrix phases.Typical measured (open circles) and fitted (red solid lines) diffraction patterns corresponding to four different positions along the scanning lines are shown in Figure 9.In order to show the different reflections of very different relative heights, the intensities are in logarithmic scale.The patterns in Figure 9a,c,e,g show the entire measured range for the positions Line-1 Run-1, Line-1 Run-8, Line-2 Run-3, and Line-3 Run-11, respectively.Figure 9b,d,f,h shows zoomed ranges from 5.5 to 7.5 degrees, showing the detailed structure of the patterns where relatively large number of diffraction peaks were in complicated sequences and overlapped.
In the Line-1 Run-1 position (in Figure 9c, the top right edge position of the blister), the Zr peaks were much stronger than the δ-ZrH peaks.This ratio was reversed in the Line-1 Run-8 pattern, where the δ-hydride peaks were much stronger than the Zr peaks (in Figure 9c, the center top position of the blister).The Line-2 Run-11 position was at a distance of about 1.8 mm from the outer edge of the blister, but the δ-hydride peaks were still clearly visible.Although the Line-3 Run-3 position, at the outer edge of the blister, was a very different position compared to Line-2 Run-11, the two patterns indicated very similar ratios of the volume fractions between ZrH and Zr, and similar dislocation and subgrain structures.The lateral and depth profiles of the dislocation density, ρ, the stacking fault density, β and the area average mean crystallite size, <x> area , in the δ-ZrH are shown as a function of position in Figure 10.The light brown shade in the figures indicates the blister.The blue and the three red dots in Figure 10a show the approximate positions of the four diffraction patterns in Figure 9.The dislocation density in the Zr matrix varied between 1(± 0.15) × and 11(± 1.5) × 10 14 m −2 , whereas in the ZrH phase, it varied between about 10(± 2) × and 350(± 50) × 10 14 m −2 .subgrain structures.The lateral and depth profiles of the dislocation density, , the stacking fault density,  and the area average mean crystallite size, <x>area, in the -ZrH are shown as a function of position in Figure 10.The light brown shade in the figures indicates the blister.The blue and the three red dots in Figure 11a show the approximate positions of the four diffraction patterns in Figure 9.The dislocation density in the Zr matrix varied between 1( 0.15)  and 11( 1.5) 10 14 m −2 , whereas in the ZrH phase, it varied between about 10( 2)  and 350( 50)  10 14 m −2 .The experimental error of the stacking fault density in the -hydride was about  20%.The volume fraction of the -hydride was determined in Reference [71].It was shown to vary from about 80% in the center region of the blister to a few percent at a distance of about 2 mm from the edge of it.For details, see Figure 2 in Reference [71].The results showed that the dislocation density and the volume fraction were in close correlation with each other in the ZrH phase.At the largest volume fractions of ZrH, the dislocation densities reached extremely large values.The dislocation density in the Zr matrix was also correlated with the volume fraction of the ZrH phase.ZrH is well-known to exert large intergranular strains on the surrounding Zr matrix [73].These strains induce large dislocation densities in both phases and cause cracking of the Zr phase, as shown in Figure 2a in Reference [71].The experimental error of the stacking fault density in the δ-hydride was about ± 20%.The volume fraction of the δ-hydride was determined in Reference [71].It was shown to vary from about 80% in the center region of the blister to a few percent at a distance of about 2 mm from the edge of it.For details, see Figure 2 in Reference [71].The results showed that the dislocation density and the volume fraction were in close correlation with each other in the ZrH phase.At the largest volume fractions of ZrH, the dislocation densities reached extremely large values.The dislocation density in the Zr matrix was also correlated with the volume fraction of the ZrH phase.ZrH is well-known to exert large intergranular strains on the surrounding Zr matrix [73].These strains induce large dislocation densities in both phases and cause cracking of the Zr phase, as shown in Figure 2a in Reference [71].

Conclusions
The CMWP line profile analysis procedure was extended by implementing a special MC method in combination with the LM algorithm to provide the global optima for dislocation densities, dislocation arrangement or dipole character parameters, planar defect densities, subgrain size, and intergranular strains in polycrystalline or single-crystal materials.The procedure has the capacity to evaluate multiphase materials where the microstructures of the different component phases are determined at the same time.Taking into account the possible large number of optimization parameters, a sectioning was introduced where the physical, peak position, and peak height values are optimized in a cycle of six subsequent steps.The six-step cycle is repeated until convergence is reached.
Profile shape plays a decisive role in determining dislocation densities from the strain profile.The dislocation arrangement or dipole character parameter is one of the more sensitive parameters in the strain profile.The Krivoglaz-Wilkens strain function was modified in order to protect it from becoming unstable in cases when a large dislocation density is accompanied by a very strong dipole character.In such cases, the effective outer cut-off radius approaches the length limit of continuum theory.An editable hard limit was introduced into the Krivoglaz-Wilkens strain function in order to prevent the effective outer cut-off radius becoming smaller than the lower length limit of continuum theory.
We developed a procedure to handle the satellite peaks appearing on the flanks of the major Bragg peaks in polycrystalline materials.Such satellites are produced by e.g., neutron-or proton-irradiated Zr alloys [58,70] or by GP zones in supersaturated solid solutions [55,56,59].
The implementation of the MC procedure and the six-step cycle procedure was tested by evaluating the coexisting hcp Zr and fcc δ-ZrH phases in a hydrated blister on a Zircaloy-4 specimen.The two-dimensional lateral and depth profiles of the dislocation densities, the stacking fault densities, and the subgrain size were determined by evaluating position-sensitive high-resolution synchrotron X-ray diffraction patterns.The results were in good correlation with the mechanical and electron microscopy analysis done on the same specimen.

Figure 2 .
Figure 2.An intuitive schematic interpretation of the dislocation dipole character parameter, M=Re√.(a) Random dislocation distribution, where M  1.(b) Dislocations arranged in strong dipole configuration where M  1.The size of the regions of the material in (a) (red circle) and (b) (gray circle) and the number of plus and minus dislocations (up and down Ts) in the two regions are the same.The red arrows are for Re and the blue double arrows are for the average dislocation distance.

2. 1 . 1 .
The Density and the Arrangement Parameter of Dislocations

Figure 2 .
Figure 2.An intuitive schematic interpretation of the dislocation dipole character parameter, M = R e √ ρ (a) Random dislocation distribution, where M 1.(b) Dislocations arranged in strong dipole configuration where M 1.The size of the regions of the material in (a) (red circle) and (b) (gray circle) and the number of plus and minus dislocations (up and down Ts) in the two regions are the same.The red arrows are for Re and the blue double arrows are for the average dislocation distance.

Figure 3 .
Figure 3. Intensity distributions of the {110} reflections of a V-5Cr-5Ti alloy deformed by highpressure torsion at 4 GPa pressure to the shear strain of  = 37 (red curve) and after heating to 800 C (blue curve) with logarithmic intensity scale normalized to the maximum intensities and the FWHM.More details about material processing and the heat treatments can be found in Reference[40].

Figure 4 .
Figure 4. (a) Half-breadth of the parameter interval,    , as a function of the number of iterations, n.Up to the first no steps (horizontal red line),    stays constant.Between the steps no and n1,    decays exponentially (blue line).Beyond n1,    stays constant (horizontal green line).(b) Schematic drawing of the parameter values    vs. the random generated number,    , in the biased cubic probability function used in the special MC algorithm.

Figure 4 .
Figure 4. (a) Half-breadth of the parameter interval, ∆ i n , as a function of the number of iterations, n.Up to the first n o steps (horizontal red line), ∆ i n stays constant.Between the steps n o and n 1 , ∆ i n decays exponentially (blue line).Beyond n 1 , ∆ i n stays constant (horizontal green line).(b) Schematic drawing of the parameter values a i n vs. the random generated number, x i n , in the biased cubic probability function used in the special MC algorithm.

Figure 5 .
Figure 5. Systematic investigation of the LM and MC optimization procedures.(a) The WSSR vs. the initial value of the dislocation density,  initial , when LM alone (open blue triangles) and the LM and MC (open red circles) optimization procedures were used.(b) The final dislocation density and dipole character parameter values,  final and  final , vs. the initial, values,  initial and  initial , in logarithmic scales, calculated using either the LM alone (open blue triangles and green triangles) or the LM and MC optimization procedures (open red circles and black dots).The vertical black lines indicate error bars.All curves in the two figures are only to guide the eye.

Figure 5 .
Figure 5. Systematic investigation of the LM and MC optimization procedures.(a) The WSSR vs. the initial value of the dislocation density, ρ initial , when LM alone (open blue triangles) and the LM and MC (open red circles) optimization procedures were used.(b) The final dislocation density and dipole character parameter values, ρ final and M final , vs. the initial, values, ρ initial and M initial , in logarithmic scales, calculated using either the LM alone (open blue triangles and green triangles) or the LM and MC optimization procedures (open red circles and black dots).The vertical black lines indicate error bars.All curves in the two figures are only to guide the eye.

Figure 6 .
Figure 6.(a) Typical satellites in the lower intensity ranges of diffraction peaks in a proton irradiated Zircaloy-2 specimen.The vertical arrows point to the satellite peaks produced either by vacancy-(on the left side of the peaks) or interstitial-type dislocation loops (on the right side of the 00.2 reflection).The interstitial-type satellite peak is blending in with the broader 10.1 Bragg peak.(b) Measured (open circles) and CMWP-calculated (red line) zoomed part of the diffraction pattern of a Zircaloy-2 cladding material specimen [58] neutron irradiated to 13.1  10 25 n/m 2 .The blue line shows the satellite peak related to vacancy-and interstitial-type dislocation loops.(c) The same as in (a) for the first three Zr peaks.The alloy contains second-phase particles (SPPs) of different precipitates, as indicated in the (b) and (c) figures.

Figure 6 . 2 . 9 .
Figure 6.(a) Typical satellites in the lower intensity ranges of diffraction peaks in a proton irradiated Zircaloy-2 specimen.The vertical arrows point to the satellite peaks produced either by vacancy-(on the left side of the peaks) or interstitial-type dislocation loops (on the right side of the 00.2 reflection).The interstitial-type satellite peak is blending in with the broader 10.1 Bragg peak.(b) Measured (open circles) and CMWP-calculated (red line) zoomed part of the diffraction pattern of a Zircaloy-2 cladding material specimen [58] neutron irradiated to 13.1 × 10 25 n/m 2 .The blue line shows the satellite peak related to vacancy-and interstitial-type dislocation loops.(c) The same as in (a) for the first three Zr peaks.The alloy contains second-phase particles (SPPs) of different precipitates, as indicated in the (b) and (c) figures.2.9.Stabilizing the Fluctuations of the Physical Parameters When the Effective Outer Cut-Off Radius, Re, Approaches the Lower Limit of Continuum Theory

Figure 7 .
Figure 7. Systematic analysis of the effect of value of Rc in the CMWP procedure.The dislocation density, , and the dipole character parameter, M, vs. Rc in the range of 0 Rc 30.The vertical black lines indicate the error margin.

Figure 7 .
Figure 7. Systematic analysis of the effect of value of R c in the CMWP procedure.The dislocation density, ρ, and the dipole character parameter, M, vs. R c in the range of 0≤ R c ≤30.The vertical black lines indicate the error margin.Neutron diffraction patterns are usually measured in time of flight (TOF) mode.In TOF patterns, the intensities are given vs. d = 1/K.The software also works in I(K) representation.The TOF patterns have to be converted to I(1/K) = I(d*) representation (where d* = 1/d), which is a very simple task.

Figure 8 .
Figure 8.(a) Schematic configuration of the hydrated platelet sectioned from a CANDU pressure tube.The conventional coordinate system with the residual stress directions, radial, axial, and hoop strain/stress directions, along with the typical basal pole directions in the textured Zr tube material, is also indicated.(b) Optical micrograph of the hydrated Zircaloy-4 platelet.The darker region is the hydrated blister.(c) Schematic drawing of the X-ray diffraction positions along the hydrated platelet.

Figure 8 .
Figure 8.(a) Schematic configuration of the hydrated platelet sectioned from a CANDU pressure tube.The conventional coordinate system with the residual stress directions, radial, axial, and hoop strain/stress directions, along with the typical basal pole directions in the textured Zr tube material, is also indicated.(b) Optical micrograph of the hydrated Zircaloy-4 platelet.The darker region is the hydrated blister.(c) Schematic drawing of the X-ray diffraction positions along the hydrated platelet.

Figure 9 .
Figure 9. Measured (open circles) and CMWP-calculated (red lines) diffraction patterns (a), (c), (e), and (g) in the entire, and (b), (d), (f), and (h) in a zoomed angular range at four different depths in the ZrH blister.(a) and (b) are along Line-1 Run-1 at the top right edge, (c) and (d) are along Line-1 Run-8 at the center, (e) and (f) are along Line-2 Run-11 at the lower center edge, and (g) and (h) are along Line-3 Run-3 at the right middle edge of the blister, respectively.

Figure 9 .
Figure 9. Measured (open circles) and CMWP-calculated (red lines) diffraction patterns (a,c,e,g) in the entire, and (b,d,f,h) in a zoomed angular range at four different depths in the ZrH blister.(a,b) are along Line-1 Run-1 at the top right edge, (c,d) are along Line-1 Run-8 at the center, (e,f) are along Line-2 Run-11 at the lower center edge, and (g,h) are along Line-3 Run-3 at the right middle edge of the blister, respectively.

Figure 9 .Figure 10 .
Figure 9. Measured (open circles) and CMWP-calculated (red lines) diffraction patterns (a), (c), (e), and (g) in the entire, and (b), (d), (f), and (h) in a zoomed angular range at four different depths in the ZrH blister.(a) and (b) are along Line-1 Run-1 at the top right edge, (c) and (d) are along Line-1 Run-8 at the center, (e) and (f) are along Line-2 Run-11 at the lower center edge, and (g) and (h) are along Line-3 Run-3 at the right middle edge of the blister, respectively.

Figure 10 .
Figure 10.(a) The dislocation density, ρ, (b) the stacking fault density, β, and (c) the area average mean crystallite size, <x> area , in δ-ZrH as a function of lateral and depth positions in the ZrH blister.
The star in    * indicates that this consecutive value has not yet been accepted as a better value than  −1