Evaluation of the Structural Deviation of Cu/Cu 2 O Nanocomposite Using the X-ray Diffraction Analysis Methods

: We successfully synthesized Cu/Cu 2 O nanocomposites using the wet chemical synthesis method. All X-ray diffraction (XRD), Reference Intensity Ratio (RIR), and Rietveld reﬁnement methods conﬁrmed that the compounds Cu and Cu 2 O are free of impurities. Scanning Electron Microscope (SEM) and Transmission electron microscopy (TEM) images show the morphology and interactions of Cu and Cu 2 O in the structure. The formation mechanism is also explained by ﬁve stages: precursor, nucleation, growth, aging, and reduction. The changes in crystallization parameters under variations in reaction temperature (T v ) and stirring speed (S v ) were conﬁrmed by agreement with the XRD database. The lattice constant in the crystal of nanocomposite increases with rising temperature in the reaction, leading to unit cell expansion, while increasing the stirring—rate leads to a random size distribution of the lattice parameter. Due to the imperfect growth of the crystal, the induced crystallite size was calculated using the Williamson-Hall model, and the precise lattice parameter values were calculated using the Nelson-Riley function.


Introduction
Metal and semiconductor oxides are used in a variety of energy applications such as batteries, supercapacitors, solar cells, and photoelectrochemical cells. Among these, Cu 2 O is attractive because of its optical characteristics, availability, abundance, and environmentallyfriendly material, which has advantages such as low synthesis cost and ease of synthesis with various methods. For these reasons, it is widely used in sensors, photo-catalysts, and solar cells [1][2][3]. Nevertheless, Cu 2 O still has drawbacks in optical applications because it has a higher e − -h + recombination rate than other materials [3]. The incorporation of Cu and Cu 2 O in hetero-junctions can enhance the photo-catalytic properties, and reaction efficiency. Moreover, it can reduce the e − -h + recombination due to the photoelectron transferability of Cu in that composite material [4].
There are various methods reported for Cu 2 O nanoparticles synthesis, such as the hydrothermal, the electro-deposition, the wet chemical method, and the thermal decomposition [5][6][7]. The wet chemical synthesis method is one of the methods that allows for adjustment of the thermodynamic and kinetic parameters to control crystal parameters such as the shape, the size, and the composition of nanoparticles [8]. Hence, the electronic, the optical, and the reactive properties of their surfaces can also be changed. Studies of controlling crystalline parameters have demonstrated their influence on the extraordinary magnetic, electrical, optical, and mechanical properties of nanomaterials [9].
Crystals should grow perfectly in all directions. No crystal, however, is perfect due to its finite size [10]. Furthermore, under different conditions of synthesis, the crystallographic parameters of the nanoparticles will grow in different directions, thus making the crystal imperfect. Deviation in crystal growth under different conditions leads to broadening of diffraction peaks compared with those from a perfect, crystalline structure. Lattice strain and crystal size are the two main parameters used to evaluate the misalignment of crystal growth [10]. Therefore, a clear understanding of the crystal size changes under different synthesis conditions will assist in the control of crystal size and shape, and further, in the control of the application of structure-based nanoparticles. In this study, the XRD technique was used to analyze the change in crystallographic parameter. In addition, we evaluated the effects of stirring speed (S v ) and reaction temperature (T v ) on crystallographic parameters of Cu/Cu 2 O nanocomposites.
Additionally, the Williamson-Hall method was used to recalculate the induced-lattice strain and induced-crystallite size; the Nelson-Riley function was applied to determine the precise lattice parameter; and the Rietveild refinement method was used to confirm the composition of our sample.

Materials and Methods
The chemicals were used in this study include copper sulfate pentahydrate (CuSO 4 ·5H 2 O, >99%), polyethylene glycol (PEG, MW 6000), L-ascorbic acid (C 6 H 8 O 6 , >99%), and sodium hydroxide (NaOH, >97%). All chemicals were from the Sigma Aldrich Company and did not need to be purified when used during the experiment. Deionized (DI) water was used to dilute all chemicals.

Synthesis of Cu/Cu 2 O PNCs
First, 50 mL precursor solution was obtained by mixing 1.0 mol/L CuSO 4 ·5H 2 O and 1.0 g/L PEG (ratio 1:1). Second, 25 mL NaOH 1M was slowly added under the stirring condition for 15 min. Third, C 6 H 8 O 6 (0.25 mmol) was added drop-by-drop into the solution, stirring for 30 min. The synthesis was performed at S v = 250 rpm and T v = 40 • C (ST). From the synthesis conditions of the ST sample, the changes in T v and S v were applied to evaluate the change in the structure of Cu/Cu 2 O PNCs. The included treatments were T1 = 60 • C, T2 = 80 • C, T3 = 100 • C, S1 = 350 rpm, S2 = 450 rpm, and S3 = 550 rpm. The experimental conditions and codes are given in Table 1. After the reaction, the powders were extracted by centrifuging the solutions at 5000 rpm for 20 min. The obtained powders were washed three times sequentially with DI water, ethanol, and then DI water again. Finally, the powders were dried on a hot plate at 95 • C until the mass remained constant.

Characterizations
The crystallite structure and chemical composition were measured as the powder XRD patterns obtained by the X-ray diffractometer under Cu-Kα radiation with λ = 1.54056 Å (Bruker D8 Focus, Model: PANalytical X' pert PRO, PANalytical, Almelo, The Netherlands). Samples (powder) were placed on a sample holder, and data were obtained using Cu Kα monochromatic radiation according to Bragg's Law convention. The data of the samples were also taken in the range of 10-80 • angles using the step size of 0.02 (2θ) per second. The surface morphology of samples was analyzed by SEM (Model JEOL JSM-7400F, Yokogushi, Tokyo, Japan) and HRTEM (FEI Company, Hillsborough, OR, USA). The diagram structure of samples was plotted by VESTA software. The Rietveld refinement method was used to scan the composition of samples by MAUD software (version 2.61, University of Trento, Trentino, Italy), and phase quantification was investigated via the reference intensity ratio (RIR) method (Chung 1974) using the MATCH! software (version 3.13, Dr. Holger Putz, Crystal Impact, Bonn, Germany). In the MATCH! software, the content of each crystalline phase was calculated as Equation (1) [11]: where C i is the percentage content of ith phase, I i is the integral intensity of ith phase, and k i is the coefficient (RIR) for ith phase obtained from the PDF database.

The Formation of Cu/Cu 2 O Powder Nanocomposite
The formation of Cu/Cu 2 O powder nanocomposite can be divided into five stages: precursor, nucleation, growth, aging, and reduction ( Figure 1). First, the precursor chemical (CuSO 4 ·5H 2 O) was dissolved in DI water; it was dissociated to yield free Cu 2+ ions in the solution. PEG was added as the capping agent to cover the ions; the capping agent increases the viscosity of the solution. As a result, the diffusion and movement of ions in the solution becomes more difficult; because there is a layer of sorbent (capping agent) on the surface of the particles, the capping agent prevents the agglomeration of the particles. Furthermore, with PEG encapsulation, it is difficult for the nanoparticles to self-align or gather together to form clusters. Accordingly, the resulting nanoparticles are smaller in size than the formed nanoparticles in solution without PEG [12]. Second, in the nucleation stage, the Cu 2+ ions react with OH − ions to form Cu(OH) 2 precipitation, the neighbor ions bond to each other and form the cluster. Over time, the increase in OH − ions causes more Cu(OH) 2 precipitates to form. During the nucleation process, the solution changes from light blue to dark blue and finally to milky blue.
Third, the growth phase of NPs has two mechanisms: the reaction on the surface and the diffusion of monomers onto the surface [13]. The capping agent not only interacts with the precursor ions but also has an important role in the growth phase. After the formation of NPs, the capping agent protects the surface of those particles by creating a stable cover layer. Moreover, the physicochemical properties of NPs mostly depend on the interaction between the NP and capping agent due to its ability to modify the interaction between the external system and NPs [14]. Additionally, in the growth phase, clusters gather together so that the size of the NPs increases. Fourth, aging is commonly used in the synthesis of NPs, and macro or micro-changes in the solid phase can occur during this process [15]. At the end of the aging period, the size of the NPs does not grow any further. Finally, L-ascorbic acid is used as a reducing agent to reduce Cu(OH) 2 so that it precipitates to form Cu 2 O particles; this process is shown in Figure 2. At first, the L-ascorbic acid reduced the Cu(OH) 2 precipitates to Cu 2 O, and then, water molecules and C 6 H 8 O 7 were released as the by-products of reaction. When excess acid is added, the surface of Cu 2 O particles were continuously reduced, which formed Cu dots on the surface of NPs, and water molecules and C 6 H 6 O 6 were liberated. Because the reduction reactions occur unevenly on the entire surface of the nanoparticle, the Cu dots are scattered on the surface of the nanoparticle forming Cu/Cu 2 O nanocomposite, but the surface is not completely reduced to form the core-shell structure. The high-resolution TEM technique was used to confirm this explanation ( Figure 3b).  Table 1 with different Tv and Sv.
Third, the growth phase of NPs has two mechanisms: the reaction on the surface and the diffusion of monomers onto the surface [13]. The capping agent not only interacts with the precursor ions but also has an important role in the growth phase. After the formation of NPs, the capping agent protects the surface of those particles by creating a stable cover layer. Moreover, the physicochemical properties of NPs mostly depend on the interaction between the NP and capping agent due to its ability to modify the interaction between the external system and NPs [14]. Additionally, in the growth phase, clusters gather together so that the size of the NPs increases. Fourth, aging is commonly used in the synthesis of NPs, and macro or micro-changes in the solid phase can occur during this process [15]. At the end of the aging period, the size of the NPs does not grow any further. Finally, L-ascorbic acid is used as a reducing agent to reduce Cu(OH)2 so that it precipitates to form Cu2O particles; this process is shown in Figure 2. At first, the L-ascorbic acid reduced the Cu(OH)2 precipitates to Cu2O, and then, water molecules and C6H8O7 were released as the by-products of reaction. When excess acid is added, the surface of Cu2O particles were continuously reduced, which formed Cu dots on the surface of NPs, and water molecules and C6H6O6 were liberated. Because the reduction reactions occur unevenly on the entire surface of the nanoparticle, the Cu dots are scattered on the surface of the nanoparticle forming Cu/Cu2O nanocomposite, but the surface is not completely reduced to form the core-shell structure. The high-resolution TEM technique was used to confirm this explanation ( Figure 3b).  Table 1 with different T v and S v .  Figure 3a shows the SEM image of Cu/Cu2O nanocomposite. We see that the composite has a cubic structure. In order to confirm that the powder is the nanocomposite, not the mixture of Cu and Cu2O particles, the high-resolution TEM technique is used and the image is shown in Figure 3b, which demonstrates the direct reduction in the Cu2O surface layer to form Cu dots to form the nanocomposite. Figure 3b also shows the posite has a cubic structure. In order to confirm that the powder is the nanocomposite, not the mixture of Cu and Cu2O particles, the high-resolution TEM technique is used and the image is shown in Figure 3b, which demonstrates the direct reduction in the Cu2O surface layer to form Cu dots to form the nanocomposite. Figure 3b also shows the (1 1 1) orientation of the crystals of both Cu and Cu2O. The measured d spacing is 0.115 nm and 0.261 nm for Cu and Cu2O, respectively, which is approximately the same as calculated from Bragg's Law (Table 1).    (2 2 2), respectively, and matched well with ICDD #01-078-2076 [16]. In addition to the peaks of Cu 2 O, the peaks of Cu also appeared at 2θ = 43.27 • , 50.41 • , and 74.14 • ; correspond to planes (1 1 1), (2 0 0), and (2 2 0), respectively; matched well with ICDD #785-1326 [17]. This indicates that the Cu/Cu 2 O nanocomposite was successfully synthesized at ambient temperature using a simple aqueous chemical solution approach. To confirm the compounds of powder samples, the RIR is very useful for elucidating the composition of sample. Figure 4b was obtained by fitting the diffraction pattern using the Match! software. It is evident that the ST sample only contains two compounds (Cu 2 O and Cu). Here, Cu 2 O accounts for 60.1%, and the rest is Cu (39.1%). No by-product (such as CuO) or residual product (Cu(OH) 2 ) remains in the sample after the completion of the synthesis. Moreover, the Rietveld refinement method was used to compare with the XRD database of the ST sample. It used MAUD software to scan the standard spectra of Cu 2 O (Pn-3 m, COD 1000063) and Cu (Fm-3 m, COD 5000216) based on the ST's XRD database. Some scanning criteria that were used include crystal structure, strain, quantitation, and texture in this method. The scanning results of Cu 2 O and Cu were shown in Figure 4c, structure, strain, quantitation, and texture in this method. The scanning results of Cu2O and Cu were shown in Figure 4c,d, respectively. It is clear that the peaks at 29.51 o , 36.37°, 42.26°, 61.36°, 73.664°, and 77.38° matched well with COD 1000063 of Cu2O, and the remaining peaks at 43.27°, 50.41°, and 74.14° agreed with COD 5000216 of Cu. Hence, good scanning refinement was achieved for both Cu and Cu2O. From the above discussion, we can conclude that the sample is nanocomposite and not a mixture of Cu and Cu2O particles.

The XRD Analysis of Cu/Cu 2 O Powder Nanocomposites
In order to investigate the lattice parameter and crystalline size of NPs under the changing conditions in the reduction phase, the powder XRD technique was used to analyze all the samples, as shown in Figure 5a,b. From these diffraction patterns, the detectable peaks were similar to the obtained peaks from the ST sample, including six peaks of Cu 2 O (ICDD #01-078-2076) and three peaks of Cu (ICDD #785-1326). The intensity versus 2θ plotted the difference under T v and S v . Normally, the crystalline size was calculated by the Scherrer equation with the XRD database (Equation (2)) [18]. The dislocation of crystal which is one of the lattice parameters considered by lattice strain and calculated by Williamson-Smallman's equation (Equation (3)) [17]: where D, K, λ, β, θ, and δ are the average crystallite size (nm), Scherrer constant, X-ray wavelength, FWHM, the Bragg's angle in degrees, and dislocation of crystal, respectively. Additionally, the lattice strain (ε) was calculated by Equation (4) [19,20]. The results are given in Table 2. The increases in T v and S v led to the increase in movement of all the ions, where D, K, λ, β, θ, and δ are the average crystallite size (nm), Scherrer constant, X-ray wavelength, FWHM, the Bragg's angle in degrees, and dislocation of crystal, respectively. Additionally, the lattice strain (ε) was calculated by Equation (4) [19,20]. The results are given in Table 2. The increases in Tv and Sv led to the increase in movement of all the ions, molecules, and particles in all processes and led to the shifting and widening of Bragg's peaks in the XRD patterns. Thus, the crystallite size values calculated using the Scherrer equation were decreased due to the increase in the dislocation, the lattice strain, and the FWHM [21,22]. Furthermore, the inter-planar spacing values of all samples are of similar lengths, reflecting similar Bragg's peak positions.   As many previous reports have shown, in the wet chemical synthesis method, the nucleation phase is supported by high temperatures, while the growth phase is supported by low temperatures [23,24]. This means that the nanoparticles tend to be larger at lower temperatures, and an increase in temperature during the reaction leads to an increase in size. Furthermore, the decrease in particle size with raising the temperature in the synthesis process is generally considered to increase the reaction rate under higher temperatures [23,25]. Furthermore, increasing the T v leads to thermal expansion of the nanoparticle's unit cell, thereby increasing the lattice constant [26]. Specifically, the lattice parameter increased Stirring speed plays an important role in nanoparticle synthesis, creating the homogeneous distribution of all compounds and giving all compounds more opportunity to react and bind together. However, several reports showed that increasing S v can decrease the particle size [27,28]. Consistent with these results, we observed in the microstructure analysis ( Table 2) that the crystallite size of the samples decreased with increasing S v . Particularly, the increased stirring speed from 50 to 550 rpm reduced the size of the nanocomposites from 3.3881 to 3.0694 nm. However, the lattice parameter and inter-planar spacing did not follow this downtrend. Further increasing the stirring speed, the size of these nanocomposites decreased significantly along with a random size distribution of the inter-planar and lattice parameter [27]. Furthermore, particles formed under the intense stirring speeds would have a higher chance to collide with other particles during the aging stage. As a result, the edges, corners, and surface of the particles will be abraded, and hence, the crystalline size will also be reduced.
To gain insight into the changes in NPs under extreme conditions during synthesis, some methods and fitting models were employed in the analysis. We investigated how increases in the width and the intensity, and how shifting the positions of the peaks were affected by the crystallite size and lattice strain using the Williamson-Hall (W-H) method. This method is an effective method for calculating the induced-lattice strain and induced-crystallite size because it uses breadth integration [18]. This method reduces strain-induced voltage and size expansion by treating the peaks of the XRD pattern as a function of 2θ. The lattice strain and crystal size were obtained from the linear regression equation between 2sinθ and βcosθ. The y-intercept was used to estimate the crystallite size, and the slope was used to fit the lattice strain. Table 2 and Figure 6 show the matching results of the XRD database with the W-H model. From these results, the increases in T v and S v led to a decrease in the induced-crystallite size. The induced-crystallite size values calculated using the W-H method were smaller than the calculated values from the Scherrer equation. Moreover, the negative y-intercept values of S3 and T3 (−0.0189 and −0.0098, respectively) indicated the shriveling of crystals [9].  Due to the influence of Tv and Sv, the crystal growth in all experiments is not perfect, so the lattice parameter is different from the ICDD (4.2696 Å). Therefore, the Nelson-Riley function (Equation (5)) [17] was used to recalculate the precise lattice parameter (PLP) from the XRD database.
f θ = cos θ sinθ + cos θ θ 2 The PLP shown in Figure 7 and Table 2 was obtained from the linear regression relationship between the function f(θ) and the lattice parameter. The PLPs extracted have shorter values compared with the values obtained from ICDD library. To explain why the PLP values have shorter values, the process of measurement and calculation method should be considered. The lattice parameters are usually measured by the direct proce- Due to the influence of T v and S v , the crystal growth in all experiments is not perfect, so the lattice parameter is different from the ICDD (4.2696 Å). Therefore, the Nelson-Riley function (Equation (5)) [17] was used to recalculate the precise lattice parameter (PLP) from the XRD database. The PLP shown in Figure 7 and Table 2 was obtained from the linear regression relationship between the function f(θ) and the lattice parameter. The PLPs extracted have shorter values compared with the values obtained from ICDD library. To explain why the PLP values have shorter values, the process of measurement and calculation method should be considered. The lattice parameters are usually measured by the direct procedure. That means the d-spacing is calculated by Bragg's derivative (λ = 2dsinθ). In this measurement, the accuracy of the crystallographic parameter (a, b, c) and the d-spacing (d) depends entirely on the accuracy of the sin(θ) calculation [29]. Moreover, the Cu/Cu 2 O nanocomposite in this research has a cubic structure (a = b = c = d h 2 + k 2 + l 2 ) Therefore, ∆a = ∆d h 2 + k 2 + l 2 . Additionally, the term ∆a/a (or ∆d/d) arises because of the error for the field in θ, which is the fractional error in a (or d). As θ approaches 90 • , the fractional error approaches 0 [27]. Therefore, ∆a/a = ∆d/a = −cotθ·∆θ.
Due to the influence of Tv and Sv, the crystal growth in all experiments is not perfect, so the lattice parameter is different from the ICDD (4.2696 Å). Therefore, the Nelson-Riley function (Equation (5)) [17] was used to recalculate the precise lattice parameter (PLP) from the XRD database.
f θ = cos θ sinθ + cos θ θ 2 (5) The PLP shown in Figure 7 and Table 2 was obtained from the linear regression relationship between the function f(θ) and the lattice parameter. The PLPs extracted have shorter values compared with the values obtained from ICDD library. To explain why the PLP values have shorter values, the process of measurement and calculation method should be considered. The lattice parameters are usually measured by the direct procedure. That means the d-spacing is calculated by Bragg's derivative (λ = 2dsinθ). In this measurement, the accuracy of the crystallographic parameter (a, b, c) and the d-spacing (d) depends entirely on the accuracy of the sin(θ) calculation [29]. Moreover, the Cu/Cu2O nanocomposite in this research has a cubic structure (a = b = c = d√h + k + l ) Therefore, Δa = Δd√h + k + l . Additionally, the term Δa/a (or Δd/d) arises because of the error for the field in θ, which is the fractional error in a (or d). As θ approaches 90°, the fractional error approaches 0 [27]. Therefore, Δa/a=Δd/a= −cotθ·Δθ.  The key to achieving high accuracy in parametric measurements is the use of a back-reflected beam with a value of 2θ as close to 180 • as possible [20,29]. When approaching 2θ = 180 • (i.e., θ = 90 • ), the value of a will approach the true value. However, because measuring the value at 2θ = 180 • will often be limited, it is necessary to plot the measured values, and then extrapolate to 2θ = 180 • against some function of θ. Moreover, there are many systematic errors, such as the absorbance of the specimen; unintentional curvedness of the sample, affecting the diffraction angles; device error; etc. [30]. The linearity of the histogram at very high angles will provide an exact value of the lattice constant. Furthermore, this function is of great significance in determining the precise lattice constant because the error function becomes minimal when f(θ) = 90 • , which corresponds to the anti-reverse reflection along the path of the incident beam [31]. Accordingly, using the Nelson-Riley function can reduce the systematic errors in the calculation process by using extrapolation from the peak positions to the magnitude 2θ.
Cu/Cu 2 O nanocomposite has many potential applications. Understanding the changes in the crystallographic parameters under the changes of synthesized conditions, scientists can enhance the efficiency of materials. First, the adsorption capacity of nanoporous coreshell Cu@Cu 2 O nanocomposite in the methyl orange removal process was investigated [4], the particle size of samples was changed (300 nm to 4 µm) by changing the precursor concentration and reaction temperature, and the physic-sorption nature of methyl orange adsorption onto the Cu@Cu 2 O nanocomposite was revealed by the activation energy (Ea) of 12.4 kJ/mol in that research. Second, Chen [32] showed the potential of Cu/Cu 2 O nanocomposite films in the solar cell, which can prevent the photo corrosion of the active layer and enhance the lifetime of perovskite solar cells. Furthermore, the combination of Cu/Cu 2 O and Pd nanosheets in glucose sensors shows excellent electrochemical activity in the analysis of glucose in real blood samples [33].

Conclusions
In this study, Cu/Cu 2 O nanocomposites were synthesized by the wet chemical method. The results of the XRD, SEM, and TEM analyses were used to confirm the composition and structure of nanocomposite. The compositions of the sample were also confirmed by the Rietveld refinement and RIR methods. Under a range of reaction times and stirring times, the crystallographic parameters changed. The Williamson-Hall method was helpful in determining that the induced-crystallite size was reduced with increasing T v and S v , and this size is also smaller than the result calculated by the Scherrer equation. Additionally, the Nelson-Riley function is used to recalculate the precise lattice parameters in less time than the results of Bragg's law. Furthermore, the control of reaction temperature and stirring speed can assist in controlling the crystal size of the nanocomposite through the reduction in lattice parameters (down from 4.12 to 4.07 Å) and inter-planar spacing (down from 1.95 to 1.88 Å).