<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving DTD v1.0 20120330//EN" "JATS-journalarchiving.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:ali="http://www.niso.org/schemas/ali/1.0">
  <front>
    <article-meta>
      <title-group>
        <article-title>Calculating thin film strain from XRD data</article-title>
      </title-group>
      <contrib-group content-type="author">
        <contrib id="person-39e7d4208113ec2cd4f16647472c1dc0" contrib-type="person" equal-contrib="no" corresp="no" deceased="no">
          <name>
            <surname>Qadr</surname>
            <given-names>Shelan Salm Sadiq</given-names>
          </name>
          <email>shelansalam699@gmail.com</email>
          <xref ref-type="aff" rid="aff-1" />
        </contrib>
      </contrib-group>
      <aff id="aff-1">
        <country>Iraq</country>
      </aff>
      <history>
        <date date-type="received" iso-8601-date="2023-07-12">
          <day>12</day>
          <month>07</month>
          <year>2023</year>
        </date>
      </history>
      <abstract />
    </article-meta>
  </front>
  <body id="body">
    <sec id="heading-2e33e82eec6703a53c7e4cddd14f69cb">
      <title>
        <bold id="bold-9d02616c09530b49147ee58f3a6fea88">Introduction</bold>
      </title>
      <p id="_paragraph-8">Different synthesis methods used to create nanostructured materials can result in a variety of structural flaws, including lattice strain. The diffraction peaks may get broader in crystals with nanometer-sized particles in the range of 10 nm to 100 nm. The crystallite size and lattice strain are determined by the widening of the diffraction peak. Scherrer's formula can be used to determine the crystallite size. In addition to these methods, the Williamson-Hall method was also used to determine the size of crystallites. The lattice strain is also calculated using this method. As a result, ZnO nanopowders were created in this work by precipitation and calcined at various temperatures. The Williamson-Hall (W-H) method was used to determine the crystallite size and lattice strain..</p>
      <p id="_paragraph-9">
        <bold id="_bold-7">Thin film </bold>
        <bold id="_bold-8">microstrain</bold>
      </p>
      <p id="_paragraph-10">Measurement of strain or deformation in the individual grains of thin films is referred to as thin film microstrain. To carry out microstrain measurements, different techniques like EBSD, -HRXRD, Raman microspectrometry, and X-Ray diffraction (XRD) can be employed. The values of microstrain measured in thin films containing multiple crystals can offer insights into various structural parameters, including crystallite size, dislocation density, and the sensitivities of the techniques employed for mapping microstrain in thin films. [1].</p>
      <p id="_paragraph-11">
        <bold id="_bold-9">Williamson-Hall (W-H) Analysis</bold>
      </p>
      <p id="_paragraph-12">X-ray profile analysis is a simple and powerful tool to estimate the crystallite size and lattice strain. This method is attributed to G. K.Williamson and his student, W. H. Hall [2]. It relies on the principle that the approximate formulae for size broadening, BL, and strain broadening, ße, vary quite differently with respect to Bragg angle, 0.</p>
      <p id="_paragraph-13">
        <bold id="_bold-10">Literature Review</bold>
      </p>
      <p id="_paragraph-14">ZnO nanopowders were created by Mote et al. using coprecipitation methods. At 450 °C, grains with an average size of roughly 50 nm were created. Using W-H analysis and several models, such as the uniform deformation model, the uniform deformation stress model, and the uniform deformation energy density mode, the physical parameters such as strain, stress, and energy density values were also computed. [3].</p>
      <p id="_paragraph-15">Fe-doped ZnO nanoparticles with a hexagonal wurtzite shape have been created by Prabhu et al. Williamson-Hall (W-H) analysis and a size-strain plot were used to examine the effects of crystallite sizes and lattice strain on the peak broadening of Fe doped ZnO nanoparticles. All of the samples' XRD peaks had their strain, stress, and energy density characteristics computed using the size-strain plot method (SSP), the uniform stress deformation model (USDM), the uniform deformation energy density model (UDEDM), and the uniform deformation model (UDM). studies from W-H analysis, SSP, and TEM studies were corroborated by the mean particle size of Fe doped ZnO-NPs[4].</p>
      <p id="_paragraph-16">Using nanocrystalline Ni and Cu powders, Brandstetter et al. utilized Williamson-Hall analysis to determine grain size and lattice strain[5].</p>
      <p id="_paragraph-17">Similar to this, Zak et al. used the sol-gel process to create ZnO nanopowders. The sample has a hexagonal wurtzite phase and is crystalline in form. XRD was used to investigate the ZnO nanoparticles' phase evolution. To investigate the individual effects of crystallite sizes and lattice strain on the peak broadening of the ZnO nanoparticles, the Williamson-Hall analysis and size-strain plot approach were utilized. ZnO-NPs that were calcined at 750 °C showed a nonuniform strain and an average particle size of roughly 20 nm in the TEM image. The TEM results and the W-H technique results were in good agreement[6].</p>
      <p id="_paragraph-18">In order to measure the lattice strain on ZnO particles produced by a straightforward, surfactant-assisted combustion synthesis, Prabhu et al. also used the W-H method[7].</p>
      <p id="_paragraph-19">ZnO flat thin film was created by Thool et al. using chemical bath deposition methods. The Scherrer method and the Williamson-Hall method were used to evaluate the crystallite size and lattice strain from X-ray line broadening[8].</p>
      <p id="_paragraph-20">After heating the precursor at 350 degrees Celsius, Sarma et al. generated ZnO nanopowders using an easy, cost-effective precipitation technique. The W-H method was used to analyze the line broadening of ZnO nanoparticles caused by the small crystallite size and strain[9].</p>
      <p id="_paragraph-21">
        <bold id="_bold-11">Aim of this study</bold>
      </p>
      <p id="_paragraph-22">The objective is to calculate strain and crystallite size from XRD data using Willium- Hall analysis.</p>
    </sec>
    <sec id="heading-a3f541d25bfb47b6e3f96947a9fc835e">
      <title>
        <bold id="_bold-13">Method</bold>
      </title>
      <p id="_paragraph-24">We will calculate crystallites size and strain from XRD data using Williamson-Hall (W-H) plot method. The broadening (<inline-formula id="inline-formula-e3ee57befb21fa4da5cee05a51bbb7d4" content-type="math/tex"><tex-math><![CDATA[B_T]]></tex-math></inline-formula>)of <inline-formula id="inline-formula-c4d039ad64e48c39a917526ef28b3ffe" content-type="math/tex"><tex-math><![CDATA[ ]]></tex-math></inline-formula>the peaks in the XRD data is caused by the combined effects of crystallite size (<inline-formula id="inline-formula-e042ce4e466149c347cb141c0f026ae7" content-type="math/tex"><tex-math><![CDATA[B_D]]></tex-math></inline-formula>) and microstrain (<inline-formula id="inline-formula-5b115069316d3e0e395b964ebe26be0d" content-type="math/tex"><tex-math><![CDATA[B_e]]></tex-math></inline-formula>).</p>
      <p id="_paragraph-25">Total broadening = Broadening due to crystallites size + Broadening due to strain</p>
      <p id="_paragraph-27">Where:</p>
      <p id="paragraph-4e612518c18418bc23a888018b515584">
        <inline-formula id="inline-formula-894c1118bfdd00d304c90effb9d70e2d" content-type="math/tex">
          <tex-math><![CDATA[B_T= B_D+ B_e  ]]></tex-math>
        </inline-formula>
      </p>
      <p id="_paragraph-28"><inline-formula id="inline-formula-3ef78e9cba4a413923ac6ae7bff1f947" content-type="math/tex"><tex-math><![CDATA[B_T]]></tex-math></inline-formula> is the total broadening</p>
      <p id="_paragraph-29"><inline-formula id="inline-formula-e7b786f2771992f93e826fe698d8efe6" content-type="math/tex"><tex-math><![CDATA[B_D]]></tex-math></inline-formula> is broadening due to crystallite size</p>
      <p id="_paragraph-30"><inline-formula id="inline-formula-61a284bcd662ccc219508b75a784345d" content-type="math/tex"><tex-math><![CDATA[B_e]]></tex-math></inline-formula>  is the broadening due to strain</p>
      <p id="_paragraph-31">From the Scherer equation, we know that,</p>
      <p id="_paragraph-32">
        <inline-formula id="inline-formula-7f6f715e31f0a29da6c36e432284d161" content-type="math/tex">
          <tex-math><![CDATA[B_D=kλ/(D cos⁡θ )]]></tex-math>
        </inline-formula>
      </p>
      <p id="_paragraph-33">Where</p>
      <p id="_paragraph-34"><inline-formula id="inline-formula-5c5c24d6014c41959da852ac222942fb" content-type="math/tex"><tex-math><![CDATA[B_D]]></tex-math></inline-formula> is FWHM (ie broadening of the peak) in radians <inline-formula id="inline-formula-e5257d6ebab59f23784d4694bb5d67e5" content-type="math/tex"><tex-math><![CDATA[k= 0.9]]></tex-math></inline-formula> is the shape factor, <inline-formula id="inline-formula-6b44122adf89a6bb34cb782d31cb5752" content-type="math/tex"><tex-math><![CDATA[λ= 0.15406  nm]]></tex-math></inline-formula> is the wavelength of X-ray source, D is the crystal size, is the peak position in radians similarly.</p>
      <p id="_paragraph-35">The XRD peak broadening due to micro strain is given by</p>
      <p id="paragraph-1aca6a628ca46634f89f0996af9a32b8">
        <inline-formula id="inline-formula-63943b0029b6beccb26a6100769deec3" content-type="math/tex">
          <tex-math><![CDATA[B_e=4 ε tan⁡〖θ                                                   (2-3)〗]]></tex-math>
        </inline-formula>
      </p>
      <p id="_paragraph-36">Where <inline-formula id="inline-formula-b3c5a8d954f30d0c96d92c4fea4b6e1c" content-type="math/tex"><tex-math><![CDATA[B_e]]></tex-math></inline-formula>. is broadening due to strain and <inline-formula id="inline-formula-0190d136ba2f182d3d8d0cb773f9963d" content-type="math/tex"><tex-math><![CDATA[θ]]></tex-math></inline-formula> is the peak position in radians.</p>
      <p id="_paragraph-37">Substituting equation (2-2) and (2-3) in equation (2-1),</p>
      <p id="_paragraph-38">we obtain:</p>
      <p id="paragraph-1c1b670d2d9e6a62d26c1bb9d5c5b403">
        <inline-formula id="inline-formula-72de9fe5e0fe9ddb9cc7120c18531dd5" content-type="math/tex">
          <tex-math><![CDATA[B_T=   kλ/(D cos⁡θ )+ 4 ε tan⁡〖θ                                   ]]></tex-math>
        </inline-formula>
      </p>
      <p id="_paragraph-39">As we know that <inline-formula id="inline-formula-a1f4798a744713351dfacae7e2a88f08" content-type="math/tex"><tex-math><![CDATA[tan⁡θ=    sin⁡θ/cos⁡θ ]]></tex-math></inline-formula></p>
      <p id="_paragraph-40">Therefore, equation (2-4) can be written as</p>
      <p id="paragraph-0c8eb6a096591c3eed285921d0ccc008">
        <inline-formula id="inline-formula-b279fd39e5da4fd471c76b8ad416ae91" content-type="math/tex">
          <tex-math><![CDATA[B_T=   kλ/(D cos⁡θ )+ 4 ε  sin⁡θ/cos⁡θ         ]]></tex-math>
        </inline-formula>
      </p>
      <p id="_paragraph-41">Multiplying both sides by <inline-formula id="inline-formula-bf8558fe4b34f24034dbd4c639e3f901" content-type="math/tex"><tex-math><![CDATA[cos⁡θ]]></tex-math></inline-formula>,</p>
      <p id="paragraph-30c5d8c13c159a58d86bf190efd6f06b">
        <inline-formula id="inline-formula-fe343c627dd0483323ea83a17b22e6c9" content-type="math/tex">
          <tex-math><![CDATA[cos⁡θ× B_T=   cos⁡θ  ×kλ/(D cos⁡θ )+cos⁡θ× 4 ε  sin⁡θ/cos⁡θ   ]]></tex-math>
        </inline-formula>
      </p>
      <p id="paragraph-0457602106a47afa20398c4f593c42ae">
        <inline-formula id="inline-formula-4223b71a53eda1552d9c5f190af37020" content-type="math/tex">
          <tex-math><![CDATA[B_T  cos⁡θ=   kλ/D+ 4 ε sin⁡θ                           ]]></tex-math>
        </inline-formula>
      </p>
      <p id="_paragraph-42">Equation (2-5) represent a straight line, in which is the gradient (slope) of the line and <inline-formula id="inline-formula-5ec204a371ef4f4b68b5c50e72026030" content-type="math/tex"><tex-math><![CDATA[kλ/D]]></tex-math></inline-formula> is the y- intercept.</p>
      <p id="_paragraph-43">Consider the standard equation of a straight line,</p>
      <p id="paragraph-532bbc4e5894b932badbbd44a2d4389c">
        <inline-formula id="inline-formula-ae2976ac9adda31f367fc12ab310770c" content-type="math/tex">
          <tex-math><![CDATA[y = mx + c ]]></tex-math>
        </inline-formula>
      </p>
      <p id="_paragraph-44">Where is the slope of line and is the y- intercept</p>
      <p id="_paragraph-45">Comparing equation (2-5) with equation (2-6), we have:</p>
      <p id="_paragraph-46"> <inline-formula id="inline-formula-f56609aceebfc4fa175c4b7d1980cd39" content-type="math/tex"><tex-math><![CDATA[y = B_T    cos⁡θ]]></tex-math></inline-formula>                                      i</p>
      <p id="_paragraph-47"><inline-formula id="inline-formula-b5125f29b08c4dbea5fd512f395e8ffb" content-type="math/tex"><tex-math><![CDATA[ m=ε]]></tex-math></inline-formula>                                                  ii</p>
      <p id="_paragraph-48"><inline-formula id="inline-formula-90d9b69666d482139d498e542fbbf4e1" content-type="math/tex"><tex-math><![CDATA[ x=4 sin⁡θ]]></tex-math></inline-formula>                                           iii</p>
      <p id="_paragraph-49"><inline-formula id="inline-formula-8e9ca7f6406065ad53f4d470984e0a39" content-type="math/tex"><tex-math><![CDATA[c=  kλ/D]]></tex-math></inline-formula>                                           iv</p>
      <p id="_paragraph-50">We will plot (<inline-formula id="inline-formula-9f67052b4bb6447d34f6540359dc7ffe" content-type="math/tex"><tex-math><![CDATA[4 sin θ]]></tex-math></inline-formula>) on x-axis and (<inline-formula id="inline-formula-889ca4ca4caa514dfb4da49a86204a16" content-type="math/tex"><tex-math><![CDATA[B_T    cos⁡θ]]></tex-math></inline-formula>) on y-axis.</p>
      <p id="_paragraph-51">The value of m which represent gradient (slope) of the line it will be the value of the strain ɛ. Finally, we will calculate crystallites size from the y-intercept</p>
      <fig id="figure-panel-104596ded9062b28c3ea92ec78468009">
        <label>Figure 1</label>
        <caption>
          <p id="paragraph-a09b063d844ff43e0229cd612a82ae39" />
        </caption>
        <graphic id="graphic-2e03ea1601a05a127cc70baaf991f980" mimetype="image" mime-subtype="png" xlink:href="1.png" />
      </fig>
      <p id="_paragraph-52">We need to follow the following steps to calculate the strain from XRD data:</p>
      <p id="_paragraph-53">1- By using originLab program we will determain 20 (degree) for each peak and FWHM, BT (in degree) from XRD data.</p>
      <p id="_paragraph-54">2- Convert and FWHM BT from degree to radians.</p>
      <p id="_paragraph-55">3- Determain the value of <inline-formula id="inline-formula-16513fb36a13c831548655c4da5dd190" content-type="math/tex"><tex-math><![CDATA[B_T    cos⁡θ]]></tex-math></inline-formula> and <inline-formula id="inline-formula-0e7861bf3e712e976b6ef99f003e3e83" content-type="math/tex"><tex-math><![CDATA[4 sin⁡θ]]></tex-math></inline-formula></p>
      <p id="_paragraph-56">4- Plot B <inline-formula id="inline-formula-2c0ffccc7ad4f4c83f1d0659f1bd22e4" content-type="math/tex"><tex-math><![CDATA[B_T    cos⁡θ]]></tex-math></inline-formula> on y-axis and <inline-formula id="inline-formula-20de6e46ce42952ce4788cc463df7158" content-type="math/tex"><tex-math><![CDATA[4 sin⁡θ]]></tex-math></inline-formula> on x-axis to obtain the final W-H plot in origin Lab program as scatter graph.</p>
      <p id="_paragraph-57">5- Preform liner fitting for the scatter graph.</p>
      <p id="_paragraph-58">6- The slope of the plot is the strain<inline-formula id="inline-formula-0a402a570c98b038cb98398306434a62" content-type="math/tex"><tex-math><![CDATA[ɛ]]></tex-math></inline-formula> .</p>
      <p id="_paragraph-59">The XRD data adoptede from the refreance [10].</p>
    </sec>
    <sec id="heading-a8931e3233792f1e531d8edac15f4afc">
      <title>
        <bold id="_bold-14">Results and Discussions</bold>
      </title>
      <p id="_paragraph-62">Figures 3-1 show XRD images of the Cu samples for various hot water treatment time.</p>
      <fig id="figure-panel-f3f7be987dbe639f1597937a4c48f9b9">
        <label>Figure 2</label>
        <caption>
          <p id="paragraph-ad9159b78a90dfc60b91af5d3a91854e" />
        </caption>
        <graphic id="graphic-51dc75dc5c9d36567d9fc2a5b89b82ce" mimetype="image" mime-subtype="png" xlink:href="2.png" />
      </fig>
      <p id="_paragraph-63">Figure 3-1: XRD profiles for the dominant oxide peaks of Cu samples for various hot water treatment time [10].</p>
      <p id="_paragraph-64">Figure 3-2 shows thr realtion between on y-axis and on x-axis and the calculatet slope represent the stanin () according to equation 2-3.</p>
      <fig id="figure-panel-925d318a5ddb5de6939356d87f4cb810">
        <label>Figure 3</label>
        <caption>
          <p id="paragraph-5c40575d62613b1e3c5f63ddb1998428" />
        </caption>
        <graphic id="graphic-dc0fc2faea7aba27990846448574e703" mimetype="image" mime-subtype="png" xlink:href="3.png" />
      </fig>
      <p id="_paragraph-65">Figure 3-2: The W-H Analysis Cu samples for various hot water treatment time.</p>
      <p id="_paragraph-66">Table 1 reveal calculated crystal size of samples for various hot water treatment time where the intercepts of the plots in figure 3-2 represent the value <inline-formula id="inline-formula-a55505b30c9446a8b7f32ca783ed5efc" content-type="math/tex"><tex-math><![CDATA[kλ/D]]></tex-math></inline-formula> in equation 2-5, where <italic id="italic-8d1cc9d4f046ba1eda999a348b3dc464">D</italic> can be obtained which is represent the crystal size by 1-2 Williamson- Hall (W-H) method. The calculate crystal size by Scherer method also included in Table 3-1 for comparison.</p>
      <table-wrap id="_table-figure-1">
        <label>Table 1</label>
        <caption>
          <title>calculate crystal size</title>
          <p id="_paragraph-68" />
        </caption>
        <table id="_table-1">
          <tbody>
            <tr id="table-row-989416ebcfc114e3d508c5f9aa0cbbd1">
              <td id="table-cell-95e8594ede4bcdc57df36a28748a90b5">Treatment time</td>
              <td id="table-cell-04536716f5b8556e7b5dda5c8996c962"><inline-formula id="inline-formula-3b66ff95d5bfb8cd032fbc93af9517ec" content-type="math/tex"><tex-math><![CDATA[D_(W-H) (nm)]]></tex-math></inline-formula>
Williamson-Half method</td>
              <td id="table-cell-c3f54924abe65032838406a9a58e5833"><inline-formula id="inline-formula-903ea98183154aadecd64fd423110370" content-type="math/tex"><tex-math><![CDATA[D_sh (nm)]]></tex-math></inline-formula>[9]
Scherer method</td>
            </tr>
            <tr id="table-row-9daffa9456c46834d65438d9f6801b7a">
              <td id="table-cell-cc48bbebdbbc831814172aec5b35cc14">0</td>
              <td id="table-cell-33c46931ad4c73f3d77189fd5dc7d374">22</td>
              <td id="table-cell-6642b35564574c56ca37c01ab81ddda9">21</td>
            </tr>
            <tr id="table-row-30576e2574304eeb7c533364a7b0bcbe">
              <td id="table-cell-f5351cdd956d760544dda4ec82f0cc9a">1</td>
              <td id="table-cell-9f6101c10230b0391dbe218782bd254b">23</td>
              <td id="table-cell-c6ca680906609efc951951efa82223aa">23</td>
            </tr>
            <tr id="table-row-507cadd7e6a433d4665b1637c8de28f1">
              <td id="table-cell-84107aee5e32e53964ed42a67d077d2b">5</td>
              <td id="table-cell-fad19fb0e99ce13d7abf08fadb282ade">26</td>
              <td id="table-cell-9680f4042b0b223e9d93b5dbc433eb35">27</td>
            </tr>
            <tr id="table-row-b2d172f544bf79fea320034a60df881e">
              <td id="table-cell-601e5cefc6faeea62e9ad536abc63271">7</td>
              <td id="table-cell-9d2e78f47564c18c9343c6ad48a69fae">21</td>
              <td id="table-cell-5ad8607c5eb0f942bda4144ef6152ccf">21</td>
            </tr>
          </tbody>
        </table>
      </table-wrap>
      <p id="_paragraph-70">Temperature can affect the mechanical properties of copper, where the dynamic softening mechanisms play a dominant role in the deformation of copper [11].</p>
    </sec>
    <sec id="heading-bd1eca16fe4c273494d43bc3cc77552f">
      <title>
        <bold id="_bold-19">Conclusion</bold>
      </title>
      <p id="_paragraph-73">1- From the calculation of strain, we can observe that increasing treatment time increases the strain of the thin films,</p>
      <p id="_paragraph-74">2- The calculated crystal size was in good agreement with the one measured by reference [10].</p>
      <p id="_paragraph-75">
        <bold id="_bold-23">Recommendation</bold>
      </p>
      <p id="_paragraph-76">Determine other relevant physical parameters such as stress and energy density.</p>
      <p id="paragraph-3e7112b5ad3d2c65f48f92287367c357" />
    </sec>
  </body>
  <back />
</article>