<?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">JMP</journal-id><journal-title-group><journal-title>Journal of Modern Physics</journal-title></journal-title-group><issn pub-type="epub">2153-1196</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jmp.2015.613197</article-id><article-id pub-id-type="publisher-id">JMP-60735</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  The Use of Phase Portraits for the Study of the Generation-Recombination Processes in Semiconductor
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>afur</surname><given-names>Gulyamov</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>A.</surname><given-names>G. Gulyamov</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>A.</surname><given-names>Q. Ergashev</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>B.</surname><given-names>T. Abdulazizov</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref></contrib></contrib-group><aff id="aff3"><addr-line>Namangan State University, Namangan, Uzbekistan</addr-line></aff><aff id="aff2"><addr-line>Physico-Technical Institute, Academy of Sciences of Uzbekistan, Tashkent, Uzbekistan</addr-line></aff><aff id="aff1"><addr-line>Namangan Engineering-Pedagogical Institute, Namangan, Uzbekistan</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>gulyamov1949@mail.ru(AG)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>15</day><month>10</month><year>2015</year></pub-date><volume>06</volume><issue>13</issue><fpage>1921</fpage><lpage>1926</lpage><history><date date-type="received"><day>6</day>	<month>July</month>	<year>2015</year></date><date date-type="rev-recd"><day>accepted</day>	<month>26</month>	<year>October</year>	</date><date date-type="accepted"><day>29</day>	<month>October</month>	<year>2015</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>
 
 
  Theoretical investigation of generation-recombination processes in silicon, which has a lifetime of charge carriers 10
  <sup>-3</sup> s and capture cross sections of 10
  <sup>-16</sup> sm
  <sup>2</sup>. For the study uses a method of phase portraits, which are widely used in the theory of vibrations. It is shown that the form of phase portraits strongly depends on the frequency of exposure to the external variable deformation.
 
</p></abstract><kwd-group><kwd>Phase Portraits</kwd><kwd> The Generation of Charge Carriers</kwd><kwd> The Recombination of Charge Carriers</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>For the analysis of all possible movements (natural vibrations) is “phase plane”. In mechanics, statistical physics and introduction of the concept of “phase space”, in the one case, it degenerates into a plane space. Analysis of the oscillating system using phase portraits has the advantage of clarity, and is a powerful tool for the study of oscillatory processes [<xref ref-type="bibr" rid="scirp.60735-ref1">1</xref>] .</p><p>Variable external influence on the semiconductor structure may cause the alternation of the generation and recombination processes. This process of generation and recombination in nature is similar to the oscillatory process. Study of generation-recombination processes using phase portraits can give a fairly complete picture of the processes in semiconductors. In [<xref ref-type="bibr" rid="scirp.60735-ref2">2</xref>] , they were obtained according to the concentration and resistance to deformation of the semiconductor film. It was shown that these relationships gave graphs pretty clear ideas ongoing processes in semiconductor films, namely the distortion of the closed curves talking about the emergence of new recombination centers. The emergence of new recombination centers is due to structural changes in semiconductor films, which are associated with films fatigue due to long exposure deformation. Thus, the above dependency graphs are allowed to diagnose conditions of semiconductor materials. However, in [<xref ref-type="bibr" rid="scirp.60735-ref2">2</xref>] , the phase portraits of the processes of the semiconductor film were not considered.</p><p>In this paper, by using the phase portraits the effect of variable frequency strain on the semiconductor film is investigated.</p></sec><sec id="s2"><title>2. The Impact of the Variable Frequency Strain Concentration of Nonequilibrium Charge Carriers</title><p>Under the effect of the deformation occurs changes the band structure leads to changes in the concentration of charge carriers, and a redistribution of charge carriers between energy valleys [<xref ref-type="bibr" rid="scirp.60735-ref3">3</xref>] . These changes in turn lead to a change in resistance of the sample [<xref ref-type="bibr" rid="scirp.60735-ref4">4</xref>] . If the external deformation will be variable, then the change in resistance of the sample at the time occurs with some lag from the applied mechanical impact. Consider the response of the excess electron concentration ne a variable deformation. Deformation can induce the generation of electrons at a rate ge moreover, excessive electrons may come from another area of the sample [<xref ref-type="bibr" rid="scirp.60735-ref5">5</xref>] [<xref ref-type="bibr" rid="scirp.60735-ref6">6</xref>] . The concentration can be determined by the equation of continuity</p><disp-formula id="scirp.60735-formula1064"><label>, (1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/20-7502333x5.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/20-7502333x6.png" xlink:type="simple"/></inline-formula> the current density of electrons, t-the electron lifetime. In a homogeneous sample is given by the equation of continuity</p><disp-formula id="scirp.60735-formula1065"><label>. (2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/20-7502333x7.png"  xlink:type="simple"/></disp-formula><p>The solution of Equation (2) in the case where the deformation of e is a random function of time, is given by [<xref ref-type="bibr" rid="scirp.60735-ref6">6</xref>]</p><disp-formula id="scirp.60735-formula1066"><label>. (3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/20-7502333x8.png"  xlink:type="simple"/></disp-formula><p>Equation (3) can be used to determine the magnitude of the concentration of non-equilibrium reacts to changing deformation. Assume that the deformation varies sinusoidally</p><disp-formula id="scirp.60735-formula1067"><label>. (4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/20-7502333x9.png"  xlink:type="simple"/></disp-formula><p>Then, the band gap varies according to the same law</p><disp-formula id="scirp.60735-formula1068"><label>. (5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/20-7502333x10.png"  xlink:type="simple"/></disp-formula><p>here <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/20-7502333x11.png" xlink:type="simple"/></inline-formula>-the deformation potential [<xref ref-type="bibr" rid="scirp.60735-ref5">5</xref>] [<xref ref-type="bibr" rid="scirp.60735-ref7">7</xref>] . The rate of heat generation has an exponential dependence of the band gap and lattice temperature [<xref ref-type="bibr" rid="scirp.60735-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.60735-ref8">8</xref>] . In this case, the general thermodynamic considerations, the rate of generation of non-equilibrium carriers can be represented as follows</p><disp-formula id="scirp.60735-formula1069"><label>. (6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/20-7502333x12.png"  xlink:type="simple"/></disp-formula><p>Then, from the expression (3) and (6) we have</p><disp-formula id="scirp.60735-formula1070"><label>. (7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/20-7502333x13.png"  xlink:type="simple"/></disp-formula><p>By specifying an explicit form of time-dependent strain ε(t) can be obtained from the time dependence of the concentration n<sub>e</sub>(t) from time to time. Equation (7) in a fairly wide range of deformation function e(t<sub>0</sub>) allows to describe the change of the carrier density.</p><p>To demonstrate the method of preparation according n<sub>e</sub>-ε, we consider a weak deformation, i.e.<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/20-7502333x14.png" xlink:type="simple"/></inline-formula>. Then,</p><p>the expression (7) can be written as</p><disp-formula id="scirp.60735-formula1071"><label>. (8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/20-7502333x15.png"  xlink:type="simple"/></disp-formula><p>Let the sample is applied a sinusoidal deformation half-period <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/20-7502333x16.png" xlink:type="simple"/></inline-formula> <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/20-7502333x17.png" xlink:type="simple"/></inline-formula> and other times there is no deformation. Then from (8) we have</p><disp-formula id="scirp.60735-formula1072"><label>. (9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/20-7502333x18.png"  xlink:type="simple"/></disp-formula><p>After the end of the deformation in the formula (8), the upper limit of the integral is to be p/w</p><disp-formula id="scirp.60735-formula1073"><label>. (10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/20-7502333x19.png"  xlink:type="simple"/></disp-formula><p>Changing the frequency of external influence is making a tangible contribution to the concentration vs. time. Since the change in the external effects of low frequencies in the non-equilibrium charge carrier concentration will increase and decrease gradually. That is, there will be some period of time between generation and recombination of charge carriers from zone to zone. And when the frequency reaches the maximum values decrease and increase of charge carriers becomes imperceptible. In other words, the generation-recombination processes do not have time to react to changes in the band spectrum of the semiconductor and the concentration of charge carriers takes a steady average.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref> and <xref ref-type="fig" rid="fig2">Figure 2</xref> present illustrative examples of the concentration of charge carriers of the time under the influence of deformation of the variable frequency. The growth rate was 3 Hz frequency in every 0.01 second. When plotting used expression (9).</p><p>Consider the case where a variable deformation influences with increasing frequency in the semiconductor structure. Assume that the initial generation and charge carrier lifetime is respectively g<sub>0</sub> = 10<sup>18</sup> and τ = 10<sup>−</sup><sup>3</sup> [<xref ref-type="bibr" rid="scirp.60735-ref9">9</xref>] . With the growth rate of deformation of the variable ω<sub>0</sub> = 3 Hz to ω = 210 Hz (see <xref ref-type="fig" rid="fig1">Figure 1</xref>(a)) a decrease in the concentration ranges. For example, if the frequency ω<sub>0</sub> = 3 Hz change in the concentration range was n ≈ 0 - 2 &#215; 10<sup>15</sup> cm<sup>−</sup><sup>3</sup>, then the frequency ω<sub>0</sub> = 210 Hz range was the change in concentration n ≈ (5 &#215; 10<sup>12</sup>) - (1.9 &#215; 10<sup>1</sup><sup>5</sup>) cm<sup>−</sup><sup>3</sup>. When the frequency of the alternating deformation reaches ω = 900 Hz (see <xref ref-type="fig" rid="fig1">Figure 1</xref>(b)) changes in the concentration range of ≈(3 &#215; 10<sup>13</sup>) - (7.5 &#215; 10<sup>1</sup><sup>4</sup>) cm<sup>−</sup><sup>3</sup>.</p><p>With increasing frequency of the alternating deformation range of the concentration of charge carriers takes some average. As can be seen in <xref ref-type="fig" rid="fig2">Figure 2</xref> when a variable frequency is 30 kHz deformation range of the charge carrier concentration varies between 10<sup>15</sup> cm<sup>−</sup><sup>3</sup> to increase the frequency up to 54 kHz. This suggests that the concentration of nonequilibrium charge carriers acquired its steady average.</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Plots of nonequilibrium charge carrier concentration of the time under the influence of variable deformation with increasing frequency: (a) increase the frequency of a variable deformation of ω<sub>0</sub> = 3 Hz to ω = 210 Hz (b) with the grown rate of deformation of the variable ω<sub>0</sub> = 3 Hz to ω = 900 Hz</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/20-7502333x20.png"/></fig><p>Apparently by prolonged exposure of the periodic pattern deformation in semiconductor recombination centers may occur in crystal lattices. As a result, the development of the processes of generation and recombination can take a different turn, and accordingly the concentration of charge carriers may have a deviation from the established mean.</p><p>The case where the frequency itself changes periodically variable deformation, ranges from 0 to 20 kHz is shown in <xref ref-type="fig" rid="fig3">Figure 3</xref>. Analysis of <xref ref-type="fig" rid="fig3">Figure 3</xref> shows that as the frequency variable deformation reaches a maximum value, set the average concentration becomes a value within 10 &#215; 15 cm<sup>−</sup><sup>3</sup>. When the frequency of the variable deformation reaches a minimum value, a range of changes in the concentration ranges from 0 - 2 &#215; 10<sup>15</sup> cm<sup>−</sup><sup>3</sup>. This is repeated for two seconds to 8 times.</p></sec><sec id="s3"><title>3. Phase Portraits of the Concentration of Charge Carriers from the Rate of Change of the Charge Carriers Concentration</title><p>To obtain the phase portrait is usually considered the dependence of some parameter x of its derivative<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/20-7502333x21.png" xlink:type="simple"/></inline-formula>. Since</p><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Graph of the nonequilibrium carrier concentration from time to time the growth rate of deformation of the variable ω<sub>0</sub> = 3 Hz to ω = 54 kHz</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/20-7502333x22.png"/></fig><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Frequency modulation effect of variable strain on the semiconductor. The concentration of charge carriers is changed periodically</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/20-7502333x23.png"/></fig><p>the deformation by changing the energy spectrum affects the generation-recombination processes, to obtain the phase portrait of such processes must be obtained from the dependence of the concentration n of its derivative</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/20-7502333x24.png" xlink:type="simple"/></inline-formula>, or the dependence of resistance <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/20-7502333x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/20-7502333x25.png" xlink:type="simple"/></inline-formula> on its derivative<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/20-7502333x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/20-7502333x25.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/20-7502333x26.png" xlink:type="simple"/></inline-formula>.</p><p>Consider the case of silicon, which has a lifetime of charge carriers 10<sup>−</sup><sup>3</sup> seconds capture cross 10<sup>−</sup><sup>16</sup> cm<sup>2</sup> [<xref ref-type="bibr" rid="scirp.60735-ref9">9</xref>] .</p><p>The phase portrait is the largest area when the frequency of the alternating deformation is 1000 Hz (see <xref ref-type="fig" rid="fig4">Figure 4</xref>(a)). That is, the phase portrait is the largest area when<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/20-7502333x27.png" xlink:type="simple"/></inline-formula>, here the generation and recombination of charge carriers takes place gradually. With the change in frequency of the alternating deformation changes the shape of the phase portraits (see Figures 4(b)-(f)). In the case where the frequency of the alternating deformation increases, and the product of frequency and time of life become<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/20-7502333x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/20-7502333x28.png" xlink:type="simple"/></inline-formula>, the shape of the phase portrait is committed to the horizontal segment (see <xref ref-type="fig" rid="fig4">Figure 4</xref>(f)), that is, generation and recombination of charge carriers is almost not felt. In the opposite case, when the frequency of the alternating deformation decreases, and the product of frequency and time of life become<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/20-7502333x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/20-7502333x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/20-7502333x29.png" xlink:type="simple"/></inline-formula>, the shape of the phase portrait is committed to the vertical segment (see <xref ref-type="fig" rid="fig4">Figure 4</xref>(f)), is the generation and recombination of charge carriers is shown sharply.</p><p>Three-dimensional phase portrait is shown in <xref ref-type="fig" rid="fig5">Figure 5</xref>. This graph is obtained for the case when the charge carrier lifetime 10<sup>−</sup><sup>2</sup> s, and capture cross sections of 10<sup>−</sup><sup>16</sup> cm<sup>2</sup>.</p><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> The phase portraits of the concentration on their derivatives</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/20-7502333x30.png"/></fig><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> The phase portraits of the concentration, their derivatives and variable- frequency distortion</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/20-7502333x31.png"/></fig></sec><sec id="s4"><title>4. Conclusions</title><p>Thus, the shape of the phase portraits allows us to determine how to respond to the frequency change of the semi- conductor strain. Namely, with increasing frequency variable strain rate and the product carrier lifetime become<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/20-7502333x32.png" xlink:type="simple"/></inline-formula>, the shape of the phase portrait is committed to the horizontal segment. If the variable frequency distortion and carrier lifetime of the product become<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/20-7502333x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/20-7502333x33.png" xlink:type="simple"/></inline-formula>, the shape of the phase portrait is committed to the vertical segment. The phase portrait is the largest area when<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/20-7502333x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/20-7502333x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/20-7502333x34.png" xlink:type="simple"/></inline-formula>.</p><p>This work was supported by government grants of Uzbekistan F2-OT-O-15494 “Improving the efficiency of the radiators, solar cells and other optoelectronic devices based on quantum dot research and hetero nanostructure exciton and bipolaron transport phenomena”, F2-FA-0-97004 “Dynamics of solitons in ultracold quantum gases, and optical systems and synergistic processes in multicomponent semiconductors” and F2-21 “Mathematical modelling of the determination of the density of surface states at the semiconductor-insulator”.</p></sec><sec id="s5"><title>Cite this paper</title><p>GafurGulyamov,A. G.Gulyamov,A. Q.Ergashev,B. T.Abdulazizov, (2015) The Use of Phase Portraits for the Study of the Generation-Recombination Processes in Semiconductor. Journal of Modern Physics,06,1921-1926. doi: 10.4236/jmp.2015.613197</p></sec></body><back><ref-list><title>References</title><ref id="scirp.60735-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Andronov, А.А., Vitt, А.А. and Haykin, C.E. (1985) Vibration Theory. Nauka, Moscow. (In Russian)</mixed-citation></ref><ref id="scirp.60735-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Shamirzaev, S.H., Gulyamov, G., Dadamirzaev, M.G., Sharibaev, N.Y. and Gulyamov, A.G. (2011) FIP (Fyzicheckaya Injeneriya Poverkhnosti), 1, 68-71.</mixed-citation></ref><ref id="scirp.60735-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Bir, G.L. and Pikus, G.E. Symmetry and Strain-Induced Effects in Semiconductors. Nauka, Moscow, 1972; Wiley, New York, 1975.</mixed-citation></ref><ref id="scirp.60735-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Bonch-Bruevich, V. and Kalashnikov, S. (1977) Semiconductor Physics. Nauka, Moscow. (In Russian)</mixed-citation></ref><ref id="scirp.60735-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Kireev, P.S. (1975) Semiconductor Physics. Vishaya shkola, Moscow. (In Russian)</mixed-citation></ref><ref id="scirp.60735-ref6"><label>6</label><mixed-citation publication-type="book" xlink:type="simple">Blakemore, J.S. (1974) The Theory of the Rate of Recombination in Semiconductors. In: Landsberg, P., Ed., Chapter 16 in the Book: Problems of Thermodynamics and Statistical Physics, Mir, Moscow, 399-430. (In Russian)</mixed-citation></ref><ref id="scirp.60735-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Gulyamov, A.G. (2015) World Journal of Condensed Matter Physics, 5, 60-65.</mixed-citation></ref><ref id="scirp.60735-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Ridley, B.K. (1999) Quantum processes in Semiconductors. Clarendon Press, Oxford.</mixed-citation></ref><ref id="scirp.60735-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Bonch-Bruevich, V.L., Zvyagin, I.P., Karpenko, I.V. and Mironov, A.G. (1987) Problems in Physics of Semiconductors. Nauka, Moscow. (In Russian)</mixed-citation></ref></ref-list></back></article>