<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">WJNST</journal-id><journal-title-group><journal-title>World Journal of Nuclear Science and Technology</journal-title></journal-title-group><issn pub-type="epub">2161-6795</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/wjnst.2013.31005</article-id><article-id pub-id-type="publisher-id">WJNST-27524</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Engineering</subject><subject> Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Dose Rate Calculation in the Vicinity of the Tunisian Gamma Irradiation
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>lassaad</surname><given-names>Jemii</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Malek</surname><given-names>Mazouz</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Lotfi</surname><given-names>Ghedira</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Département de Physique, Faculté des Sciences de Monastir, Université de Monastir, Monastir, Tunisia</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>elassaad_fsm@yahoo.fr(LJ)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>29</day><month>01</month><year>2013</year></pub-date><volume>03</volume><issue>01</issue><fpage>28</fpage><lpage>32</lpage><history><date date-type="received"><day>November</day>	<month>4,</month>	<year>2012</year></date><date date-type="rev-recd"><day>December</day>	<month>6,</month>	<year>2012</year>	</date><date date-type="accepted"><day>December</day>	<month>27,</month>	<year>2012</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
   Dose rates calculations, in PMMA dosimeters, placed in the vicinity of the Tunisian <sup>60</sup>Co gamma ray irradiator, have been achieved using a pencil like model. The obtained results are in good agreement with recent experimental data. Moreover, in this work we determine also the conversion factor between the dose rate deposited in a PMMA dosimeter and the one deposited in a reel medium. This factor is used to determine the dose deposited in a real irradiated medium such as foodstuff products.  
 
</p></abstract><kwd-group><kwd>Dose Rate; &lt;sup&gt;60&lt;/sup&gt;Co Gamma Source; GEANT4</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Many authors have calculated dose rates for various sources geometries such as a cylindrical source [<xref ref-type="bibr" rid="scirp.27524-ref1">1</xref>] and a linear source [2-4]. These authors calculated the dose rate, at an arbitrary point in a given infinite medium, essentially air. In these calculations, the dose rate is proportional to the photon flux. The proportionality constant is equal to the gamma rate exposure and the energy loss is equal to the photon energy. In the present work, we will not calculate the dose rate in any point on the infinite medium but we calculate the dose rate deposited in a PMMA dosimeter [<xref ref-type="bibr" rid="scirp.27524-ref5">5</xref>] used in the experience and placed in a given direction from the source. However, the obtained results have been compared with recent experimental data. The advantage of this work is the possibility, to obtain a dose rate in each dosimeter position, using a straightforward calculation.</p></sec><sec id="s2"><title>2. Materials and Methods</title><p>The irradiator of the Tunisian National Center for Nuclear Sciences and Technologies (CNSTN) is designed for medical devices sterilisation and foodstuff preservation [<xref ref-type="bibr" rid="scirp.27524-ref6">6</xref>]. It is built of eight <sup>60</sup>Co pencils, each of 45 cm height and 1 cm diameter. The pencils are arranged around the Z-axis in two levels. The overall height of the irradiator is 90 cm, with a mean radial extension of about 6 cm. In a previous work [<xref ref-type="bibr" rid="scirp.27524-ref7">7</xref>], we showed that this irradiator can be simulated by a single pencil of equal height and equal activity placed along the Z-axis.</p><p>Gharbi et al. [<xref ref-type="bibr" rid="scirp.27524-ref8">8</xref>] selected 29 positions, with a 10 cm step, parallel to the Z axis, at a constant distance X = 150 cm from this axis. Dose measurements were carried out using PMMA dosimeters which are oriented perpendicularly to the X-axis (<xref ref-type="fig" rid="fig1">Figure 1</xref>). In order to increase the dosimeter thickness, the authors superposed three films at each position. The films have a parallelepipedic shape: 3 cm long by 1 cm large and 3 mm thick.</p></sec><sec id="s3"><title>3. Dose Rate Calculation</title><p>The gamma rate generated by a pencil-like gamma source of a length 2L and activity A, intercepting a dosimeter placed at a point M such us OM = r, is given by [<xref ref-type="bibr" rid="scirp.27524-ref7">7</xref>]:</p><disp-formula id="scirp.27524-formula104539"><label>(3-1)</label><graphic position="anchor" xlink:href="5-1090101\f8aa241a-4d42-4e36-bf3d-12c9aa785111.jpg"  xlink:type="simple"/></disp-formula><p>The air attenuation factor of the photon flux is equal to<img src="5-1090101\9c2a5c82-6cb7-43fc-977e-8ff62edf8749.jpg" />. Given that, X = 150 cm and μ = 6.8 &#215; 10<sup>−5 cm−1</sup> [<xref ref-type="bibr" rid="scirp.27524-ref9">9</xref>], the air attenuation contribution is clearly negligible (about 1%). Therefore, the number of free electrons per unit time generated in the dosimeter is given by:</p><disp-formula id="scirp.27524-formula104540"><label>(3-2)</label><graphic position="anchor" xlink:href="5-1090101\a16fafac-a690-4013-97d0-ae58c0284fa8.jpg"  xlink:type="simple"/></disp-formula><p>where ξ = 0.061 cm<sup>2</sup>/g [<xref ref-type="bibr" rid="scirp.27524-ref9">9</xref>] is the photon mass attenuation coefficient, ρ = 1.19 g/cm<sup>3</sup> is the PMMA density [<xref ref-type="bibr" rid="scirp.27524-ref9">9</xref>] and e = 3 &#215; 0.3 cm is the dosimeter thickness which corresponding to a mass thickness m = 1.07 g/cm<sup>2</sup>.</p><p>Since<img src="5-1090101\3d6967a7-51aa-4926-a6cb-94f6d7e9f41d.jpg" />, the photon interaction probability ξ, ρe is small (0.065), the development of Equation (3-2) to the first order gives:</p><disp-formula id="scirp.27524-formula104541"><label>. (3-3)</label><graphic position="anchor" xlink:href="5-1090101\aff1e5de-a5ad-429c-b28e-b4f0805ce756.jpg"  xlink:type="simple"/></disp-formula><p>The dose rate is equal to the electrons rate multiplied by the mean energy deposited by each electron in the dosimeter. At <img src="5-1090101\dc1a7487-6734-4d8b-b2bb-b9659048969f.jpg" /> (the mean energy of the two <sup>60</sup>Co photons), Compton scattering is predominant and the forward electron energy is E<sub>0</sub> = 1.03 MeV having a range R = 0.45 g/cm<sup>2</sup> [<xref ref-type="bibr" rid="scirp.27524-ref9">9</xref>] in the PMMA. In order to calculate the energy deposited in the dosimeter, two cases will be treated.</p><sec id="s3_1"><title>3.1. Front End Region</title><p>If a Compton electron is produced in the dosimeter front end region, equal to<img src="5-1090101\8c77d9fd-9af4-4f87-b8bd-3385578d436a.jpg" />, the electron energy loss is equal to<img src="5-1090101\768853b4-3df3-4251-87a8-9a45d55083a2.jpg" />, where <img src="5-1090101\62bb7286-f5a2-436b-9c6b-affcc253f344.jpg" /> is the mean Compton electron energy:</p><disp-formula id="scirp.27524-formula104542"><label>(3-1-1)</label><graphic position="anchor" xlink:href="5-1090101\b674dce7-a77b-4f13-8eea-c2165932f92f.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="5-1090101\a1c8a6da-3b68-420f-a909-5066ee75b8d8.jpg" /> is the electron kinetic energy at a given angle <img src="5-1090101\ed09093b-c4a2-45d4-adf4-f2175fe1b4d6.jpg" /> and <img src="5-1090101\e7e93cc5-2f53-4141-9275-153cd0192a7e.jpg" /> is the Compton differential cross section.</p><p>After integration, we obtain<img src="5-1090101\e3bf133a-a390-4ad3-b8a7-b119e872e204.jpg" />.</p><p>Consequently, the dose rate in the dosimeter is given by:</p><disp-formula id="scirp.27524-formula104543"><label>, (3-1-2)</label><graphic position="anchor" xlink:href="5-1090101\e5f07d83-6bc0-475a-81a1-db984bd3a895.jpg"  xlink:type="simple"/></disp-formula><p>where the factor 2 corresponds to the two <sup>60</sup>Co photons and m<sub>d</sub> = ρeds is the dosimeter mass.</p><p>Thus:</p><disp-formula id="scirp.27524-formula104544"><label>. (3-1-3)</label><graphic position="anchor" xlink:href="5-1090101\e185b418-eb13-436b-b897-869888a2c68d.jpg"  xlink:type="simple"/></disp-formula><p>After integration, we have:</p><disp-formula id="scirp.27524-formula104545"><label>, (3-1-4)</label><graphic position="anchor" xlink:href="5-1090101\0eaae3e2-8098-4209-9b1c-da5e8cd1694a.jpg"  xlink:type="simple"/></disp-formula><p>where: <img src="5-1090101\2a66abbf-7ab9-42be-a5d9-1beaa42fca9a.jpg" />is a constant. In our case A = 1.81 &#177;</p><p>0.18 PBq [<xref ref-type="bibr" rid="scirp.27524-ref8">8</xref>] and L = 45 cm, we obtain K = 660 Gy m/h.</p><p>We obtain the same expression found by Huttlin [<xref ref-type="bibr" rid="scirp.27524-ref4">4</xref>], for a linear source of length 2L. However the two constants K and f are not equal, precisely</p><p><img src="5-1090101\233acf73-8f6c-4d0b-b727-4c5502d650cc.jpg" />.</p></sec><sec id="s3_2"><title>3.2. Back End Region</title><p>If a Compton electron is produced in the dosimeter back end of a mass thickness equal to R, this electron will lose only a fraction of its kinetic energy. Therefore, the electron mean energy loss in the dosimeter, is no more equal to<img src="5-1090101\e2b3d012-5e8c-49b1-91da-d22db95474c5.jpg" />, but to <img src="5-1090101\cea7ddc6-3673-4ae3-bc5e-4a97fdce165c.jpg" /> given by:</p><p><img src="5-1090101\8d74f130-1aea-4090-b679-f868826ec4e2.jpg" />which is equivalent to:</p><disp-formula id="scirp.27524-formula104546"><label>, (3-2-1)</label><graphic position="anchor" xlink:href="5-1090101\69759a61-e150-4b50-a717-18c0a4da4b02.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="5-1090101\60d69332-7d01-4ebe-aa26-cd2cdd404986.jpg" /> is the mean energy loss of escaping electrons.</p><p>The dose rate is then:</p><disp-formula id="scirp.27524-formula104547"><label>, (3-2-2)</label><graphic position="anchor" xlink:href="5-1090101\cbf5b6f0-9c81-4c5a-ab33-8bd1bc4155bd.jpg"  xlink:type="simple"/></disp-formula><p>After integration, we obtain:</p><disp-formula id="scirp.27524-formula104548"><label>(3-2-3)</label><graphic position="anchor" xlink:href="5-1090101\0e5647f8-dbde-40c8-bc79-65b9352b7907.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="5-1090101\8efec1df-85ab-426e-ac26-8a811f11dda7.jpg" /> is given by:</p><disp-formula id="scirp.27524-formula104549"><label>, (3-2-4)</label><graphic position="anchor" xlink:href="5-1090101\ee49220b-aafe-44d1-8da7-01465a914bfa.jpg"  xlink:type="simple"/></disp-formula><p>is the correction factor due to escaping electrons.</p><p>The constant <img src="5-1090101\fd15f39c-2ce1-4655-923a-a463f5bf48e1.jpg" /> is determined by estimating<img src="5-1090101\450e8480-ede5-4649-aa55-e3ceca9b9c23.jpg" />.</p><sec id="s3_2_1"><title>3.2.1. Energy Loss of Escaping Electron</title><p>The mean energy loss of escaping electrons is given by:</p><p><img src="5-1090101\a37b54b6-3899-4483-a87d-76dbb7b06d7a.jpg" /></p><p>where <img src="5-1090101\32d98854-b5b6-497f-9534-ec68d3746799.jpg" /> is the probability distribution. <img src="5-1090101\135d55ab-60b7-4fb6-af65-3422ffdbea0a.jpg" />is the electrons rate per unit mass such as: <img src="5-1090101\df860819-f493-4ccc-921d-8fa53dff8fdb.jpg" />and <img src="5-1090101\50a6c50c-1a02-49df-b5ad-1b8de3212d55.jpg" /> is the energy loss of escaping electrons per unit mass.</p><p>Since in our case<img src="5-1090101\dc97f445-66f1-4170-9ce0-13a632f33bb8.jpg" />, we can write<img src="5-1090101\ae3a46a2-0703-46db-a8b9-e50ed49111b0.jpg" />, thus<img src="5-1090101\44a62b17-1450-441c-8cdd-a471488c3ddc.jpg" />.</p><p>Therefore, we can write:</p><disp-formula id="scirp.27524-formula104550"><label>. (3-2-5)</label><graphic position="anchor" xlink:href="5-1090101\657fed69-f8ce-46a2-8871-5d95133e9d10.jpg"  xlink:type="simple"/></disp-formula><p>Now, we use the approximation:<img src="5-1090101\41d33d2e-b256-4d2e-bded-f326e6a393b0.jpg" />, where a is a constant and E is the electron energy.</p><p>Given that<img src="5-1090101\bec4e4cd-b80d-4240-b541-6deab0d52998.jpg" />, where <img src="5-1090101\6a861e39-565a-4113-b220-652c31356c5e.jpg" /> is the forward Compton electron energy and <img src="5-1090101\f249995b-b160-46f6-b287-9cd375595bbb.jpg" /> is the energy loss of escaping electron, thus we have:</p><disp-formula id="scirp.27524-formula104551"><label>. (3-2-6)</label><graphic position="anchor" xlink:href="5-1090101\2c31c92a-2a76-41c9-af4e-ef505ad37758.jpg"  xlink:type="simple"/></disp-formula><p>Finally, we obtain:</p><disp-formula id="scirp.27524-formula104552"><label>. (3-2-7)</label><graphic position="anchor" xlink:href="5-1090101\b4569a20-9da1-47f0-9b5e-674cd90489c1.jpg"  xlink:type="simple"/></disp-formula><p>By injecting Equation (3-2-14) in Equation (3-2-8), we find:</p><disp-formula id="scirp.27524-formula104553"><label>. (3-2-8)</label><graphic position="anchor" xlink:href="5-1090101\1e8ee850-edf3-4ae3-8a18-61db6800fd8c.jpg"  xlink:type="simple"/></disp-formula></sec><sec id="s3_2_2"><title>3.2.2. Simulation</title><p>To check the validity of Equation (3-2-13), it is equivalent to show Equation (3-2-14). To determinate the constant a, we integrate Equation (3-2-14) and we find</p><p><img src="5-1090101\085eed1d-8d4d-4511-9601-3e0212f2d9ec.jpg" />.</p><p>Finally Equation (3-2-14) can be written as:</p><disp-formula id="scirp.27524-formula104554"><label>. (3-2-9)</label><graphic position="anchor" xlink:href="5-1090101\7874cec6-cd10-4bad-bf85-b08ded26562c.jpg"  xlink:type="simple"/></disp-formula><p>To validate this equation, we used the GEANT4 Monte Carlo code. The irradiator is simulated by a single pencil of a 90 cm length (<xref ref-type="fig" rid="fig1">Figure 1</xref>). A PMMA film with a mass thickness equal to the range R of the forward electron is placed at 150 cm from the center of the source. <xref ref-type="fig" rid="fig2">Figure 2</xref></p><p>shows the obtained probability distribution <img src="5-1090101\0d2db0d6-c077-4dc4-bf0b-30bd8f86c338.jpg" /> as a function of the escaping electron energy loss. The fit of the obtained distribution by Equation (3-2-9) gives a good mean squared error<img src="5-1090101\80066c29-d457-428e-b1e9-1f800ffa9432.jpg" />.</p></sec></sec><sec id="s3_3"><title>3.3. Conversion Factor</title><p>To determine the dose rate d<sub>med</sub> in a very thick medium from the dose rate d<sub>dos</sub> measured using a PMMA dosimeter we write:</p><disp-formula id="scirp.27524-formula104555"><label>, (3-3-1)</label><graphic position="anchor" xlink:href="5-1090101\dc75450e-2182-4518-a3f9-1d5ca1e18290.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="5-1090101\604148fe-24d0-44a4-81b8-a8d304f3c64e.jpg" /> is the conversion factor.</p><p>Usually, the length L is small relatively to the distance r, in our case L = 45 cm and X = 150 cm.</p><p>Therefore, the correction factor becomes:</p><disp-formula id="scirp.27524-formula104556"><label>. (3-3-2)</label><graphic position="anchor" xlink:href="5-1090101\1145c1b1-0ea2-4e4d-b760-6556230524ba.jpg"  xlink:type="simple"/></disp-formula><p>Given that:</p><p><img src="5-1090101\26cf7d11-2879-4e0d-b433-b5d21898bbda.jpg" />, <img src="5-1090101\881b00ba-411e-4436-8887-a3bcb3952c1c.jpg" />, <img src="5-1090101\08c97d83-ca12-4366-9482-3e6853293189.jpg" />, R = 0.45 g/cm<sup>2</sup> and m = 1.07 g/cm<sup>2</sup>.</p><p>so k = 0.178, then for Z = 0, we obtain: C<sub>f</sub> = 2.62.</p><p>This means that the measured dose rate is only about 38% of the dose rate that will be deposited in an infinite medium. <xref ref-type="fig" rid="fig3">Figure 3</xref> shows the curve of the conversion factor as a function of dosimeter position.</p></sec></sec><sec id="s4"><title>4. Results</title><p>The comparison between the calculated and the experimental dose rates is shown in <xref ref-type="fig" rid="fig4">Figure 4</xref>. Given that the uncertainty on the source activity is about 10% [<xref ref-type="bibr" rid="scirp.27524-ref8">8</xref>], the distribution of the experimental data and the theoretical curve show a very good agreement which is a proof of the validity of the calculation. The conversion factor between</p><p>the dose rates in a given thick medium using the one measured by a PMMA dosimeter can be very useful, for instance, when commissioning a similar irradiation facility.</p></sec><sec id="s5"><title>5. Discussion</title><p>The authors of reference [<xref ref-type="bibr" rid="scirp.27524-ref8">8</xref>] have determined the dose rate by simulating the dosimeters with spheres of 8 cm diameter filled with water. The mass thickness of the sphere is more than 20 times the range of the forward Compton electrons. Then, the multiple scattering is not negligible and the electron energy loss in such spheres is higher than in the dosimeter case. Consequently, the escaping effect has then been neglected, which is not the case, in PMMA dosimeters. The same authors find a too good agreement between their simulation and the data, which is questionable. It seems that, their simulation is not absolute but relative to the dosimeter placed at Z = 0. Finally, if the simulation is necessary to check or validate some calculations, we see no reason for not using the real dosimeter geometry in order to obtain an absolute result.</p></sec><sec id="s6"><title>6. Conclusion</title><p>In conclusion, using a relatively simple calculation, we are able, to determine the dose rates in the vicinity of the CNSTN <sup>60</sup>Co irradiator. The agreement between our calculation and the experimental data is satisfactory. We have also shown that the escaping effect, in a PMMA dosimeter, is not negligible relatively to the spheres used by some references authors. Finally, the conversion factor and the dose rates calculated in such dosimeters can be used to obtain the dose rates deposited in a real medium.</p></sec><sec id="s7"><title>REFERENCES</title></sec><sec id="s8"><title>NOTES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.27524-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">M. Contineanu, S. Perisanu and A. Neacsu, “Calculation of the Dose Rate in an External Point of a Cylindrical Gamma Radioactive Source,” Anatele Universitatti di Bucuresti—Chimi (Serie Noua), Vol. 19, No. 1, 2010, pp. 69-77.</mixed-citation></ref><ref id="scirp.27524-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">J. J. Fletcher, G. Emi-Renolds and E. T. Glover, “Shielding Calculations for Changing from Circular to a Rectangular Source Configuration for a Pilot Scale Irradiator,” Journal of Applied Science and Technology, Vol. 5, No. 1-2, 2000, pp. 39-40.</mixed-citation></ref><ref id="scirp.27524-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">G. Z. Gochaliev, “Calculation of the Average Dose Rate for Linear Gamma-Radiation Sources,” High Energy Chemistry, Vol. 44, No. 6, 2010, p. 539.  
doi:10.1134/S0018143910060147</mixed-citation></ref><ref id="scirp.27524-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">G. A. Huttlin, “Analytical 60Co Dose-Rate Contours Maps,” Army Research Laboratory. Adelphi, MD. ARLMR-324, 1996.</mixed-citation></ref><ref id="scirp.27524-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">www.harwell-dosimeters.co.uk/specifications</mixed-citation></ref><ref id="scirp.27524-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">A. Loussaief, A. Trabelsi and B. Baccari, “Extended Gamma Sources Modelling Using Multipole Expansion: Application to the Tunisian Gamma Source Load Planning,” Radiation Physics and Chemistry, Vol. 75, No. 4, 2006, pp. 463-472.  
doi:10.1016/j.radphyschem.2005.12.024</mixed-citation></ref><ref id="scirp.27524-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">E. Jemii, M. Mazouz, A. B. Fredj, et al., “Modeling of the Tunisian 60Co Gamma Irradiator by a Coaxial Equal Height and Equal Activity Single Pencil,” Radiation Physics and Chemistry, Vol. 80, No. 11, 2011, pp. 11581161. doi:10.1016/j.radphyschem.2011.05.003</mixed-citation></ref><ref id="scirp.27524-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">F. Gharbi, O. Kadri, K. Farah, et al., “Validation of Geant Code as Predictive Tool of Dose Rate Measurement in the Tunisian Gamma Irradiation Facility,” Radiation Physics and Chemistry, Vol. 74, No. 2, 2005, pp. 102-110.  
doi:10.1016/j.radphyschem.2005.02.006</mixed-citation></ref><ref id="scirp.27524-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">www.physics.nist.gov/pml/data</mixed-citation></ref></ref-list></back></article>