Experimental Methods in Chemical Engineering: X-Ray Diffraction spectroscopy-XRD
Experimental Methods in Chemical Engineering: X-Ray Diffraction spectroscopy-XRD
DOI: xxx/xxxx
MINI-REVIEW
KEYWORDS:
XRD, nanoparticle, crystallinity, Debye-Scherrer method, limit of quantification
† This
is an example for title footnote.
This article has been accepted for publication and undergone full peer review but has not been
through the copyediting, typesetting, pagination and proofreading process which may lead to
differences between this version and the Version of Record. Please cite this article as doi:
10.1002/cjce.23747
This article is protected by copyright. All rights reserved.
2 Khan ET AL
TABLE 1 Comparing XRD with other analytical techniques; G= gas, L= liquid, S = solid
() = plane
A (101)
A = anatase
B = brookite
Intensity
A (200)
Crystal R= rutile
A (004)
A (211)
A (105)
A (204)
X-rays
R (121)
B (110)
A (220)
A (215)
A (116)
reflections Fourier e- density Atomic
transform map model
20 25 30 35 40 45 50 55 60 65 70 75 8
2θ (degree)
FIGURE 2 Single crystal XRD reflections. [20] A crystal
diffracts an x-ray beam and a film or electronic detector record
reflections on a 2-D surface. A Fourier transform of the pattern FIGURE 3 XRD diffractogram of TiO2
and intensity of the reflections produces an electron density
map from which we derive an atomic model This approach requires that the component has the essential
absorption coefficient when a standard is unavailable:
detector or film records the reflections as the x-ray source com-
𝑘B 𝐼B
pletes an arc over the sample. The position and intensity of the 𝑋B = = 0.14, brookite (3)
reflection relate to the identity and position of the atoms in the 𝑘A 𝐼A + 𝑘B 𝐼A + 𝐼R
unit cell [21] (Figure 2).
𝐼R
Powder x-ray diffraction has a wider applicability com- 𝑋R = = 0.05, rutile (4)
𝑘A 𝐼A + 𝑘B 𝐼A + 𝐼R
pared to single crystal analysis (SCXRD), but peaks overlap
substantially more in the diffractogram compared to SCXRD Although the peak intensity of the anatase is almost 20× higher
of multicomponent samples, which obscures the position and than rutile, its mass fraction is only 6× higher. The peak inten-
intensities of the diffraction maxima (Figure 3). [22,23] In pow- sities correlate with composition, the peak position identifies
der diffraction XRD, a large number of individual crystallites the phase, unit cell parameters, group spacing, stress-strain
intercept the incident x-rays and the individual beams of inten- analysis, and crystal system. The software completes the find
sity become cones; the detector individual spots generate rings and match procedure and identifies crystallographic planes—
of varying intensity. We compare the peak intensity and d- ((101), for example). The peak width (FWHM: full width
spacing against the JCPDS library to identify crystals. at half maximum) correlates with crystal size, lattice distor-
The mass fraction of each phase, 𝑋𝑖 , in the mixture is based tion, and structural dislocation. The Williamson-Hill plot of
on the Zhang and Benfield equation that considers the peak 𝛽 cos 𝜃∕𝜆 versus sin 𝜃 gives a straight line with the slope equal
intensity, 𝐼𝑖 , and a coefficient derived for each component, 𝑘𝑖 to 𝜂∕𝜆 and the intercept is 1∕𝑑p :
(𝑘A = 0.884, 𝑘B = 2.721, 𝑘R = 1) [24] : 𝛽 cos 𝜃 1 𝜂 sin 𝜃
= + (5)
𝑘A 𝐼A 𝜆 𝑑p 𝜆
𝑋A = = 0.81, anatase (2)
𝑘A 𝐼A + 𝑘B 𝐼A + 𝐼R where 𝛽 = FWHM (in radians), 𝑑p is the effective crystallite
size, and 𝜂 is the effective strain. We derive the crystallite size
TABLE 2 TiO2 crystallographic parameters derived from the x-ray diffractogram (Figure 3). We assign the peak with the highest
intensity 1000 and other peak heights are ratios of this value † Full width at half maximum
from Scherrer’s equation: extremely bright and monochromatic and produces diffrac-
0.9𝜆 tograms with minimal peak overlap. Electrons enter a stor-
𝑑p = (6)
𝛽 cos 𝜃 age ring several kilometers in circumference at high vacuum
The anatase crystallite size (Figure 3) is 20 nm, which is a lit- (10 Torr to 12 Torr) to minimize particle loss by collision
tle smaller than the rutile (23 nm) but larger that the brookite with residual gas atoms. In the storage ring, a magnetic
(17 nm) based on a copper metal target with 𝜆 = 0.154 18 nm. field changes the direction of the high speed electrons so
Under reaction conditions, catalysts calcine and become they emit electromagnetic radiation because of the angular
more crystalline, which increases the 𝑑p and thereby the acceleration—synchrotron radiation (SR)—that ranges from
intensity as the peaks become narrower (Figure 4). [25,26] In from microwaves (𝜆 > 1 m) to hard x-rays (0.05 nm). [30] A
situ XRD measures changes in catalyst and other materials synchrotron beam line incorporates collimating mirrors (col-
structure during synthesis and reaction to identify structural limates the beam in vertical direction into parallel light to
changes associated with phase transition and chemical reac- improve the energy resolution [31] ), slits, focusing mirrors, and
tions Figure 5. Si (111) or germanium (220) double crystal monochroma-
It has become straightforward to follow changes as a func- tors to create tunable monochromatic x-rays and focus the
tion of temperature, pressure, and gaseous environment while beam in a horizontal direction. [32] The characteristics of syn-
maintaining the signal quality throughout the duration of the chrotron radiation that make it suitable for XRD include [33] :
test. [27,28] he following example shows the usefulness of the (a) bright x-rays that are 100× to 1000× more intense than con-
powder diffractometer with a hot stage that allows for in-situ ventional laboratory analyzers; (b) highly collimated beams
heating and controlled ambient. In this case, we anneal a thin to increase resolution; (c) a wide energy spectrum; and (d) a
bilayer of Al (≈ 10 nm) and Ag (≈ 200 nm) that resides on short pulsed time structure. These characteristics yield high
a SiO2 in an oxygen environment (for safety, Ar is the carrier signal/noise ratios and high angular resolution, while minimiz-
gas). ing peak overlap and improving peak positioning (Figure 7).
At the onset of annealing, Al diffuses into the Ag to form Furthermore, the 𝐾𝛼2 and 𝐾𝛽 diffraction peaks are absent.
a solid-solution that changes the lattice constant (Figure 6) At
725 ◦C, the Ag peak gradually shifts to the original direction
but less at 500 ◦C. This confirms that the high temperature 3 APPLICATIONS
anneal accelerates Al diffusion through Ag films. Also, the
(111) peaks are stronger at 725 ◦C than 500 ◦C at all time peri- Researchers apply x-ray diffraction across a wide range of
ods, which confirms that a higher temperature enhances the materials in dozens of scientific fields (Figure 8): metals and
texture of Ag films. [29] alloys, clays and minerals, catalysts, cement, ceramics, plas-
Synchrotron radiation, generated as charged particles accel- tics, composites, corrosion products, fly ash, asbestos, solar
erate in a straight-line path or travelling in a curved path, is cells, films, and semiconductors. Carbon-based materials—
graphene, diamond, carbon nanotubes, carbon nanobuds, car-
bon nanofoams, and activated carbon—are finding application
VHP
Spent
Calcined
VPP
Precursor
VPP
VPP
VHP
VHP
VHP
VHP
VHP
*
*
20 22 24 26 28 30 32 34
2 (degrees)
20 30 40 50 60 70 80
2 (degrees)
FIGURE 4 XRD diffractograms of vanadyl hemihydrate (precursor), calcined vanadium pyrophosphate (VPP), and equilibrated
VPP after over 12 000 h of continuous operation at 400 ◦C [25,26]
Al Al
AlOx
SiO2
O2
SiO2
Ag Ag(Al)
Forensic sciences
Corrosion • Trace analysis
• Solids phase composition • Qualitative phase analysis &
• Identifies, quantifies compounds composition
• Recognize samples, diluents &
Pharmaceutical adulterant
XRD • Identification of torn specimens
• Resolves complex structures
• Physicochemical properties analysis
• Identifies drug polymorphism
Microelectronics
• Excipient analysis
• Materials characterization
• Temperature gradient & moisture
• Topography defects and imaging
analysis
ŶĂŶŽƐŚĞĞƚƐ
>ŝŝŽŶ
ŝƌƌĂĚŝĂƚŝŽŶ
ƉŚŽƚŽĚĞŐƌĂĚĂƚŝŽŶ
ĨĂďƌŝĐĂƚŝŽŶ
ƉŚŽƚĐĂƚĂůLJƚŝĐ
,Ϯ ĂĐƚŝǀŝƚLJ
ƉƌŽĚƵĐƚŝŽŶ
ĐŽŵƉŽƐŝƚĞƐ
ŐƌĂƉŚĞŶĞŽdžŝĚĞ
ĚĞƉŽƐŝƚŝŽŶ
ĞĨĨŝĐŝĞŶĐLJĨŝůŵƐ
ĐĂƚĂůLJƐƚƐ
ZĂŵĂŶ
ŵĞĐŚĂŶŝƐŵ ƐŽůĂƌ
ĐĞůůƐ
ĐŽŶǀĞƌƐŝŽŶ ŶĂŶŽĐƌLJƐƚĂůƐ
ĚŝĞůĞĐƚƌŝĐ
ƉƌŽƉĞƌƚŝĞƐ
ĂƋ ƐŽů͛Ŷ ŚĞƚĞƌŽŐĞŶĞŽƵƐ
ĐĂƚĂůLJƐƚƐ ĚĞƐŝŐŶ ŵŝĐƌŽƐƚƌƵĐƚƵƌĞ
ĂĐŝĚ ĞŵŝƐƐŝŽŶ
ĞŶĞƌŐLJ
ƚƌĂŶƐĨĞƌ
FIGURE 9 XRD bibliometric map generated by VOSviewer. [36,37] The data base consists of the top 107 keywords of the 10 000
most cited articles in WoS (2016-2017). The size of the circle represents the number of occurrences of the keyword in the 10 000
articles, while related research is grouped into clusters of the same colour. nanoparticles appears in 1493 articles followed by
performance (1054), oxidation (416), degradation (597), and adsorption (659). The least frequent of the top 100 keywords appear
fewer than 130 times: sol gel method (117), supercapacitor (121), heterogeneous catalysts (115), sensitized solar cells (116),
and sorption (126)
! "
FIGURE 10 XRD scans (𝜃 − 2𝜃) of Ag on SiO2 and PEN FIGURE 11 XRD data of encapsulated Ag layers annealed at
substrates 625 ◦C and 725 ◦C
Standards are 100 % crystalline to minimize order/dis- (in radians). For a displacement of 𝑠 = 0.15 mm, the
order variability of amorphous materials. [48] Only the peak shift is 0.08◦ at 2𝜃 = 28.4◦ and 𝑅 = 200 mm.
(111) and (200) reflections appear on the Ag thin film Algorithms, which compensate for sample displacement
diffractograms over SiO2 and the polymer substrate PEN errors, require internal calibration standards. Zero back-
(Figure 10). The major reflections of PEN substrates ground sample holders and parallel-beam optics further
appear at 37◦ and 41◦ (Figure10). Relative intensities of minimize displacement errors.
(200) that are normalized to the intensities of the (111) iv Sample transparency error: X-ray penetration depth
peak were calculated for Ag on SiO2 as I200/I111 = 4.69 depends on the specimen mass absorption coefficient,
and PEN as I200/I111 = 13.19. These values indicated 𝜇, and the incident angle of the x-ray beam. Conse-
that Ag on SiO2 has more (111) planes parallel to the sur- quently, diffracted x-rays arrive from different points,
face than Ag on PEN. The full width at half maximum (dotted lines versus solid lines in Figure 12(B)) and
(FWHM) is inversely proportional to the grain size of the introduce peak position errors and peak asymmetry—
thin film. This implies that the average grain size of Ag the most common error for organic and low absorbing
on SiO2 is higher than on PEN (Figure 10, FWHMPEN (low atomic number) specimens. Thin specimens and
= 0.5341◦ and FWHMSiO2 = 0.5092◦ ). parallel beam optics minimize transparency error.
Figure 11 is another example of the preferred orienta-
v Plate sample error: When the entire surface of the flat
tion during sample processing. X-ray diffraction spectra
sample cannot lie on the focusing circle (Figure 12(C))
(𝜃 − 2𝜃) of Ag (100 nm)/Al (30 nm) bi-layers annealed
the peaks are broader and asymmetric (at low 2𝜃 angles).
at 625 ◦C and 725 ◦C in NH4 form an aluminium-
Narrow divergence slits with a shorter beam and parallel
encapsulation layer on the surface. The as-deposited
beam optics reduce this error.
sample shows a strong (111) texture that increases for
the high temperatures anneal at 725 ◦C. vi Counting statistics: Statistical characterization of pow-
ders is best when the particle size is less than 10 µm.
ii Peak shape: Peaks of amorphous materials are typically
Diffractograms of large crystallite sizes and non-random
broad and asymmetrical, which compromises quantify-
orientations lead to peak intensity variations that are
ing how much is present.
incongruent with an ideal powder (measured for many
iii Sample displacement: If the sample moves (dashed line crystallites randomly oriented) and thus are poorly iden-
above green block Figure 12(A)), the diffracted lines tifiable with reference patterns in the Powder Diffraction
miss the detector aperature (solid lines versus dashed File (PDF) database.
lines). As a result, the instrument incorrectly reports the
vii Axial/vertical divergence error: Soller slits, capillary
peak position. The displacement error, 𝑠, is a function of
lenses, and decreasing the vertical opening of the
the 𝜃 and the goniometer radius, 𝑅:
counter slit maximize the diffracted intensity and reduce
2𝑠 cos 𝜃
Δ2𝜃 = − (7) peak assymetric broadening due the divergence of the
𝑅
x-ray beam in the plane with the sample. [49]
viii Microabsorption: This error stems from differences in interact with materials (clays and zeolite minerals, for
the interactions of each material with the x-ray beam, example) and alter their structure.
volume fractions of the components (large particles not
crystallites), and on die geometrical peculiarities of their
distribution. [50] Complex composites such as cements/- 4.3 Detection limits
concretes, coal combustion by-products (CCBs), and Following the rules of a Poisson distribution, for 𝑁 counts at
geologic materials are among the materials that are 2𝜃 the absolute and relative standard deviation 𝜎 and 𝜎rel are:
susceptible to microabsorption errors. Each material √
absorbs x-ray radiation depending on its linear absorp- 𝜎= 𝑋 (8)
tion coefficient for the particular wavelength (energy). √
For example, CCBs have linear absorption coefficients 𝑋 1
𝜎rel = =√ (9)
ranging from 81 cm−1 for quartz to 1153 cm−1 for mag- 𝑋 𝑋
netite using copper radiation. Grinding and milling to a The limit of detection (LOD) of a particular reflection is:
smaller particle size reduces this error.
𝑋ref lection > 𝑋background > 3𝜎background (10)
For example, measuring an 𝑋max of 10 000 counts, the 𝜎rel =
4.2.2 Instrumental errors 0.01, corresponding to a relative error of 1 %, the counting sta-
i Sample fluorescence: Instrumental error in which inci- tistical error. (Note: Suppose we have a background of 100
dent x-ray radiation excite electrons of certain elements counts and a small hump of 120 counts. This cannot be clas-
(V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Ga). The excited elec- sified as a reflection for the reason that 3𝜎background = 30 is
tron emits a characteristic x-ray and the detector will obviously higher than 20. The only solution is to increase the
record signals that are close to Cu K𝛼—increasing back- measurement time to improve the peak background ratio.)
ground noise but not position or intensity. Monochro-
mators, 𝛽 filters, energy sensitive detectors, and pulse
4.4 Limit of quantification
height distribution levels (PHD) reduce the fluorescence
contribution to the signal. In the analysis of multiphase samples, an important considera-
tion is the question of the lower limit of quantification (LOQ).
ii Equipment misalignment: The error of the “zero” 2𝜃 What is the smallest amount of a given phase that can be iden-
position, ie, the offset of the instrument, introduces a tified by XRD? The LOQ depends on preferred orientation,
systematic peak position error that is proportional to matrix effect crystal symmetry, peak overlap, and amorphous
2𝜃. [51] The source, sample, and detector must all be content. Generally, the LOQ of a phase must be determined by
perfectly in line at 0◦ 𝜃. a calibration curve.
Practitioners cite a lower limit of 5 %, but this value
4.2.3 Compositional variations errors varies with the composition of the constituents. Catalysts, for
example, are often deposited on high surface area substrates to
i Grinding: Excessive grinding induces changes amor-
minimize the mass of costly metals. To demonstrate the effect
phism, strain, decomposition via local heating, and loss
of the substrate on the LOQ, we prepared five samples of Ru on
of volatile components. Some clays, zeolites, and engi-
TiO2 , activated carbon, and zeolite (ZSM-5). Ruthenium(III)
neered materials are sensitive to low-temperature dam-
chloride (45 % to 55 %) and activated charcoal (DARCO R —
age. Shatter boxes or ball-mills produce a tail’ of very
100 mesh particle size, powder) were purchased from Sigma-
fine particles that broaden the diffraction peaks. Tung-
Aldrich. ZSM-5, ammonium (powder, S.A 400 m2 g−1 , 30:1
sten carbide and other brittle grinding media contam-
mole ratio SiO2 :Al2 O3 ) were obtained from Alfa Aesar. Vena-
inate the sample and introduce additional peaks. Non-
tor forwarded the titanium oxide (TP Hombikat, 100 µm).
percussive techniques, like a mortar and pestle, reduce
We followed a wetness impregnation method in which we
these errors.
first added water to determine the total pore volume. We dis-
ii Irradiation: Incident x-rays can interact with samples solved the ruthenium(III) chloride in the calculated volume.
and change the composition, particularly for organic For a more intimate contact between the active phase and the
compounds. X-rays change the colour or cloud inorganic catalyst, we mixed the liquid and the solid in a rotary evap-
compounds without changing the diffractogram. orator for 3 h at 70 rpm. After increasing the temperature by
5 ◦C every 30 min until 85 ◦C, we decreased the pressure to
iii Environment: Elevated temperatures expand solids,
300 mbar for 2 h to dry the catalyst. A furnace calcined the
inducing stress and strain. Water and other liquids
Tube Tube 00 00 00 00 00 00 00
Receiving slits Receiving slits 00 00 00 00 00 00 00
00 00 00 00 00 00 00
00 00 00 00 00 00 00
0000000
00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00
00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00
000 000 000 000 000 000 000 000 000 000 000 000 000 000 000 000 000 000 000
(c) 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00
Detector Sample
Sample
Tube
Receiving slits
00 00 00 00 00 00
000 000 000 000 000 000
00 00 00 00 00 00
000 000 000 000 000 000
00 00 00 00 00 00
00 00 00 00 00 00
000 000 000 000 000 000
000000
Sample
FIGURE 12 XRD instrument error: A, sample displacement error; B, sample transparency error; and, C, flat sample error
TABLE 3 Ru mass fraction on supports: The second column to see more RuO2 on the surface. The Ru shoulder on the TiO2
summarizes concentration based on the recipe and the third is is more evident for the sample with a mass fraction of 5 %
based on nuclear activation analysis. The last column is the (Figure 14). The surface area of ZSM-5 is 240 m2 g−1 so per-
uncertainty in the nuclear activation. The Ru/AC concentration haps more of the Ru enters the structure. [52] AC is composed
is high since the carbon is undetected. of carbon (graphite) layers of fused hexagons held by weak
van de Waals forces and surface area can exceed 1000 m2 g−1 .
Catalyst Wet deposition Nuclear Δ Interestingly, XRD readily detects the 1.5 % Ru and the peaks
g g−1 , % g g−1 , % g g−1 , % grow for the sample with a mass fraction of 3 % (Figure 15).
At 2𝜃 between 65◦ to 70◦ the two peaks are better resovled at a
Ru/ZSM-5 1.5 1.38 0.06
mass fraction of 3 % compared to 1.5 %. SEM-EDX analysis of
Ru/AC 1.5 32.2 1.3
these samples showed that the RuO2 interacted little with the
Ru /AC 3 46.2 1.8
AC and rather formed large independent crystals. XRD anal-
Ru/TiO2 5 4.76 0.19
ysis alone was incapable of differentiating between a species
Ru/TiO2 1.5 2.01 0.08
deposited on a support and the same species present as a seg-
regated phase in the same sample. SEM-EDX analysis locates
the different phase compositions in a small specimen and is
solid: first it maintained the temperature constant at 120 ◦C for
complementary to XRD.
4 h to remove traces of water followed by a 2.5 ◦C min−1 ramp
to 600 ◦C, then a 4 h hold in air.
To ensure that the Ru remained on the supports after syn- 4.5 Crystallinity
thesis, we measured the concentration with a SLOWPOKE
nuclear reactor (Table 3). Nuclear bombardment detects all ele- Sample crystallinity is the ratio of the total intensities in the
ments from F to higher molar mass (except Pb because of its diffraction pattern, 𝑁net , and the sum of all the measured inten-
very short half life). For this reason, the analysis reports high sities, 𝑁tot , including the amorphous part and air scatter, 𝑁scat .
mass fractions of Ru for the AC samples. Otherwise, agree- The latter must be determined separately (measuring a zero
ment between the expected value based on the recipe agrees background holder) and then subtracted:
with that measured value by nuclear bombardment. 𝛴𝑁net
𝐶 = 100 (11)
The minimum concentration that XRD detects depends on 𝛴𝑁tot − 𝛴𝑁scat
the nature of the support and distribution (Figure 13): a mass The ratio between the pure crystalline phase and the amor-
fraction of 1.5 % RuO2 on ZSM-5 is undetectable, while a phous phase of the sample, 𝐶, depends on the peak position of
slight shoulder is evident at 38◦ on the TiO2 sample. TiO2 has
a very compact tetragonal crystalline structure and we expect
6 CONCLUSIONS
[3] R. Das, E. M. Ali, S. B. A. Hamid, Rev. on Adv. Mater. [23] D. Kenneth, E. Cheung, Chem. Soc. Rev. 2004, 33, 526.
Sci. 2014, 38, 95.
[24] H. Zhang, J. Banfield, The J. Phys. Chem. B 2000, 104,
[4] T. L. Alford, L. C. Feldman, J. W. Mayer, Fundamentals 3481.
of Nanoscale Film Analysis, Springer 2007.
[25] G. S. Patience, R. E. Bockrath, J. D. Sullivan, H. S.
[5] G. S. Patience, Can. J. Chem. Eng. 2018, 96, 2312. Horowitz, Ind. & Eng. Chem. Res. 2010, 46, 4374.
[6] S. Strandman, R. Vachon, M. Dini, S. Giasson, X. Zhu, [26] G. S. Patience, R. E. Bockrath, Appl. Catal. A: Gen. 2010,
The Can. J. Chem. Eng. 2018, 96, 20. 376, 4.
[7] C. Han, H. Liu, L. Zhang, J. Deng, Y. Luo, Can. J. Chem. [27] K. Harris, M. Tremayne, B. Kariuki, Angewandte Chemie
Eng. 2017, 95, 543. Int. Ed. 2001, 40, 1626.
[8] F. Tari, M. Shekarriz, S. Zarrinpashne, A. Ruzbehani, [28] B. Kariuki, K. Psallidas, K. Harris, R. Johnston, R. Lan-
Can. J. Chem. Eng. 2017, 95, 1632. caster, S. Staniforth, S. Cooper, Chem. Commun. 1999,
17, 1677.
[9] Y. Dai, Z. Fei, X. Xu, X. Chen, J. Tang, M. Cui, X. Q.
and, Can. J. Chem. Eng. 2016, 94, 1140. [29] Y. Wang, Diffusion Barrier Formation and Characteri-
zation for Silver Metallization, PhD thesis, Arizona State
[10] G. Fan, R. Lin, Z. Su, X. Lin, R. Xu, W. Chen, Can. J. University, Tempe, USA 2001.
Chem. Eng. 2017, 95, 717.
[30] Z. Dauter, J. Mol. Struct. 1996, 374, 29.
[11] H. Nguyen, T. Dang, N. Nguyen, H. Nguyen, N. Dinh,
The Can. J. Chem. Eng. 2018, 96, 832. [31] M. Isshiki, Y. Ohishi, S. Goto, K. Takeshita, T. Ishikawa,
Nucl. Instruments Methods Phys. Res. Sect. A: Accel.
[12] J. Hu, X. Yang, J. Yu, G. Dai, Can. J. Chem. Eng. 2016, Spectrometers, Detect. Assoc. Equip. 2001, 467–468,
94, 340. 663.
[13] A. Farjoo, S. M. Kuznicki, Can. J. Chem. Eng. 2016, 94, [32] T.-Y. Yang, W. Wen, G.-Z. Yin, X.-L. Li, M. Gao, Y.-L.
2219. Gu, L. Li, Y. Liu, H. Lin, X.-M. Zhang, B. Zhao, T.-K.
[14] S. Mondal, A. Elkamel, D. Reinalda, K. Wang, Can. J. Liu, Y.-G. Yang, Z. Li, X.-T. Zhou, X.-Y. Gao, Nucl. Sci.
Chem. Eng. 2017, 95, 1993. Tech. 2015, 26, 1.
[15] M. Geramian, D. Ivey, Q. Liu, T. Etsell, The Can. J. [33] H. Cheng, C. Lu, J. Liu, Y. Yan, X. Han, H. Jin, Y. Wang,
Chem. Eng. 2018, 96, 49. Y. Liu, C. Wu, Prog. Nat. Sci. Mater. Int. 2017, 27, 66.
[16] R. Multani, K. Waters, The Can. J. Chem. Eng. 2018, 96, [34] E. Cobas, M. Cabrera, L. Pardo, L. Montero, Mater.
1185. 2016, 9, 1.
[17] M. A. Rahimi, A. S. A. Ramazani, H. A. Alijanvand, [35] P. Fewster, Acta Crystallogr. Sect. A, Foundations Adv-
M. H. Ghazanfari, M. Ghanavati, Can. J. Chem. Eng. naces 2014, A70, 257.
2017, 95, 83.
[36] Clarivate Analytics, Web of Science Core [52] M. Menéndez, J. Herguido, A. Bérard, G. S. Patience,
Collection, accessed on 25 June 2018, Can. J. Chem. Eng. 2019, 97, 2383.
http://apps.webofknowledge.com 2018.